燕山绒山羊体尺、体重、绒厚与产绒量的相关回归分析

葛丛丛，刘华云，李寒妹，付艳芳，付志新，刘铮铸，巩元芳*

(1.河北科技师范学院 动物科技学院，河北省特色动物种质资源挖掘与创新重点实验室，河北 秦皇岛 066004；2.河北省畜牧总站，河北 石家庄 050035)

家畜的体重和体尺指标与其生产性能、繁殖性能等密切相关，并随生态环境和选育方式的变更而发生改变[1]。绒山羊是经过长期自然选择和人工选育而形成的一类绒毛用山羊品种，因其采食消化能力强，纤维消化率高，可在贫瘠的山区放牧下生存[2]，受到养殖者的喜爱。我国绒山羊主要分布在华北、西北、东北及青藏高原等地区，主要包括辽宁绒山羊、内蒙古绒山羊和河西绒山羊等品种[3]。体重、产绒量是直接影响绒山羊经济效益的重要指标[4-5]。燕山绒山羊主产区在河北燕山山区，多在秦皇岛的青龙县、承德的宽城县等地饲养[6]，是当地农民以本地羊为基础，经过数十年不太系统的选育和杂交，形成的一个绒肉兼用型的绒山羊新类群[7]，截止到目前为止，成年羊群体内仍存在体型大小、产绒量及绒品质等性能参差不齐的问题，严重影响其经济效益。因此，加快燕山山区优良绒山羊的选育，提高产绒量和羊绒品质势在必行。本研究拟对燕山绒山羊的体尺指标、体重、绒厚与产绒量等指标进行测定与分析，建立最优回归方程，旨在找出燕山绒山羊各体尺指标、体重、绒厚与产绒量等的相互关系，为燕山绒山羊的选育和进一步开发、利用提供参考依据。

1 材料与方法

1.1 试验材料

2019年3～4月份，在河北省秦皇岛市青龙满族自治县某绒山羊养殖场选取年龄在1～5岁健康无病的燕山绒山羊母羊162只，采集体尺、抓绒后体重、绒厚和产绒量等相关数据。

1.2 测定项目

利用皮尺、钢尺、测杖、电子秤等工具对受试山羊的产绒量(Y)、体长(X1)、体高(X2)、胸围(X3)、管围(X4)、绒厚(X5)和抓绒后的体重(X6)进行了测定，具体操作方法见文献[8]。

1.3 数据的统计与分析

利用Excel和SPSS20.0软件进行简单归纳处理，剔除无效数据后，对各指标进行相关系数、回归方程、通径系数等参数的分析。

2 结果与分析

2.1 燕山绒山羊各指标的描述性统计

