乌骨山羊体尺体重性状的相关及回归分析

索效军，张 年，杨前平，陶 虎，熊 琪，李晓锋，陈明新

（1.湖北省农业科学院畜牧兽医研究所，武汉 430064；2.动物胚胎工程及分子育种湖北省重点实验室，武汉 430064）

乌骨山羊体尺体重性状的相关及回归分析

索效军1，2，张 年1，杨前平1，陶 虎1，熊 琪1，李晓锋1，陈明新1

（1.湖北省农业科学院畜牧兽医研究所，武汉 430064；2.动物胚胎工程及分子育种湖北省重点实验室，武汉 430064）

为乌骨山羊品种选择和育种提供理论依据，利用SPSS软件对64只乌骨山羊成年母羊体尺体重性状进行相关分析，并采用逐步回归法，建立最优回归方程。结果表明：乌骨山羊体高（X1）、体长（X2）、胸围（X3）、管围（X4）与体重（Y）均呈极显著正相关（P＜0.01），胸围、体长及体高的直接作用和间接作用对体重影响极大（P＜0.01）；乌骨山羊体尺指标（Xi）对体重（Y）的最优回归方程为Y=-69.923 19 +0.326 27X1+0.138 46X2+0.961 72X3。

乌骨山羊；体重；体尺；相关性分析；逐步回归分析

通山乌骨山羊是在独特的自然生态环境条件下，经过长期的闭锁繁育而形成，体型与一般山羊基本相同，毛色以黑色、花色、黄色、黑白花色为主，唇、眼睑、舌、粘膜、肌肉、骨骼、内脏呈乌黑色，是集科研、药用、食用于一身的珍稀肉羊品种资源。体重是乌骨山羊最重要的遗传改良性状之一，体尺不仅反映乌骨山羊的体格大小、体躯结构和器官发育，且与肉羊的繁育、抗逆性等有一定的相关性［1］，但目前尚未有乌骨山羊体重与体尺关系的研究报道。本研究利用相关及回归分析等方法，分析了乌骨山羊体重与体尺性状间的内在关系，旨在为乌骨山羊的选育和改良提供科技支撑。

1 材料与方法

1.1 材料

乌骨山羊选自通山阿淼养殖专业合作社，供试羊均为成年母羊，健康状况良好。

1.2 测定指标与方法

测定指标为体高、体长、胸围、管围和体重，按文献提供方法测定［2］。

运用SPSS 19.0软件分析各体尺性状对体重的直接作用和间接作用，对体重进行正态性检验，利用逐步回归法建立体尺性状对体重的最优回归方程，并对方程进行显著性检验［3-5］。

2 结果与分析

2.1 体重与体尺

由表1可以看出，乌骨山羊体重和体长的变异系数较大，其中体重的变异系数最大（22.13%），说明乌骨山羊体重有较大的选择潜力。

2.2 体重与体尺性状间的相关分析

乌骨山羊体尺性状指标作为自变量，体重作为依变量，相关系数见表2，可以看出，乌骨山羊体重与体尺均呈极显著正相关（P＜0.01），相关程度最高的是胸围（r= 0.870 60），其次是管围、体长（r=0.571 38，r=0.517 27）；各体尺性状间（除体高与体长外）均呈较强正相关（P＜0.05），其中胸围与管围的相关程度最大（r=0.537 37）。

表1 体重与体尺性状的平均数、标准差及变异系数

表2 体重与体尺性状的相关系数

2.3 体重与体尺性状间相关系数的分解

根据相关系数的组成效应，可将相关系数（rij）分解为体尺对体重的直接作用（Pi）和某一体尺性状（Xi）通过其他体尺性状（Xj）对体重（Y）的间接作用（∑rijPj）两部分［6-7］，即rij=Pi+∑rijPj。从表3可知，胸围对体重的直接影响作用最大，为0.754 86；而体长对体重的间接作用最大，为0.355 81。

表3 简单相关系数的分解

2.4 体重正态性检验

3.1.2 暂时大比分领先。调查中发现，有4个队员认为比赛的胜利已经稳握在自己手中，对现在的领先已经满足，因此比赛紧张的心理状态放松了下来，慢慢不敢进攻，害怕失误，对教练在赛前布置战术和打法抛于脑后，没有继续发挥比赛的后半程优势，失误越来越多。而另外6名运动员发挥比较稳定，对这种短暂领先没有放在心上，抓住对方进攻和防守的弱点，发挥自己的战术，继续拉大领先优势，最后取得比赛胜利。（见表2）

对体重进行正态性检验，输出结果见表4，经检验在Kolmogorov-Smirnov Test和Shapiro-Wilk Test状态下的P值分别为0.200 00和0.223 92，均大于0.05，因此，依变量（体重）服从正态分布，即体重是正态变量，可以进行回归分析。

表4 正态性检验输出结果

