Rate matrix estimation from site frequency data

详细信息    查看全文

• 作者：Conrad J. Burden^a ; ^b ; ^{conrad.burden@anu.edu.au} ; Yurong Tang^a ; ^{yurong.tang@anu.edu.au}
• 关键词：Multi-allele Wright&ndash ; Fisher ; Decoupled Moran ; Neutral evolution ; Evolutionary rate matrices
• 刊名：Theoretical Population Biology
• 出版年：2017
• 出版时间：February 2017
• 年：2017
• 卷：113
• 期：Complete
• 页码：23-33
• 全文大小：889 K
• 卷排序：113

文摘

A procedure is described for estimating evolutionary rate matrices from observed site frequency data. The procedure assumes (1) that the data are obtained from a constant size population evolving according to a stationary Wright–Fisher or decoupled Moran model; (2) that the data consist of a multiple alignment of a moderate number of sequenced genomes drawn randomly from the population; and (3) that within the genome a large number of independent, neutral sites evolving with a common mutation rate matrix can be identified. No restrictions are imposed on the scaled rate matrix other than that the off-diagonal elements are positive, their sum is ≪1≪1, and that the rows of the matrix sum to zero. In particular the rate matrix is not assumed to be reversible. The key to the method is an approximate stationary solution to the diffusion limit, forward Kolmogorov equation for neutral evolution in the limit of low mutation rates.