成年酉州烏羊體尺長與體質量的關系

王高富　劉良佳　任航行　蔣婧　余琴　黎年富　夏元友　周鵬

摘要：為研究成年酉州烏羊體尺長與體質量的關系，測定16只酉州烏羊成年公羊和74只成年母羊的體斜長、體高、胸圍、管圍和體質量等指標，運用SPSS軟件分別分析公母羊體尺長與體質量間的相關系數、體尺長對體質量的直接和間接影響以及對體質量的決策程度，分別建立公羊和母羊體尺長與體質量的最優回歸模型。結果表明，胸圍和體斜長是酉州烏羊母羊體質量的最主要決策變量，胸圍是公羊體質量的最主要決策變量。母羊和公羊體尺對體質量的最優回歸模型分別是Y=-37569+0262X1+0702X3，Y=-44157+1035X3，其中X1表示體斜長，X3表示胸圍。

關鍵詞：酉州烏羊；體斜長；胸圍；體質量；回歸分析

中圖分類號： S8268+92文獻標志碼： A

文章編號：1002-1302（2017）21-0165-03

收稿日期：2016-06-15

基金項目：重慶市農發資金項目（編號：15404）；重慶市基本科研業務費計劃項目（編號：15435）。

作者簡介：王高富（1980—），男，重慶人，碩士，副研究員，研究方向為羊遺傳育種。E-mail：wanggaofs20031216@163com。

通信作者：周鵬，高級畜牧師，研究方向為動物遺傳育種。E-mail：cqzp2006@163com。

酉州烏羊主要分布于武陵山區腹地的酉陽土家族苗族自治縣，是重慶市特有的本地山羊品種，2009年正式被評定為國家級遺傳資源1]，該羊肉質細嫩，全身皮膚、眼、鼻、嘴、肛門、陰門等可視黏膜均為烏色2]，肉、骨、內臟均可入藥3-4]，被當地群眾稱為“藥羊”3]。但酉州烏羊群體數量小，現存欄數量不足5 000只，生長速度慢，成年個體體質量小，對其進行選育提高已經迫在眉睫。

體質量性狀是肉羊的重要經濟性狀，是肉羊育種過程中主要選擇性狀之一。而體尺性狀反映肉羊體格大小、體軀結構和組織器官發育情況，它與肉羊的繁殖、生長、抗病力以及對外界生活條件的適應能力等均有一定的相關性。因為體尺性狀數量較多，研究這些性狀與體質量的關系，找出較為主要的體尺性狀作為育種指標具有一定的實用價值。為了使酉州烏羊的保種選育工作更加高效快捷，本研究擬測定成年酉州烏羊的4個體尺和體質量性狀，分析它們之間的關系，以期建立體質量與體尺之間的回歸模型。

1材料與方法

11試驗對象、測定指標及方法

試驗在重慶市酉州烏羊資源保護場進行，隨機選擇常規飼養條件下的成年公羊16只，成年母羊74只進行體尺和體質量測定，試驗羊體況中上等，健康無病。測定指標包括體斜長（X1）、體高（X2）、胸圍（X3）、管圍（X4）和成年體質量（Y）；測定方法按照《畜禽遺傳資源調查技術手冊》執行5]，主要測定工具是測杖、卷尺、電子秤等。

12統計方法

試驗數據錄入Excel數據庫中，用SPSS統計軟件分析體質量與體尺間的相關系數、通徑系數、體尺對體質量的決策程度，并建立體質量與體尺的最優回歸模型6]。

2結果與分析

21各性狀的表型參數

由表1所測各性狀的平均數、標準差及變異系數可知，酉州烏羊成年公母羊體質量的變異系數在各性狀中是最大的，分別達到3153%和2403%，這表明酉州烏羊的體質量變化較大，具有較大的選育提高空間。另外，成年公羊的體斜長、體高、胸圍和管圍等體尺指標的變異系數均大于成年母羊，并且變異系數均較大，說明酉州烏羊體尺性狀也有較大的選育提高空間。

22體尺、體質量表型相關分析

由表2、表3可以看出，成年公母羊體尺與體質量都呈極顯著正相關（P<001），且公羊各性狀間的相關系數均大于母羊。其中公羊胸圍與體質量的相關系數最大，達到0972，體斜長次之，為0920，其他依次為管圍和體高；母羊胸圍與體質量的相關系數最大，為0893，其他依次為體斜長、體高和管圍。而公母羊在各體尺之間亦有較強的正相關，除母羊體斜長與管圍呈顯著正相關（P<005），其他性狀間都達到了極顯著相關（P<001），這說明成年酉州烏羊的體尺、體質量是相互依存的不同性狀， 對其中任何性狀進行選擇都應該綜合考慮。

23不同性別酉州烏羊體質量的正態性檢驗

