青海無角歐拉羊羔羊體重與體尺多元回歸分析

莊蕾，閻明毅，余忠祥，張強龍，唐燕花，多杰才讓，唐燕青，李永欽，吳森

(1.青海大學畜牧獸醫(yī)科學院 青海省高原家畜遺傳資源保護與創(chuàng)新利用重點實驗室，西寧810016;2.青海省河南蒙古族自治縣生態(tài)有機畜牧業(yè)技術(shù)服務(wù)中心，河南縣811500;3.青海省河南蒙古族自治縣有機產(chǎn)業(yè)科技園區(qū)管理委員會，河南縣811500;4.青海省河南蒙古族自治縣托葉瑪鄉(xiāng)農(nóng)牧業(yè)技術(shù)服務(wù)站，河南縣811500;5.青海省河南蒙古族自治縣獸醫(yī)站，河南縣811500)

歐拉羊是青藏高原藏系綿羊的一個重要類群，曾被稱為“歐拉型藏羊”，2018年被列入國家畜禽資源，更名為“歐拉羊”，因其主產(chǎn)區(qū)位于歐拉山而得名。歐拉羊具有發(fā)育快、體格結(jié)實、肢體高大、生長快、繁殖性能好、適應(yīng)高寒氣候、耐粗飼、肉質(zhì)好等特點［1，2］，其成年羊體重在60 kg以上，遠大于其他藏羊。另外，歐拉羊主產(chǎn)區(qū)位于“中華水塔”三江源腹地的河南蒙古族自治縣，擁有國家認證的有機牧場，其肉質(zhì)獨特，被世人所關(guān)注［3］。歐拉羊具備較高的開發(fā)利用價值，近年來，在本品種不斷選育的基礎(chǔ)上，現(xiàn)選育出大群無角型歐拉羊，性格溫順，易于管理，且符合現(xiàn)代規(guī)?；⒓s化養(yǎng)殖趨勢及高效畜牧業(yè)生產(chǎn)理念，歐拉羊是良好的現(xiàn)代化藏羊品種育成資源。

羔羊期是綿羊生長的關(guān)鍵時期，容易受代謝印記影響后期發(fā)育。羔羊期生長發(fā)育相關(guān)研究有助于群體早期選種及生長發(fā)育評估，特別是公羔羊。公羔羊的體況評估研究有利于種公羊的早期選留，減少不必要的養(yǎng)殖負擔，實現(xiàn)種群的快速優(yōu)化。本研究對無角型歐拉羊公羔羊的體重與體尺性狀進行了線性回歸分析、通徑分析、逐步回歸分析，建立了最優(yōu)回歸方程，旨在為無角型歐拉羊的早期選種選育、優(yōu)良種公羊培育等工作提供技術(shù)和理論參考。

1 材料與方法

1.1 試驗動物

本研究以無角型歐拉羊為研究對象，測定收集青海省黃南藏族自治州河南蒙古族自治縣優(yōu)干寧鎮(zhèn)經(jīng)多年本品種選育的無角型歐拉羊公羔羊(8～10月齡)653只，所有測定羊均發(fā)育正常、生長健康，實行常規(guī)免疫，自然放牧。

1.2 測定指標與方法

剪毛后對其體重(Y)和主要體尺性狀，即體高(X1)、體長(X2)、胸寬(X3)、胸深(X4)、尻寬(X5)、胸圍(X6)進行測定，生產(chǎn)性能測定主要參照張沅主編《家畜育種學》的第三章第七節(jié)體型鑒定［4］及《綿山羊生產(chǎn)性能測定技術(shù)規(guī)范》行業(yè)標準［5］進行，由固定測量專員對體高、體長、胸寬、胸深、尻寬、胸圍進行測定，每個部位測定3次，取平均值。

1.3 逐步回歸分析

