澳洲白羊体质量与体尺指标的相关及回归分析

施海娜， 刘雨田， 李世恩， 张金霞， 高 钰， 梁万鹏， 朱正生， 陈兆勇， 李 凯， 徐振飞*

(1.庆阳市农业科学研究院，甘肃 庆阳 745000；2.甘肃中盛农牧发展有限公司，甘肃 庆阳 745000)

“澳洲白”是澳大利亚第1个利用现代基因测定手段培育的品种。该品种集成了白杜泊绵羊、万瑞绵羊、无解道赛特羊和特克赛尔羊等品种基因，通过对多个品种羊特定肌肉生长基因标记和抗寄生虫基因标记的选择(MyoMAX，LoinMAX，WormSTAR)，培育而成的专门用于与杜泊绵羊配套的，属于粗毛型的中、大型肉羊品种，于2009年10月在澳大利亚注册。其特点是体型大、生长快、成熟早、全年发情，有很好的自动换毛能力。在放牧条件下5～6月龄可达到23 kg胴体，舍饲条件下，该品种6月龄胴体重可达26 kg，且脂肪覆盖均匀，板皮质量具佳。澳洲白在我国已多次引入，截至2017年甘肃省已引入发展数量达到2 000只以上[1]。此品种使养殖者能够在各种养殖条件下用作三元配套的终端父本，可以产出在生长速率、个体重量、出肉率和出栏周期短等方面理想的商品羔羊。

家畜体尺数据直接反映了家畜的体格大小和体躯结构、发育等情况，与家畜的生理机能、生产性能、抗病力以及对外界环境条件的适应能力等密切相关。国内外学者针对有关澳洲白羊种羊饲养管理技术和杂交改良地方品种试验等方面的研究较多[2-11]，但对其体质量与体尺指标相关性及回归分析研究尚缺乏报道。本文通过对以湖羊为受体采用胚胎移植技术扩繁的澳洲白羊体质量与其体尺指标进行测定，初步掌握澳洲白羊的生长发育情况，并进行了相关性分析、通径分析、逐步回归分析以及对其表型相关系数分解，探究其体尺指标与体质量的关系，建立体质量与其体尺指标之间的最优回归方程，以期能为将来的肉羊新品种选育工作提供参考数据。

1 材料与方法

1.1 试验动物

试验动物澳洲白羊为2017年9月底在甘肃中盛华美羊业西峰羊场(以湖羊为受体，从天津奥群公司购买的纯种澳洲白A级冷冻胚胎)采用胚胎移植技术所产，年龄为6月龄左右，健康状况良好。

1.2 测定项目

随机抽取32只(公母羊各16只)健康羊只对其进行体质量(y)、体尺(xi)测定。体质量称量用电子秤，精度为0.01 kg。体尺指标测量参考家畜育种学[12]，测量工具为测杖、卷尺，测量指标包括体高(x1)、体斜长(x2)、胸围(x3)、胸深(x4)、胸宽(x5)和管围(x6)。

1.3 数据处理

采用Excel 2007对所测数据进行初步整理，用SPSS 19.0统计分析软件进行相关性分析、通径分析、逐步回归分析，并对相关系数进行分解，建立体质量与体尺指标之间的最优回归方程。试验数据均以平均数±标准差表示。

2 结果与分析

2.1 体质量与体尺指标的表型统计量

表1为澳洲白羊体质量与体尺指标的表型统计量。由表1可见，体质量(y)的变异系数最大，为10.0%；其次为胸宽(x5)、胸深(x4)、管围(x6)、胸围(x3)、体斜长(x2)，其变异系数分别为7.0%，6.1%，5.3%，4.8%，4.6%；体高(x1)的变异系数最小，为4.2%。

表1 澳洲白羊体质量与体尺指标的表型统计量

2.2 体质量与体尺指标的表型相关性分析

由表2可见，澳洲白羊体质量(y)与胸围(x3)和管围(x6)相关性极显著(P<0.01)，与体高(x1)相关性显著(0.010.05)。体质量(y)与管围(x6)的表型相关系数最大，为0.680；与胸围(x3)、体高(x1)、体斜长(x2)、胸深(x4)、胸宽(x5)的表型相关系数为0.501，0.393，0.314，0.182，0.154。其中体斜长(x2)与管围(x6)、胸围(x3)与胸宽(x5)相关性极显著(P<0.01)，其表型相关系数分别为0.493，0.450。

表2 澳洲白羊体质量与体尺指标的表型相关系数

2.3 体质量与体尺指标之间的通径分析

为了进一步揭示体质量与体尺指标间的相关关系，对体质量、体尺指标进行了通径分析，用体尺指标与体质量的表型相关系数建立以通径系数为变量的多元一次方程组，方程组为：

