烏骨山羊體尺體重性狀的相關(guān)及回歸分析

索效軍，張 年，楊前平，陶 虎，熊 琪，李曉鋒，陳明新

（1.湖北省農(nóng)業(yè)科學(xué)院畜牧獸醫(yī)研究所，武漢 430064；2.動(dòng)物胚胎工程及分子育種湖北省重點(diǎn)實(shí)驗(yàn)室，武漢 430064）

烏骨山羊體尺體重性狀的相關(guān)及回歸分析

索效軍1，2，張 年1，楊前平1，陶 虎1，熊 琪1，李曉鋒1，陳明新1

（1.湖北省農(nóng)業(yè)科學(xué)院畜牧獸醫(yī)研究所，武漢 430064；2.動(dòng)物胚胎工程及分子育種湖北省重點(diǎn)實(shí)驗(yàn)室，武漢 430064）

為烏骨山羊品種選擇和育種提供理論依據(jù)，利用SPSS軟件對(duì)64只烏骨山羊成年母羊體尺體重性狀進(jìn)行相關(guān)分析，并采用逐步回歸法，建立最優(yōu)回歸方程。結(jié)果表明：烏骨山羊體高（X1）、體長(zhǎng)（X2）、胸圍（X3）、管圍（X4）與體重（Y）均呈極顯著正相關(guān)（P＜0.01），胸圍、體長(zhǎng)及體高的直接作用和間接作用對(duì)體重影響極大（P＜0.01）；烏骨山羊體尺指標(biāo)（Xi）對(duì)體重（Y）的最優(yōu)回歸方程為Y=-69.923 19 +0.326 27X1+0.138 46X2+0.961 72X3。

烏骨山羊；體重；體尺；相關(guān)性分析；逐步回歸分析

通山烏骨山羊是在獨(dú)特的自然生態(tài)環(huán)境條件下，經(jīng)過(guò)長(zhǎng)期的閉鎖繁育而形成，體型與一般山羊基本相同，毛色以黑色、花色、黃色、黑白花色為主，唇、眼瞼、舌、粘膜、肌肉、骨骼、內(nèi)臟呈烏黑色，是集科研、藥用、食用于一身的珍稀肉羊品種資源。體重是烏骨山羊最重要的遺傳改良性狀之一，體尺不僅反映烏骨山羊的體格大小、體軀結(jié)構(gòu)和器官發(fā)育，且與肉羊的繁育、抗逆性等有一定的相關(guān)性［1］，但目前尚未有烏骨山羊體重與體尺關(guān)系的研究報(bào)道。本研究利用相關(guān)及回歸分析等方法，分析了烏骨山羊體重與體尺性狀間的內(nèi)在關(guān)系，旨在為烏骨山羊的選育和改良提供科技支撐。

1 材料與方法

1.1 材料

烏骨山羊選自通山阿淼養(yǎng)殖專業(yè)合作社，供試羊均為成年母羊，健康狀況良好。

1.2 測(cè)定指標(biāo)與方法

測(cè)定指標(biāo)為體高、體長(zhǎng)、胸圍、管圍和體重，按文獻(xiàn)提供方法測(cè)定［2］。

運(yùn)用SPSS 19.0軟件分析各體尺性狀對(duì)體重的直接作用和間接作用，對(duì)體重進(jìn)行正態(tài)性檢驗(yàn)，利用逐步回歸法建立體尺性狀對(duì)體重的最優(yōu)回歸方程，并對(duì)方程進(jìn)行顯著性檢驗(yàn)［3-5］。

2 結(jié)果與分析

2.1 體重與體尺

由表1可以看出，烏骨山羊體重和體長(zhǎng)的變異系數(shù)較大，其中體重的變異系數(shù)最大（22.13%），說(shuō)明烏骨山羊體重有較大的選擇潛力。

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

