Maximum likelihood estimation of a migration matrix and effective population sizes in n subpopulations by using a coalescent approach

AUTOR(ES)

Beerli, Peter

FONTE

The National Academy of Sciences

RESUMO

A maximum likelihood estimator based on the coalescent for unequal migration rates and different subpopulation sizes is developed. The method uses a Markov chain Monte Carlo approach to investigate possible genealogies with branch lengths and with migration events. Properties of the new method are shown by using simulated data from a four-population n-island model and a source–sink population model. Our estimation method as coded in migrate is tested against genetree; both programs deliver a very similar likelihood surface. The algorithm converges to the estimates fairly quickly, even when the Markov chain is started from unfavorable parameters. The method was used to estimate gene flow in the Nile valley by using mtDNA data from three human populations.

Documentos Relacionados

• Maximum Likelihood Estimation of Population Parameters
• Using rare mutations to estimate population divergence times: A maximum likelihood approach
• The Maximum Likelihood Estimation of the Number of Self-Sterility Alleles in a Population
• Methodology and Accuracy of Estimation of Quantitative Trait Loci Parameters in a Half-Sib Design Using Maximum Likelihood
• Alu Evolution in Human Populations: Using the Coalescent to Estimate Effective Population Size