SPSS對一組數據進行正態性檢驗有2種方法，即JP3]Kolmogorov-Smirnov Test和Shapiro-Wilk Test。Kolmogorov-Smimov Test檢驗結果較精確，適用于大樣本的檢測。而Shapiro-Wilk Test適用于小樣本的檢驗。本研究母羊樣本數量73只，公羊樣本數量16只，均屬于小樣本，因此對因變量Y進行正態性檢驗后利用Shapiro-Wilk Test的輸出結果。母羊Shapiro-Wilk統計量0968，P值=0061>005（表4），公羊統計量0953，P值=0535>005（表5），所以公母羊體質量均服從正態分布，可以進行回歸分析。

24通徑分析

各體尺性狀對體質量的相關僅反映兩性狀間的表型相關，而不能反映各體尺性狀與體質量的直接關系大小。因此，必須深入分析各因變量（體尺）對自變量（體質量）的直接影響，即通徑分析，進一步揭示各性狀間的相關原因。根據通徑分析原理，利用各性狀間的表型相關系數建立體尺對體質量的通徑系數正規方程組。

241母羊正規方程組方程組如下：

P1+0563P2+0728P3+0240P4=0755；

0563P1+P2+0626P3+0329P4=0624；

0728P1+0626P2+P3+0503P4=0893；

0240P1+0329P2+0503P3+P4=0402。

解方程得：P1=017；P2=0085；P3=0709；P4=-0025。endprint

通過對通徑系數進行顯著性檢驗可知：P3達到了極顯著的水平（P<001），P1達到顯著水平（P<005），P2、P4均未達到顯著水平，表明胸圍和體斜長對母羊體質量的直接影響較大，其中以胸圍對母羊體質量的直接影響作用最大（P3=0709），其次是胸體斜長（P1=017），體高和管圍對體質量的直接影響較小，這與表型相關的分析結果一致。

242公羊正規方程組方程組如下：

P1+0871P2+0928P3+0866P4=0920；

0871P1+P2+0873P3+0824P4=0874；

0928P1+0873P2+P3+0890P4=0972；

0866P1+0824P2+0890P3+P4=0912。

解方程得：P1=0046；P2=0055；P3=0697；P4=-0206。

通過對通徑系數進行顯著性檢驗可知：P3達到了極顯著的水平（P<001），P1、P2、P4均未達到顯著水平，表明胸圍對公羊體質量的直接影響較大，體斜長、體高和管圍對體質量的直接影響較小。

25相關系數分解

通徑分析在多元回歸的基礎上將相關系數riy分解為直接通徑系數（某一自變量對因變量的直接作用）和間接通徑系數（該自變量通過其他自變量對因變量的間接作用）。通徑分析的理論已證明，任一自變量（x）與因變量（y）之間的簡單相關系數（riy）=xi與y之間的直接通徑系數（Piy）和所有xi與Y的間接通徑系數，任一自變量xi對y的間接通徑系數=相關系數（riy）×通徑系數（Pjy）。

251母羊各體尺對體質量的相關系數分解

由表6可知，在4個自變量對體質量的直接影響方面，胸圍（X3）的直接作用最大，體斜長次之，管圍直接作用最小。通過分析各個間接通徑系數發現，體斜長通過胸圍對體質量的間接作用最大，其間接通徑系數達到0516 2，其次是體高和管圍通過胸圍對體質量的間接作用，其間接通徑系數分別達到了0443 8和0356 6，而體高和管圍對體質量的直接影響卻較小。由此可見，體高（X2）和管圍（X4）對體質量的影響主要取決于間接作用；胸圍（X3）對體質量不但有較強的直接影響，而且輔助其他性狀對體質量產生較大的間接作用。因此，胸圍和體斜長對酉州烏羊母羊體質量的增加具有重要作用，而體高和管圍對體質量改變影響不大，在選育時可以不必過多考慮。

252公羊各體尺對體質量的相關系數分解

由表7可知，4個自變量對體質量的直接影響中，胸圍X3的直接作用最大。通過分析各個間接通徑系數發現，體斜長通過胸圍對體質量的間接作用最大，其間接通徑系數達到0646 8，其次是管圍和體高通過胸圍對體質量的間接作用，其間接通徑系數分別達到了0620 3和0608 5，而體斜長、體高和管圍對體質量的直接影響卻較小。由此可見，體斜長（X1）、體高（X2）和管圍（X4）對體質量的影響主要取決于間接作用；胸圍（X3）對體質量不但有較強的直接影響，而且還輔助其他性狀對體質量產生較大的間接作用。因此，胸圍對酉州烏羊公羊體質量的增加具有重要作用。

26決策程度分析

袁志發等提出決策系數的概念，指出在復雜的路徑信息中，選擇什么樣的路徑對因變量Y最好，可由決策系數來決定7]。決策系數計算公式為R2（i）=2Piriy-Pi2。

由決策系數公式計算出酉州烏羊母羊各變量決策系數：

R2（1）=2×017×0755-0172=0227 8；

R2（2）=2×0085×0624-00852=0098 9；

R2（3）=2×0709×0893-07092=0763 6；

R2（4）=2×（-0025）×0402-（-0025）2=-0207 3。

決策系數排序為R2（3）﹥R2（1）﹥R2（2）﹥R2（4），且R2（4）﹤0。由決策系數可知，酉州烏羊母羊胸圍（X3）和體斜長（X1）是體質量的最主要決策變量，其次是體高（X2）和管圍（X4）。

