新疆也木勒白羊公羊体尺指标对体重影响力的分析

梁小鹏，张振良，阿曼开迪·莫哈麦提汗，姚力丹，包敖敦格日乐，决肯·阿尼瓦什

(1. 新疆农业大学动物科学学院，乌鲁木齐830052；2. 新疆塔城地区畜牧科技推广研究中心，新疆塔城834700)

0 引 言

【研究意义】也木勒白羊，主要分布在额敏县周边。也木勒白羊作为新疆塔城地区优良的肉毛兼用地方品种，它具有适应性强、抗病力高、抗寒冷、耐炎热、耐粗饲、生长发育迅速的特点[1]。2012年统计显示，也木勒白羊存栏仅万余只，2013年额敏县也木勒白羊选育推广中心养殖基地成立，负责也木勒白羊的保种和繁育工作，但是，由于各种原因，一直未能得到足够的重视，在饲养和选种、选育方面存在很多问题[2]。在选种过程中，由于体重的测量工作较大，而用几个体尺指标来估计体重相对简单，建立一个合理的线性方程尤为重要。【前人研究进展】对也木勒白羊生长发育规律、生产性能的研究，阿曼开迪·莫哈买提汗等[3,4]研究显示也木勒白羊属于早熟型绵羊品种，羔羊在哺乳期生长发育最快，4个半月龄就可达到出栏体重，此时屠宰率可达53.5%，净肉率41.2%。【本研究切入点】国内外用体尺或其他指标来评估体重估计体重的研究很多[5-7]，方法也多种多样，但对也木勒白羊体尺指标与体重的相关关系的研究还未有报道。对木勒白羊公羊体尺指标对体重的影响力进行关联性分析。【拟解决的关键问题】使用Excel和SPSS19.0软件对50只也木勒白羊公羊体尺与体重数据进行整理和分析，研究也木勒白羊公羊体重与体尺指标的相关关系，分析影响也木勒白羊公羊体重的较主要体尺指标。建立体重与体尺指标间最优的线性回归方程，为也木勒白羊公羊以后的选育留种提供数据参考。

1 材料与方法

1.1 材 料

试验于额敏县也木勒白羊选育推广中心养殖基地进行。

随机抽取的50只也木勒白羊种公羊来自也木勒白羊选育推广中心养殖基地，年龄在2周岁以上，健康无病，处于半放牧状态。

1.2 方 法

试验测定项目包括体重(y)、体高(x1)、体长(x2)、胸宽(x3)、胸深(x4)、胸围(x5)、管围(x6)、脂臀宽(x7)、脂臀长(x8)、脂臀厚(x9)。使用家用型直视式电子吊秤来称量体重，精度为0.01 kg，体尺的测量依照《家畜育种学》[8]，使用测杖和1.5 m皮软尺测量。

1.3 数据处理

使用Excel 初步整理数据，然后用SPSS 19.0进行统计分析、相关性分析、通径分析、逐步回归分析，方法依照罗凤明等[9]，最后建立体重与体尺指标最优的线性回归方程。

2 结果与分析

2.1 体重与体尺表型

研究表明，体重的平均值为63.4 kg，极大值为83.95 kg，极小值为50.55 kg。脂臀厚的平均值为16.4 cm，极大值为22 cm，极小值为12 cm。胸围的极大值为107 cm，极小值为87 cm，平均值为97.24 cm。在所有指标中，脂臀厚的变异系数最大，为13.32%，接下来是脂臀长和体重、变异系数为13.21%、13.08%，然后是脂臀宽、胸宽、体长、体高、胸深、胸围、管围，其变异系数为11.32%、8.25%、7.05%、6.57%、6.21%、5.6%、5.09%。表1

表1 也木勒白羊体重与体尺表型统计
Table 1 Quantities of body weight and body size phenotype

数量Quantity极小值Minimum极大值maximum均值Mean标准差Standard deviation变异系数(%)Coefficient of variation体重 y (kg)Body weight (y)5050.5583.9563.48.29413.08体高 x1 (cm)body height 50618168.964.5316.57体长 x2 (cm)body length 50557966.364.6767.05胸围 x3 (cm)chest circumference 508710997.245.4495.6胸宽 x4 (cm)chest width 50202723.521.948.25胸深 x5 (cm)chest depth 50273732.762.0366.21管围 x6 (cm)tube circumference 50898.220.4185.09脂臀长 x7 (cm)fat hip length 50112013.921.83913.21脂臀宽 x8 (cm)fat hip width 50172821.522.43511.32脂臀厚 x9 (cm)fat hip thickness 50122216.42.18513.32

