Estimating effective population size and mutation rate from sequence data using Metropolis-Hastings sampling - PubMed (original) (raw)

Comparative Study

Estimating effective population size and mutation rate from sequence data using Metropolis-Hastings sampling

M K Kuhner et al. Genetics. 1995 Aug.

Abstract

We present a new way to make a maximum likelihood estimate of the parameter 4N mu (effective population size times mutation rate per site, or theta) based on a population sample of molecular sequences. We use a Metropolis-Hastings Markov chain Monte Carlo method to sample genealogies in proportion to the product of their likelihood with respect to the data and their prior probability with respect to a coalescent distribution. A specific value of theta must be chosen to generate the coalescent distribution, but the resulting trees can be used to evaluate the likelihood at other values of theta, generating a likelihood curve. This procedure concentrates sampling on those genealogies that contribute most of the likelihood, allowing estimation of meaningful likelihood curves based on relatively small samples. The method can potentially be extended to cases involving varying population size, recombination, and migration.

PubMed Disclaimer

Similar articles

Cited by

References

    1. Theor Popul Biol. 1975 Apr;7(2):256-76 - PubMed
    1. Genet Res. 1992 Apr;59(2):139-47 - PubMed
    1. Proc Natl Acad Sci U S A. 1991 Oct 1;88(19):8720-4 - PubMed
    1. J Mol Evol. 1980 Dec;16(2):111-20 - PubMed

Publication types

MeSH terms

Substances

LinkOut - more resources