利用SPSS20.0软件对所得燕山绒山羊体尺、抓绒后体重、绒厚和产绒量的数据进行了描述统计，结果见表1。由表1可知，燕山绒山羊体长和体高的变异系数较小，为弱变异(CV<10%)[9]，而胸围、管围、绒厚、抓绒后体重和产绒量为中等变异(10%

表1 燕山绒山羊各体尺指标、抓绒后体重、绒厚与产绒量的描述性统计Table 1 Descriptive statistics of body size,body weight after fleece,cashmere thickness and cashmere yield of Yanshan Cashmere Goat

2.2 燕山绒山羊各指标的相关分析

相关分析可以确定两性状之间的相关程度，用相关系数r来表示，r为正且值越接近1，表示两性状之间的正相关程度越强。利用SPSS20.0软件统计分析了燕山绒山羊各体尺指标、抓绒后体重、绒厚与产绒量之间的相关程度，获得的相关系数见表2。

由表2可知，体长、体高、胸围、绒厚、抓绒后体重与产绒量均呈极显著正相关(P<0.01)，其中绒厚与产绒量的相关系数最大，为0.728；体长、胸围、抓绒后体重与产绒量相关系数较大，分别为0.449、0.639和0.558；体高与产绒量相关系数较小，为0.274。管围与产绒量相关不显著(P>0.05)，相关系数仅为0.071。体高、管围与绒厚相关性不显著(P>0.05)，其余各性状之间均为极显著正相关(P<0.01)。

表2 燕山绒山羊体尺、抓绒后体重、绒厚和产绒量的相关分析结果Table 2 Correlation analysis results of body size,body weight after fleece,cashmere thickness and cashmere yield of Yanshan Cashmere Goat

2.3 燕山绒山羊各指标的通径分析

根据以上各指标相关系数的组成效应，可将燕山绒山羊各体尺指标、抓绒后体重、绒厚中某一性状(Xi)与产绒量(Y)的相关系数剖分为某一性状(Xi)对产绒量的直接影响(Piy)和某一性状(Xi)通过另一性状(Xj)对产绒量的间接影响(∑rijPj)两部分，即相关系数riy=Xi与Y的直接通径系数(Piy)+所有Xi与Y的间接通径系数(∑rijPj)，Xi与Y的间接通径系数=相关系数(rij)×通径系数(Pjy)[10]。利用SPSS20.0软件对各指标与产绒量之间的关系进行建模，获得相关系数和通径系数，所得结果见表3。

由表3各指标相关系数的分析可知，体长、体高、胸围、管围、抓绒后体重、绒厚对产绒量的直接影响中，绒厚对产绒量的直接作用最大，为0.547，其次是胸围、体长对产绒量的直接作用，分别为0.194和0.138，而体高、管围、抓绒后体重对产绒量的直接作用均较小。间接影响中，抓绒后体重对产绒量的间接作用最大，为0.461，其次是胸围，为0.447，管围对产绒量的间接作用最小，为0.150。综上，绒厚、胸围、体长主要通过直接作用影响产绒量，体高、管围、抓绒后体重主要通过胸围对产绒量产生间接影响，其中抓绒后体重产生的间接影响最大。

表3 燕山绒山羊各体尺指标、抓绒后体重、绒厚与产绒量相关系数的分析Table 3 Analysis of correlation coefficients of body size,body weight after fleece,cashmere thickness and cashmere yield of Yanshan Cashmere Goat

2.4 燕山绒山羊各指标与产绒量的回归分析

利用逐步法进行回归分析，得到燕山绒山羊各体尺指标、抓绒后体重、绒厚与产绒量关系的回归方程。根据显著性检验，选择回归关系极显著(P<0.01)的最优回归方程为：Y=-598.036+4.527X1+4.81X3+66.282X5，将此回归方程标准化后为：Y=0.161X1+0.244X3+0.555X5。

标准化回归方程内部各指标经显著性分析可知，胸围(X3)、绒厚(X5)与产绒量(Y)的回归关系呈极显著水平(P<0.01)，体长(X1)与产绒量(Y)的回归关系呈显著水平(P<0.05)，说明这些性状对产绒量的关系是真实存在的。

3 讨 论

绒山羊的产绒量和体重是衡量其经济效益的重要指标，也是养殖者选育的首要指标，二者受遗传、营养和环境等多方面的影响[11]。在选育过程中，既要考虑产绒量和体重等重要经济性状，更要考虑各体尺指标对产绒量的影响，达到精确选育的目的。《中国畜禽遗传资源志-羊志》记录辽宁绒山羊母羊的产绒量是641 g[12]，明显高于本研究所得燕山绒山羊的产绒量(548.58 g)，说明燕山绒山羊的产绒量仍有提高的空间。另外，本研究测得燕山绒山羊抓绒后体重和产绒量的变异系数分别为25.20 %和33.30 %，说明采集的个体间指标数据变动幅度较大，群体内绒山羊个体之间生长发育差异较大。其中，产绒量的变异系数大于各体尺指标的变异系数，这与依明·苏来曼等[13]对青格里绒山羊的研究和李景玉等[14]对辽宁绒山羊的研究结果基本一致，说明燕山绒山羊产绒量的群体增长趋势与其他绒山羊一致。燕山绒山羊作为一种绒肉兼用型品种，抓绒后体重、产绒量变异系数较大，也为其群体内选育提供可能。

通过相关分析和通径分析发现，燕山绒山羊的体长、体高、胸围、绒厚、抓绒后体重与产绒量均呈极显著正相关(P<0.01)，这与前人对陕北绒山羊[15-16]和辽宁绒山羊[17]的研究结果一致。本研究发现,燕山绒山羊绒厚(0.547)、胸围(0.194)、体长(0.138)主要通过直接作用对产绒量产生影响；间接作用最大的是抓绒后体重(0.461)，抓绒后体重主要通过影响胸围(0.158)和绒厚(0.222)等指标间接影响产绒量。管围对产绒量的直接作用和间接作用均较小。综上，燕山绒山羊的体长、胸围和绒厚对产绒量的直接作用和间接作用都较强。虽然尚未有研究证实绒厚是产绒量的第一影响因子，但胸围对产绒量的直接作用在藏山羊[18-19]、青格里绒山羊等[13]多个绒山羊品种上均得到证实。说明在选育高产绒量的燕山绒山羊时，除了参考绒厚这一指标外，胸围也应作为一个重要的参考指标。

由回归分析可知，虽然体高、管围、抓绒后体重也通过绒厚、胸围和体长对产绒量产生了间接影响，但在建立最优线性回归方程时发现，仅绒厚、胸围与产绒量的回归关系呈极显著水平(P<0.01)，体长与产绒量的回归关系呈显著水平(P<0.05)，体高、管围、抓绒后体重与产绒量关系不显著(P>0.05)，因此并未计入回归方程中。体高未计入回归方程，与贺晓宏等[17]对辽宁绒山羊的研究结果一致；抓绒后体重、管围未计入回归方程，与李景玉等[14]和闫海龙等[20]对辽宁绒山羊和陕北白绒山羊上的研究结果不一致，绒山羊的体尺指标只有胸围和体长进入了回归方程，导致这一结果可能是绒山羊品种之间存在差异，或者是与各指标测量的样本量不足够多有关。后续将继续扩大样本量，多方面寻找影响燕山绒山羊经济性状的影响因素，加快其选育的进程。

4 结 论

本研究结果表明，燕山绒山羊体长、体高、胸围、绒厚和抓绒后体重与产绒量均呈极显著正相关。体长、胸围和绒厚对燕山绒山羊的直接作用和间接作用都较强。燕山绒山羊体尺与产绒量的最优线性回归方程标准化后为Y=0.161X1(体长)+0.244X3(胸围)+0.555X5(绒厚)。