Maximum likelihood estimation of a migration matrix and effective population sizes in n subpopulations by using a coalescent approach
Maximum likelihood estimation of a migration matrix and effective population sizes in n subpopulations by using a coalescent approach
复制标题
DOI：
10.1073/pnas.081068098
复制
发表时间：
2001-04-10
影响因子：
11.1
通讯作者：
Felsenstein, J
中科院分区：
文献类型：
--
作者：
Beerli, P;Felsenstein, J
对于不同的迁移率和不同的子总体大小，提出了一种基于合并的极大似然估计。该方法使用马尔可夫链蒙特卡罗方法来调查可能的家谱与分支长度和迁移事件。用四种群n岛模型和源-汇种群模型的模拟数据说明了新方法的性质。我们的估计方法编码在迁移测试对GENETREE，两个程序提供了一个非常相似的可能性表面。该算法收敛到估计相当快，即使马尔可夫链是从不利的参数。该方法被用来估计基因流在尼罗河流域使用mtDNA数据从三个人群。
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.