注：体重(y)体高(x1)体长(x2)胸围(x3)胸宽(x4)胸深(x5)管围(x6)脂臀长(x7)脂臀宽(x8)脂臀厚(x9)

Note: Body weight (y) body height (x1) body length (x2) chest circumference (x3) chest width (x4) chest depth (x5) tube circumference (x6) fat hip length (x7) fat hip width (x8) fat hip thickness (x9)

2.2 体重与体尺相关性

研究表明，体重数据正态分布。体高、体长、胸围、胸宽、胸深、脂臀长、脂臀宽、脂臀厚与体重相关性均呈极显著正相关(P< 0.01),体重与管围呈显著正相关(P< 0.05)。其中体重与胸围的相关系数最大，为0.886；体重与管围的相关系数最小，为0.336；其他体尺指标与体重相关性由大到小依次为胸宽、脂臀宽、胸深、体长、体长、体高、脂臀长、脂臀厚，相关系数依次为0.82、0.793、0.783、0.646、0.571、0.559、0.391。也木勒白羊各体尺指标间也有着极显著的相关性，其中体长、胸围、胸宽、胸深、脂臀宽、脂臀厚与体高相关性极显著(P< 0.01)；胸围、胸宽、胸深、脂臀宽与体长相关性极显著(P< 0.01)；胸宽、胸深、管围、脂臀长、脂臀宽与胸围相关性极显著(P< 0.01)；胸深、脂臀长、脂臀宽、脂臀厚与胸宽相关性极显著(P< 0.01)；管围、脂臀长、脂臀宽、脂臀厚与胸深相关性极显著(P< 0.01)；脂臀厚与管围相关性极显著；脂臀宽、脂臀厚与脂臀长相关性极显著(P< 0.01)。表2，图1

表2 也木勒白羊体重与体尺数据相关性
Table 2 Results of correlation analysis between body weight and body size data of Amalair White Sheep

项目yx1x2x3x4x5x6x7x8x9y10.571∗∗0.646∗∗0.886∗∗0.820∗∗0.783∗∗0.336∗0.559∗∗0.793∗∗0.391∗∗x10.571∗∗10.658∗∗0.511∗∗0.395∗∗0.621∗∗0.220.294∗0.433∗∗0.369∗∗x20.646∗∗0.658∗∗10.543∗∗0.393∗∗0.599∗∗0.282∗0.343∗0.573∗∗0.261x30.886∗∗0.511∗∗0.543∗∗10.795∗∗0.804∗∗0.379∗∗0.576∗∗0.741∗∗0.321∗x40.820∗∗0.395∗∗0.393∗∗0.795∗∗10.668∗∗0.2080.412∗∗0.698∗∗0.369∗∗x50.783∗∗0.621∗∗0.599∗∗0.804∗∗0.668∗∗10.399∗∗0.540∗∗0.614∗∗0.467∗∗x60.336∗0.220.282∗0.379∗∗0.2080.399∗∗10.289∗0.306∗0.370∗∗x70.559∗∗0.294∗0.343∗0.576∗∗0.412∗∗0.540∗∗0.289∗10.401∗∗0.394∗∗x80.793∗∗0.433∗∗0.573∗∗0.741∗∗0.698∗∗0.614∗∗0.306∗0.401∗∗10.301∗x90.391∗∗0.369∗∗0.2610.321∗0.369∗∗0.467∗∗0.370∗∗0.394∗∗0.301∗1

注：**表示相关性极显著(P< 0.01),*表示相关性显著(P< 0.05)。体重(y)体高(x1)体长(x2)胸围(x3)胸宽(x4)胸深(x5)管围(x6)脂臀长(x7)脂臀宽(x8)脂臀厚(x9)

Note:**indicatesthatthecorrelationisextremelysignificant(P < 0.01),and*indicatesasignificantcorrelation(P < 0.05).Bodyweight(y)bodyheight(x1)bodylength(x2)chestcircumference(x3)chestwidth(x4)chestdepth(x5)tubecircumference(x6)fathiplength(x7)fathipwidth(x8)fathipthickness(x9)

图1 体重的正太分布Q-QFig.1NormaldistributionQ-Qdiagramofbodyweight

2.3 体重与体尺的通径分析

体重与体尺指标的相关分析结果仅能反映两个指标之间表型的相关作用程度，为了进一步探索体尺指标与体重的因果关系程度，需要对也木勒白羊公羊体重与体尺指标进行通径系数分析，通径系数是表示两个变量间因果关系程度的指标。根据体尺指标与体重的相关性分析结果，建立了以通径系数为变量的多元一次方程组，方程组如下：

