鲁中肉羊体重与体尺指标的相关性及回归分析

杨存明，张梦华，何军敏，张晓雪，刘桂芬，魏晨，任一帆，毛静艺，李雪，张文静，田可川，黄锡霞

(1. 新疆农业大学动物科学学院，新疆 乌鲁木齐 830052；2. 山东省农业科学院畜牧兽医研究所，山东 济南 250100)

绵羊育种以提高其重要经济性状的生产水平为主要目的[1]。 鲁中肉羊是以白头杜泊羊为父本，湖羊为母本，杂交繁育十余年育成的肉用绵羊新品种，具有成长速度快、繁殖率高、产肉性能好、适应性强等优良特性[2-3]。 生长性状是肉羊育种的重要考察因素，体尺和体重都是衡量生长发育的重要指标，体尺数据反映畜体的生长发育状况及各部位间的相对发育关系，体现家畜个体或群体的生产性能，而体重则是检验饲养管理水平的主要依据[4]。 对体重与体尺指标进行相关性分析、通径系数分析及逐步回归分析，有助于提高相关性状的选择精度，为保障绵羊生产系统的整体效率和盈利能力提供依据[5-6]。 前人采用这些方法对不同品种羊，如杜泊羊、萨福克母羊、陶赛特母羊[7]、青海无角欧拉羊羔羊[8]、大尾寒羊[9]、科罗基羊[10]、乌达羊[11]、杜湖杂种羔羊[12]、云上黑山羊[13]、澳洲白羊[14]、公滩羊[15]等进行了研究。

本研究以周岁鲁中肉羊为对象，测定体重和体尺指标参数，对相关数据进行了初步统计，并对体重和体尺指标进行了相关性分析、通径分析以及多元回归分析，建立了最优回归方程，以期为优化鲁中肉羊的饲养管理结构奠定基础，并为后期生产实践以及品种选育工作提供参考。

1 材料与方法

1.1 试验动物

试验于山东嬴泰农牧科技有限公司进行。 随机选取周岁鲁中肉羊共计1 289 只，其中公羊330只(平均体重83.49 kg)，母羊959 只(平均体重60.73 kg)，所有羊只均发育良好，健康无病，实行常规免疫，舍饲圈养。

1.2 测定指标及方法

参照《家畜育种学》[1]和《绵、山羊生产性能测定技术规范》[16]对周岁鲁中肉羊公羊和母羊的体重(Y)、体高(X1)、体长(X2)、胸围(X3)、胸宽(X4)、胸深(X5)、管围(X6)进行测定，每个指标2次重复。

1.3 数据处理与分析

通过Microsoft Excel 2019 对试验数据进行初步整理，采用拉依达准则(平均值±3 倍标准差)将指标异常值进行剔除。 采用SPSS 25.0 软件对数据进行描述性统计，以各体尺指标为自变量，体重为因变量进行相关性分析，并采用多元回归分析中的逐步回归法建立周岁鲁中肉羊体重和体尺指标的最优回归方程，最后通径分析各体尺性状对体重的直接和间接作用。

2 结果与分析

2.1 鲁中肉羊体重和体尺指标的测定

周岁鲁中肉羊公羊和母羊的体重和体尺指标测定结果见表1、表2。 结果显示，周岁鲁中肉羊公羊体重的变异系数较大，为9.17%，母羊体重变异系数为2.67%，公羊和母羊体尺性状中变异系数最大的分别是胸宽和体高，分别为5.01%和6.91%。

表1 周岁鲁中肉羊公羊体重、体尺指标测定结果(n＝330)

表2 周岁鲁中肉羊母羊体重、体尺指标测定结果(n＝959)

2.2 鲁中肉羊体重与体尺指标的相关性分析

从表3 可知，鲁中肉羊公羊体重与胸深呈显著负相关(P＜0.05)，胸深与胸围呈显著正相关(P＜0.05)，体高与胸深呈极显著负相关(P＜0.01)，体高与管围间无显著相关性(P＞0.05)，其余各性状间均呈极显著正相关(P＜0.01)。 各体尺指标与体重间的相关系数依次为体高＞胸围＞胸宽＞体长＞管围＞胸深。 母羊体长与各体尺指标及体重之间均呈极显著负相关(P＜0.01)，其余各体尺指标之间以及各体尺指标与体重之间均呈现极显著正相关(P＜0.01)，各体尺指标与体重的相关系数依次为胸宽＞胸围＞体高＞胸深＞管围＞体长。

表3 鲁中肉羊各性状间相关系数分析

2.3 鲁中肉羊体重与体尺指标的通径分析

由表4 可知，在鲁中肉羊公羊各体尺指标中，体高对体重的直接作用最大，通径系数为0.325，其他指标对体重的直接作用由高到低依次为胸围、胸宽、体长和胸深，通径系数分别为0.289、0.228、0.218 和-0.157。 胸围对体重的间接作用最大，通径系数为0.615，其他指标对体重的间接作用由高到低依次为体高、胸宽、体长、胸深，通径系数分别为0.552、0.494、0.375 和-0.011。 由此可见，体高、胸围和胸宽对鲁中肉羊公羊体重的影响无论是直接作用还是间接作用均较大，说明其体高、胸围和胸宽与体重密切相关。

表4 鲁中肉羊公羊体尺对体重的通径分析

由表5 可知，在鲁中肉羊母羊各体尺指标中，胸宽对体重的直接作用最大，通径系数为0.689，其次为胸围、管围、体长和体高，通径系数分别为0.342、-0.151、-0.279 和-0.294。 胸宽对体重的间接作用最大，通径系数为0.776，其他指标由高到低依次为体长、胸围、管围和体高。 说明鲁中肉羊母羊胸围和胸宽同样与体重密切相关。