2.5 最优回归方程的建立

2.5.1 体尺指标对体重回归方程 依据体尺性状对体重贡献率大小和标准回归系数的显著性，通过逐步回归分析，剔除对体重影响不显著的体尺性状，建立最优的多元线性回归方程。本研究中随着体尺性状被逐步引入，回归方程的复相关系数和决定系数在逐渐增大，说明引入的体尺性状对体重的贡献率在增加。详见表5。

表5 体尺指标对体重回归模型

2.5.2 方差分析 由表6可知，自变量X3引入回归方程后，F=194.147 28；自变量X3、X2引入回归方程后，F= 107.979 66；自变量 X3、X2、X1引入回归方程后，F= 79.817 04。回归方程1、2、3方差分析均达到极显著水平（P<0.01），而X4的偏回归系数无统计学意义，即对回归方程的影响不大，故未被引入回归方程。

2.5.3 偏回归系数及其t检验 各步引入对回归方程影响最大的变量后，有关的偏回归系数及t检验结果见表7，回归方程3中自变量X3、X2、X1所对应的P值均小于0.05，而没有引入方程的 X4在回归方程中其P值均大于0.05，无显著统计学意义，故为不重要变量。

表6 乌骨山羊多元回归方差分析

表7 体重与体尺的多元回归系数的参数估计值检验

3 讨论

3.1 体重与体尺的变异系数

体重与体尺是全面、综合分析羊生长发育的重要指标，对肉羊育种工作及生产实践有重要作用。变异系数是标志变动度相对程度的指标，又称离散系数。在生产实际中，体重可能会由于各种原因变化较大，因此可利用体尺指标进行肉羊的品种选育。本研究中体重（Y）、体长（X2）二者变异系数均大于10%，这表明乌骨山羊体重、体长变异幅度较大；而体高（X1）、胸围（X3）、管围（X4）变异系数小于10%，说明乌骨山羊体高、胸围、管围的变动程度相对较小。

3.2 体重与体尺性状间的相关关系

相关系数可反映体重与体尺两个随机变量间的表型相关程度，由于没有排除其他变量的影响，因此并不能反映各体尺性状与体重的直接关系。本研究中胸围与体重的相关程度最高（r=0.870 60），其次是管围（r= 0.571 38）、体长（r=0.517 27）；除体高与体长外，各体尺性状间亦呈较强正相关（P＜0.05），其中胸围与管围的相关系数最大（r=0.537 37）。体重与胸围之间的相关是由于家畜体重与骨骼、肌肉、内脏的生长发育密切相关，同时体重与体长、管围的相关也是与家畜骨骼、肌肉的生长发育密不可分。前人已将体尺性状与体重性状的相关关系成功应用于羊、鸡等品种的改良和新品系选育中［8-10］。因此，今后在乌骨山羊的品种选育方面应以胸围为主并兼顾体长、管围。

3.3 体重与体尺指标间相关及其相关系数的分解

在分析一个目标性状（依变量）和若干个不同性状（自变量）之间相关关系时，常采用多元回归分析方法［11］。由相关系数剖分可知，乌骨山羊体重（Y）和胸围（X3）的相关主要是由两者直接作用而形成的，此结果与买买提明巴拉提等［12］、徐铁山等［13］研究结果一致，与王同军等［14］的研究结果不一致，可能是品种差异所致；体高（X1）、体长（X2）对体重的影响直接效应较小，主要通过间接作用产生，这表明肉羊是一个交互作用的有机整体，其决定体重的因素较多，影响途径也有所不同，不仅有单一体尺性状的影响，也存在多个体尺性状的同时作用。因此，在实际生产中在加强对胸围选择的同时，也应注意关注间接作用的选择，进行综合考虑，从而有一定的目的性，避免盲目选择。

3.4 最优回归方程的建立

通过对体尺指标与体重的逐步回归分析，得到体尺对体重的最优回归方程：Y=-69.923 19+0.326 27X1+ 0.138 46X2+0.961 72X3，入选回归方程的体尺性状胸围（X3）、体长（X2）和体高（X1）回归关系都达到显著水平（P<0.05）。本研究中胸围、体长是影响体重的主要因素，这一结果与高志英等［15］的研究结果相吻合。由线性回归分析得出胸围、体长和体高每增加1 cm，体重将分别增加0.961 72、0.138 46和0.326 27 kg，说明这3个体尺指标对乌骨山羊体重增长影响比较大。本回归方程拟合度的决定系数（R2）值为0.799 63，校正R2为0.789 62，剩余因子e=0.600 49，该值较大，说明还有对体重有影响的因素尚未考虑进去，有待进一步加大样本进行研究，以完善乌骨山羊体尺、体重的最优回归方程。

4 结论

胸围、体长和体高是影响乌骨山羊母羊体重的主要因素，乌骨山羊成年母羊体尺与体重的最优回归方程为Y=-69.923 19+0.326 27X1+0.138 46X2+0.961 72X3。

