青海无角欧拉羊羔羊体重与体尺多元回归分析

庄蕾，阎明毅，余忠祥，张强龙，唐燕花，多杰才让，唐燕青，李永钦，吴森

(1.青海大学畜牧兽医科学院 青海省高原家畜遗传资源保护与创新利用重点实验室，西宁810016;2.青海省河南蒙古族自治县生态有机畜牧业技术服务中心，河南县811500;3.青海省河南蒙古族自治县有机产业科技园区管理委员会，河南县811500;4.青海省河南蒙古族自治县托叶玛乡农牧业技术服务站，河南县811500;5.青海省河南蒙古族自治县兽医站，河南县811500)

欧拉羊是青藏高原藏系绵羊的一个重要类群，曾被称为“欧拉型藏羊”，2018年被列入国家畜禽资源，更名为“欧拉羊”，因其主产区位于欧拉山而得名。欧拉羊具有发育快、体格结实、肢体高大、生长快、繁殖性能好、适应高寒气候、耐粗饲、肉质好等特点［1，2］，其成年羊体重在60 kg以上，远大于其他藏羊。另外，欧拉羊主产区位于“中华水塔”三江源腹地的河南蒙古族自治县，拥有国家认证的有机牧场，其肉质独特，被世人所关注［3］。欧拉羊具备较高的开发利用价值，近年来，在本品种不断选育的基础上，现选育出大群无角型欧拉羊，性格温顺，易于管理，且符合现代规模化、集约化养殖趋势及高效畜牧业生产理念，欧拉羊是良好的现代化藏羊品种育成资源。

羔羊期是绵羊生长的关键时期，容易受代谢印记影响后期发育。羔羊期生长发育相关研究有助于群体早期选种及生长发育评估，特别是公羔羊。公羔羊的体况评估研究有利于种公羊的早期选留，减少不必要的养殖负担，实现种群的快速优化。本研究对无角型欧拉羊公羔羊的体重与体尺性状进行了线性回归分析、通径分析、逐步回归分析，建立了最优回归方程，旨在为无角型欧拉羊的早期选种选育、优良种公羊培育等工作提供技术和理论参考。

1 材料与方法

1.1 试验动物

本研究以无角型欧拉羊为研究对象，测定收集青海省黄南藏族自治州河南蒙古族自治县优干宁镇经多年本品种选育的无角型欧拉羊公羔羊(8～10月龄)653只，所有测定羊均发育正常、生长健康，实行常规免疫，自然放牧。

1.2 测定指标与方法

剪毛后对其体重(Y)和主要体尺性状，即体高(X1)、体长(X2)、胸宽(X3)、胸深(X4)、尻宽(X5)、胸围(X6)进行测定，生产性能测定主要参照张沅主编《家畜育种学》的第三章第七节体型鉴定［4］及《绵山羊生产性能测定技术规范》行业标准［5］进行，由固定测量专员对体高、体长、胸宽、胸深、尻宽、胸围进行测定，每个部位测定3次，取平均值。

1.3 逐步回归分析

建立最优回归方程的方法很多，常用的是逐步回归法(Stepwise)［6］，其原理是将自变量逐个引入，对全部的自变量X1，X2，...，Xp，按其对因变量Y影响的大小进行比较，通过F检验，选择偏回归平方和显著的变量进入回归方程。同时，每引入一个新变量后，对前面已引入回归模型的变量逐个进行检验，将经检验后偏回归平方和不显著的变量删除，以保证所得自变量子集中每一个变量都是显著的，经过多次循环，直到不能再引入新变量为止。此时的多元线性回归方程即为最优多元线性回归方程［7，8］。

1.4 数据处理

用Excel进行数据基本整理，用SPSS 23.0统计软件进行统计分析，对相关系数进行分解，建立体重与各体尺性状之间的最优回归方程。

2 结果与分析

2.1 无角型欧拉羊体重与体尺性状的基本统计量

由表1可知，无角型欧拉羊体高变异系数是7.87%，为弱变异(CV＜10%)。其余表型均为中等变异(10%＜CV＜100%)，胸宽变异系数为57.36%，接下来依次是体长43.66%、尻宽36.63%、体重23.93%、胸深12.36%、胸围11.98%。

表1 无角型欧拉羊体重与体尺性状的基本统计量

2.2 无角型欧拉羊体重与体尺性状的相关性分析

由图1可知，体重与体高、胸深、胸围之间存在极显著正相关(P＜0.01)，相关程度由大到小为体高＞胸围＞胸深;体重与体长、胸宽、尻宽之间存在显著正相关(P＜0.05)，相关程度由大到小为体长=胸宽＞尻宽;其余各性状两两之间也存在正相关。

图1 无角型欧拉羊体重与体尺性状的相关性分析(Pearson)

2.3 无角型欧拉羊体尺指标和体重相关系数的分解