表5 鲁中肉羊母羊体尺对体重的通径分析

2.4 鲁中肉羊体尺指标与体重的回归分析

以周岁鲁中肉羊体重为因变量，各体尺指标为自变量，采用多元回归分析中的逐步回归分析法建立多元线性回归方程。 由表6 可知，模型1～5 中鲁中肉羊公羊体重分别与体高、胸围、体长、胸深、胸宽呈极显著线性回归关系(P＜0.01)。 由此建立以下5 个线性回归方程:模型1:Y ＝2.166X1-67.778，校正R2＝0.693；模型2:Y ＝1.381X1+1.181X3-132.499，校正R2＝0.776；模型3:Y ＝1.273X1+0.482X2+0.990X3-144.109，校正R2＝0.801；模型4:Y ＝1.073X1+0.577X2+1.125X3-0.640X5-129.685，校正R2＝0.810；模型5:Y ＝0.845X1+0.554X2+0.818X3+1.245X4-0.868X5-107.903，校正R2＝0.826。 模型5 的校正R2值最高，因此其建立的方程为最优回归方程。

表6 鲁中肉羊公羊体重与体尺的回归分析

由表7 可知，模型1～5 中鲁中肉羊母羊体重分别与胸宽、管围、胸围、体长、体高呈极显著线性回归关系(P＜0.01)。 由此建立以下5 个线性回归方程:模型1:Y ＝1.537X4+22.071，校正R2＝0.800；模型2:Y＝1.616X4-0.998X6+28.625，校正R2＝0.819；模型3:Y ＝0.345X3+1.299X4-1.792X6+11.829，校正R2＝0.836；模型4:Y ＝-0.032X2+0.342X3+1.220X4-1.770X6+15.911，校正R2＝0.839；模型5:Y ＝-0.104X1-0.125X2+0.420X3+1.184X4-1.033X6+16.088，校正R2＝0.844。 模型5 的校正R2值最高，为最优回归方程。

3 讨论

3.1 鲁中肉羊体重与体尺性状之间的相关性分析

本试验分析结果表明，鲁中肉羊的体尺指标与体重有极显著的相关性(P＜0.01)，其中公羊体高与体重的相关性最强，其次为胸围；母羊胸宽与体重的相关性最强，其次为胸围。 可能由于本研究所选羊只均为周岁鲁中肉羊，体型基本稳定，所以其体高、胸宽和胸围等指标与羊只的肌肉、骨骼、内脏的生长发育以及膘情等情况密切相关。在青海无角欧拉羊羔羊[8]中，体高对体重的影响最大，其次为胸围；大尾寒羊[9]体高和胸围不仅对体重有较大的直接作用，还辅助其他性状有较大的间接作用，这均与本研究结果类似。 而对科罗基羊[10]、乌达羊[11]、杜湖杂种羔羊[12]、云上黑山羊[13]、澳洲白羊[14]、滩羊[15]、杜泊羊、特克塞尔羊[17]等研究中，均显示胸围与体重的相关系数最大，其差异可能是由于品种不同和羊只年龄不一致所导致的。 因此，在鲁中肉羊的品种选育过程中，应以体重为主，兼顾体高、胸围、胸宽等其他体尺性状。

3.2 鲁中肉羊体重与体尺指标的通径分析

本研究通过对周岁鲁中肉羊体尺体重间的通径分析发现，公羊和母羊体尺指标中对体重的直接作用最大的分别为体高和胸宽，其次均为胸围，间接作用最大的均为胸宽，其次分别为胸围和体高。 对成年欧拉羊[18]、云上黑山羊[13]、酉州乌羊[19]、海南黑山羊[20]的分析均表明胸围是影响体重的主要体尺指标，与本研究中对鲁中肉羊公羊的分析结果类似，但存在一定差异，可能是由所选品种和饲养水平不同造成的。 在青海无角欧拉羊羔羊[8]的研究中，体高对体重的直接影响最大，与本研究中对鲁中肉羊公羊的分析结果一致。本研究中，通径分析结果和相关性分析结果均表明鲁中肉羊公羊主要通过体高、胸围，母羊主要通过胸宽、胸围直接或间接影响体重。 因此，在进行鲁中肉羊体重性状选育时，应区分性别，公羊、母羊分别以体高、胸围和胸宽、胸围指标为主要参考。

3.3 鲁中肉羊体尺指标与体重的回归分析

本研究通过多元回归分析建立了周岁鲁中肉羊公羊与母羊体尺指标与体重的最优回归方程，与对青海无角欧拉羊[8]、杜湖杂种羔羊[12]、巴特肯羊[21]、巴美肉羊[22]、布拖黑绵羊育成母羊[23]等品种羊的研究结果相似。 剩余因子(e ＝1-R2)分别为0.174 和0.156，说明入选回归方程的体尺性状对体重有较大的决定作用。 但还有其他因素未考虑，因此下一步需要对影响鲁中肉羊体重的其他因素进行进一步分析。

4 结论

周岁鲁中肉羊公羊体高与体重的相关性最强，母羊胸宽与体重的相关性最强，其次均为胸围。 体高和胸宽主要通过直接作用影响体重，而胸围、体长、胸深和管围主要通过间接作用影响体重。 因此，在鲁中肉羊品种选育过程中，公羊应重点关注体高、胸围，母羊应重点关注胸宽，从而提高选育速度。 鲁中肉羊公羊体尺指标与体重的最优方程为Y ＝0.845X1+0.554X2+0.818X3+1.245X4-0.868X5-107.903(P＜0.01，R2＝0.826)，母羊为Y＝-0.104X1- 0. 125X2+ 0. 420X3+ 1. 184X4-1.033X6+16.088(P＜0.01，R2＝0.844)，可作为利用体尺指标评估体重的参考依据。