烏骨山羊體尺性狀指標(biāo)作為自變量，體重作為依變量，相關(guān)系數(shù)見(jiàn)表2，可以看出，烏骨山羊體重與體尺均呈極顯著正相關(guān)（P＜0.01），相關(guān)程度最高的是胸圍（r= 0.870 60），其次是管圍、體長(zhǎng)（r=0.571 38，r=0.517 27）；各體尺性狀間（除體高與體長(zhǎng)外）均呈較強(qiáng)正相關(guān)（P＜0.05），其中胸圍與管圍的相關(guān)程度最大（r=0.537 37）。

表1 體重與體尺性狀的平均數(shù)、標(biāo)準(zhǔn)差及變異系數(shù)

表2 體重與體尺性狀的相關(guān)系數(shù)

2.3 體重與體尺性狀間相關(guān)系數(shù)的分解

根據(jù)相關(guān)系數(shù)的組成效應(yīng)，可將相關(guān)系數(shù)（rij）分解為體尺對(duì)體重的直接作用（Pi）和某一體尺性狀（Xi）通過(guò)其他體尺性狀（Xj）對(duì)體重（Y）的間接作用（∑rijPj）兩部分［6-7］，即rij=Pi+∑rijPj。從表3可知，胸圍對(duì)體重的直接影響作用最大，為0.754 86；而體長(zhǎng)對(duì)體重的間接作用最大，為0.355 81。

表3 簡(jiǎn)單相關(guān)系數(shù)的分解

2.4 體重正態(tài)性檢驗(yàn)

3.1.2 暫時(shí)大比分領(lǐng)先。調(diào)查中發(fā)現(xiàn)，有4個(gè)隊(duì)員認(rèn)為比賽的勝利已經(jīng)穩(wěn)握在自己手中，對(duì)現(xiàn)在的領(lǐng)先已經(jīng)滿足，因此比賽緊張的心理狀態(tài)放松了下來(lái)，慢慢不敢進(jìn)攻，害怕失誤，對(duì)教練在賽前布置戰(zhàn)術(shù)和打法拋于腦后，沒(méi)有繼續(xù)發(fā)揮比賽的后半程優(yōu)勢(shì)，失誤越來(lái)越多。而另外6名運(yùn)動(dòng)員發(fā)揮比較穩(wěn)定，對(duì)這種短暫領(lǐng)先沒(méi)有放在心上，抓住對(duì)方進(jìn)攻和防守的弱點(diǎn)，發(fā)揮自己的戰(zhàn)術(shù)，繼續(xù)拉大領(lǐng)先優(yōu)勢(shì)，最后取得比賽勝利。（見(jiàn)表2）

對(duì)體重進(jìn)行正態(tài)性檢驗(yàn)，輸出結(jié)果見(jiàn)表4，經(jīng)檢驗(yàn)在Kolmogorov-Smirnov Test和Shapiro-Wilk Test狀態(tài)下的P值分別為0.200 00和0.223 92，均大于0.05，因此，依變量（體重）服從正態(tài)分布，即體重是正態(tài)變量，可以進(jìn)行回歸分析。

表4 正態(tài)性檢驗(yàn)輸出結(jié)果

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

2.5.1 體尺指標(biāo)對(duì)體重回歸方程 依據(jù)體尺性狀對(duì)體重貢獻(xiàn)率大小和標(biāo)準(zhǔn)回歸系數(shù)的顯著性，通過(guò)逐步回歸分析，剔除對(duì)體重影響不顯著的體尺性狀，建立最優(yōu)的多元線性回歸方程。本研究中隨著體尺性狀被逐步引入，回歸方程的復(fù)相關(guān)系數(shù)和決定系數(shù)在逐漸增大，說(shuō)明引入的體尺性狀對(duì)體重的貢獻(xiàn)率在增加。詳見(jiàn)表5。

表5 體尺指標(biāo)對(duì)體重回歸模型

