甘肅高山細毛羊體重與體尺指標的相關性研究

王志明，王婕姝，王 毅，韓向敏，王 華

(1. 甘肅農業大學 動物科學技術學院，甘肅 蘭州 730070; 2. 蘭州市畜牧獸醫研究所，甘肅 蘭州 730050)

(1. College of Animal Science and Technology， Gansu Agricultural University， Lanzhou 730070， China; 2. Research Institute of Animal Husbandry and Veterinary of Lanzhou， Lanzhou 730050， China)

?

甘肅高山細毛羊體重與體尺指標的相關性研究

王志明，王婕姝，王 毅，韓向敏，王 華

(1. 甘肅農業大學 動物科學技術學院，甘肅 蘭州 730070; 2. 蘭州市畜牧獸醫研究所，甘肅 蘭州 730050)

采用SPSS 19.0軟件對甘肅高山細毛羊體重與體尺指標進行相關性分析、通徑分析、逐步回歸分析，并對其相關系數進行分解。通過對體重與體尺指標的逐步回歸分析，得出了體重與體尺指標之間的最優回歸模型。結果表明：體高(X1)、體長(X2)、胸圍(X3)、胸深(X4)均對甘肅高山細毛羊體重產生極顯著影響(P<0.01)，對體重(Y)的相關系數分別為：0.766，0.837，0.834，0.801。體長(X2)、胸圍(X3)指標均通過直接作用和間接作用對體重(Y)產生影響，而體高(X1)、胸深(X4)則主要通過輔助其他性狀對體重(Y)產生間接作用。甘肅高山細毛羊體尺指標(Xi)對體重(Y)的最優回歸方程：Y=-43.851+0.425X2+0.587X3。估測值與實測值之間的相關程度為0.914，決定系數R2為0.836，擬合程度較高。

甘肅高山細毛羊；體重；體尺；相關性分析；通徑分析；逐步回歸分析

(1.CollegeofAnimalScienceandTechnology，GansuAgriculturalUniversity，Lanzhou730070，China; 2.ResearchInstituteofAnimalHusbandryandVeterinaryofLanzhou，Lanzhou730050，China)

甘肅高山細毛羊育成于甘肅皇城綿陽育種試驗場和天祝藏族自治縣內的場、社，屬于毛肉兼用型品種，它具有肉質鮮嫩、膻味輕、產肉和沉積脂肪能力強、對高寒牧區(海拔>2 600 m)適應能力強等特點。1981年被批準為新品種，命名為甘肅高山細毛羊。它是由當地藏羊、蒙古羊為母本，新疆細毛羊、高加索細毛羊為父本，并導入少量莎力斯羊、斯達夫羊和多毛品系(含澳血)羊的血液育成的[1-2]。有關甘肅高山細毛羊生產性能、生長發育及飼養管理等方面的研究較多[3-9]，而對其表型性狀中體重與體尺指標的相關性研究尚缺乏報道。本文通過對甘肅高山細毛羊體重、體尺指標的相關性分析、通徑分析、逐步回歸分析以及對其表型相關系數分解，探究體尺指標與體重的相關關系，包括直接相關關系、間接相關關系，并建立甘肅高山細毛羊體重、體尺間的最優回歸方程，為后續的選種、選育工作提供參考。

1 材料與方法

1.1 試驗動物甘肅高山細毛羊取自甘肅省天祝藏族自治縣牧興婦女牛羊養殖專業合作社，年齡為6月齡左右，且供試羊均為母羊，健康狀況良好。

1.2 測定項目

隨機抽取50只健康羊只對其進行體重(Y)、體尺(Xi)測定。體重稱量用耀華牌XK3190A1型電子秤，精度為0.01 kg。體尺指標測量參考家畜育種學[10]，測量工具為測杖、卷尺、圓形測定器，測量指標包括體高(X1)、體長(X2)、胸圍(X3)、胸深(X4)和管圍(X5)。

1.3 數據處理

采用Excel 2013對所測數據進行初步整理，用SPSS 19.0統計分析軟件進行相關性分析、通徑分析、逐步回歸分析，并對相關系數進行分解，建立體重、體尺間的最優回歸方程。

2 結果與分析

2.1 體重和體尺指標的表型統計量