P1+0.291P2+0.272P3+0.079P4+0.042P5-0.050P6=0.393

0.291P1+P2+0.177P3+0.006P4+0.208P5+0.493P6=0.314

0.272P1+0.177P2+P3+0.240P4+0.450P5+0.286P6=0.501

0.079P1+0.006P2+0.240P3+P4+0.110P5+0.313P6=0.182

-0.050P1+0.208P2+0.450P3+0.110P4+P5+0.740P6=0.154

0.345P1+0.493P2+0.286P3+0.313P4+0.740P5+P6=0.680

用矩阵法对多元一次方程组进行求解，得到澳洲白羊体尺指标对体质量的通径系数分别为：P1=0.121；P2=-0.081；P3=0.334；P4=-0.099；P5=-0.008；P6=0.614。

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

根据相关系数的组成效应，可将澳洲白羊各体尺指标(xi)与体质量(y)的表型相关系数(皮尔森相关系数)分解为体尺对体质量的直接作用和该体尺指标通过其他体尺性状对体质量的间接作用两部分。澳洲白羊体尺对体质量相关系数的分解见表3。

由表3可知，澳洲白羊的管围(x6)对体质量的直接作用最大，为0.614；胸深(x4)对体质量的直接作用最小，为-0.099；体斜长(x2)对体重(y)的综合间接作用最大，为0.394；其次为体高(x1)、胸深(x4)、胸围(x3)、胸宽(x5)，分别为：0.271，0.209，0.167，0.161；管围(x6)的综合间接作用最小，为0.066。各体尺指标对体质量的直接作用大小为：管围(x6)>胸围(x3)>体高(x1)>胸宽(x5)>体斜长(x2)>胸深(x4)；各体尺指标对体质量的间接作用大小为：体斜长(x2)>体高(x1)>胸深(x4)>胸围(x3)>胸宽(x5)>管围(x6)。可以看出，体高(x1)对体质量既有较强的直接作用，又有较强的间接作用；管围(x6)、胸围(x3)主要是通过直接作用对体质量产生影响，对体质量的间接作用较弱；体斜长(x2)、胸深(x4)主要是通过间接作用对体质量产生影响，对体质量的直接作用较弱；胸宽(x5)对体质量的直接作用、间接作用均较弱。体高(x1)、体斜长(x2)、胸围(x3)胸深(x4)主要是通过管围(x6)对体质量产生间接作用，胸宽(x5)和管围(x6)主要是通过胸围(x3)对体质量产生间接作用。

表3 澳洲白羊体尺对体质量相关系数的分解

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

2.5.1 体质量和各体尺指标整体回归分析模型 从表4～6中可以看出，体高(x1)、体斜长(x2)、胸深(x4)、胸宽(x5)、对体质量(y)的影响不显著，胸围(x3)、管围(x6)对体质量(y)的影响显著，因此应该剔除自变量体高(x1)、体斜长(x2)、胸深(x4)、胸宽(x5)后，对因变量体质量(y)与胸围(x3)、管围(x6)进行多元逐步回归分析。

表4 体质量和各体尺指标整体回归分析模型

表5 体质量和各体尺指标整体回归分析的多元回归方差分析

表6 体质量和各体尺指标总体回归分析的多元回归系数参数估计值

2.5.2 体质量和各体尺指标逐步回归分析模型 逐步回归分析即按自变量对因变量作用的显著程度，从大到小逐个引入回归方程，即将自变量的t值从大到小逐个引入回归方程。由表6可知，各个自变量的显著性由大到小依次为管围(x6)>胸围(x3)>体高(x1)(不显著)>胸深(x4)(不显著)>体斜长(x2)(不显著)>胸宽(x5)(不显著)，因此建立了体质量(y)和管围(x6)、胸围(x3)之间的回归分析模型2，结果见表7。由模型2所得到的方程为体质量(y)和体尺指标之间的最优回归方程为y=-39.698+0.355x3+5.064x6(R=0.752，P<0.05，表7)，标准化后的回归方程为y=0.334x3+0.585x6(表9)，方差分析的F值为18.815(表8)，达到极显著水平(P<0.01)。

表7 澳洲白羊体质量和胸围、管围指标逐步回归分析模型

表8 体质量和胸围、管围指标逐步回归分析的多元回归方差分析

表9 体质量和胸围、管围指标逐步回归分析的多元回归系数参数估计值

3 讨 论

3.1 体质量和体尺指标的表型统计量