P1+0.658P2+0.511P3+0.395P4+0.621P5+0.22P6+0.294P7+0.433P8+0.369P9=0.571.

0.658P1+P2+0.5437P3+0.393P4+0.599P5+0.282P6+0.343P7+0.573P8+0.261P9=0.646.

0.511P1+0.543P2+P3+0.795P4+0.804P5+0.379P6+0.576P7+0.741P8+0.321P9=0.886.

0.395P1+0.393P2+0.795P3+P4+0.668P5+0.208P6+0.412P7+0.698P8+0.369P9=0.82.

0.621P1+0.599P2+0.804P3+0.668P4+P5+0.399P6+0.54P7+0.614P8+0.467P9=0.783.

0.22P1+0.282P2+0.379P3+0.208P4+0.399P5+P6+0.289P7+0.306P8+0.37P9=0.336.

0.294P1+0.343P2+0.576P3+0.412P4+0.54P5+0.289P6+P7+0.306P8+0.37P9=0.559.

0.433P1+0.573P2+0.741P3+0.698P4+0.614P5+0.306P6+0.401P7+P8+0.301P9=0.793.

0.369P1+0.261P2+0.321P3+0.369P4+0.467P5+0.37P6+0.394P7+0.301P8+P9=0.391.

在Excel表格中，使用逆矩阵法对上多元一次方程组求解，得到也木勒白羊公羊各体尺指标与体重的直接通径系数为：P1=0.05、P2=0.17、P3=0.34、P4=0.29、P5=0.03、P6=-0.01、P7=0.08、P8=0.16、P9=0.02。

研究表明，胸围对体重(0.34)的直接影响作用最大，管围对体重(-0.01)的直接影响作用最小。其他体尺指标对体重的直接影响作用按大小顺序为胸宽(0.29)、体长(0.17)、脂臀宽(0.16)、脂臀长(0.08)、体高(0.05)，胸深(0.03)、脂臀厚(0.02)。其中胸围、体长、胸宽、脂臀宽是直接影响体重较主要的指标，相对重要。间接影响作用最大的是胸深(0.751)，接着是脂臀宽(0.627)、胸围(0.545)、胸宽(0.535)、体高(0.521)、脂臀长(0.477)、体长(0.476)、脂臀厚(0.373)，间接影响作用最小的是管围(0.341)。其中胸围和胸宽对体重的影响除了有较强的直接影响作用，还有较强的间接影响作用。表3

表3 也木勒白羊公羊不同体尺指标下对体重的变化
Table3EffectofbodysizeindexofGongmulawhitesheepramonbodyweight

体高x1 (cm)体长x2 (cm)胸围x3(cm)胸宽x4(cm)胸深x5(cm)管围x6(cm)脂臀长x7(cm)脂臀宽x8(cm)脂臀厚x9(cm)相关系数0.5710.6460.8860.820.7830.3360.5590.7930.391直接作用0.050.170.340.290.03-0.010.080.160.02间接作用总和0.5210.4760.5450.5350.7510.3410.4770.6270.373体高0.0340.0260.020.0320.0110.0150.0220.019体长0.1120.0920.0670.1020.0480.0580.0970.044胸围0.1740.1850.2710.2740.1290.1970.2530.11胸宽0.1130.1120.2270.1910.0590.1180.1990.105胸深0.0210.020.0270.0220.0130.0180.020.015管围-0.001-0.002-0.002-0.001-0.002-0.002-0.002-0.002脂臀长0.0240.0280.0470.0340.0450.0240.0330.032脂臀宽0.0710.0940.1220.1150.1010.050.0660.05脂臀厚0.0070.0050.0060.0070.0080.0070.0070.005

注：体高(height)，体长(Body length)，胸围(bust)，胸宽(Chest width)，胸深(Thoracic depth)，管围(Pipe circumference)，脂臀长(Fat hip length)，脂臀宽(Fat hip width)，脂臀厚(Fat hip thick)，相关系数(Correlation coefficient)，直接作用(Direct action)，间接作用总和(Indirect effect sum)

2.4 体重与体尺的逐步回归分析

研究表明，各模型的体尺指标与体重的拟合度和决定系数，其中模型4的决定系数最大，所以模型4为最优模型，负相关系数R为0.935，校正的决定系数R2为0.874。所以可以剔除模型1、2和3，只考虑模型4建立方程。模型4所选取的对体重影响最显著的体尺指标为胸围、脂臀宽、胸宽、体长，各指标偏回归系数分别为0.647、0.518、1.248、0.378。标准偏回归系数为0.425、0.152、0.292、0.213。经T检验，其中胸围、胸宽和体长达到极显著水平(P< 0.01)，脂臀宽达到显著水平(P< 0.05)。F值为78.069，可以看出，达到了极显著水平(P< 0.01)，模型有效。表4～6