表1為甘肅高山細毛羊體重和體尺指標的表型統計量，由此可知，體重(Y)的變異系數最大，為24.26%，其次為胸深(X4)、體長(X2)，分別為11.39%，10.55%，三者均大于10%。體高(X1)、胸圍(X3)、管圍(X5)變異系數依次為7.27%，7.33%，6.74%，均小于10%。

表1 甘肅高山細毛羊體重和體尺指標的表型統計量

Table 1 Phenotypic statistic of the body weight and body size in Gansu alpine fine-wool sheep

項目體重(Y)/kg體高(X1)/cm體長(X2)/cm胸圍(X3)/cm胸深(X4)/cm管圍(X5)/cm平均值25.5667.2230.8271.4464.656.89標準差6.2004.8883.2525.2347.3660.465CV/%24.267.2710.557.3311.396.74

注：CV(%)=標準差/平均值×100。

2.2 體重與體尺指標的相關性分析

用SPSS 19.0對體重、體尺指標進行相關性分析，得出體重、體尺間的相關系數如表2。甘肅高山細毛羊體重(Y)與體尺指標中體高(X1)、體長(X2)、胸圍(X3)、胸深(X4)呈極顯著正相關(P<0.01)，與管圍(X5)相關性不顯著。體重(Y)與體長(X2)的表型相關系數最大，為0.837，與胸圍(X3)、胸深(X4)、體高(X1)的表型相關系數依次為0.834，0.801，0.766。除管圍(X5)外，各個體尺指標之間也表現出極強的相關性。其中，體高(X1)與體長(X2)、胸圍(X3)、胸深(X4)之間呈極顯著正相關(P<0.01)，表型相關系數分別為0.633，0.862，0.696；體長(X2)與胸圍(X3)、胸深(X4)呈極顯著正相關(P<0.01)，相關系數分別為0.671，0.739；胸圍(X3)與胸深(X4)之間呈極顯著正相關(P<0.01)，相關系數為0.717。

2.3 體重與體尺指標間的通徑分析

為了進一步揭示體重與體尺指標間的相關關系，對體重、體尺指標進行了通徑分析，用體尺指標與體重的表型相關系數建立以通徑系數為變量的多元一次方程組，方程組如下：

P1+0.633P2+0.862P3+0.696P4+0.150P5=0.766;

0.633P1+P2+0.671P3+0.739P4+0.164P5=0.837;

0.862P1+0.671P2+P3+0.717P4-0.003P5=0.834;

0.696P1+0.739P2+0.717P3+P4+0.138P5=0.801;

0.150P1+0.164P2-0.003P3+0.138P4+P5=0.126。

用矩陣法對多元一次方程組進行求解，得到甘肅高山細毛羊體尺指標對體重的通徑系數分別為：P1=0.021；P2=0.403；P3=0.408；P4=0.191；P5=0.032。

表2 甘肅高山細毛羊體重與體尺指標的表型相關系數

Table 2 Phenotypic correlation between the body weight and body size of Gansu alpine fine-wool sheep

注：*表示相關性顯著(P<0.05)，**表示相關性極顯著(P<0.01)。

2.4 體尺指標對體重相關系數的分解

根據相關系數的組成效應，可將甘肅高山細毛羊各體尺指標(Xi)與體重(Y)的表型相關系數(皮爾森相關系數)分解為體尺對體重的直接作用和該體尺指標通過其他體尺性狀對體重的間接作用兩部分。甘肅高山細毛羊體尺對體重相關系數的分解見表3。

由表3可知，甘肅高山細毛羊體長(X2)對體重(Y)的綜合間接作用最大，為0.676；其次為胸深(X4)、體高(X1)、胸圍(X3)，分別為：0.611，0.502，0.425；管圍(X5)對體重(Y)的綜合間接作用最小，為0.094。各體尺指標對體重的直接作用大小為：胸圍(X3)>體長(X2)>胸深(X4)>管圍(X5)>體高(X1)；各體尺指標對體重的間接作用大小為：體長(X2)>胸深(X4)>體高(X1)>胸圍(X3)>管圍(X5)。可以看出，體長(X2)、胸圍(X3)對體重既有較強的直接作用，又有較強的間接作用；而體高(X1)、胸深(X4)主要是通過間接作用對體重產生影響，對體重的直接作用較弱；管圍(X5)對體重的直接作用、間接作用均較弱。除管圍(X5)外，各體尺指標通過胸圍(X3)對體重均有較強的間接作用。

