对于可靠性我写过很多次，但是随着使用模型的复杂，也会逐渐增加难度。
这里参考一篇老文献对多性状多胎次的随机回归测定日模型的可靠性进行讲解（可能还没理解很清楚）。
Interbull对于遗传估计（genetic evaluation）的权重因子给出的参考程序：
第1步骤： 计算可靠性根据动物自己的行为记录。
第2步骤：使用后代和伴侣的可靠性计算每头公牛的可靠性。
在MTMTRRM中一般很难转换第1步骤，因为其模型中一般含有除了遗传和永久环境外的其他随机因子，如herd-test day与Resiudal， 这两个的矩阵大小与性状的矩阵不直接相对。
所以该作者主要对第1步骤进行研究和提出解决方法。
众所周知，可靠性与PEV直接相关，公式为：
可靠性= 1 PEV/var(G)
但是大的矩阵，不能直接对系数求逆，得到PEV。所以作者提出近似PEV（approximation PEV）。
首先LRH系数分为：
image.png
其中：
image.png
分两步骤计算某头牛的近似PEV：
A. 将HTD消解（obsorbed）到Ci中：
image.png
B. 将PE的方程消解到动物遗传部分：
image.png
发送给Interbull的育种值为：
image.png
近似PEV就为：
image.png
image.png
参考文献：Stranden, I., M. Lidauer, E. A. Mäntysaari, and J. Pösö. 2000. Calculation of Interbull weighting factors for the Finnish test day model. Interbull Bull., 26:78–79.
image.png
网友评论