甘肅高山細(xì)毛羊體重與體尺指標(biāo)的相關(guān)性研究

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

(1. 甘肅農(nóng)業(yè)大學(xué) 動物科學(xué)技術(shù)學(xué)院，甘肅 蘭州 730070; 2. 蘭州市畜牧獸醫(yī)研究所，甘肅 蘭州 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)

?

甘肅高山細(xì)毛羊體重與體尺指標(biāo)的相關(guān)性研究

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

(1. 甘肅農(nóng)業(yè)大學(xué) 動物科學(xué)技術(shù)學(xué)院，甘肅 蘭州 730070; 2. 蘭州市畜牧獸醫(yī)研究所，甘肅 蘭州 730050)

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

甘肅高山細(xì)毛羊；體重；體尺；相關(guān)性分析；通徑分析；逐步回歸分析

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

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

1 材料與方法

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

1.2 測定項目

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

1.3 數(shù)據(jù)處理

采用Excel 2013對所測數(shù)據(jù)進行初步整理，用SPSS 19.0統(tǒng)計分析軟件進行相關(guān)性分析、通徑分析、逐步回歸分析，并對相關(guān)系數(shù)進行分解，建立體重、體尺間的最優(yōu)回歸方程。

2 結(jié)果與分析

2.1 體重和體尺指標(biāo)的表型統(tǒng)計量

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

表1 甘肅高山細(xì)毛羊體重和體尺指標(biāo)的表型統(tǒng)計量

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標(biāo)準(zhǔn)差6.2004.8883.2525.2347.3660.465CV/%24.267.2710.557.3311.396.74

注：CV(%)=標(biāo)準(zhǔn)差/平均值×100。

2.2 體重與體尺指標(biāo)的相關(guān)性分析

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

2.3 體重與體尺指標(biāo)間的通徑分析

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

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。

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

表2 甘肅高山細(xì)毛羊體重與體尺指標(biāo)的表型相關(guān)系數(shù)

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

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

2.4 體尺指標(biāo)對體重相關(guān)系數(shù)的分解

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

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

表3 甘肅高山細(xì)毛羊體尺對體重相關(guān)系數(shù)的分解

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

自變量對Y的相關(guān)系數(shù)直接作用間接作用體高(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 體尺指標(biāo)對體重的逐步回歸分析

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

3 討論

3.1 體重和體尺指標(biāo)的表型統(tǒng)計量

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

表4 體尺指標(biāo)對體重回歸模型匯總

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

模型RR2調(diào)整R2估計標(biāo)準(zhǔn)誤更改統(tǒng)計量R2F自由度1自由度2顯著性10.8370.7010.6923.4380.70179.7931340.00020.9140.8360.8262.5850.13527.1791330.000

表5 甘肅高山細(xì)毛羊多元回歸系數(shù)的參數(shù)估計值檢驗

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

模型B值標(biāo)準(zhǔn)誤差標(biāo)準(zhǔn)系數(shù)t值P值B的95.0%置信區(qū)間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 甘肅高山細(xì)毛羊多元回歸方差分析

Table 6 ANOVA of multiple regression of rams

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

3.2 體重與體尺指標(biāo)的相關(guān)性分析

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

3.3 體重與體尺指標(biāo)間的通徑分析及其相關(guān)系數(shù)的分解

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

3.4 體尺指標(biāo)對體重的逐步回歸分析

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

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

4 結(jié)論

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

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

[1] 楊正謙. 甘肅高山細(xì)毛羊育種工作進展及發(fā)展建議[J]. 甘肅農(nóng)業(yè)，2004(11)：81.

[2] 馬乃祥. 甘肅高山細(xì)毛羊育種現(xiàn)狀及發(fā)展途徑[J]. 甘肅畜牧獸醫(yī)，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] 高芳山，張永堂. 導(dǎo)入澳美提高甘肅高山細(xì)毛羊生產(chǎn)性能研究[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] 孫曉萍，肖西山. 提高甘肅高山細(xì)毛羊生產(chǎn)性能的研究[J]. 草與畜雜志，1993(2)：22.

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

[8] 王喜軍. 甘肅高山細(xì)毛羊羔羊飼養(yǎng)管理關(guān)鍵技術(shù)[J]. 畜牧獸醫(yī)雜志，2014，33(4)：110-111.

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

[10] 張沅. 家畜育種學(xué)[M]. 北京：中國農(nóng)業(yè)出版社，2001.

[11] 李林杰，董正信. 經(jīng)濟應(yīng)用統(tǒng)計學(xué)[M]. 北京：中國物資出版社，2001.

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

[13] 杜榮騫. 生物統(tǒng)計學(xué)[M]. 北京：高等教育出版社，2003.

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

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

[16] 明道緒. 通徑分析[M]. 雅安：四川農(nóng)業(yè)大學(xué)學(xué)報編輯部，1990.

[17] 王欣榮，吳建平，楊聯(lián)，等. 甘南草地型藏羊體質(zhì)量與體尺指標(biāo)的相關(guān)性研究[J]. 甘肅農(nóng)業(yè)大學(xué)學(xué)報，2011，46(5)：7-11.

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

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

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

(責(zé)任編輯 盧福莊)

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

甘肅省農(nóng)業(yè)科技創(chuàng)新項目(GNCX—2013-34)

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

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

S826.8

A

1004-1524(2016)01-0028-05

浙江農(nóng)業(yè)學(xué)報ActaAgriculturaeZhejiangensis， 2016，28(1):28-32

王志明，王婕姝，王毅，等. 甘肅高山細(xì)毛羊體重與體尺指標(biāo)的相關(guān)性研究[J].浙江農(nóng)業(yè)學(xué)報，2016，28(1):28-32.