［1］ Semakula J，Mutetikka D，Kugonza R D，et a1.Variability in body morphometric measurements and their application in predicting live body weightof Mubende and Small East African goatbreeds in Uganda［J］.Middle-East Journal of Scientific Research，2010，5（2）：98-105．

［2］ 陈伟生，徐桂芳.中国家畜地方品种资源图谱（下）［M］.北京：中国农业出版社，2004.

［3］ 苏理云，陈彩霞，高红霞.SPSS 19统计分析基础与案例应用教程［M］.北京：北京希望电子出版社，2012.

［4］ 明道绪.生物统计附试验设计［M］.北京：中国农业出版社，2002.

［5］ 索效军，张年，熊琪，等.麻城黑山羊母羊体质量与体尺的回归分析［J］.东北农业大学学报，2013，44（12）：63-67.

［6］ 杨佩燕，杨姣，刘志学.崇明白山羊和徐淮山羊体重与体尺的相关和通径分析［J］.上海农业学报，2010，26（1）：24-27.

［7］ 蔡惠芬，陈志，罗卫星，等.黔北麻羊体尺与体重的相关分析［J］.畜牧与兽医，2012，44（7）：53-55.

［8］ 韩廷义.青海高原半细毛羊选育技术的研究［J］.中国畜牧兽医，2012，39（12）：230-232.

［9］ 张力，江宵兵，王长康.福建闽清毛脚鸡体尺性状与屠体性状的典型相关分析［J］.西北农林科技大学学报，2010，38（2）：8-12.

［10］马猛，王克华，曲亮，等.13周龄母鸡体尺性状主基因+多基因混合遗传模型分析［J］.东北农业大学学报，2015，46（9）：75-82.

［11］Yücel C.Correlation and path coefficientanalyses ofseed yield components in the narbon bean（Vicia narbonensis L.）［J］.Turkish Journalof Agriculture and Forestry，2004，28（5）：371-376.

［12］买买提明·巴拉提，怕尔哈提，买合木提江，等.周岁麦盖提羊体重与体尺指标相关性的研究［J］.四川畜牧兽医，2001，28（7）：20-21.

［13］徐铁山，王东劲，刘小林，等.海南黑山羊体尺与体重的通径分析及最优回归模型的建立［J］.家畜生态学报，2005，26（1）：50-53.

［14］王同军.建昌黑山羊体尺体重相关性研究［J］.四川畜牧兽医，2003，30（9）：22-26.

［15］高志英，乌云毕力克，艾尔肯江，等.德国美利奴成年母羊体尺与体重相关回归分析［J］.草食家畜，2012，157（4）：31-36.

信息专栏

《中国草食动物科学》面向全国读者开展免费阅读

为进一步提高本刊杂志论文的扩散力，进而提高杂志的影响因子和被引频次，更为广大科研人员创造良好的文献查阅平台，《中国草食动物科学》于2013年在教育部科技发展中心的“中国科技论文在线”网站上面向全国读者开展免费阅读，这将大大提高论文的引用率（3倍）。同时，本刊于2013年开始应用DOI，意为“数字物品身份证符”。由于DOI的全球唯一性，在科技文献中如果给出被引参考文献的DOI，通过DOI可以直接定位到被引参考文献上，这样就提高了被引文献的表现度，以及论文的扩散力和影响力。

Correlation and Regression Analysis of Body Weight and Body Sizes in Black-bone Goats

Suo Xiaojun1，2，Zhang Nian1，Chen Mingxin1，etal
（1.Institute ofAnimalScience，HubeiAcademy ofAgriculturalSciences，Wuhan 430064，China；2.HubeiKey Laboratory ofAnimalEmbryo and Molecular Breeding，Wuhan 430064，China）

To provide a theoreticalbasis for species selection and breeding ofBlack-bone goat，the data ofbody weights and body sizes from 64 adult Black-bone does were analyzed by SPSS software with regression analysis and correlation analysis，stepwise regression analysis method，the optimalregression equation was established.The results showed that the body height（X1），body length（X2），chestcircumference（X3），cannon circumference（X4）showed very significantcorrelations with body weight（Y）atthe 0.01 level．The body weightwasextremely affected directly and indirectly by the chestcircumference，body length，body height.The bestregression equation ofbody weightwith body sizeswasobtained as：Y=-69.923 19+0.326 27X1+0.138 46X2+0.961 72X3.

Black-bone goat；body weight；body size；correlation analysis；stepwise analysis

S827.2

A

2095-3887（2017）05-0010-04

10.3969/j.issn.2095-3887.2017.05.003

2017-06-01

湖北省农业科技创新中心资助项目（2016-620-004-002）

索效军（1977-），男，副研究员。

陈明新，男，研究员，主要从事牛羊育种研究工作。