表4 体重与体尺逐步回归分析模型拟合汇总
Table 4 Summary of model fit analysis of body weight and body size stepwise regression analysis

模型ModelRR2调整 R2Adjustment R2标准估计的误差Standard estimated error10.886a0.7840.783.891 7320.909b0.8260.8183.535 0830.920c0.8460.8363.360 8840.935d0.8740.8633.071 65

注：a. 预测变量: (常量), 胸围。b. 预测变量: (常量), 胸围, 脂臀宽。c. 预测变量: (常量), 胸围, 脂臀宽, 胸宽。d. 预测变量: (常量), 胸围, 脂臀宽, 胸宽, 体长。c. 预测变量: (常量), 胸围, 脂臀宽, 体高

Note: a. Predictors: (constant), bust. b. Predictors: (constant), bust, fat hip width. c. Predictors: (constant), bust, fat hip width, chest width. d. Predictors: (constant), bust, fat hip width, chest width, body length. c. Predictors: (constant), bust, fat hip width, height

表5 体重与体尺逐步回归分析后模型3回归系数参数
Table 5 Regression coefficient parameters of model 3 after stepwise regression analysis of body weight and body size

模型model非标准化系数Non-standardized coefficient标准系数Standard coefficientB标准 误差Standard coefficient(试用版)trial versiontSig4(常量)-65.1139.009-7.2270胸围0.6470.1530.4254.2430脂臀宽0.5180.2960.1521.7530.086胸宽1.2480.3920.2923.1810.003体长0.3780.1190.2133.1730.003

表6 体重与体尺逐步回归分析后模型3方差分析
Table 6 Model 3 Variance Analysis Table after Stepwise Regression Analysis of Body Weight and Body Size

模型model平方和sum of squaredf均方Mean squareFSig4回归(return)2 946.3464736.58778.0690.000d残差(Residual)424.578459.435总计(total)3 370.92449

注：d. 预测变量: (常量), 胸围, 脂臀宽, 体高

Note: d. Predictors: (constant), bust, fat hip width, height

3 讨 论

3.1 体重与体尺表型统计

试验中，所测得的体尺与体重数据相比较于其他新疆地方绵羊品种公羊体尺数据值相对较小，可能与测量时间和羊场本身饲养条件较差有关。试验数据测量于2018年6月，为喀拉也木勒乡牧区牧草刚刚生长的时候，所有羊都刚刚越过冬，加上羊场重视力度不够，饲草料供应不足，导致所有羊膘情和体重都不是很理想。一组数据变异系数的大小直接反映这一指标数据的离散程度。试验中也木勒白羊公羊体重的变异系数为13.08%，脂臀长、脂臀宽和脂臀厚的变异系数为13.21%、11.21%和13.32%，说明也木勒白羊公羊体重和脂尾数据的变化幅度较大，群体生长发育一致性较差。此结果与许鑫等[10]对特克塞尔羊公羊的相关性研究结果和索效军等[11]对麻城黑山羊母羊的相关性研究结果基本一致，体重的变异系数大于体尺指标的变异系数。说明也木勒白羊公羊体重群体增长趋势与其他品种一样，品种间差异不是很明显。试验结果与买买提明巴拉提等[12]对麦盖提羊周岁公羊的研究结果相同，麦盖提羊周岁公羊体重变异系数13.05%，臀长变异系数15.16%，臀宽变异系数24.37%，臀厚变异系数23.20%，体重和脂尾指标变异系数都较大，但试验中也木勒白羊公羊脂臀的变异系数较周岁麦盖提羊要小，也木勒白羊公羊在脂尾发育上较麦盖提羊更加稳定。变异系数值较大也彰显着体重和脂臀在表型选育上有着较高的可选择性。体重与肉羊的产肉性能息息相关，脂臀指标的大小直接表明了肉羊的脂臀的重量。现代人，对肥油的需求不是很高，更多的是喜欢吃低脂肪的瘦肉，而且也木勒白羊脂臀重量有甚者能占胴体重的1/8，加上肥肉的价格十分廉价，所以脂臀的大小可以影响了整只羊的价格和口感，直接影响了牧民的经济收入。对于也木勒白羊肉羊养殖来说，在育种过程中加强对脂臀的选育是十分有必要的。

