原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
第七章个体遗传评定-blup法
对于y，总可以找到一个线性函数 ，它满足 ，使得： 亦即 不受固定效应影响。 若y服从正态分布，则 也服从正态分布，即： 的对数似然函数，也即约束似然函数的对数为： 其中，C为一常数，且有： n 为有记录的个体数， 为固定效应关联矩阵的秩。 求解上述对数约束似然函数关于 和 的最大值，便可得到 和 的REML估值。 于是： REML和ML的区别 思路上：ML是求观察值的似然函数的最大值；而REML是求观察值的某一特定线性函数的似然函数的最大值。 结果上：由于 不考虑估计固定效应，因此，REML法克服了ML法中产生有偏估值的一个来源，即克服了ML法中因估计固定效应而造成的自由度的损失。 REML和ML的相同点 都具有大样本特性，即小样本数据估值有偏。 非求导的REML—DFREML DFREML、MTDFREML、DMU 求一阶偏导数的REML—EMREML VCE、REMLF90 利用一阶和二阶偏导数的REML—AIREML DFREML、AIREMLF90、VCE、DMU REML的有关算法与相关软件 约束与假设 所有犊牛都来自同一个品种； 母亲年龄对犊牛体重无影响； 犊牛的性别相同或性别对体重无影响； 除日龄组外的其他环境条件相同。 对每一观察值建立方程式 y11=198= ? + a1 +e11 y12=204= ? + a1 +e12 y13=201= ? + a1 +e13 y21=203= ? + a2 +e21 y22=206= ? + a2 +e22 y23=210= ? + a2 +e23 y31=205= ? + a3 +e31 y32=212= ? + a3 +e32 y33=216= ? + a3 +e33 y41=225= ? + a4 +e41 y42=220= ? + a4 +e42 日龄组 1 2 3 4 观察值 残 差 个 体 1 2 3 4 日 龄 组 y11 =198= ? + a1 + e11 y12 =204= ? + a1 + e12 y13 =201= ? + a1 + e13 y21 =203= ? + a2 + e21 y22 =206= ? + a2 + e22 y23 =210= ? + a2 + e23 y31 =205= ? + a3
文档评论（0）