白萨福克羊体重与体尺指标的相关及回归分析

王伟峰

（甘肃省会宁县世界银行贷款畜牧综合发展肉羊项目办公室，会宁 730799）

白萨福克羊体重与体尺指标的相关及回归分析

王伟峰

（甘肃省会宁县世界银行贷款畜牧综合发展肉羊项目办公室，会宁 730799）

采用SPSS 17.0对白萨福克羊的体重与体尺指标进行相关性分析、通径分析、逐步回归分析，并对其相关系数进行分解，得出了体重与体尺指标之间的最优回归方程。结果表明：体高（x1）和管围（x4）均对白萨福克羊体重影响不显著（P＞0.05），其相关系数分别为0.382和0.525；胸围（x3）和体斜长（x2）对白萨福克羊的体重影响极显著（P＜0.01），其相关系数分别为0.822和0.519；体高（x1）、体斜长（x2）和管围（x4）主要是通过胸围（x3）对体重产生间接作用，胸围（x3）主要是通过体斜长（x2）和管围（x4）对体重产生间接作用；白萨福克羊体重与体尺指标之间的最优回归方程为y=-105.123+1.217x3+ 0.619x2，标准化后的方程为y=0.738x3+0.186x2（R2=0.704，P＜0.01）。

白萨福克羊；体重；相关性分析；通径分析；回归模型

白萨福克羊原产于英国英格兰东南部的萨福克、诺福克、剑桥和艾塞克斯等地。该品种羊是以南丘羊为父本，当地体型较大、瘦肉率高的旧型黑头有角诺福克羊为母本进行杂交培育的，于1958年育成［1］。白萨福克羊是目前世界上体格、体重最大的肉用羊品种之一，早熟，生长发育快，产肉性能好，瘦肉率高，是生产大胴体和优质羔羊肉的理想品种［1］。国内外学者针对有关白萨福克羊生产性能、生长发育及饲养管理等方面的研究较多［2-3］，但对其体重与体尺指标相关性研究尚缺乏报道。本研究通过对饲养在甘肃河西走廊地区的永昌县肉用种羊繁育技术推广站的白萨福克羊体重与体尺指标进行相关性分析、通径分析、逐步回归分析以及对其表型相关系数分解，探究其体尺指标与体重的相关关系，并建立体重与体尺指标之间的最优回归方程，为将来肉羊种质利用提供参考。

1 材料与方法

1.1 动物

113只白萨福克母羊均来自于甘肃省永昌县肉用种羊繁育技术推广站，年龄24月龄以上，体质健康，发育良好。

1.2 测定项目及方法

对每只绵羊均进行体重（y）和体尺指标［体高（x1）、体斜长（x2）、胸围（x3）、管围（x4）］测定，体尺指标的测量参考张沅［4］测定方法，体重用电子秤测量，精度0.01 kg，体尺指标用测杖、卷尺测量。

1.3 数据处理

采用Excel2010对所测数据初步整理后，用SPSS 17.0统计软件进行相关性分析、通径分析、逐步回归分析，并对相关系数进行分解，建立体重与体尺指标之间的最优回归方程。

2 结果与分析

2.1 体重和体尺指标的表型分析

由表1可知，体重的变异系数最大，为19.8%；其次分别为管围、胸围、体斜长，其变异系数分别为9.7%、8.7%、6.2%；体高的变异系数最小，为4.0%。

2.2 体重和体尺指标的表型相关性分析

体重和体尺指标之间的相关系数如表2所示：白萨福克羊的体重与体高、体斜长、胸围、管围间相关性均极显著（P＜0.01）。体重与胸围表型相关系数最大，为0.822，与管围、体斜长、体高的表型系数分别为0.525、0.519、0.382。每两个体尺指标之间相关性也均极显著（P＜0.01）。

表1 白萨福克羊体重和体尺指标的表型统计量

表2 白萨福克羊体重和体尺指标的表型相关系数

2.3 体重和体尺指标之间的通径分析

为了进一步明确体重和各体尺指标之间的关系，用体尺指标与体重的表型相关系数建立以通径系数为变量的多元一次方程组，方程组如下：

p1+0.495p2+0.331p3+0.306p4=0.382

0.495 p1+p2+0.451p3+0.338p4=0.519

0.331 p1+0.451p2+p3+0.544p4=0.822

0.306 p1+0.338p2+0.544p3+p4=0.525