建立最優(yōu)回歸方程的方法很多，常用的是逐步回歸法(Stepwise)［6］，其原理是將自變量逐個引入，對全部的自變量X1，X2，...，Xp，按其對因變量Y影響的大小進行比較，通過F檢驗，選擇偏回歸平方和顯著的變量進入回歸方程。同時，每引入一個新變量后，對前面已引入回歸模型的變量逐個進行檢驗，將經(jīng)檢驗后偏回歸平方和不顯著的變量刪除，以保證所得自變量子集中每一個變量都是顯著的，經(jīng)過多次循環(huán)，直到不能再引入新變量為止。此時的多元線性回歸方程即為最優(yōu)多元線性回歸方程［7，8］。

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

用Excel進行數(shù)據(jù)基本整理，用SPSS 23.0統(tǒng)計軟件進行統(tǒng)計分析，對相關(guān)系數(shù)進行分解，建立體重與各體尺性狀之間的最優(yōu)回歸方程。

2 結(jié)果與分析

2.1 無角型歐拉羊體重與體尺性狀的基本統(tǒng)計量

由表1可知，無角型歐拉羊體高變異系數(shù)是7.87%，為弱變異(CV＜10%)。其余表型均為中等變異(10%＜CV＜100%)，胸寬變異系數(shù)為57.36%，接下來依次是體長43.66%、尻寬36.63%、體重23.93%、胸深12.36%、胸圍11.98%。

表1 無角型歐拉羊體重與體尺性狀的基本統(tǒng)計量

2.2 無角型歐拉羊體重與體尺性狀的相關(guān)性分析

由圖1可知，體重與體高、胸深、胸圍之間存在極顯著正相關(guān)(P＜0.01)，相關(guān)程度由大到小為體高＞胸圍＞胸深;體重與體長、胸寬、尻寬之間存在顯著正相關(guān)(P＜0.05)，相關(guān)程度由大到小為體長=胸寬＞尻寬;其余各性狀兩兩之間也存在正相關(guān)。

圖1 無角型歐拉羊體重與體尺性狀的相關(guān)性分析(Pearson)

2.3 無角型歐拉羊體尺指標和體重相關(guān)系數(shù)的分解

對無角型歐拉羊的各體尺性狀與體重數(shù)據(jù)進行逐步回歸分析。根據(jù)相關(guān)系數(shù)的組成效應(yīng)，可將無角型歐拉羊羔羊的各體尺性狀(Xi)對體重(Y)的相關(guān)系數(shù)劃分為直接作用(Pi)和通過性狀(Xj)的間接作用(∑RijPj)兩部分，即rxiy=Pi+∑RijPj。由表2可知，體高對體重的直接作用最大，為0.421;胸深對體重的間接作用最大，為0.2034。

表2 無角型歐拉羊體高、體長、胸寬、胸深、尻寬、胸圍和體重相關(guān)系數(shù)的分解

2.4 最優(yōu)回歸方程的建立

回歸方程1引入的體尺性狀是體高，回歸方程2引入的體尺性狀是體高和胸圍，回歸方程3引入的體尺性狀是體高、胸圍、胸深。由表3可知，體高、胸圍、胸深被依次引入回歸方程后，其R值逐漸變大，估計標準誤逐漸變小。

表3 無角型歐拉羊體重與體尺的逐步回歸模型與參數(shù)估計值檢驗

由上述分析可知，回歸方程3的決定系數(shù)最大，所以回歸方程3為最優(yōu)回歸方程，復(fù)相關(guān)系數(shù)R為0.616，校正的決定系數(shù)R2為0.376?；貧w方程3中對體重影響最顯著的體尺性狀是體高、胸圍、胸深，各性狀標準偏回歸系數(shù)為0.425、0.253、0.122。經(jīng)T檢驗，體高、胸圍、胸深均達到極顯著水平(P＜0.01)。綜上所述，逐步回歸方程Y=0.425X1+0.253X6+0.122X4為最優(yōu)回歸方程，但Y=0.544X1可作為野外無稱重條件下快速預(yù)測歐拉羊體重的最優(yōu)回歸方程。

3 討論

3.1 體重與體尺的表型分析

