《混合线性模型软件包的介绍.docx》由会员分享,可在线阅读,更多相关《混合线性模型软件包的介绍.docx(2页珍藏版)》请在三一文库上搜索。
1、混合线性模型软件包的介绍混合线性模型是一般线性模型的延伸和拓展,在作物、林木、动物、水产育种和科研中应用广泛。相对于一般线性模型,它能处理缺失值和不平衡数据,可以支持数据的方差不齐次和不独立,使得数据分析更准确和高效。分析混合线性模型的软件很多,这里将其分为R包和非R包。R包:nlme:是一个在S-Plus应用广泛的混线性模型包,后来转换到R平台上,对于镶嵌结构(nested)的随机因子定义简单,但对于交叉的随机因子(crossed)定义困难。具有多个功能,比如lme应用于线性混合模型,nlme应用于非线性混合模型。可以定义复杂的方差结构,不支持广义线性混合模型(GLMM)。lme4:是nlm
2、e的进一步发展,比nlme运行速度要快,支持GLMM,但很难处理交叉的随机因子。, 比nlme运行速度要快MCMCglmm:用马尔可夫链蒙特卡罗(Markov Chain Monte Carlo,MCMC)的方法拟合模型,贝叶斯先验分布,可以定义一些复杂的方差结构(heterogeneous yes, AR1 no)。asreml:是ASReml软件的R版本,运算速度快,支持复杂的模型(G矩阵和R矩阵),支持系谱信息和多性状分析,在动物、作物、林木、水产育种和科研中应用广范。glmmADMB:是ADMB软件的R版本,很灵活,但是运行速度很慢。非R包:l ASReml 商业软件:对于第三世界的科
3、研单位是免费的,有单机版(Win、Linux、Mac)和R版(ASReml-R),也有窗口化的版本(GenStat),应用稀疏矩阵和Ai算法,速度很快,广泛应用于植物和动物育种,支持随机因子的矩阵定义(G)和残差的矩阵定义(R),Splines也被很好的整合,对于广义线性模型,应用的是PQL方法。l ADMB:自动模型微分(Automatic Differentiation Model Builder),主要应用在森林、水产和野生动物中,开始时是个商业软件,现在开源了,支持非线性混合模型。l SAS 商业软件:n PROC MIXED:一般线性混合模型(LMM),应用广泛,但是速度比较慢。n
4、PROC GLIMMIX:增加了广义线性模型(GLMM),它现在支持了Laplace approximation和adaptive Gaussian quadrature方法,但对于复杂的模型,用的还是PQL方法。n HPMIXED:是MIXED的改进版,在速度上有明显的提升,但是支持的模型较少。n PROC NLMIXED:支持非线性混合模型功能介绍:软件函数方法置信区间随机(G矩阵)残差(R矩阵)其它nlmelmeMLREML标准差(SD)Nestedblock相关联(corStruct)异质性(varStruct)速度较慢没有给出方差组分lme4lmerMLREML标准误(se)Nest
5、edDiagonalblockNO不支持系谱不支持多性状分析ASremlasremlAIREML标准误(se)NestedCrossBlockSplinesARCorStrUnstructedMaternFactor Analysis速度快支持系谱支持多性状分析SASMixedHPMIXEDGLIMMIXMLREML标准误(se)NestedCrossBlockNO速度慢支持系谱支持多性状分析MCMCglmmMCMCglmmMCMCCredible intervalsNestedCrossBlock速度慢支持系谱支持多性状分析StataXtmixedXtregLMREMLNestedCrossBlockUnstructed功能方法函数NlmeLme4MCMCglmmASRemlSummaryEstimateTDfPEstimateStd errTDfPPost.meanCIEff.sampleLoglikeVarcompStd.errZ-ratioConstraintCoef效应值Fixed(BLUE)Ranef(BLUP)LoglikPlotPredictPedIntervalsUpdateLRT