用矩阵法在Excel 2010中对多元一次方程组进行求解，得出白萨福克母羊各体尺指标对体重的通径系数分别为：p1=0.052；p2=0.154；p3=0.691；p4=0.081。

2.4 体尺指标对体重相关系数的分解

由表3可知，白萨福克羊的胸围对体重的直接作用最大，为0.691；体高对体重的直接作用最小，为0.052；管围对体重的间接作用最大，为0.444；胸围对体重的间接作用最小，为0.131。在各体尺指标中对体重的直接作用由大到小依次为胸围＞体斜长＞管围＞体高，对体重的间接作用由大到小依次为管围＞体斜长＞体高＞胸围。由此可以看出，体高对体重有一定的间接作用，直接作用较弱；体斜长对体重既有一定的直接作用，又有一定的间接作用；胸围对体重既有明显的直接作用，又有一定的间接作用；管围对体重有明显的间接作用，直接作用较弱。体高、体斜长和管围主要是通过胸围对体重产生间接作用，胸围主要是通过体斜长和管围对体重产生间接作用。

2.5 体尺指标对体重的逐步回归分析

2.5.1 白萨福克羊体重和各体尺指标整体回归分析模型 从表4中可以看出，各个自变量的显著性由大到小依次为胸围＞体斜长＞管围（不显著）＞体高（不显著）。因此剔除自变量体高、管围后，对因变量体重与体斜长、胸围进行多元逐步回归分析。

2.5.2 白萨福克羊体重和各体尺指标逐步回归分析模型汇总 由表5所得，体重和体尺指标之间的最优回归方程为y=-105.123+1.217x3+0.619x2，标准化后的回归方程为y=0.738x3+0.186x2（R2=0.704，P＜0.01），方差分析的F值为130.804（见表6），达到极显著水平（P＜0.01）。

表3 白萨福克羊体尺指标对体重相关系数的分解

表4 白萨福克羊体重和各体尺指标总体回归分析后多元回归系数参数估计值

表5 白萨福克羊体重和胸围、体斜长指标逐步回归分析后多元回归系数参数估计值

表6 白萨福克羊体重和胸围、体斜长指标逐步回归分析后多元回归方差分析

3 讨论

本研究通过对体重和体尺指标的相关性分析，可以看出体重与体高、体斜长、胸围、管围的相关关系均极显著，其表型相关系数分别为0.382、0.519、0.822和0.525。说明体高、体斜长、胸围和管围对体重的综合作用较大，其表现为直接作用或辅助其他性状对体重的间接作用。在体高、体斜长、胸围、管围4个自变量中，每两个自变量均表现出极显著相关性。这与白俊艳等［5］对大尾寒羊、高志英等［6］对多浪羊、高晔等［7］对陕北白绒山羊、顾亚玲等［8］对白萨福克羊、李少斌等［9］对甘肃高山细毛羊杂交一代、喇永富等［10］对湖羊等体重与其体尺指标分析结果有类似的地方。说明白萨福克羊的体重和各体尺指标是相互联系的，选种时对各个性状的选择都需要综合考虑，以达到最优选育效果。

通径系数表示所研究的自变量对依变量的直接作用大小，它能真实地反映自变量对依变量的直接作用程度，它是标准化的偏回归系数［11-12］，它是一个在生物实验数据分析时非常重要的统计量。通过对体重和各体尺指标之间的通径分析以及其相关系数的分解可以看出：体重与体高、体斜长、胸围、管围相关关系均极显著（其表型相关系数分别为0.382、0.519、0.822、0.525）。从表3中相关系数的分解中可以看出，体高对体重有一定的间接作用，直接作用较弱；体斜长对体重既有一定的直接作用，又有一定的间接作用；胸围对体重既有明显的直接作用，又有一定的间接作用；管围对体重有明显的间接作用，直接作用较弱。在各体尺指标中对体重的直接作用由大到小依次为胸围＞体斜长＞管围＞体高，对体重的间接作用由大到小依次为管围＞体斜长＞体高＞胸围。通过上述分析说明，生物体是一个有机作用的综合体，其体重的决定因素不仅较多，而且影响途径也不同。因此，在实际选种过程中不仅要加强对体重直接作用强的因素体斜长、胸围的选择，也应关注间接作用的选择，进行综合考虑，在选择时可参考本研究建立的最优回归模型，根据生物统计学原理进行选育，从而有一定的目的性，避免盲目选择。