2.5.2 方差分析 由表6可知，自變量X3引入回歸方程后，F(xiàn)=194.147 28；自變量X3、X2引入回歸方程后，F(xiàn)= 107.979 66；自變量 X3、X2、X1引入回歸方程后，F(xiàn)= 79.817 04?；貧w方程1、2、3方差分析均達(dá)到極顯著水平（P<0.01），而X4的偏回歸系數(shù)無(wú)統(tǒng)計(jì)學(xué)意義，即對(duì)回歸方程的影響不大，故未被引入回歸方程。

2.5.3 偏回歸系數(shù)及其t檢驗(yàn) 各步引入對(duì)回歸方程影響最大的變量后，有關(guān)的偏回歸系數(shù)及t檢驗(yàn)結(jié)果見(jiàn)表7，回歸方程3中自變量X3、X2、X1所對(duì)應(yīng)的P值均小于0.05，而沒(méi)有引入方程的 X4在回歸方程中其P值均大于0.05，無(wú)顯著統(tǒng)計(jì)學(xué)意義，故為不重要變量。

表6 烏骨山羊多元回歸方差分析

表7 體重與體尺的多元回歸系數(shù)的參數(shù)估計(jì)值檢驗(yàn)

3 討論

3.1 體重與體尺的變異系數(shù)

體重與體尺是全面、綜合分析羊生長(zhǎng)發(fā)育的重要指標(biāo)，對(duì)肉羊育種工作及生產(chǎn)實(shí)踐有重要作用。變異系數(shù)是標(biāo)志變動(dòng)度相對(duì)程度的指標(biāo)，又稱離散系數(shù)。在生產(chǎn)實(shí)際中，體重可能會(huì)由于各種原因變化較大，因此可利用體尺指標(biāo)進(jìn)行肉羊的品種選育。本研究中體重（Y）、體長(zhǎng)（X2）二者變異系數(shù)均大于10%，這表明烏骨山羊體重、體長(zhǎng)變異幅度較大；而體高（X1）、胸圍（X3）、管圍（X4）變異系數(shù)小于10%，說(shuō)明烏骨山羊體高、胸圍、管圍的變動(dòng)程度相對(duì)較小。

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

相關(guān)系數(shù)可反映體重與體尺兩個(gè)隨機(jī)變量間的表型相關(guān)程度，由于沒(méi)有排除其他變量的影響，因此并不能反映各體尺性狀與體重的直接關(guān)系。本研究中胸圍與體重的相關(guān)程度最高（r=0.870 60），其次是管圍（r= 0.571 38）、體長(zhǎng)（r=0.517 27）；除體高與體長(zhǎng)外，各體尺性狀間亦呈較強(qiáng)正相關(guān)（P＜0.05），其中胸圍與管圍的相關(guān)系數(shù)最大（r=0.537 37）。體重與胸圍之間的相關(guān)是由于家畜體重與骨骼、肌肉、內(nèi)臟的生長(zhǎng)發(fā)育密切相關(guān)，同時(shí)體重與體長(zhǎng)、管圍的相關(guān)也是與家畜骨骼、肌肉的生長(zhǎng)發(fā)育密不可分。前人已將體尺性狀與體重性狀的相關(guān)關(guān)系成功應(yīng)用于羊、雞等品種的改良和新品系選育中［8-10］。因此，今后在烏骨山羊的品種選育方面應(yīng)以胸圍為主并兼顧體長(zhǎng)、管圍。

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

