Welcome, Guest. Please Login
YaBB - Yet another Bulletin Board
 
  HomeHelpSearchLogin  
 
Statistics/info print outs (Read 2533 times)
jh26
YaBB Newbies
*
Offline


Feed your monkey!

Posts: 2
Statistics/info print outs
Jan 23rd, 2009 at 5:37am
 
Hi

I have a couple of (perhaps blindingly) obvious questions, but after numerous analyses using Hyphy/Datamoney and looking at the guides and boards, I still haven't found a way to get a print out of:

1) The number of synonymous and nonsynonymous substitutions per codon site, and
2) The overall codon tree length (as opposed to subs/nucleotide).

What are the necessary steps to achieve this? (Dataset of 35 sequences and 210 codon sites).

Also, Hyphy tends to crash after I try and save results to file...is this anything to do with an urgent need to upgrade my (Windows)computer?

Thank you very much in advance,

Jo  Huh
Back to top
 
 
IP Logged
 
Sergei
YaBB Administrator
*****
Offline


Datamonkeys are forever...

Posts: 1658
UCSD
Gender: male
Re: Statistics/info print outs
Reply #1 - Jan 24th, 2009 at 9:44am
 
Dear Jo,

Try using Standard Analyses>Miscellaneuos>Phylohandbook.bf>dSdN.bf
and following the prompts.

When the analysis is done, you should get a printout containing this list (obviously with different numbers for your data):

Code:
Total nucleotide sites :1698
Synonymous  sites	:373.961
Nonsynonymous  sites   :1252.14

Total tree lengths:
	E[syn]	  : 0.0556991
	E[non-syn]    : 0.0295723
	dS		: 0.252906
	dN		: 0.0401025

 



Sounds like you are looking for E[syn] (syn subs/codon site) and E[non-syn] (non-syn subs/codon site). The total length of the tree is the sum of the two.

S.

Back to top
 

Associate Professor
Division of Infectious Diseases
Division of Biomedical Informatics
School of Medicine
University of California San Diego
WWW WWW  
IP Logged
 
jh26
YaBB Newbies
*
Offline


Feed your monkey!

Posts: 2
Re: Statistics/info print outs
Reply #2 - Jan 25th, 2009 at 7:40am
 
Hi Sergei

That's great, thank you. However is there a way to get a print out of the number of nonsynonymous and synonymous subistitutions for each codon site? I think you can do this using the SLAC method in Datamonkey, e.g:

Codon      Observed S       Observed NS
                Changes            Changes      
1             1                      0
2             3                      1
3             1                      0
etc

Also, will this depend upon the selection method used (SLAC,REL,FEL),  or is there a quick and general way of calculating the above?

Thanks!
Jo
Back to top
 
 
IP Logged
 
Sergei
YaBB Administrator
*****
Offline


Datamonkeys are forever...

Posts: 1658
UCSD
Gender: male
Re: Statistics/info print outs
Reply #3 - Jan 25th, 2009 at 5:34pm
 
Dear Jo,

SLAC will indeed count substitutions for you.
There is no general way to do it, however. One must fit an evolutionary model to the data, infer ancestral sequences and then map mutations to the tree to obtain the two numbers -- they will therefore depend on the model chosen. SLAC is a reasonably good default model - it will perform well unless you have very long tree branches (>1 sub long), when SLAC can undercount the number of substitutions per site.

Best,
Sergei
Back to top
 

Associate Professor
Division of Infectious Diseases
Division of Biomedical Informatics
School of Medicine
University of California San Diego
WWW WWW  
IP Logged