3.2 体重与体尺相关性

根据也木勒白羊公羊体尺指标与体重的相关性分析结果，也木勒白羊公羊各体尺指标与体重的相关性都达到显著相关水平(P< 0.05)，且除开管围指标，都呈极显著正相关(P< 0.01)。各相关系数中，胸围与体重的相关系数最大，为0.886，胸宽次之，为0.82，相差不大，说明胸围和胸宽对体重的影响作用较明显。此次试验结果与田亚磊等[13]对139只豫西脂尾羊母羊的体重与体尺的相关性分析结果不同，豫西脂尾羊所以体尺指标与体重相关性极显著，但胸深与体重的相关系数最大。研究与韩学平[14]对欧拉型藏羊公羊体重与体尺相关性研究结果也有差别，拉型藏羊公羊管围与体重相关系数最大，尾长和尾宽呈负相关性，而研究中脂臀长和脂臀宽与体重呈极显著正相关，这可能是因为品种差异和年龄不同造成的。试验结果与田亚磊等[15]对太行黑山羊和郭万春等[16]对青海藏体尺及体重性状相关分析结果相同，胸围与体重的相关系数最大，体重与所有体尺指标的相关性极显著，说明在体尺指标相关性上，也木勒白羊公羊与河南小尾寒羊有一定的相似性，胸围是影响体重最主要因素。试验中各体尺指标间也有极显著的相关性，其中胸围与胸深、胸宽和脂臀宽的相关系数较大，依次为0.804、0.795、0.741。

根据相关系数的组成成分关系，可将也木勒白羊公羊体尺指标与体重的相关系数拆分为直接影响体重的直接影响作用和通过影响其它体尺指标间接影响体重的间接影响作用两部分。通径系数的大小就直接反映了体尺指标与体重的直接相关程度，从通径分析的结果可知道也木勒白羊公羊各体尺指标对体重的直接作用程度。胸围(0.34)、胸宽(0.29)、体长(0.17)、脂臀宽(0.16)对体重的直接影响作用最为明显，这4个体尺指标直接影响体重的大小。间接作用最大的是胸深(0.751)，胸深主要通过影响胸围((0.274)、胸宽0.191)等体尺指标间接影响体重。管围的直接作用为-0.01，间接作用为0.341，主要通过影响胸围(0.191)间接影响体重，但直接作用和间接作用都较小，对体重的影响不大。试验中胸围和胸宽的直接影响作用和间接影响作用都较强。此结果与对豫西脂尾羊的研究结果不相同，豫西脂尾羊的体长和胸围的直接作用为负数，主要是通过影响其他体尺指标来间接影响豫西脂尾羊的体重。白俊艳等[17]通过对大尾寒羊的研究结果表明，胸围和体高有较强的直接和间接影响作用，体长直接影响作用较弱，尾宽指标的直接影响作用和间接影响作用都较弱。此结果与研究不相同，也木勒白羊脂臀宽无论是直接作用还是间接作用都较大，这可能是品种差异或者本文只是单一考虑公羊造成的。试验结果与王志明等[18]对甘肃省高山细毛羊体尺与体重相关性分析结果基本相同，体长和胸围对体重的影响有较强的直接作用，而且还有较强的间接作用，间接作用较为明显的是胸深。在对不同品种羊的研究中，胸围基本都有较强的直接作用和间接作用。

3.3 体重与体尺数据的逐步回归分析

研究中影响也木勒白羊公羊体重最主要的体尺指标为胸围、体长、胸宽和脂臀宽。此结果与杨燕等[19]对沂蒙山黑山羊母羊的分析结果不同，管围在研究中没有较明显的影响作用，这可能是所选择羊年龄、性别不同造成的。试验结果田亚磊等[20]对鲁山牛腿山羊的相关性研究结果较为相同，但鲁山牛腿山羊没有考虑尾巴指标对体重的影响，当然，也不需要,也木勒白羊属于脂尾型绵羊，脂臀的重量较大，在研究中不可忽视。

4 结 论

也木勒白羊公羊体尺指标与体重呈显著正相关(P< 0.05),除管围外，均呈极显著正相关(P< 0.01)。胸围和胸宽对也木勒白羊公羊体重的影响有较强的直接影响作用和较强的间接影响作用。筛选出来的影响也木勒白羊公羊体重的主要体尺指标为胸围、胸宽、体长、脂臀宽。也木勒白羊公羊体尺指标与体重的最优线性回归方程为y=-65.113+0.647x3(胸围)+0.518x8(脂臀宽)+1.248x1(胸宽)+0.378x2(体长)。