在分析一個(gè)目標(biāo)性狀（依變量）和若干個(gè)不同性狀（自變量）之間相關(guān)關(guān)系時(shí)，常采用多元回歸分析方法［11］。由相關(guān)系數(shù)剖分可知，烏骨山羊體重（Y）和胸圍（X3）的相關(guān)主要是由兩者直接作用而形成的，此結(jié)果與買買提明巴拉提等［12］、徐鐵山等［13］研究結(jié)果一致，與王同軍等［14］的研究結(jié)果不一致，可能是品種差異所致；體高（X1）、體長(zhǎng)（X2）對(duì)體重的影響直接效應(yīng)較小，主要通過(guò)間接作用產(chǎn)生，這表明肉羊是一個(gè)交互作用的有機(jī)整體，其決定體重的因素較多，影響途徑也有所不同，不僅有單一體尺性狀的影響，也存在多個(gè)體尺性狀的同時(shí)作用。因此，在實(shí)際生產(chǎn)中在加強(qiáng)對(duì)胸圍選擇的同時(shí)，也應(yīng)注意關(guān)注間接作用的選擇，進(jìn)行綜合考慮，從而有一定的目的性，避免盲目選擇。

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

通過(guò)對(duì)體尺指標(biāo)與體重的逐步回歸分析，得到體尺對(duì)體重的最優(yōu)回歸方程：Y=-69.923 19+0.326 27X1+ 0.138 46X2+0.961 72X3，入選回歸方程的體尺性狀胸圍（X3）、體長(zhǎng)（X2）和體高（X1）回歸關(guān)系都達(dá)到顯著水平（P<0.05）。本研究中胸圍、體長(zhǎng)是影響體重的主要因素，這一結(jié)果與高志英等［15］的研究結(jié)果相吻合。由線性回歸分析得出胸圍、體長(zhǎng)和體高每增加1 cm，體重將分別增加0.961 72、0.138 46和0.326 27 kg，說(shuō)明這3個(gè)體尺指標(biāo)對(duì)烏骨山羊體重增長(zhǎng)影響比較大。本回歸方程擬合度的決定系數(shù)（R2）值為0.799 63，校正R2為0.789 62，剩余因子e=0.600 49，該值較大，說(shuō)明還有對(duì)體重有影響的因素尚未考慮進(jìn)去，有待進(jìn)一步加大樣本進(jìn)行研究，以完善烏骨山羊體尺、體重的最優(yōu)回歸方程。

4 結(jié)論

胸圍、體長(zhǎng)和體高是影響烏骨山羊母羊體重的主要因素，烏骨山羊成年母羊體尺與體重的最優(yōu)回歸方程為Y=-69.923 19+0.326 27X1+0.138 46X2+0.961 72X3。

［1］ Semakula J，Mutetikka D，Kugonza R D，et a1.Variability in body morphometric measurements and their application in predicting live body weightof Mubende and Small East African goatbreeds in Uganda［J］.Middle-East Journal of Scientific Research，2010，5（2）：98-105．

［2］ 陳偉生，徐桂芳.中國(guó)家畜地方品種資源圖譜（下）［M］.北京：中國(guó)農(nóng)業(yè)出版社，2004.

［3］ 蘇理云，陳彩霞，高紅霞.SPSS 19統(tǒng)計(jì)分析基礎(chǔ)與案例應(yīng)用教程［M］.北京：北京希望電子出版社，2012.

［4］ 明道緒.生物統(tǒng)計(jì)附試驗(yàn)設(shè)計(jì)［M］.北京：中國(guó)農(nóng)業(yè)出版社，2002.

［5］ 索效軍，張年，熊琪，等.麻城黑山羊母羊體質(zhì)量與體尺的回歸分析［J］.東北農(nóng)業(yè)大學(xué)學(xué)報(bào)，2013，44（12）：63-67.

［6］ 楊佩燕，楊姣，劉志學(xué).崇明白山羊和徐淮山羊體重與體尺的相關(guān)和通徑分析［J］.上海農(nóng)業(yè)學(xué)報(bào)，2010，26（1）：24-27.

［7］ 蔡惠芬，陳志，羅衛(wèi)星，等.黔北麻羊體尺與體重的相關(guān)分析［J］.畜牧與獸醫(yī)，2012，44（7）：53-55.

［8］ 韓廷義.青海高原半細(xì)毛羊選育技術(shù)的研究［J］.中國(guó)畜牧獸醫(yī)，2012，39（12）：230-232.

