Welcome, Guest. Please Login
YaBB - Yet another Bulletin Board
 
  HomeHelpSearchLogin  
 
dN/dS versus rateN/rateS (Read 2014 times)
Danny
YaBB Newbies
*
Offline



Posts: 31
dN/dS versus rateN/rateS
Mar 13th, 2009 at 12:11am
 
Dear Sergei,

Page 13 of the hyphybook says that these two ratios are close but not equal and this is true for my data using the MG model. I would like to hear or see a reference to a more detailed explanation of the difference between these two ratios.  My intuition is that the nt substitution and equilibrium frequency values in the Q matrix attempt to account for (normalize for) changes that are due simply to background nt substitution rates.  However, dS is based partly on the codon frequencies.  In contrast, in the MG Q matrix the frequencies used are the nt freqs at each codon position.  Is this part of the reason for the difference in the ratios?  The nt substitution rates do not seem to be part of the dS calculation; is this another reason for the difference?  And then there's the T/S which seems rough.  I guess one way to ask is, is there, in principle, a way to calculate dN and dS such that dN/dS = rateN/rateS.  For the purposes of calculating omega in HYPHY though, I take it that rateN/rateS is the answer and dN/dS is some kind of ad hoc approximation?

-Danny
Back to top
 
 
IP Logged
 
Sergei
YaBB Administrator
*****
Offline


Datamonkeys are forever...

Posts: 1658
UCSD
Gender: male
Re: dN/dS versus rateN/rateS
Reply #1 - Mar 13th, 2009 at 12:24am
 
Dear Danny,

Spencer Muse's paper on this (Multimedia File Viewing and Clickable Links are available for Registered Members only!!  You need to Login Login) is still probably the best reference to read. I am not sure I can give you a better intuition into why dN/dS and omega are different, other than to note that you could actually (for MG94 or GY94) write out exactly what the ratio between dN/dS and omega is -- it will depend entirely on base frequencies and nucleotide bias parameters (e.g. Ts/Tv), but not on non-synonymous or synonymous rates. If I remember correctly, the ratio becomes one when all nucleotides are equally frequent and the nucleotide model is F81.

Cheers,
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