錦江母牛體質量與體尺指標的相關與回歸關系

許 超， 高雨飛， 彭志鵬， 歐陽克蕙， 瞿明仁， 黎力之， 熊小文， 溫慶琪， 許蘭嬌

(江西農業大學江西省動物營養重點實驗室/營養飼料開發工程研究中心，江西南昌 330045)

錦江母牛體質量與體尺指標的相關與回歸關系

許 超， 高雨飛， 彭志鵬， 歐陽克蕙， 瞿明仁， 黎力之， 熊小文， 溫慶琪， 許蘭嬌

(江西農業大學江西省動物營養重點實驗室/營養飼料開發工程研究中心，江西南昌 330045)

以江西省高安市錦江牛保種區2015年測量獲得的132條相關數據為基礎材料，對錦江牛體質量與年齡、體高、十字部高、體斜長、胸圍、腹圍、管圍的相關系數進行分析，并進一步分析了估測錦江牛成年母牛體質量的回歸模型。結果表明：錦江牛體質量與年齡、體高、十字部高、體斜長、胸圍、腹圍、管圍的相關系數分別為0.496、0.704、0.699、0.661、0.893、0.838、0.605；得到了3個估測體質量的回歸模型，估測值與實測值之間的相關程度分別為 0.893、0.907、0.913。

錦江牛；體質量；相關系數；回歸模型；體尺指標

錦江牛是優良的地方品種之一，產于我國南方亞熱帶丘陵地區，主要分布于江西省高安市、上高縣。該品種具有體型偏小、耐粗飼、抗病力強、性成熟早、適應性強的遺傳特點，是進行黃牛雜交改良的優秀母本，也是培育當地肉牛新品種不可多得的育種素材。本研究旨在探討錦江牛成牛母牛體質量與體尺指標的相關性，進一步估測體質量的回歸分析，以期為該品種牛的實際生產及選育工作提供依據。

1 材料與方法

1.1 材料來源與樣本數量

試驗采用江西省高安市錦江牛產區2015年8月測量獲得的132條相關數據。體質量(kg)為第2天清晨測量的?？崭官|量，體斜長(cm)為肩端到坐骨端的距離，體高(cm)為鬐甲最高點至地面的垂直距離，十字部高(cm)為兩腰角中央到地面的垂直距離，胸圍(cm)為沿肩胛骨后角處量取的體軀垂直周徑，管圍(cm)為牛左前肢管骨上1/3最細處測量的水平周徑。

1.2 統計方法

采用SPSS 19.0軟件的Pearson和Linear過程[1]對試驗數據進行處理。

1.2.1 Pearson相關系數分析計算模型

式中：R為Pearson相關系數，X為自變量，Y為因變量。R值介于-1～1之間，當R值為正數時為正相關，表示因變量隨自變量的增大而增大；當R值為負數時為負相關，表示因變量隨自變量的增大而減小；當R值等于0時，表示因變量與自變量之間沒有相關性。

1.2.2 多元線性回歸分析計算模型

Y=b0+b1x1+b2x2+…+bnxn。

式中：Y為因變量，b0為常數，b1、b2…bn為回歸系數，x1、x2…xn為回歸系數對應的自變量。

2 結果與分析

2.1 相關性分析結果

以年齡、體高、十字部高、體斜長、胸圍、腹圍、管圍作為自變量，以體質量作為因變量。各指標之間的相關系數見表1。

由表1可知，年齡、體高、十字部高、體斜長、胸圍、腹圍、管圍與體質量之間的相關系數分別為0.496、0.704、0.699、0.661、0.893、0.838、0.605。經檢驗，年齡、體尺指標與體質量之間的相關性均極顯著(P<0.01)。其中，胸圍與體質量之間的相關性最大；在體尺指標中，體高與十字部高的相關性最大，管圍與其他體尺指標的相關性較小。

2.2 通徑分析結果

各體尺指標與體質量的相關系數只能反映2個性狀間的表型相關，如果進一步分析各性狀指標間的直接或間接相互作用程度還需進行通徑分析[2]。

由通徑分析(表2)可知，錦江牛胸圍通過體斜長和腹圍對體質量的間接通徑系數分別為0.177×0.596=0.105、0.206×0.867=0.179，總間接作用為0.105+0.179=0.284。由通徑分析的理論可知，其相關系數為0.609+0.177×0.596+0.206×0.867=0.893，該結果與表1一致。直接作用中，胸圍(0.609)最大，腹圍(0.206)次之，體斜長(0.177)最小；間接作用中，胸圍(0.891)最大，腹圍(0.300)次之，體高(0.209)最小。

2.3 體質量與體尺指標的回歸分析結果

分別采用Linear過程“Enter”法和“Stepwise”法建立多元線性回歸方程，統計結果見表3、表4。 由表3可知，3個回歸模型分別為：體質量=-291.738+3.560×胸圍；體質量=-353.458+3.085×胸圍+1.186×體斜長；體質量=-359.525+2.428×胸圍+1.052×體斜長+0.707×腹圍。

表1 體質量與體尺指標的相關性分析結果

注：“*”“**”分別表示縱橫指標間在0.05、0.01水平下相關性顯著。

表2 錦江牛體重與體尺指標通徑分析結果

表3 回歸模型系數