表3 甘肅高山細毛羊體尺對體重相關系數的分解

Table 3 Subdivision of phenotypic correlations between body size and body weight in Gansu alpine fine-wool sheep

自變量對Y的相關系數直接作用間接作用體高(X1)體長(X2)胸圍(X3)胸深(X4)管圍(X5)合計體高(X1)0.7660.021—0.0130.3520.1330.0040.502體長(X2)0.8370.4040.256—0.2740.1410.0050.676胸圍(X3)0.8340.4080.0180.271—0.137-0.0010.425胸深(X4)0.8010.1910.0150.2990.293—0.0040.611管圍(X5)0.1260.0320.0030.066-0.0010.026—0.094

2.5 體尺指標對體重的逐步回歸分析

運用SPSS 19.0軟件對體重、體尺數據進行多元逐步回歸分析。結果，排除了因變量體高(X1)、胸深(X4)、管圍(X5)，得到體重、體尺指標的最優回歸方程(模型2)(表4，表5)，方程為Y=-43.851+0.587X3+0.425X2(R=0.914，P<0.01)，方差分析的F值為84.204，達到極顯著水平(P<0.01)(表6)。經標準化后的方程為Y=0.505X2+0.495X3(R=0.914，P<0.01)(表5)。

3 討論

3.1 體重和體尺指標的表型統計量

變異系數又稱離散系數，它反映的是標志變動相對程度的指標[11]。體重(Y)、胸深(X4)、體長(X2)三者變異系數均大于10%，且體重最大為24.26%，說明甘肅高山細毛羊體重、胸深、體長變異幅度較大，其中體重指標具有最大的選擇與提高潛力[12]。體高(X1)、胸圍(X3)、管圍(X5)變異系數均小于10%，說明甘肅高山細毛羊體高、胸圍、管圍的變異幅度較小。

表4 體尺指標對體重回歸模型匯總

Table 4 The model of stepwise analysis between the body weight and body size

模型RR2調整R2估計標準誤更改統計量R2F自由度1自由度2顯著性10.8370.7010.6923.4380.70179.7931340.00020.9140.8360.8262.5850.13527.1791330.000

表5 甘肅高山細毛羊多元回歸系數的參數估計值檢驗

Table 5 Multiple regression correlation of Gansu alpine fine-wool sheep

模型B值標準誤差標準系數t值P值B的95.0%置信區間1(常量)-20.0155.134-3.8990.000-30.448-9.582體長0.7050.0790.8378.9330.0000.5440.8652(常量)-43.8515.983-7.3290.000-56.023-31.679體長0.4250.0800.5055.3120.0000.2620.588胸圍0.5870.1130.4955.2130.0000.3580.816

表6 甘肅高山細毛羊多元回歸方差分析

Table 6 ANOVA of multiple regression of rams

模型平方和自由度均方F值P值1943.4011943.40179.7930.00021124.9522562.47684.2040.000

3.2 體重與體尺指標的相關性分析

相關系數是衡量兩個隨機變量間的表型相關程度的指標，由于沒有消除其他變量對這兩個變量的影響，所以表型相關系數不能真實地反映自變量、因變量的相關關系[13]。只能反映該自變量對因變量的綜合作用程度。本文通過對體重、體尺的相關性分析，可以看出體重與體高(X1)、體長(X2)、胸圍(X3)、胸深(X4)均有極顯著的相關性(P<0.01)，其表型相關系數分別為：0.837，0.834，0.801，0.766。說明體尺指標中體高、體長、胸圍、胸深對體重的綜合作用較大，其表現為直接作用和輔助其他性狀對體重的間接作用中的其一或者二者皆有。除管圍(X5)外，體高(X1)、體長(X2)、胸圍(X3)、胸深(X4)相互之間均表現出了極強的正相關，這與白俊艷等[14]對大尾寒羊、陳碧紅等[15]對戴云山羊體重、體尺相關性分析結果一致，說明動物體的體重、體尺指標是相互聯系的，選育中對各個性狀的選擇都需要綜合考慮，以達到最優的選育效果。