體重和體尺性狀是綿羊外貌體況評定的重要指標，不僅反映綿羊品種自身的體型特征，同時與綿羊的生產(chǎn)性能有密切聯(lián)系［9，10］。變異系數(shù)可以將不同單位觀測數(shù)據(jù)的相對變化程度直觀的表現(xiàn)出來，變異系數(shù)低于10%時稱為弱變異，在10%～100%之間稱為中等變異［7］。在本研究中，無角型歐拉羊的體高變異系數(shù)是7.87%，為弱變異;體重及其余體尺性狀變異系數(shù)均處于10%～100%之間，為中等變異，且各性狀變異系數(shù)均較大，說明無角型歐拉羊公羔羊在體重、體尺性狀方面具有較大的選育潛力。

3.2 體重與體尺性狀間的相關(guān)性分析

體重是重要的經(jīng)濟與選種指標，主要受到遺傳和環(huán)境的共同影響，因此在同品種內(nèi)和同環(huán)境下，研究體重與各體尺性狀之間的關(guān)系具有一定的代表性［11］。無角型歐拉羊的體高、胸深、胸圍與體重之間存在極顯著正相關(guān)，體長、胸寬、尻寬與體重之間存在顯著正相關(guān)，且體高、胸深、胸圍對體重均有直接或間接作用，說明這3個體尺性狀對無角型歐拉羊在體重增長方面影響較大。在對無角型歐拉羊體重增長影響最大的3個體尺性狀中，由相關(guān)系數(shù)可知，體高對體重的影響最大，為0.544;其次是胸圍，為0.429;胸深對體重的影響最小，為0.324。體重與體高、胸圍、胸深之間的相關(guān)是由于體重與動物的肌肉、骨骼、內(nèi)臟的生長發(fā)育緊密相關(guān)［12］。在梁小鵬等［13］對新疆也木勒白羊、許鑫等［14］對特克塞爾羊的研究中發(fā)現(xiàn)，不僅胸圍對體重會有較大影響，還有體長等體尺性狀也會對體重有較大影響，結(jié)果與本研究不同，可能是本文測定的無角型歐拉羊的月齡較小和品種差異造成的。白俊艷等［15］對大尾寒羊的研究中發(fā)現(xiàn)，體高和胸圍不僅對體重有較強的直接作用，還輔助其他性狀產(chǎn)生較大的間接作用，與本研究結(jié)果相似。在李江文等［16］對呼倫貝爾短尾羊、海占忠等［17］對加什科羊、楊卓等［18］對青海加什科羊等的研究中也有相似的結(jié)果。因此，今后在無角型歐拉羊的品種選育、提高生產(chǎn)性能等研究中，應(yīng)以體重為主，兼顧體高、胸圍、胸深，以取得較好的選育效果。

3.3 最優(yōu)回歸方程的建立

體尺性狀不僅反映其外貌的特征，同時也可以為體重的預(yù)估提供參考［17］。本研究以體重為因變量，體尺性狀為自變量建立逐步回歸方程:Y=0.425X1+0.253X6+0.122X4，因為體高(X1)、胸深(X4)、胸圍(X6)對體重影響較大，所以被選入回歸方程，擬合優(yōu)度的決定系數(shù)為0.379，校正R2值為0.376，回歸關(guān)系達到極顯著水平(P＜0.01)，表明應(yīng)用該回歸方程能夠較準確地評估無角型歐拉羊的體重變化。

綜上所述，體高、胸深、胸圍是影響無角型歐拉羊公羔羊體重的主要因素，無角型歐拉羊公羔羊的體重-體尺最優(yōu)回歸方程為Y=0.425X1+0.253X6+0.122X4為最優(yōu)回歸方程，Y=0.544X1可作為野外無稱重條件下快速預(yù)測歐拉羊體重的最優(yōu)回歸方程。

4 結(jié)論

無角型歐拉羊公羔羊各體尺性狀兩兩之間存在正相關(guān)，體高、胸深、胸圍與體重之間存在極顯著正相關(guān)，這3個體尺性狀是影響體重的主要因素，因此在對無角型歐拉羊進行選育時，要重點關(guān)注體高、胸深、胸圍3個性狀，以獲得更好的選育結(jié)果。