澳洲白羊體質(zhì)量與體尺指標的相關(guān)及回歸分析

施海娜， 劉雨田， 李世恩， 張金霞， 高 鈺， 梁萬鵬， 朱正生， 陳兆勇， 李 凱， 徐振飛*

(1.慶陽市農(nóng)業(yè)科學(xué)研究院，甘肅 慶陽 745000；2.甘肅中盛農(nóng)牧發(fā)展有限公司，甘肅 慶陽 745000)

“澳洲白”是澳大利亞第1個利用現(xiàn)代基因測定手段培育的品種。該品種集成了白杜泊綿羊、萬瑞綿羊、無解道賽特羊和特克賽爾羊等品種基因，通過對多個品種羊特定肌肉生長基因標記和抗寄生蟲基因標記的選擇(MyoMAX，LoinMAX，WormSTAR)，培育而成的專門用于與杜泊綿羊配套的，屬于粗毛型的中、大型肉羊品種，于2009年10月在澳大利亞注冊。其特點是體型大、生長快、成熟早、全年發(fā)情，有很好的自動換毛能力。在放牧條件下5～6月齡可達到23 kg胴體，舍飼條件下，該品種6月齡胴體重可達26 kg，且脂肪覆蓋均勻，板皮質(zhì)量具佳。澳洲白在我國已多次引入，截至2017年甘肅省已引入發(fā)展數(shù)量達到2 000只以上[1]。此品種使養(yǎng)殖者能夠在各種養(yǎng)殖條件下用作三元配套的終端父本，可以產(chǎn)出在生長速率、個體重量、出肉率和出欄周期短等方面理想的商品羔羊。

家畜體尺數(shù)據(jù)直接反映了家畜的體格大小和體軀結(jié)構(gòu)、發(fā)育等情況，與家畜的生理機能、生產(chǎn)性能、抗病力以及對外界環(huán)境條件的適應(yīng)能力等密切相關(guān)。國內(nèi)外學(xué)者針對有關(guān)澳洲白羊種羊飼養(yǎng)管理技術(shù)和雜交改良地方品種試驗等方面的研究較多[2-11]，但對其體質(zhì)量與體尺指標相關(guān)性及回歸分析研究尚缺乏報道。本文通過對以湖羊為受體采用胚胎移植技術(shù)擴繁的澳洲白羊體質(zhì)量與其體尺指標進行測定，初步掌握澳洲白羊的生長發(fā)育情況，并進行了相關(guān)性分析、通徑分析、逐步回歸分析以及對其表型相關(guān)系數(shù)分解，探究其體尺指標與體質(zhì)量的關(guān)系，建立體質(zhì)量與其體尺指標之間的最優(yōu)回歸方程，以期能為將來的肉羊新品種選育工作提供參考數(shù)據(jù)。

1 材料與方法

1.1 試驗動物

試驗動物澳洲白羊為2017年9月底在甘肅中盛華美羊業(yè)西峰羊場(以湖羊為受體，從天津奧群公司購買的純種澳洲白A級冷凍胚胎)采用胚胎移植技術(shù)所產(chǎn)，年齡為6月齡左右，健康狀況良好。

1.2 測定項目

隨機抽取32只(公母羊各16只)健康羊只對其進行體質(zhì)量(y)、體尺(xi)測定。體質(zhì)量稱量用電子秤，精度為0.01 kg。體尺指標測量參考家畜育種學(xué)[12]，測量工具為測杖、卷尺，測量指標包括體高(x1)、體斜長(x2)、胸圍(x3)、胸深(x4)、胸寬(x5)和管圍(x6)。

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

采用Excel 2007對所測數(shù)據(jù)進行初步整理，用SPSS 19.0統(tǒng)計分析軟件進行相關(guān)性分析、通徑分析、逐步回歸分析，并對相關(guān)系數(shù)進行分解，建立體質(zhì)量與體尺指標之間的最優(yōu)回歸方程。試驗數(shù)據(jù)均以平均數(shù)±標準差表示。

2 結(jié)果與分析

2.1 體質(zhì)量與體尺指標的表型統(tǒng)計量