3.3 體重與體尺指標間的通徑分析及其相關系數的分解

通徑系數表示自變量對依變量的直接作用大小，其真實地反映了兩個變量間的直接相關程度，它是標準化的偏回歸系數[16]。通過對體重、體尺的通徑分析以及對其相關系數的剖分可以看出，體重與體長(X2)、胸圍(X3)有極強的相關性(P<0.01)，與體高(X1)、胸深(X4)相關性不顯著(P>0.05)。體高(X1)、胸深(X4)雖然在相關性分析時對體重表現出了極強的表型相關性(P<0.01)，但通徑分析中卻對體重相關性不顯著(P>0.05)，由此說明：體高(X1)、胸深(X4)對體重的作用主要是通過輔助其他性狀對體重產生間接影響，而本身對體重的直接作用較弱，這點也可從表4相關系數的分解中看出。體高(X1)、胸深(X4)對體重的間接作用合計為0.502，0.611。體高(X1)主要是通過胸圍(X3)、胸深(X4)對體重產生作用，胸深(X4)主要通過體長(X2)、胸圍(X3)對體重產生作用。在相關性分析中，體長(X2)、胸圍(X3)對體重Y的相關系數r2y=0.837>r3y=0.834，而通徑分析中，體長(X2)、胸圍(X3)對體重Y的通徑系數P2=0.404

3.4 體尺指標對體重的逐步回歸分析

通過對體尺指標與體重的逐步回歸分析，得到體尺對體重的最優回歸方程：Y=-43.851+0.425X2+0.587X3。方程組排除了體高、胸深、管圍指標，因對逐步回歸系數的顯著性檢驗中，體高(X1)P=0.593>0.05，胸深(X4)P=0.090>0.05，管圍(X5)P=0.522>0.05，回歸系數顯著程度不高，所以予以排除。說明體高(X1)、胸深(X4)、管圍(X5)對甘肅高山細毛羊體重的直接影響不明顯。回歸結果入選指標與王欣榮等[17]、徐鐵山等[18]、田亞磊等[19]研究結果類似，而他們的回歸方程中除了體長、胸圍被入選外，也分別入選了體高、胸寬、腰角寬。可能是由于試驗羊的品種、生活環境、飼養管理水平等因素的不同造成的。

回歸方程中R2=0.836，說明該回歸方程的回歸平方和占總變異量的83.6%，回歸方程顯著存在，且所選體尺指標對體重具有較大的決定作用，方程擬合度較高，具有一定的參考價值。剩余因子 e=0.405，說明對甘肅高山細毛羊體重產生影響的自變量不止體長、胸圍，還有一些影響因素沒考慮到，與田亞磊等[19]、王瑞芳等[20]對河南奶山羊、黃淮山羊的分析結果一致。因此，對甘肅高山細毛羊體重產生影響的因素需進行全面分析。

4 結論

(1)甘肅高山細毛羊體尺指標中體長、胸圍對體重既有較強的直接作用，又通過輔助其他性狀對體重產生較強的間接作用。而體高、胸深對體重的作用主要是間接作用；管圍對體重的直接作用和間接作用都不明顯。

(2)甘肅高山細毛羊體尺指標對體重的最優回歸方程為：Y=-43.851+0.425X2+0.587X3(R=0.914，P<0.01)。

[1] 楊正謙. 甘肅高山細毛羊育種工作進展及發展建議[J]. 甘肅農業，2004(11)：81.

[2] 馬乃祥. 甘肅高山細毛羊育種現狀及發展途徑[J]. 甘肅畜牧獸醫，1992，22(1)：1-3.

[3] CHEN K K， YANG B H， GUO J. Correlation analysis between microsatellite marker and procreative traits for mutton-wool strain of Gansu alpine fine wool sheep[J].GuizhouAgriculturalSciences， 2009(6)：116-122.

[4] 高芳山，張永堂. 導入澳美提高甘肅高山細毛羊生產性能研究[J]. 中國草食動物，2004，24(6)：20-21.

[5] WANG W M， HU T X， LI F D， et al. Study on differential expression of ovine RARG in ovary at different estrous cycles[J].ScientiaAgriculturaSinica， 2014， 47(15)：3069-3076.