通过对白萨福克羊体重与其体尺指标之间的逐步回归分析，说明体高、管围对白萨福克羊的体重影响不明显。回归方程中R2=0.704，说明该回归方程的回归平方和占总变异的70.4%，回归方程显著存在，所选体尺指标对体重有较大的决定作用，方程拟合度较高，具有一定的参考价值。剩余因子e=1-R2=0.296，说明对白萨福克羊体重产生影响的自变量不止胸围、体斜长，还有一些影响因素没有考虑到，因此对影响白萨福克羊体重的因素需要进行全面分析。

4 结论

本研究条件下，白萨福克羊的体高对体重有一定的间接作用，直接作用极弱；体斜长对体重既有一定的直接作用，又有一定间接作用；胸围对体重具有明显的直接作用，又有一定的间接作用；管围对体重有明显的间接作用，直接作用较弱。体高、体斜长和管围主要是通过胸围对体重产生间接作用，胸围主要是通过体斜长和管围对体重产生间接作用。白萨福克羊的体重和体尺指标之间的最优回归方程为y=-105.123+1.217x3+0.619x2，标准化后的回归方程为y=0.738x3+0.186x2（R2=0.704，P<0.01）。

［1］ 赵有璋.中国养羊学［M］.北京：中国农业出版社，2013.

［2］ 杨宇泽，王大愚，赵有璋.白萨福克种羊在河西走廊地区行为学研究初探［J］.家畜生态报，2009(2)：78-80.

［3］ 赵有璋，王大愚，杨宇泽.白萨福克羊和特克塞尔羊在河西走廊地区适应性比较研究［J］.甘肃农业大学学报，2008（6）：38-44.

［4］ 张沅.家畜育种学［M］.北京：中国农业出版社，2001：138.

［5］ 白俊艳，庞有志，王永伟.大尾寒羊体重与体尺的回归分析［J］.安徽农业科学，2007，35（15）：4537-4538.

［6］ 高志英，於建国，阿曼姑阿巴斯，等.多浪羊体重及其体尺指标相关性研究［J］.草食家畜，2009，12（4）：29-31.

［7］ 高晔，闫海龙，冯平，等.陕北白绒山羊体尺与体重相关回归分析［J］.中国草食动物，2010，30（5）：30-32.

［8］ 顾亚玲，马吉锋.萨福克羔羊断奶体重体尺指标及其相关性的研究［J］.内蒙古畜牧科学，2003，24（5）：28-29.

［9］ 李少斌，王继卿，成述儒，等.甘肃高山细毛羊杂交一代生长模型及杂种优势分析［J］.草业学报，2011，20（5）：177-182.

［10］喇永富，席锐，李发弟，等.西北寒旱地区湖羊及其杂交后代生长发育模型分析［J］.草业学报，2016，25（9）：182-188.

［11］王志明，王婕姝，王毅，等.甘肃高山细毛羊体重与体尺指标的相关性研究［J］.浙江农业学报，2016，28（1）：28-32.

［12］宋小园，朱仲元，刘艳伟，等.通径分析在SPSS逐步线性回归中的实现［J］.干旱区研究，2016，33（1）：108-113.

Correlation and Regression Analysis ofBody Weightand Body Sizes in White Suffolk Sheep

Wang Weifeng
（Sheep ProjectOffice ofWorld Bank Loan Comprehensive Livestock DevelopmentProject，Huining 730799，Gansu，China）

The phenotypic correlations，path coefficients，directand indirectcorrelations between the body weightand body sizes in White Suffolk sheep were analyzed by SPSS 17.0 software.Then the optimum regression equation between the body weightand body sizes was established.The results indicated that the body height（x1），cannon-bone circumference（x4）did not show significantcorrelation with the body weightatthe 0.05 leveland the phenotypic correlation were 0.382 and 0.525，respectively.The chest circumference（x3）and body length（x2）showed significant correlations with the body weight at the 0.01 level and the phenotypic correlation coefficientwere 0.822 and 0.519，respectively.The body height，body length and cannon bone circumference had indirectinfluence by chestcircumference mainly.The chestcircumference had indirectinfluence by body length and cannon bone circumference.The optimum regression equation between the body weight and body size in White Suffolk sheep was：y=-105.123+1.217x3+0.619x2，and the normalized formofitwas：y=0.738x3+0.186x2（R2=0.704，P＜0.01).

White Suffolk sheep；body weight；correlation analysis；path analysis；regression model

S826.2

A

2095-3887（2017）05-0014-04

10.3969/j.issn.2095-3887.2017.05.004

2017-06-23

王伟峰（1981-），女，畜牧师，主要从事畜牧技术推广工作。