一些育种的基本概念，需要看教科书理解一下。这次我学习教科书，加深对相关概念的理解，最后用实际数据演示一下如何计算遗传相关及显著性。

参考书： Raphael A.Mrode编著，于向春 张豪 主译《线性模型在动物育种值预测中的应用》 第三版 陈国宏, 张勤. 动物遗传原理与育种方法[M]. 中国农业出版社, 2009.

在育种过程中，经常对多个性状进行选择，这些性状可能有遗传相关。

多性状模型可以充分利用性状的表型相关和遗传相关等信息，针对两性状或多性状对个体进行评估。

多性状模型最佳线性无偏预测（MBLUP）的优点之一是增加了评估的精确性。精确度的增加取决于遗传相关和残差相关大小。这些相关的差异越大，精确性增加越多。

低遗传力和高遗传力同时进行多性状分析时，低遗传力性状收益更多。而且多性状模型分析时，因为考虑了性状间的残差协方差，所以增加了评定的准确性。

另外，多性状模型可以避免淘汰偏差（culling bias），比如断奶体重和周岁体重这两个性状，很多个体都会根据断奶体重进行淘汰，因此很多个体都没有周岁体重。如果只针对周岁体重的个体进行遗传评估（选择后的个体），结果会产生偏差，因为它没有利用断奶体重的信息。所以，对断奶体重和周岁体重进行双性状模型评估，可以消除这种偏差。

遗传相关的分子原因，可以分为两类：

遗传相关可用于确定间接选择的依据和预测间接选择反应大小。间接选择是指当一个性状（如X）不能直接选择或者选择效果很差时，借助与之相关的另一个性状（如Y）的选择来达到对性状X的选择目的

间接选择在育种实践中具有很重要的意义，比如有些性状只有屠宰后才能度量，有些性状只能在个体的生命晚期才能度量。有些性状受性别限制无法度量，还有些性状直接选择效果不理想，在这些情况下都可以考虑采用间接选择

遗传相关可用于比较不同环境条件下的选择效果。实际上，不但不同性状可以来估计遗传相关，而且同一性状在不同环境下的表现也可以作为不同的性状来估计遗传相关。

不同环境下的遗传相关，为解决育种工作中的一个重要实际问题提供了理论依据，即在条件优良的种畜场选育的优良品种，推广到条件较差的其它条件生产厂是否能保持其优良特性。

实际上，遗传相关可以推断同一性状在不同环境下的反应是否一致。

一般而言，只要涉及两个性状以上的选择问题，无不需要用到遗传相关这一参数，用于指定相关形状的选择指数，这也是遗传相关主要的用途之一。

数据来源：

数据来源我编写的learnasreml包，原数据是一篇论文中的示例数据，我将其放到包中，方便调用。

构建模型：

这里使用个体动物模型，双性状模型：BWT和TARSUS两个性状。固定因子为SEX，随机因子为加性效应。

查看模型迭代情况和方差组分：

可以看到，模型收敛，根据方差组分，计算遗传相关和表型相关：

遗传相关和表型相关的显著性

这里使用似然比检验（LRT），另外构建两个模型：

moda：

结果：可以看出遗传相关和表型相关都达到了极显著。

注意，原始的blup值是这个样子的：

我这里用tidyverse对其进行了清洗，将ID和性状提取出来，并且重新进行了整理，结果为三列：ID，trait1-blup, trait2-blup，然后将其输出到csv格式。