燕山絨山羊體尺、體重、絨厚與產絨量的相關回歸分析

葛叢叢，劉華云，李寒妹，付艷芳，付志新，劉錚鑄，鞏元芳*

(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(絨厚)。