Beerli P, Felsenstein J
Department of Genetics, University of Washington, Seattle, WA 98195-7360, USA.
Proc Natl Acad Sci U S A. 2001 Apr 10;98(8):4563-8. doi: 10.1073/pnas.081068098. Epub 2001 Apr 3.
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.
开发了一种基于合并的最大似然估计器，用于估计不等迁移率和不同亚种群大小的情况。该方法采用马尔可夫链蒙特卡罗方法来研究具有分支长度和迁移事件的可能系谱。通过使用来自四群体n岛模型和源 - 汇种群模型的模拟数据，展示了新方法的特性。我们在migrate中编码的估计方法与genetree进行了测试比较；两个程序产生的似然表面非常相似。即使马尔可夫链从不利参数开始，该算法也能相当快地收敛到估计值。该方法被用于通过使用来自三个人类群体的线粒体DNA数据来估计尼罗河流域的基因流。