变异系数又称离散系数，它反映的是相关研究变量变化的相对程度指标。澳洲白羊的体质量与其体尺指标的变异系数由大到小依次为：体质量(y)>胸宽(x5)>胸深(x4)>管围(x6)>胸围(x3)>体斜长(x2)>体高(x1)，体质量的变异系数为10%，说明澳洲白羊的体质量变异幅度较大，体质量具有较大的选择与提高潜力。体高的变异系数最小，为4.2%，说明体高的选择与提高潜力很小。

3.2 体质量与体尺指标的相关性分析

相关系数是用来衡量2个自变量或1个自变量与1个因变量之间相关程度的指标，由于在数据分析过程中笔者没有消除其他研究变量对这两个变量的影响，所以表型相关系数并不能够真实地反映2个自变量或自变量与因变量之间的相关关系，只能反映所研究的2个自变量或自变量与因变量之间的综合作用程度。本文通过对澳洲白羊体质量与其体尺指标进行相关性分析，可以看出澳洲白羊体质量(y)与胸围(x3)和管围(x6)相关性极显著(P<0.01)，与体高(x1)相关性显著(0.01

3.3 体质量与其体尺指标间的通径分析及其相关系数的分解

通径系数表示自变量对因变量的直接作用大小，其真实地反映了2个变量间的直接相关程度，它是标准化的偏回归系数[18-19]。通过对体质量和各体尺指标之间的通径分析以及其相关系数的分解可以看出：6月龄澳洲白羊体质量(y)与胸围(x3)和管围(x6)均表现出极显著相关性(P<0.01)。从表3中相关系数的分解中可以看出，体高(x1)对体质量的直接作用和间接作用均较强；管围(x6)、胸围(x3)主要是通过直接作用对体质量产生影响，对体质量的间接作用较弱；体斜长(x2)、胸深(x4)主要是通过间接作用对体质量产生影响，对体质量的直接作用较弱；胸宽(x5)对体质量的直接作用和间接作用均较弱。通过上述分析说明澳洲白羊的生长发育是一个复杂的有机整体，各指标间有着种种密切而复杂的联系。影响其体质量的因素不但多，而且影响途径亦不同。所以在选种过程中除了加强对体质量影响直接作用强的因素如澳洲白羊的管围、胸围选择的同时，也应该关注其他因素对体质量的作用。在选种时，可以参考本实验数据分析结果，即本实验所建立的澳洲白羊体质量与其体尺指标之间的最优回归方程，进行综合全面考虑，避免盲目选择，从而能达到最优的选种效果。

3.4 体尺指标对体质量的逐步回归分析

通过对澳洲白羊体质量与其体尺指标之间的逐步回归分析，得到它们之间的最优回归方程：y=-39.698+5.064x6+0.355x3(R=0.752，P<0.05)，标准化后的回归方程：y=0.334x3+0.585x6，方差分析的F值为18.815，达极显著水平(P<0.01)。体高(x1)、体斜长(x2)、胸深(x4)和胸宽(x5)回归系数显著程度不高(P值均大于0.05)，所以予以排除，说明这4个体尺指标对澳洲白的体质量影响不明显。这与杨燕等[15]对沂蒙黑山羊、李胜开等[16]对隆林山羊以及王欣荣等[17]对甘南草地型藏羊的研究结果类似。而他们的回归方程中除了管围、胸围入选外，也入选了体高、体斜长、胸深等因素，这主要是由于羊的品种、饲养管理水平、生活环境等诸多因素不同而造成的。从方程可知，当胸围保持不变时，管围每增加1 cm，澳洲白羊的体质量增加0.585 kg；当管围保持不变时，胸围每增加1 cm，杜泊羊的体质量增加0.334 kg。回归方程中R2=0.565，说明该回归方程的回归平方和占总变异的56.5%，回归方程显著存在，所选体尺指标对体质量有较大的决定作用，方程拟合度较高，具有一定的参考价值。剩余因子e=1-R2=0.435，说明对澳洲白羊体质量产生影响的自变量不止管围、胸围，还有一些其他影响因素没有考虑到，因此下一步需要对影响澳洲白羊体质量的因素(包括澳洲白羊的体尺指标、营养水平、生存环境)进行全面分析。

4 结 论

在本试验条件下，6月龄澳洲白羊管围(x6)和胸围(x3)是影响其体质量的主要因素，体高(x1)、体斜长(x2)、胸围(x3)、胸深(x4)、主要是通过管围(x6)对体质量产生间接作用，胸宽(x5)和管围(x6)，主要是通过胸围(x3)对体质量产生间接作用。其体质量和体尺指标之间的最优回归方程y=-39.698+5.064x6+0.355x3(R=0.752，P<0.05)，标准化后的回归方程：y=0.334x3+0.585x6(F=18.815，P<0.01)。