REML法和MCMC法在数量性状核心家系遗传方差分量模型中参数估计的比较

摘要

目的比较限制性最大似然估计(REML)和马尔可夫链蒙特卡罗(MCMC)方法在数量性状遗传方差分量模型中参数估计的偏差和精度.方法计算机模拟50个数量性状核心家系数据集,运用SAS软件中的PROC MIXED程序和WinBUGS软件进行参数估计,比较不同样本含量时两法参数估计结果.结果在本试验参数指定条件下,对固定效应参数,两法估计的相对偏差低于2%;对随机效应,估计偏差较高.小样本时,MCMC法估计结果比REML法更接近真值.结论在模拟数量性状方差分量模型时,大样本家系资料可根据实际情况选用两法;小样本家系资料,MCMC法可得到比REML法更稳健的估计,推荐使用.