錦江母牛體質(zhì)量與體尺指標(biāo)的相關(guān)與回歸關(guān)系

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

(江西農(nóng)業(yè)大學(xué)江西省動(dòng)物營養(yǎng)重點(diǎn)實(shí)驗(yàn)室/營養(yǎng)飼料開發(fā)工程研究中心，江西南昌 330045)

錦江母牛體質(zhì)量與體尺指標(biāo)的相關(guān)與回歸關(guān)系

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

(江西農(nóng)業(yè)大學(xué)江西省動(dòng)物營養(yǎng)重點(diǎn)實(shí)驗(yàn)室/營養(yǎng)飼料開發(fā)工程研究中心，江西南昌 330045)

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

錦江牛；體質(zhì)量；相關(guān)系數(shù)；回歸模型；體尺指標(biāo)

錦江牛是優(yōu)良的地方品種之一，產(chǎn)于我國南方亞熱帶丘陵地區(qū)，主要分布于江西省高安市、上高縣。該品種具有體型偏小、耐粗飼、抗病力強(qiáng)、性成熟早、適應(yīng)性強(qiáng)的遺傳特點(diǎn)，是進(jìn)行黃牛雜交改良的優(yōu)秀母本，也是培育當(dāng)?shù)厝馀Ｐ缕贩N不可多得的育種素材。本研究旨在探討錦江牛成牛母牛體質(zhì)量與體尺指標(biāo)的相關(guān)性，進(jìn)一步估測體質(zhì)量的回歸分析，以期為該品種牛的實(shí)際生產(chǎn)及選育工作提供依據(jù)。

1 材料與方法

1.1 材料來源與樣本數(shù)量

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

1.2 統(tǒng)計(jì)方法

采用SPSS 19.0軟件的Pearson和Linear過程[1]對(duì)試驗(yàn)數(shù)據(jù)進(jìn)行處理。

1.2.1 Pearson相關(guān)系數(shù)分析計(jì)算模型

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

1.2.2 多元線性回歸分析計(jì)算模型

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

式中：Y為因變量，b0為常數(shù)，b1、b2…bn為回歸系數(shù)，x1、x2…xn為回歸系數(shù)對(duì)應(yīng)的自變量。

2 結(jié)果與分析

2.1 相關(guān)性分析結(jié)果

以年齡、體高、十字部高、體斜長、胸圍、腹圍、管圍作為自變量，以體質(zhì)量作為因變量。各指標(biāo)之間的相關(guān)系數(shù)見表1。

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

2.2 通徑分析結(jié)果

各體尺指標(biāo)與體質(zhì)量的相關(guān)系數(shù)只能反映2個(gè)性狀間的表型相關(guān)，如果進(jìn)一步分析各性狀指標(biāo)間的直接或間接相互作用程度還需進(jìn)行通徑分析[2]。

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

2.3 體質(zhì)量與體尺指標(biāo)的回歸分析結(jié)果

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

表1 體質(zhì)量與體尺指標(biāo)的相關(guān)性分析結(jié)果

注：“*”“**”分別表示縱橫指標(biāo)間在0.05、0.01水平下相關(guān)性顯著。

表2 錦江牛體重與體尺指標(biāo)通徑分析結(jié)果

表3 回歸模型系數(shù)

表4 回歸模型擬合度

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

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

2.4 回歸系數(shù)的顯著性檢驗(yàn)

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

3 結(jié)論與討論

3.1 錦江牛母牛體質(zhì)量與體尺指標(biāo)的相關(guān)性

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

表5 回歸系數(shù)顯著性檢驗(yàn)結(jié)果

注：顯著性檢驗(yàn)概率值為雙尾概率，99%置信區(qū)間。

3.2 通徑分析

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

3.3 錦江牛母牛體質(zhì)量與體尺指標(biāo)的回歸分析

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

3.4 結(jié)論

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

[1]蘇金明，傅榮華，周建斌，等. 統(tǒng)計(jì)軟件SPSS for Windows實(shí)用指南[M]. 北京：電子工業(yè)出版社，2000：418-439.

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

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

[4]李永祿，邱 懷，馮仰廉，等. 養(yǎng)牛學(xué)[M]. 北京：中國農(nóng)業(yè)出版社，1987：235-236.

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

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

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

10.15889/j.issn.1002-1302.2017.02.045

2015-11-23

公益性行業(yè)(農(nóng)業(yè))科研專項(xiàng)(編號(hào)：201303143)；國家現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術(shù)體系(編號(hào)：nycytx-38)。

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

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

S823.1

A

1002-1302(2017)02-0152-02

許 超，高雨飛，彭志鵬，等. 錦江母牛體質(zhì)量與體尺指標(biāo)的相關(guān)與回歸關(guān)系[J]. 江蘇農(nóng)業(yè)科學(xué)，2017，45(2)：152-154.