LAMARC 2.0: maximum likelihood and Bayesian estimation of population parameters

UNLABELLED We present a Markov chain Monte Carlo coalescent genealogy sampler, LAMARC 2.0, which estimates population genetic parameters from genetic data. LAMARC can co-estimate subpopulation Theta = 4N(e)mu, immigration rates, subpopulation exponential growth rates and overall recombination rate, or a user-specified subset of these parameters. It can perform either maximum-likelihood or Bayesian analysis, and accomodates nucleotide sequence, SNP, microsatellite or elecrophoretic data, with resolved or unresolved haplotypes. It is available as portable source code and executables for all three major platforms. AVAILABILITY LAMARC 2.0 is freely available at http://evolution.gs.washington.edu/lamarc

[1]  Jon A Yamato,et al.  Usefulness of single nucleotide polymorphism data for estimating population parameters. , 2000, Genetics.

[2]  J. Felsenstein,et al.  A Hidden Markov Model approach to variation among sites in rate of evolution. , 1996, Molecular biology and evolution.

[3]  Motoo Kimura,et al.  A model of mutation appropriate to estimate the number of electrophoretically detectable alleles in a finite population*. , 1973, Genetical research.

[4]  J. Felsenstein,et al.  Sampling among haplotype resolutions in a coalescent‐based genealogy sampler , 2000, Genetic Epidemiology.

[5]  C. Geyer.  Estimating Normalizing Constants and Reweighting Mixtures , 1994 .

[6]  J. Felsenstein,et al.  Maximum-likelihood estimation of migration rates and effective population numbers in two populations using a coalescent approach. , 1999, Genetics.

[7]  Jon A Yamato,et al.  Maximum likelihood estimation of recombination rates from population data. , 2000, Genetics.

[8]  Peter Beerli,et al.  Maximum likelihood estimation of a migration matrix and effective population sizes in n subpopulations by using a coalescent approach , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[9]  Jon A Yamato,et al.  Maximum likelihood estimation of population growth rates based on the coalescent. , 1998, Genetics.

[10]  Jon A Yamato,et al.  Estimating effective population size and mutation rate from sequence data using Metropolis-Hastings sampling. , 1995, Genetics.

[11]  C. Geyer.  Markov Chain Monte Carlo Maximum Likelihood , 1991 .