Penalized estimation of haplotype frequencies
Penalized estimation of haplotype frequencies
复制标题
DOI:
10.1093/bioinformatics/btn236
复制
发表时间:
2008-07-15
期刊:
影响因子:
5.8
通讯作者:
Lange, Kenneth
中科院分区:
文献类型:
--
作者:
Ayers, Kristin L.;Lange, Kenneth
Motivation: Low haplotype diversity and linkage disequilibrium are the rule in short genomic segments. This fact suggests that parsimony should be enforced in estimation of haplotype frequencies. The current article introduces a diversity penalty that automatically discards potential haplotypes with low explanatory power. The standard EM algorithm for haplotype frequency estimation can accommodate the penalty if one passes over to a more general minorizemaximize (MM) scheme for estimation.Results: Our new MM algorithm converges in fewer iterations, eliminates marginal haplotypes from further consideration and reduces the computational complexity of each iteration. Estimation by the MM algorithm also improves haplotyping and genotype imputation compared to naive application of the EM algorithm. Thus, the MM algorithm is a useful substitute for the EM algorithm. Compared to the most sophisticated current methods of haplotyping and genotype imputation, the MM algorithm is slightly less accurate but at least an order of magnitude faster.