表1為澳洲白羊體質(zhì)量與體尺指標的表型統(tǒng)計量。由表1可見，體質(zhì)量(y)的變異系數(shù)最大，為10.0%；其次為胸寬(x5)、胸深(x4)、管圍(x6)、胸圍(x3)、體斜長(x2)，其變異系數(shù)分別為7.0%，6.1%，5.3%，4.8%，4.6%；體高(x1)的變異系數(shù)最小，為4.2%。

表1 澳洲白羊體質(zhì)量與體尺指標的表型統(tǒng)計量

2.2 體質(zhì)量與體尺指標的表型相關(guān)性分析

由表2可見，澳洲白羊體質(zhì)量(y)與胸圍(x3)和管圍(x6)相關(guān)性極顯著(P<0.01)，與體高(x1)相關(guān)性顯著(0.010.05)。體質(zhì)量(y)與管圍(x6)的表型相關(guān)系數(shù)最大，為0.680；與胸圍(x3)、體高(x1)、體斜長(x2)、胸深(x4)、胸寬(x5)的表型相關(guān)系數(shù)為0.501，0.393，0.314，0.182，0.154。其中體斜長(x2)與管圍(x6)、胸圍(x3)與胸寬(x5)相關(guān)性極顯著(P<0.01)，其表型相關(guān)系數(shù)分別為0.493，0.450。

表2 澳洲白羊體質(zhì)量與體尺指標的表型相關(guān)系數(shù)

2.3 體質(zhì)量與體尺指標之間的通徑分析

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

P1+0.291P2+0.272P3+0.079P4+0.042P5-0.050P6=0.393

0.291P1+P2+0.177P3+0.006P4+0.208P5+0.493P6=0.314

0.272P1+0.177P2+P3+0.240P4+0.450P5+0.286P6=0.501

0.079P1+0.006P2+0.240P3+P4+0.110P5+0.313P6=0.182

-0.050P1+0.208P2+0.450P3+0.110P4+P5+0.740P6=0.154

0.345P1+0.493P2+0.286P3+0.313P4+0.740P5+P6=0.680

用矩陣法對多元一次方程組進行求解，得到澳洲白羊體尺指標對體質(zhì)量的通徑系數(shù)分別為：P1=0.121；P2=-0.081；P3=0.334；P4=-0.099；P5=-0.008；P6=0.614。

2.4 體尺指標對體質(zhì)量的相關(guān)系數(shù)分解

根據(jù)相關(guān)系數(shù)的組成效應(yīng)，可將澳洲白羊各體尺指標(xi)與體質(zhì)量(y)的表型相關(guān)系數(shù)(皮爾森相關(guān)系數(shù))分解為體尺對體質(zhì)量的直接作用和該體尺指標通過其他體尺性狀對體質(zhì)量的間接作用兩部分。澳洲白羊體尺對體質(zhì)量相關(guān)系數(shù)的分解見表3。

由表3可知，澳洲白羊的管圍(x6)對體質(zhì)量的直接作用最大，為0.614；胸深(x4)對體質(zhì)量的直接作用最小，為-0.099；體斜長(x2)對體重(y)的綜合間接作用最大，為0.394；其次為體高(x1)、胸深(x4)、胸圍(x3)、胸寬(x5)，分別為：0.271，0.209，0.167，0.161；管圍(x6)的綜合間接作用最小，為0.066。各體尺指標對體質(zhì)量的直接作用大小為：管圍(x6)>胸圍(x3)>體高(x1)>胸寬(x5)>體斜長(x2)>胸深(x4)；各體尺指標對體質(zhì)量的間接作用大小為：體斜長(x2)>體高(x1)>胸深(x4)>胸圍(x3)>胸寬(x5)>管圍(x6)?？梢钥闯?，體高(x1)對體質(zhì)量既有較強的直接作用，又有較強的間接作用；管圍(x6)、胸圍(x3)主要是通過直接作用對體質(zhì)量產(chǎn)生影響，對體質(zhì)量的間接作用較弱；體斜長(x2)、胸深(x4)主要是通過間接作用對體質(zhì)量產(chǎn)生影響，對體質(zhì)量的直接作用較弱；胸寬(x5)對體質(zhì)量的直接作用、間接作用均較弱。體高(x1)、體斜長(x2)、胸圍(x3)胸深(x4)主要是通過管圍(x6)對體質(zhì)量產(chǎn)生間接作用，胸寬(x5)和管圍(x6)主要是通過胸圍(x3)對體質(zhì)量產(chǎn)生間接作用。