由決策系數公式計算出酉州烏羊公羊各變量決策系數：

R2（1）=2×0046×0920-0046 2=0082 5；

R2（2）=2×0055×0874-0055 2=0093 1；

R2（3）=2×0697×0972-0697 2=0869 2；

R2（4）=2×（-0206）×0912-（-0206）2=-0418 2。

決策系數排序為R2（3）﹥R2（2）﹥R2（1）﹥R2（4），且R2（4）﹤0。由決策系數可知，酉州烏羊公羊胸圍（X3）是體質量的最主要決策變量，其次是體斜長（X1）、體高（X2）和管圍（X4）。

27最優線性回歸方程的建立

根據逐步線性回歸原理，從所有可供選擇的自變量中逐步選擇加入或剔除某個自變量，直到建立最優的回歸方程。運用SPSS分析程序分別獲得酉州烏羊母羊和公羊體質量和體尺之間的最優線性回歸方程，回歸系數等見表8。

由表8可知母羊各自變量的偏回歸系數、方程截距、標準回歸系數（即通徑系數）、標準誤差以及相對應的顯著性檢驗結果，從而可得線性回歸方程：

Y=-37569+0262 X1+0702 X3。

對成年母羊回歸模型的截距（a）和2個偏回歸系數（b1，b3）進行顯著性檢驗，均達到極顯著的水平（P<001）；同時，入選的2個體尺性狀對體質量的決定系數R2=0821。說明回歸方程極也顯著存在且所選性狀對體質量有較大的決定作用，擬合的方程也具有較強的參考價值。

由表9可知，公羊各自變量的偏回歸系數、方程截距、標準回歸系數（即通徑系數）、標準誤差以及相對應的顯著性檢驗結果，從而可得線性回歸方程如下：endprint

Y=-44157+1035X3。

對成年公羊回歸模型的截距（a）和一個偏回歸系數（b3）進行顯著性檢驗，表明達到了極顯著的水平（P<001）；同時，入選的此體尺性狀對體質量的決定系數R2=0945。說明回歸方程極顯著存在且所選性狀對體質量有較大的決定作用，擬合的方程具有較強的參考價值。

3討論

31主要選育指標的確定

本研究選擇了酉州烏羊的4個成年體尺指標作為研究對象，從相關分析、通徑分析及決定程度分析表明，無論是成年公羊還是母羊，其胸圍均是影響體質量的一個主要指標，與湖北烏羊和大足黑山羊等研究結果8-9]有所不一致，湖北烏羊公羊體質量最主要影響因素是體長，大足黑山羊體質量最主要影響因素是胸寬和管圍，這可能是由于品種和生長環境不同所造成的。因此，在酉州烏羊本品種選育過程中應該以胸圍為主并兼顧體斜長和管圍作為主要選擇指標。

32體質量與體尺性狀的作用關系

由相關系數剖分可知，不管母羊或者公羊的體質量和胸圍的相關主要由兩者直接作用而形成的，間接作用影響較小。而與體質量的相關性次之的為體斜長、體高和管圍，對體質量的影響主要通過間接作用產生，直接效應較小。這表明畜體是一個交互作用的有機整體，選擇時應綜合考慮才能取得較好的效果。

33需要進一步研究的問題

本研究只選擇了幾個重要的體尺性狀，其他體尺性狀如胸深、胸寬和腰角寬等是否存在與體質量更強的相關性，還有待進一步探討。此外，由于本研究公羊樣本量較小，變異系數較大，在以后的研究中要擴大公羊樣本量，進一步研究體尺與體質量的相關關系。

HS223]參考文獻：

1]ZK（#]趙金紅，王高富，陳靜，等 酉州烏羊胚胎移植效果初報J] 上海畜牧獸醫通訊，2012（1）：014

2]周鵬，王高富，任航行，等 酉州烏羊遺傳資源的保護對策J] 上海畜牧獸醫通訊，2013（5）：75-75

3]夏元友，彭春江，鄭義 酉州烏羊品種特性及生產性能初探J] 畜禽業，2005（7）：44-45

4]周鵬，黃勇富，王高富，等 酉州烏羊生理生化指標的研究J] 上海畜牧獸醫通訊，2012（4）：32-33

5]陳偉生，徐桂芳 中國家畜地方品種資源圖譜（下）M] 北京：中國農業出版社，2004

6]杜家菊，陳志偉 使用SPSS線性回歸實現通徑分析的方法J] 生物學通報，2010，45（2）：4-6

7]袁志發，周靜芋，郭滿才，等 決策系數——通徑分析中的決策指標J] 西北農林科技大學學報（自然科學版），2001，29（5）：131-133

8]韓燕國，王黨偉，何春波，等 湖北烏羊體尺與體重的關系研究J] 中國草食動物科學，2012（專輯）：366-369

9]陳永軍，趙中權，張家驊，等 成年大足黑山羊體重與體尺通徑分析及最優回歸模型的建立J] 草食家畜，2008（3）：71-74endprint