LAMARC 2.0 Maximum Likelihood and Bayesian Estimation of
LAMARC 2.0 Maximum Likelihood and Bayesian Estimation of
LAMARC 2.0 Maximum Likelihood and Bayesian Estimation of
The Author 2006. Published by Oxford University Press. All rights reserved. For Permissions, please email: [email protected]
The online version of this article has been published under an open access model. Users are entitled to use, reproduce, disseminate, or display the open access
version of this article for non-commercial purposes provided that: the original authorship is properly and fully attributed; the Journal and Oxford University
Press are attributed as the original place of publication with the correct citation details given; if an article is subsequently reproduced or disseminated not in its
entirety but only in part or as a derivative work this must be clearly indicated. For commercial re-use, please contact [email protected]
LAMARC 2.0
3 IMPLEMENTATION 4 DISCUSSION
3.1 Search strategy 4.1 Model assumptions
LAMARC 2.0 provides several mechanisms to improve its search LAMARC 2.0 assumes that individuals are drawn from panmictic
efficiency. Metropolis-Coupled MCMC or ‘heating’ allows auxil- subpopulations and that the subpopulation structure has been con-
liary searches with more permissive acceptance criteria to act as stant throughout the lifespan of the underlying coalescent tree. It
‘scouts’ for the main analysis (Geyer, 1991a). For likelihood is not suitable for populations which have recently diverged from
analyses, multiple replicated searches can be combined using a common ancestor. It assumes that the rate at which a lineage
reverse logistic regression (Geyer, 1991b). For Bayesian analyses, immigrates into a population is independent of the size of both
the Bayesian priors and the ratio of parameter change steps to source and recipient populations. It also assumes that exponential
genealogy change steps can be set by the user. growth rates and immigration rates have been constant throughout
the lifespan of the coalescent tree and that recombination rate
3.2 Mutational models does not vary by position, subpopulation or with time. Finally, it
LAMARC 2.0 offers the Felsenstein 84 (F84) and General Time- assumes that the variation being observed is neutral, though puri-
Reversible (GTR) models for DNA or RNA data, and for SNP data fying selection removing harmful mutations does not disrupt the
769
M.K.Kuhner
Jon Yamato and Wang Yi. Development of this program was sup- Geyer,C.J. (1991a) Markov chain Monte Carlo maximum likelihood. In Keramidas
ported by NIH grant 5R01GM51929-11 to M.K. Funding to pay the (ed.), Computing Science and Statistics: Proceedings of 23rd Symposium on the
Interface, Interface Foundation, Fairfax Station, pp. 156–163.
Open Access publication charges was provided by the National
Geyer,C.J. (1991b) Estimating normalizing constants and reweighting mixtures
Institutes of health grant 5R01CM51929-11 to M.K. in Markov chain Monte Carlo. Technical Report No. 568, School of Statistics,
University of Minnesota, MN revised 1994.
Conflict of Interest: none declared. Kuhner,M.K. and Felsenstein,J. (2000) Sampling among haplotype resolutions in
a coalescent-based genealogy sampler. Genet. Epidemiol., 19 (Suppl. 1), S15–S21.
Kuhner,M.K. et al. (1995) Estimating effective population size and mutation rate from
sequence data using Metropolis-Hastings sampling. Genetics, 140, 1421–1430.
REFERENCES Kuhner,M.K. et al. (1998) Maximum likelihood estimation of population growth rates
Beerli,P. and Felsenstein,J. (1999) Maximum-likelihood estimation of migration rates based on the coalescent. Genetics, 149, 429–434.
and effective population numbers in two populations using a coalescent approach. Kuhner,M.K. et al. (2000a) Usefulness of single nucleotide polymorphism data for
Genetics, 152, 763–773. estimating population parameters. Genetics, 156, 439–447.
Beerli,P. and Felsestein,J. Maximum likelihood estimation of a migration matrix and Kuhner,M.K. et al. (2000b) Maximum likelihood estimation of recombination rates
effective population sizes in n subpopulations using a coalescent approach. Proc. from population data. Genetics, 156, 1393–1401.
Natl Acad. Sci. USA, 98, 4563–4568. Ohta,T. and Kimura,M. (1973) A model of mutation appropriate to estimate the number
Felsenstein,J. and Churchill,G.A. (1996) A hidden Markov model approach to variation of electrophoretically detectable alleles in a finite population. Genet. Res., 22,
770