努比亚山羊体尺体重的聚类和主成分分析

朱俊红，刘兴能，和晓明，任洪辉，陆颖，邓卫东，席冬梅

（ 云南农业大学动物科学技术学院，云南 昆明 650201 ）

努比亚山羊原产自非洲东北部等地区，属乳肉皮兼用型，毛色多为黑色或棕色，1995年引入我国饲养培育[1]。努比亚山羊具有生产性能优良、遗传性能稳定、适应性强、肉质细嫩、膻味小等优点[2-4]。国内学者针对努比亚山羊进行了大量研究。占今舜等[5]对努比亚山羊的肉品质、血清指标等进行了研究；王坤等[6]和周波[7]对努比亚山羊的养分消化和瘤胃微生物等进行了分析研究。国外学者对努比亚山羊遗传参数估计、遗传多样性评价和生长性能、瘤胃发酵、血液化学等[8-10]也进行了大量研究。家畜体尺数据直接反映了家畜体躯的结构、发育等情况，与家畜的生理功能、抗病力以及对外界环境条件的适应能力等密切相关。目前对努比亚山羊体重和体尺指标的分析研究较少，本文对云南省红龙养殖场的努比亚山羊体重和体尺指标进行分析研究，为日后更好地开发和保护努比亚山羊遗传资源提供参考。

1 材料与方法

1.1 试验材料

随机选取云南省昭通市彝良县红龙养殖场3月龄大小的羔羊32只，公、母各半，测量其体重、体高、体斜长、胸围、胸宽、胸深、管围、最大额宽，测量方法参考家畜育种学和羊生产学测定方法[11-12]。

1.2 数据统计与分析

运用Excel软件整理努比亚山羊的体重和体尺指标数据，SPSS 25.0软件进行主成分分析与聚类分析。

聚类分析是指物理或抽象对象的集合，分组为由类似对象组成的多个类的分析过程[13]，具体操作如下：计算性状指标的相关矩阵、相关矩阵的特征根和特征向量，选取累计贡献率＞85%的特征根，进行主成分分析。聚集8个指标性质相同或相近的对象，并按照这些对象的定性或定量特征数值将其分组归类。

2 结果与分析

2.1 测定指标的平均数与变异性指标（见表1）

表1 测定指标的平均值与变异系数Tab.1 Mean and coefficient of variation of determination index

由表1可知，体重的变异系数为33.34%，变异幅度较大；其他性状指标的变异系数均小于7%，变异幅度小。结果表明，努比亚山羊体尺指标变异比较稳定，但体重指标变异幅度稍大，选育潜力也相对较大。

2.2 主成分分析

2.2.1 体重和体尺指标间的相关性（见表2、表3）

表2 体重和体尺指标的相关系数矩阵Tab.2 Correlation coefficient matrix of body weight and body conformation indexes

表3 相关系数矩阵的特征根及累计贡献率Tab.3 Characteristic roots and cumulative contribution rate of correlation coefficient matrix

由表2可知，努比亚山羊各指标之间均具有显著相关性，其中仅最大额宽和胸深显著性相关（P＜0.05），其余各指标之间均极显著相关（P＜0.01），表明各指标间关联性很强，可利用显著性较强的指标进行努比亚山羊的选育。

由表3可知，选取前3个主成分，使累计贡献率大于85%（达到87.798%），损失信息仅有12.202%，相对较小。其中第1主成分贡献率为71.091%，反映的综合信息量最大；第2、第3主成分贡献率分别为10.501%、6.206%。

2.2.2 主成分的特征根和贡献率（见表4）

表4 主成分的特征根和贡献率Tab.4 Characteristic roots and contribution rates of principal components

由表4可知，贡献率较大的3个主成分代表分别为Y1、Y2、Y3。根据各主成分的特征向量，列出各主成分的关系表达式如下：

第1主成分（Y1）特征根为5.687，贡献率为71.091%，其中体重系数（X1）最大，其次是体高（X2）、体斜长（X3）、胸围（X4），可将其命名为体重因子，反映体重和体型外貌特征，主要决定努比亚山羊体形的长短高矮。第2主成分（Y2）特征根为0.840，贡献率为10.501%，其中胸宽（X5）、胸深（X6）、管围（X7）系数较大，可将其命名为躯体因子，反映躯体宽度特征，主要决定努比亚山羊的宽度大小。第3主成分（Y3）特征根为0.496，贡献率为6.206%，最大额宽（X8）系数最大，其余较小，可将其命名为额宽因子，反映努比亚山羊的额宽特征。

