Journal
GENETICS
Volume 170, Issue 2, Pages 955-967Publisher
GENETICS
DOI: 10.1534/genetics.104.038349
Keywords
-
Categories
Funding
- NIGMS NIH HHS [GM-40282, R01 GM040282] Funding Source: Medline
Ask authors/readers for more resources
This article presents an efficient importance-sampling method for computing the likelihood of the effective size of a population under the coalescent model of Berthier et al. Previous computational approaches, using Markov chain Monte Carlo, required many minutes to several hours to analyze small data sets. The approach presented here is orders of magnitude faster and can provide an approximation to the likelihood curve, even for large data sets, in a matter of seconds. Additionally, confidence intervals on the estimated likelihood curve provide a useful estimate of the Monte Carlo error. Simulations show the importance sampling to be stable across a wide range of scenarios and show that the N-e estimator itself performs well. Further simulations show that the 95% confidence intervals around the N-e estimate are accurate. User-friendly software implementing the algorithm for Mac, Windows, and Unix/Linux is available for download. Applications of this computational framework to other problems are discussed.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available