[6] 孫曉萍，肖西山. 提高甘肅高山細毛羊生產性能的研究[J]. 草與畜雜志，1993(2)：22.

[7] 楊瑞基，魏玉明，王凱，等. 高寒牧區德國肉用美利奴羊、邊區萊斯特羊和甘肅高山細毛羊雜交F1生長發育及肉用性能觀測試驗報告[J]. 中國草食動物，2006，26(5)：32-33.

[8] 王喜軍. 甘肅高山細毛羊羔羊飼養管理關鍵技術[J]. 畜牧獸醫雜志，2014，33(4)：110-111.

[9] 李少斌，王繼卿，成述儒，等. 甘肅高山細毛羊雜交一代生長模型及雜種優勢分析[J]. 草業學報，2011，20(5)：177-182.

[10] 張沅. 家畜育種學[M]. 北京：中國農業出版社，2001.

[11] 李林杰，董正信. 經濟應用統計學[M]. 北京：中國物資出版社，2001.

[12] 高曄，閆海龍，馮平，等. 陜北白絨山羊體尺與體重相關回歸分析[J]. 中國草食動物，2010，30(5)：30-32.

[13] 杜榮騫. 生物統計學[M]. 北京：高等教育出版社，2003.

[14] 白俊艷，龐有志，王永偉. 大尾寒羊體重與體尺的回歸分析[J]. 安徽農業科學，2007，35(15)：4537-4538.

[15] 陳碧紅，劉慶華. 戴云山羊體尺與體重的回歸分析[J]. 中國草食動物，2001，31(5)：34-36.

[16] 明道緒. 通徑分析[M]. 雅安：四川農業大學學報編輯部，1990.

[17] 王欣榮，吳建平，楊聯，等. 甘南草地型藏羊體質量與體尺指標的相關性研究[J]. 甘肅農業大學學報，2011，46(5)：7-11.

[18] 徐鐵山，王東勁，劉小林，等. 海南黑山羊體尺與體重的通徑分析及最優回歸模型的建立[J]. 家畜生態學報，2005，26(1)：49-53.

[19] 田亞磊，高騰云，宋云青，等. 河南奶山羊體尺與體重的相關性分析[J]. 山東農業科學，2009(6)：96-98.

[20] 王瑞芳，龐訓勝，王鋒，等. 黃淮山羊體尺體重性狀的相關及回歸分析[J]. 畜牧與獸醫，2008，40(8)：47-49.

(責任編輯 盧福莊)

Correlation analyses between the body weight and body sizes of Gansu alpine fine-wool sheep

WANG Zhi-ming1， WANG Jie-shu1， WANG Yi2， HAN Xiang-min1，*， WANG Hua1

The phenotypic correlations，path coefficients， direct and indirect correlations between the body weight and several body sizes of Gansu alpine fine-wool sheep were analyzed by SPSS 19.0 software. Then the optimum regression equation between the body weight and body sizes was established. The results indicated that the body height (X1)， body length (X2)， chest circumference (X3) and chest depth (X4) showed significant correlations with the body weight at the 0.01 level. The phenotypic correlations were 0.766， 0.837， 0.834， 0.801， respectively. The body length (X2) and chest circumference (X3) had significant influences on the body weight directly and indirectly. However， the body height (X1) and chest depth (X4) showed effects on the body weight indirectly by assisting other characters. The optimal regression model of the body sizes (Xi) and body weight (Y) was established by the stepwise regression asY=-43.851+0.425X2+0.587X3. (R2=0.836，P<0.01).

Gansu alpine fine-wool sheep; body weight; body size; correlation analysis; path analysis; stepwise analysis

http://www.zjnyxb.cn

10.3969/j.issn.1004-1524.2016.01.05

2015-05-14

甘肅省農業科技創新項目(GNCX—2013-34)

王志明(1991—)，男，碩士研究生，研究方向為反芻動物營養及動物生產系統與工程。E-mail：282046956@qq.com

*通信作者，韓向敏，E-mail：min65@sina.com

S826.8

A

1004-1524(2016)01-0028-05

浙江農業學報ActaAgriculturaeZhejiangensis， 2016，28(1):28-32

王志明，王婕姝，王毅，等. 甘肅高山細毛羊體重與體尺指標的相關性研究[J].浙江農業學報，2016，28(1):28-32.