Estimating Effective Population Size and Mutation Rate from Sequence Data Using Metropolis-Hastings Sampling

AUTOR(ES)

Kuhner, M. K.

RESUMO

We present a new way to make a maximum likelihood estimate of the parameter 4N(e)μ (effective population size times mutation rate per site, or θ) 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 θ must be chosen to generate the coalescent distribution, but the resulting trees can be used to evaluate the likelihood at other values of θ, 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.

Documentos Relacionados

• Description and validation of a method for simultaneous estimation of effective population size and mutation rate from human population data.
• Estimating Effective Population Size or Mutation Rate Using the Frequencies of Mutations of Various Classes in a Sample of DNA Sequences
• A Phylogenetic Estimator of Effective Population Size or Mutation Rate
• A New Method for Estimating the Effective Population Size from Allele Frequency Changes
• A Generalized Approach for Estimating Effective Population Size from Temporal Changes in Allele Frequency