Question: Effective Population Size (Ne) From Nucleotide Diversity (Pi) And Mutation Rate (Mu)
gravatar for Stingery
7.9 years ago by
Stingery20 wrote:

I am a bit confused and would be grateful for some advice. I know that the effective population size can be calculated using the following relation (for diploid organisms, under neutrality): Ne = pi / (4 * mu)

By copying the work of someone else, I mistakenly used Ne = 4 * (pi / mu) This, however, resulted in much more realistic estimates, despite being larger by a factor of 16. Estimates obtained from the first equation do not make sense, especially when compared to the results of another estimation method of Ne (extended Bayesian skyline plots), which are very close to the estimates obtained from the second equation. Because I only found the first equation in published work, am I (and my colleague) wrong to use the second equation, or is the latter one indeed valid?


nucleotide population mutation • 5.7k views
ADD COMMENTlink written 7.9 years ago by Stingery20

Would you mind sharing the reference for that paper where they used the other equation? As I am sure you are aware, the pi=4mu*Ne calculation is dependent on an idealized diploid population. I am not sure why both calculations give such different values. What are your expectations based on anyway?

ADD REPLYlink modified 7.9 years ago • written 7.9 years ago by Whetting1.5k

More precisely, it is theta, not pi. Pi is just an estimator of theta. And you cannot take skyline plot as the answer. I forgot how skyline plot works, but it may be more sensitive to sequencing errors. You need to give more details, including the skyline plot.

ADD REPLYlink written 7.9 years ago by lh332k

In my own limited experience, it can be difficult to get the units correct for pi and mu. Maybe you have an inconsistency in the computation of pi or mu that accounts for the factor of 16?

ADD REPLYlink written 7.9 years ago by bdemarest460

I think the factor 16 is easy enough. 4x = simply 16 times as large as x/4. And that is what he writes, where x=pi/mu.

ADD REPLYlink written 7.9 years ago by Chris Evelo10k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 880 users visited in the last hour