对无角型欧拉羊的各体尺性状与体重数据进行逐步回归分析。根据相关系数的组成效应，可将无角型欧拉羊羔羊的各体尺性状(Xi)对体重(Y)的相关系数划分为直接作用(Pi)和通过性状(Xj)的间接作用(∑RijPj)两部分，即rxiy=Pi+∑RijPj。由表2可知，体高对体重的直接作用最大，为0.421;胸深对体重的间接作用最大，为0.2034。

表2 无角型欧拉羊体高、体长、胸宽、胸深、尻宽、胸围和体重相关系数的分解

2.4 最优回归方程的建立

回归方程1引入的体尺性状是体高，回归方程2引入的体尺性状是体高和胸围，回归方程3引入的体尺性状是体高、胸围、胸深。由表3可知，体高、胸围、胸深被依次引入回归方程后，其R值逐渐变大，估计标准误逐渐变小。

表3 无角型欧拉羊体重与体尺的逐步回归模型与参数估计值检验

由上述分析可知，回归方程3的决定系数最大，所以回归方程3为最优回归方程，复相关系数R为0.616，校正的决定系数R2为0.376。回归方程3中对体重影响最显著的体尺性状是体高、胸围、胸深，各性状标准偏回归系数为0.425、0.253、0.122。经T检验，体高、胸围、胸深均达到极显著水平(P＜0.01)。综上所述，逐步回归方程Y=0.425X1+0.253X6+0.122X4为最优回归方程，但Y=0.544X1可作为野外无称重条件下快速预测欧拉羊体重的最优回归方程。

3 讨论

3.1 体重与体尺的表型分析

体重和体尺性状是绵羊外貌体况评定的重要指标，不仅反映绵羊品种自身的体型特征，同时与绵羊的生产性能有密切联系［9，10］。变异系数可以将不同单位观测数据的相对变化程度直观的表现出来，变异系数低于10%时称为弱变异，在10%～100%之间称为中等变异［7］。在本研究中，无角型欧拉羊的体高变异系数是7.87%，为弱变异;体重及其余体尺性状变异系数均处于10%～100%之间，为中等变异，且各性状变异系数均较大，说明无角型欧拉羊公羔羊在体重、体尺性状方面具有较大的选育潜力。

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

体重是重要的经济与选种指标，主要受到遗传和环境的共同影响，因此在同品种内和同环境下，研究体重与各体尺性状之间的关系具有一定的代表性［11］。无角型欧拉羊的体高、胸深、胸围与体重之间存在极显著正相关，体长、胸宽、尻宽与体重之间存在显著正相关，且体高、胸深、胸围对体重均有直接或间接作用，说明这3个体尺性状对无角型欧拉羊在体重增长方面影响较大。在对无角型欧拉羊体重增长影响最大的3个体尺性状中，由相关系数可知，体高对体重的影响最大，为0.544;其次是胸围，为0.429;胸深对体重的影响最小，为0.324。体重与体高、胸围、胸深之间的相关是由于体重与动物的肌肉、骨骼、内脏的生长发育紧密相关［12］。在梁小鹏等［13］对新疆也木勒白羊、许鑫等［14］对特克塞尔羊的研究中发现，不仅胸围对体重会有较大影响，还有体长等体尺性状也会对体重有较大影响，结果与本研究不同，可能是本文测定的无角型欧拉羊的月龄较小和品种差异造成的。白俊艳等［15］对大尾寒羊的研究中发现，体高和胸围不仅对体重有较强的直接作用，还辅助其他性状产生较大的间接作用，与本研究结果相似。在李江文等［16］对呼伦贝尔短尾羊、海占忠等［17］对加什科羊、杨卓等［18］对青海加什科羊等的研究中也有相似的结果。因此，今后在无角型欧拉羊的品种选育、提高生产性能等研究中，应以体重为主，兼顾体高、胸围、胸深，以取得较好的选育效果。

3.3 最优回归方程的建立

体尺性状不仅反映其外貌的特征，同时也可以为体重的预估提供参考［17］。本研究以体重为因变量，体尺性状为自变量建立逐步回归方程:Y=0.425X1+0.253X6+0.122X4，因为体高(X1)、胸深(X4)、胸围(X6)对体重影响较大，所以被选入回归方程，拟合优度的决定系数为0.379，校正R2值为0.376，回归关系达到极显著水平(P＜0.01)，表明应用该回归方程能够较准确地评估无角型欧拉羊的体重变化。

综上所述，体高、胸深、胸围是影响无角型欧拉羊公羔羊体重的主要因素，无角型欧拉羊公羔羊的体重-体尺最优回归方程为Y=0.425X1+0.253X6+0.122X4为最优回归方程，Y=0.544X1可作为野外无称重条件下快速预测欧拉羊体重的最优回归方程。

4 结论

无角型欧拉羊公羔羊各体尺性状两两之间存在正相关，体高、胸深、胸围与体重之间存在极显著正相关，这3个体尺性状是影响体重的主要因素，因此在对无角型欧拉羊进行选育时，要重点关注体高、胸深、胸围3个性状，以获得更好的选育结果。