表4 回歸模型擬合度

注：因變量為體質量。類別1的回歸模型組成：常量、胸圍；類別2的回歸模型組成：常量、胸圍、體斜長；類別3的回歸模型組成：常量、胸圍、體斜長、腹圍。

由表4可知，3個回歸模型的復相關系數(R)分別為 0.893、0.907、0.913，擬合度的確定系數(R2)分別為0.797、0.823、0.833。可見，回歸模型的線性度較好，尤其是第2、3個回歸模型擬合度的確定系數(R2)更高，線性度更好。

2.4 回歸系數的顯著性檢驗

由回歸系數顯著性檢驗結果(表5)可知，3個回歸模型中回歸系數的F值分別為510.223、299.316、212.626，其相應的顯著性檢驗值均小于0.001，差異極顯著(P<0.01)，表明3個回歸模型中回歸系數是合理的。

3 結論與討論

3.1 錦江牛母牛體質量與體尺指標的相關性

胸圍是對中國黃牛及其雜交后代影響最大的體尺因素[3]，本研究對體高、十字部高、體斜長、胸圍、腹圍、管圍6個體尺指標進行相關分析并得到了相似結果。管圍是評估公牛體質量回歸模型中是常見的自變量，與公牛體質量呈明顯正相關關系，但對于母牛其相關性卻很弱[4-6]。本研究中，錦江牛母牛管圍與體質量的相關性很弱，這與王占紅等的研究結論[7]一致。

表5 回歸系數顯著性檢驗結果

注：顯著性檢驗概率值為雙尾概率，99%置信區間。

3.2 通徑分析

通徑分析表明，胸圍、腹圍、體斜長對體質量的直接作用力和間接作用力均很大，且三者與體質量的相關性極顯著，因此對錦江牛母牛選育時應首先考慮胸圍、腹圍、體斜長3個因素。本研究所測體尺指標對體質量回歸方程的最大決定系數R2=0.833，表明方程中包含了估計體質量的主要體尺因素；1-R2=0.167，表明還有一些對體質量有影響的因素尚未考慮，仍需進一步研究。

3.3 錦江牛母牛體質量與體尺指標的回歸分析

本研究共得到3個估測錦江成年母牛體質量的回歸模型，表明錦江牛成年母牛體質量與胸圍、體斜長、腹圍存在明顯的線性關系，此結論與國內外相關研究相吻合。判斷回歸模型優劣的標準是多方面的，就數理統計角度而言，根據回歸模型計算擬合值與實測值之間的決定系數(R2)，通過比較R2值的大小來擇其優劣更為適宜。由表3可知，第2個回歸模型的R2值0.823與第3個回歸模型的R2值0.833大致相等，兩者均大于第1個回歸模型的R2值0.793。但如果同時考慮體尺測量的繁瑣程度可選取第1個回歸模型，即：體質量=-291.738+3.560×胸圍。

3.4 結論

錦江牛母牛體高、十字部高、體斜長、胸圍、腹圍、管圍6個體尺指標與體質量的相關系數分別為0.704、0.699、0.661、0.893、0.838、0.605，其中胸圍對體質量的影響最大。錦江牛母牛體尺指標與體質量的多元線性回歸模型為：(1)體質量=-291.738+3.560×胸圍(R=0.893，P<0.01)；(2)體質量=-353.458+3.085×胸圍+1.186×體斜長(R=0.907，P<0.01)；(3)體質量=-359.525+2.428×胸圍+1.052×體斜長+0.707×腹圍(R=0.913，P<0.01)。

[1]蘇金明，傅榮華，周建斌，等. 統計軟件SPSS for Windows實用指南[M]. 北京：電子工業出版社，2000：418-439.

[2]田亞磊，高騰云，張花菊，等. 魯山牛腿山羊體重與體尺性狀相關性分析[J]. 江蘇農業科學，2009(6)：288，291.

[3]焦平林. 秦川淮北及秦淮F1黃牛體尺體重回歸分析[J]. 黃牛雜志，1998，24(4)：6-7.

[4]李永祿，邱 懷，馮仰廉，等. 養牛學[M]. 北京：中國農業出版社，1987：235-236.

[5]Faroppa V. Standard of the breed and regulations[J]. La Razza Bovina Piermontese，1993(4)：9-19.

[6]曹紅鶴. 意大利皮埃蒙特肉牛線性體型評分方法[J]. 黃牛雜志，1999，25(4)：17-19.

[7]王占紅，何永濤，郭維春，等. 沿江牛成年母牛體重與體尺指標的相關與回歸分析[J]. 中國畜牧雜志，2007，43(1)：12-13.

10.15889/j.issn.1002-1302.2017.02.045

2015-11-23

公益性行業(農業)科研專項(編號：201303143)；國家現代農業產業技術體系(編號：nycytx-38)。

許 超(1988—)，男，江蘇宿遷人，碩士研究生，主要從事反芻動物營養研究。E-mail：229462924@qq.com。

歐陽克蕙。E-mail：ouyangkehui@sina.com。

S823.1

A

1002-1302(2017)02-0152-02

許 超，高雨飛，彭志鵬，等. 錦江母牛體質量與體尺指標的相關與回歸關系[J]. 江蘇農業科學，2017，45(2)：152-154.