［9］ 張力，江宵兵，王長(zhǎng)康.福建閩清毛腳雞體尺性狀與屠體性狀的典型相關(guān)分析［J］.西北農(nóng)林科技大學(xué)學(xué)報(bào)，2010，38（2）：8-12.

［10］馬猛，王克華，曲亮，等.13周齡母雞體尺性狀主基因+多基因混合遺傳模型分析［J］.東北農(nóng)業(yè)大學(xué)學(xué)報(bào)，2015，46（9）：75-82.

［11］Yücel C.Correlation and path coefficientanalyses ofseed yield components in the narbon bean（Vicia narbonensis L.）［J］.Turkish Journalof Agriculture and Forestry，2004，28（5）：371-376.

［12］買買提明·巴拉提，怕爾哈提，買合木提江，等.周歲麥蓋提羊體重與體尺指標(biāo)相關(guān)性的研究［J］.四川畜牧獸醫(yī)，2001，28（7）：20-21.

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

［14］王同軍.建昌黑山羊體尺體重相關(guān)性研究［J］.四川畜牧獸醫(yī)，2003，30（9）：22-26.

［15］高志英，烏云畢力克，艾爾肯江，等.德國(guó)美利奴成年母羊體尺與體重相關(guān)回歸分析［J］.草食家畜，2012，157（4）：31-36.

信息專欄

《中國(guó)草食動(dòng)物科學(xué)》面向全國(guó)讀者開展免費(fèi)閱讀

為進(jìn)一步提高本刊雜志論文的擴(kuò)散力，進(jìn)而提高雜志的影響因子和被引頻次，更為廣大科研人員創(chuàng)造良好的文獻(xiàn)查閱平臺(tái)，《中國(guó)草食動(dòng)物科學(xué)》于2013年在教育部科技發(fā)展中心的“中國(guó)科技論文在線”網(wǎng)站上面向全國(guó)讀者開展免費(fèi)閱讀，這將大大提高論文的引用率（3倍）。同時(shí)，本刊于2013年開始應(yīng)用DOI，意為“數(shù)字物品身份證符”。由于DOI的全球唯一性，在科技文獻(xiàn)中如果給出被引參考文獻(xiàn)的DOI，通過(guò)DOI可以直接定位到被引參考文獻(xiàn)上，這樣就提高了被引文獻(xiàn)的表現(xiàn)度，以及論文的擴(kuò)散力和影響力。

Correlation and Regression Analysis of Body Weight and Body Sizes in Black-bone Goats

Suo Xiaojun1，2，Zhang Nian1，Chen Mingxin1，etal
（1.Institute ofAnimalScience，HubeiAcademy ofAgriculturalSciences，Wuhan 430064，China；2.HubeiKey Laboratory ofAnimalEmbryo and Molecular Breeding，Wuhan 430064，China）

To provide a theoreticalbasis for species selection and breeding ofBlack-bone goat，the data ofbody weights and body sizes from 64 adult Black-bone does were analyzed by SPSS software with regression analysis and correlation analysis，stepwise regression analysis method，the optimalregression equation was established.The results showed that the body height（X1），body length（X2），chestcircumference（X3），cannon circumference（X4）showed very significantcorrelations with body weight（Y）atthe 0.01 level．The body weightwasextremely affected directly and indirectly by the chestcircumference，body length，body height.The bestregression equation ofbody weightwith body sizeswasobtained as：Y=-69.923 19+0.326 27X1+0.138 46X2+0.961 72X3.

Black-bone goat；body weight；body size；correlation analysis；stepwise analysis

S827.2

A

2095-3887（2017）05-0010-04

10.3969/j.issn.2095-3887.2017.05.003

2017-06-01

湖北省農(nóng)業(yè)科技創(chuàng)新中心資助項(xiàng)目（2016-620-004-002）

索效軍（1977-），男，副研究員。

陳明新，男，研究員，主要從事牛羊育種研究工作。