论文部分内容阅读
在全基因组混合模型关联分析（GWMMAS）中,利用剩余多基因遗传力代替方差比值（剩余多基因方差/误差方差）,将多基因遗传力求解限制在（0,1）区间内。引入R语言Rcpp Armadillo程序包中的极速线性模型拟合函数（fast Lm Pure函数）快速估计单核苷酸多态性（SNP）效应和完整LMM最大似然值。从由GBLUP估计的性状基因组遗传力出发,逐个高通量SNP的GWMMAS约需4次全基因组回归扫描。当仅关注EMMAX法估计的大效应或高显著水准标记时,GWMMAS运行时间缩短在两次扫描之内。与采用lm