2.2.3 努比亚山羊8个性状的碎石图（见图1）

图1 努比亚山羊8个性状的碎石图Fig.1 Gravage map of eight traits of Nubian goat

由图1可知，努比亚山羊的8个性状中，因子1、2、3之间的连线最为陡峭，表明这3个因子均为主要因子，占累计贡献率约90%。较陡的直线断点也表明，相应因子之间具有较大的特征值差异，进一步证实了所选取的3个主成分的准确性。

2.3 努比亚山羊8个性状的聚类分析（见图2）

图2 努比亚山羊8个性状的聚类分析Fig.2 Cluster analysis of eight traits of Nubian goat

由图2可知，相似程度较高的指标聚为一类，随着相似程度降低，聚类的类数呈增加趋势。类间距为15时，8个性状指标聚类为2个亚类：第1类包括体重（X1）、胸深（X6）、胸宽（X5）、管围（X7）和最大额宽（X8），第2类包括体高（X2）、体斜长（X3）和胸围（X4）。

3 讨论

体重和体尺指标是家畜选育中重要的经济性状，亦是遗传效应与环境效应共同作用的结果。本试验对努比亚山羊的体重和体尺指标进行研究，结果显示，各指标之间呈显著或极显著相关。这与努尔比亚·吾布力等[14]对中国美利奴羊以及侯洪梅[15]对加什科羊的主要性状的研究结果一致。杨存明等[16]通过对和田公羊进行相关性分析揭示了影响其体重的主要体尺指标是胸围、体高、胸深，与本研究结果一致。谭晓山等[17]对贵州黑山羊的体尺指标和体重进行了相关性研究，发现体长、胸围、腰角宽可作为选育的重要参考指标。本研究结果表明，努比亚山羊的体重、体尺指标间具有显著相关性，可作为重要选育指标。

主成分分析是一种统计方法，可用于描述家畜形态结构并说明表型性状对整体形态结构的影响[18]。主成分分析能够从原始数据中提取出最具代表性的信息，并将其表示为一组相互独立的主成分。这些主成分能够综合反映原始数据中的多个指标，并使用较少的变量解释数据的变异情况[19]。本研究主成分分析将8个指标区分为3个独立的主成分，其中第1主成分（Y1）反映体重和体型外貌特征，第2主成分（Y2）反映躯体特征，第3主成分（Y3）反映额宽特征。这与叶昌辉等[20]将雷州山羊的8个性状指标区分为相对独立的3个主成分研究和徐晓莉等[21]对绵羊体尺等指标的分析研究结果相似。有研究对不同杂交组合羔羊6月龄的体重和体尺指标进行聚类分析，将体尺性状分为2类，体高、体长、胸围为一类，管围为一类[22]。本试验也是将8个性状指标聚类为2个亚类，且聚类的体尺指标也非常一致。李明明等[23]对藏羊羔羊的研究结果也与本试验结果相似，进一步说明本研究的结果在不同品种山羊之间具有一定的适用性。目前针对牛体尺指标已有大量研究分析[24-26]，但针对羊体尺指标的研究较少，因此本研究可为养殖场日后羊的选育规划提供参考。

4 结论

本文对努比亚山羊的体重、体高、体斜长、胸围、胸宽、胸深、管围、最大额宽等8个性状指标进行分析，根据贡献率选取了3个独立的主成分，其中第1主成分贡献率为71.091%，是努比亚山羊体尺等指标变异的主要来源，主要反映体重和体型外貌特征；第2主成分贡献率为10.501%，主要反映躯体特征；第3主成分贡献率为6.206%，主要反映额宽特征。聚类分析将8个指标聚类为2个亚类，第1类包括体重、胸深、胸宽、管围和最大额宽，第2类包括体高、体斜长和胸围。本试验结果可为今后对努比亚山羊选育和生产提供一定的参考。