表3 澳洲白羊體尺對體質(zhì)量相關(guān)系數(shù)的分解

2.5 體尺指標對體質(zhì)量的逐步回歸分析

2.5.1 體質(zhì)量和各體尺指標整體回歸分析模型 從表4～6中可以看出，體高(x1)、體斜長(x2)、胸深(x4)、胸寬(x5)、對體質(zhì)量(y)的影響不顯著，胸圍(x3)、管圍(x6)對體質(zhì)量(y)的影響顯著，因此應(yīng)該剔除自變量體高(x1)、體斜長(x2)、胸深(x4)、胸寬(x5)后，對因變量體質(zhì)量(y)與胸圍(x3)、管圍(x6)進行多元逐步回歸分析。

表4 體質(zhì)量和各體尺指標整體回歸分析模型

表5 體質(zhì)量和各體尺指標整體回歸分析的多元回歸方差分析

表6 體質(zhì)量和各體尺指標總體回歸分析的多元回歸系數(shù)參數(shù)估計值

2.5.2 體質(zhì)量和各體尺指標逐步回歸分析模型 逐步回歸分析即按自變量對因變量作用的顯著程度，從大到小逐個引入回歸方程，即將自變量的t值從大到小逐個引入回歸方程。由表6可知，各個自變量的顯著性由大到小依次為管圍(x6)>胸圍(x3)>體高(x1)(不顯著)>胸深(x4)(不顯著)>體斜長(x2)(不顯著)>胸寬(x5)(不顯著)，因此建立了體質(zhì)量(y)和管圍(x6)、胸圍(x3)之間的回歸分析模型2，結(jié)果見表7。由模型2所得到的方程為體質(zhì)量(y)和體尺指標之間的最優(yōu)回歸方程為y=-39.698+0.355x3+5.064x6(R=0.752，P<0.05，表7)，標準化后的回歸方程為y=0.334x3+0.585x6(表9)，方差分析的F值為18.815(表8)，達到極顯著水平(P<0.01)。

表7 澳洲白羊體質(zhì)量和胸圍、管圍指標逐步回歸分析模型

表8 體質(zhì)量和胸圍、管圍指標逐步回歸分析的多元回歸方差分析

表9 體質(zhì)量和胸圍、管圍指標逐步回歸分析的多元回歸系數(shù)參數(shù)估計值

3 討 論

3.1 體質(zhì)量和體尺指標的表型統(tǒng)計量

變異系數(shù)又稱離散系數(shù)，它反映的是相關(guān)研究變量變化的相對程度指標。澳洲白羊的體質(zhì)量與其體尺指標的變異系數(shù)由大到小依次為：體質(zhì)量(y)>胸寬(x5)>胸深(x4)>管圍(x6)>胸圍(x3)>體斜長(x2)>體高(x1)，體質(zhì)量的變異系數(shù)為10%，說明澳洲白羊的體質(zhì)量變異幅度較大，體質(zhì)量具有較大的選擇與提高潛力。體高的變異系數(shù)最小，為4.2%，說明體高的選擇與提高潛力很小。

3.2 體質(zhì)量與體尺指標的相關(guān)性分析

