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

被引：1418

作者：

Beerli, P ^{[1
]}

Felsenstein, J ^{[1
]}

机构：

[1] Univ Washington, Dept Genet, Seattle, WA 98195 USA

来源：

PROCEEDINGS OF THE NATIONAL ACADEMY OF SCIENCES OF THE UNITED STATES OF AMERICA | 2001年 / 98卷 / 08期

关键词：

coalescence theory; population genetics;

D O I：

10.1073/pnas.081068098

中图分类号：

O [数理科学和化学]; P [天文学、地球科学]; Q [生物科学]; N [自然科学总论];

学科分类号：

07 ; 0710 ; 09 ;

摘要：

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.

引用

页码：4563 / 4568

页数：6

共 42 条

[1]

[Anonymous], 1979, Monte Carlo Methods, DOI DOI 10.1007/978-94-009-5819-7

[2] Inference from gene trees in a subdivided population [J].