Dear Denise,
Quote:It is unclear to me however, how, in PAML, a single dS rate is used. Where does this single dS rate come from and how is it used in the calculations when estimating the omega category for individual sites? Denise
In a codon model (like MG94 or GY94) one can generally estimate two rates: alpha (~dS) and beta (~dN). In PAML, the synonymous rate is effectively fixed at 1 (for technical reasons you can't estimate dS and tree length separately, only their product, so dS = 1 is a convenient relative scale). In datamonkey (and HyPhy), the synonymous rate can be set to one on average, but allowed to vary from site to site. This is the main difference. I strongly encourage you to read Multimedia File Viewing and Clickable Links are available for Registered Members only!! You need to
where all these issues are touched upon.
The paper may be a bit technical, but it thorough. I'll be happy to assist you with further questions, if something in the manuscript is confusing.
Cheers,
Sergei