相關(guān)系數(shù)是用來衡量2個自變量或1個自變量與1個因變量之間相關(guān)程度的指標，由于在數(shù)據(jù)分析過程中筆者沒有消除其他研究變量對這兩個變量的影響，所以表型相關(guān)系數(shù)并不能夠真實地反映2個自變量或自變量與因變量之間的相關(guān)關(guān)系，只能反映所研究的2個自變量或自變量與因變量之間的綜合作用程度。本文通過對澳洲白羊體質(zhì)量與其體尺指標進行相關(guān)性分析，可以看出澳洲白羊體質(zhì)量(y)與胸圍(x3)和管圍(x6)相關(guān)性極顯著(P<0.01)，與體高(x1)相關(guān)性顯著(0.01

3.3 體質(zhì)量與其體尺指標間的通徑分析及其相關(guān)系數(shù)的分解

通徑系數(shù)表示自變量對因變量的直接作用大小，其真實地反映了2個變量間的直接相關(guān)程度，它是標準化的偏回歸系數(shù)[18-19]。通過對體質(zhì)量和各體尺指標之間的通徑分析以及其相關(guān)系數(shù)的分解可以看出：6月齡澳洲白羊體質(zhì)量(y)與胸圍(x3)和管圍(x6)均表現(xiàn)出極顯著相關(guān)性(P<0.01)。從表3中相關(guān)系數(shù)的分解中可以看出，體高(x1)對體質(zhì)量的直接作用和間接作用均較強；管圍(x6)、胸圍(x3)主要是通過直接作用對體質(zhì)量產(chǎn)生影響，對體質(zhì)量的間接作用較弱；體斜長(x2)、胸深(x4)主要是通過間接作用對體質(zhì)量產(chǎn)生影響，對體質(zhì)量的直接作用較弱；胸寬(x5)對體質(zhì)量的直接作用和間接作用均較弱。通過上述分析說明澳洲白羊的生長發(fā)育是一個復(fù)雜的有機整體，各指標間有著種種密切而復(fù)雜的聯(lián)系。影響其體質(zhì)量的因素不但多，而且影響途徑亦不同。所以在選種過程中除了加強對體質(zhì)量影響直接作用強的因素如澳洲白羊的管圍、胸圍選擇的同時，也應(yīng)該關(guān)注其他因素對體質(zhì)量的作用。在選種時，可以參考本實驗數(shù)據(jù)分析結(jié)果，即本實驗所建立的澳洲白羊體質(zhì)量與其體尺指標之間的最優(yōu)回歸方程，進行綜合全面考慮，避免盲目選擇，從而能達到最優(yōu)的選種效果。

3.4 體尺指標對體質(zhì)量的逐步回歸分析

通過對澳洲白羊體質(zhì)量與其體尺指標之間的逐步回歸分析，得到它們之間的最優(yōu)回歸方程：y=-39.698+5.064x6+0.355x3(R=0.752，P<0.05)，標準化后的回歸方程：y=0.334x3+0.585x6，方差分析的F值為18.815，達極顯著水平(P<0.01)。體高(x1)、體斜長(x2)、胸深(x4)和胸寬(x5)回歸系數(shù)顯著程度不高(P值均大于0.05)，所以予以排除，說明這4個體尺指標對澳洲白的體質(zhì)量影響不明顯。這與楊燕等[15]對沂蒙黑山羊、李勝開等[16]對隆林山羊以及王欣榮等[17]對甘南草地型藏羊的研究結(jié)果類似。而他們的回歸方程中除了管圍、胸圍入選外，也入選了體高、體斜長、胸深等因素，這主要是由于羊的品種、飼養(yǎng)管理水平、生活環(huán)境等諸多因素不同而造成的。從方程可知，當胸圍保持不變時，管圍每增加1 cm，澳洲白羊的體質(zhì)量增加0.585 kg；當管圍保持不變時，胸圍每增加1 cm，杜泊羊的體質(zhì)量增加0.334 kg?；貧w方程中R2=0.565，說明該回歸方程的回歸平方和占總變異的56.5%，回歸方程顯著存在，所選體尺指標對體質(zhì)量有較大的決定作用，方程擬合度較高，具有一定的參考價值。剩余因子e=1-R2=0.435，說明對澳洲白羊體質(zhì)量產(chǎn)生影響的自變量不止管圍、胸圍，還有一些其他影響因素沒有考慮到，因此下一步需要對影響澳洲白羊體質(zhì)量的因素(包括澳洲白羊的體尺指標、營養(yǎng)水平、生存環(huán)境)進行全面分析。

4 結(jié) 論

在本試驗條件下，6月齡澳洲白羊管圍(x6)和胸圍(x3)是影響其體質(zhì)量的主要因素，體高(x1)、體斜長(x2)、胸圍(x3)、胸深(x4)、主要是通過管圍(x6)對體質(zhì)量產(chǎn)生間接作用，胸寬(x5)和管圍(x6)，主要是通過胸圍(x3)對體質(zhì)量產(chǎn)生間接作用。其體質(zhì)量和體尺指標之間的最優(yōu)回歸方程y=-39.698+5.064x6+0.355x3(R=0.752，P<0.05)，標準化后的回歸方程：y=0.334x3+0.585x6(F=18.815，P<0.01)。