新疆地區(qū)西門塔爾犢牛體重與體尺指標(biāo)的相關(guān)及回歸分析

張曉雪，葛建軍，魏 趁，張夢(mèng)華，王 丹,尤震晨，黃錫霞*，馬光輝*

(1.新疆農(nóng)業(yè)大學(xué) 動(dòng)物科學(xué)學(xué)院，新疆 烏魯木齊 830052；2.新疆呼圖壁種牛場(chǎng)，新疆 呼圖壁 831203)

西門塔爾牛是原產(chǎn)于瑞士阿爾卑斯山區(qū)的乳肉兼用品種，被毛顏色為黃白花或紅白花。西門塔爾牛具有繁殖率高、泌乳性能好、耐粗飼、適應(yīng)性強(qiáng)、個(gè)體大、肉質(zhì)好、屠宰率高等優(yōu)良特性。各個(gè)國(guó)家都引進(jìn)西門塔爾牛，根據(jù)不同的需求，培育出了肉用西門塔爾、乳用西門塔爾和肉乳兼用西門塔爾。中國(guó)1990年在山東引進(jìn)該品種，對(duì)本地的黃牛進(jìn)行雜交改良，效果顯著。西門塔爾牛分布在中國(guó)東北的森林草原和科爾沁草原，中南的南嶺山脈和其山區(qū)，還有新疆的廣大草原和青藏高原等地。根據(jù)培育和改良品種的地點(diǎn)及環(huán)境不同，分為平原、草原、山區(qū)三個(gè)類群。

目前，對(duì)牛羊的生長(zhǎng)發(fā)育規(guī)律已有廣泛的研究[1-3]，而新疆地區(qū)西門塔爾犢牛生長(zhǎng)發(fā)育規(guī)律的研究鮮有報(bào)道。生長(zhǎng)發(fā)育規(guī)律研究是開(kāi)展遺傳育種、品種資源開(kāi)發(fā)利用的基礎(chǔ)工作，本文旨在分析新疆地區(qū)西門塔爾犢牛體尺與體重的相關(guān)性，總結(jié)犢牛生長(zhǎng)發(fā)育規(guī)律，為根據(jù)體尺指標(biāo)估測(cè)新疆地區(qū)西門塔爾犢牛的體重提供依據(jù)。

1 材料與方法

1.1 材料來(lái)源

數(shù)據(jù)來(lái)自新疆呼圖壁種牛場(chǎng)牧三場(chǎng)，對(duì)2011～2015年1～6月齡西門塔爾犢牛體尺和體重跟蹤測(cè)定6個(gè)月，共計(jì)4460條體尺、體重記錄。

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

測(cè)試指標(biāo)有體重、體高、體長(zhǎng)、胸圍、腿圍、管圍等。體重為清晨測(cè)量的牛只空腹重；體高為鬐甲最高點(diǎn)到地面的垂直距離，用測(cè)杖量取；胸圍是肩胛后角處體軀的垂直周徑，用卷尺量??；體長(zhǎng)是從肩端至坐骨端的距離，用卷尺量??；管圍是左前肢管部上1/3處的周徑，在前管的最細(xì)處用卷尺量取[4-5]；年齡(月)以呼圖壁種牛場(chǎng)牧三場(chǎng)西門塔爾牛產(chǎn)犢記錄為準(zhǔn)；性別以現(xiàn)場(chǎng)測(cè)量時(shí)記錄為準(zhǔn)。

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

把所有數(shù)據(jù)用Excel進(jìn)行整理，采用SPSS 13.0程序Linear和Bivariate過(guò)程進(jìn)行處理，并求Pearson相關(guān)系數(shù)。

1.3.1 多元線性回歸模型Y=b0+b1x1+b2x2+……+bnxn

其中：Y為依變量，

b0為常數(shù)；b1、b2、……bn為回歸系數(shù)；

x1、x2、……xn為回歸系數(shù)對(duì)應(yīng)的自變量。

2 結(jié)果與分析

2.1 新疆地區(qū)西門塔爾犢牛體尺體重測(cè)定結(jié)果

由表1可以看出，隨著月齡的增加新疆地區(qū)西門塔爾1～6月齡犢牛的體重、體尺指標(biāo)也不斷增加，同一性狀的增長(zhǎng)趨勢(shì)基本一致，相對(duì)于其他幾個(gè)性狀，犢牛更趨向于胸圍和體重的快速增長(zhǎng)，但性別不同各性狀的生長(zhǎng)快慢是明顯不同的。這可能和公犢初生重一般高于母犢有關(guān)[6-8]。由統(tǒng)計(jì)結(jié)果可以看出，1～2月齡公犢體重和體尺指標(biāo)均高于母犢。

對(duì)于體斜長(zhǎng)，其中4～5月齡時(shí)生長(zhǎng)速率最大，比上月增長(zhǎng)了8.19 cm，而此時(shí)公犢生長(zhǎng)速率最慢，2～3月齡公犢生長(zhǎng)速率最大，增長(zhǎng)了8.27 cm；對(duì)于胸圍，整體上生長(zhǎng)速度較快，3～4月齡增長(zhǎng)最快，5月齡后生長(zhǎng)速率下降；對(duì)于腿圍，整體上生長(zhǎng)速度較快，其中公犢在2～3月齡間生長(zhǎng)速率最快，母犢在4～5月齡增長(zhǎng)最快；對(duì)于管圍，隨著月齡的增加犢牛的管圍增加緩慢；對(duì)于體重，1～6月齡犢牛體重增長(zhǎng)的速率最大，其中1～3月齡公、母犢生長(zhǎng)速率基本一致，4月齡后公犢生長(zhǎng)速率明顯高于母犢。

表1 新疆地區(qū)西門塔爾犢牛的體重與體尺Table 1 Body weight and body size of Simmental calf in Xinjiang area cm/kg

注：同列數(shù)據(jù)(相同年齡公母牛之間比較)肩標(biāo)為相同字母或無(wú)字母表示差異不顯著(P>0.05)，不同小寫字母表示差異顯著(P<0.05)，不同大寫字母表示差異極顯著(P<0.01)。下表同。

Notes: In the same column(comparing between male and female of the same age), values with the same or no superscripts mean insignificant difference (P>0.05), different lowercase letter superscripts mean significant difference (P<0.05), different capital letter superscripts mean extremely significant difference (P<0.01). The same below.

2.2 相關(guān)性分析

由表2可以看出，對(duì)實(shí)測(cè)的新疆地區(qū)西門塔爾犢牛體重、體尺指標(biāo)進(jìn)行相關(guān)分析，結(jié)果體重與體高、體長(zhǎng)、胸圍、腿圍、管圍之間的相關(guān)系數(shù)分別為0.924、0.759、0.945、0.853、0.763。經(jīng)檢驗(yàn)這四個(gè)體尺指標(biāo)與體重的相關(guān)系數(shù)均達(dá)到極顯著水平(P<0.01)。

2.3 回歸分析

依體高(X1)、體斜長(zhǎng)(X2)、胸圍(X3)、腿圍(X4)、管圍(X5)作自變量，體重(Y)作依變量，分別采用Linear過(guò)程“Enter”法和“Stepwise”法建立多元線性回歸模型，從而得到了2個(gè)估測(cè)新疆地區(qū)西門塔爾犢牛體重的回歸方程(見(jiàn)表3)，但他們的擬合效果不同。

對(duì)于公犢牛，2個(gè)回歸模型的相關(guān)系數(shù)(R)是0.953和0.949，決定系數(shù)(R2)值為0.909和0.901；母犢牛的相關(guān)系數(shù)(R)是0.957和0.955，決定系數(shù)(R2)分別為0.915和0.911。R與1都十分接近，且R2>0，所以Y與Xi高度正相關(guān)[9]，說(shuō)明線性度很好，或者說(shuō)，當(dāng)采用四元或二元線性回歸方程來(lái)估測(cè)新疆西門塔爾犢牛體重時(shí)，它們能夠“解釋”公犢體重變異性的90.9%和90.1%、母犢體重變異性的91.5%和91.1%，可見(jiàn)應(yīng)用四元回歸模型更準(zhǔn)確一點(diǎn)。

表2 體重與體尺指標(biāo)的相關(guān)系數(shù)Table 2 Correlation coefficient of body weight with body size

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

Notes:** Probability valne of significance test is 2-tailed, loith 90% confidence inter vals.

表3 估測(cè)犢牛體重的回歸模型Table 3 Regression model of predicting body weight on calf

3 討 論

本研究發(fā)現(xiàn)舍飼新疆地區(qū)西門塔爾犢牛各性狀生長(zhǎng)最快集中在4～5月齡，1～5月齡間，母犢的體尺、體重平均增長(zhǎng)速度高于公犢。這與王鵬等[12]哺乳期短角犢牛的長(zhǎng)發(fā)育規(guī)律探討的結(jié)果基本相符。王文奇[10]研究表明：公母犢牛隨著月齡的增加，其相對(duì)于體高的發(fā)育而言，犢牛更趨向于體長(zhǎng)和胸圍的快速發(fā)育，趨向于橫向(胸圍向變粗的方向)發(fā)展。本文得到的新疆地區(qū)西門塔爾牛生長(zhǎng)規(guī)律與其他國(guó)內(nèi)外研究者的結(jié)果類似[7,11-13]。

本研究共得到了新疆地區(qū)西門塔爾犢牛體重與體高、體長(zhǎng)、胸圍、腿圍、管圍五個(gè)指標(biāo)之間的相關(guān)系數(shù)，均達(dá)到極顯著水平(P<0.01)。該結(jié)果與國(guó)內(nèi)鄭林等[14]和國(guó)外研究者[11]統(tǒng)計(jì)結(jié)果相似。

本研究建立的2個(gè)估測(cè)新疆地區(qū)西門塔爾犢牛體重的回歸模型表明了新疆地區(qū)西門塔爾犢牛體重與體高、體斜長(zhǎng)、胸圍、腿圍、管圍存在極顯著的多元線性回歸關(guān)系，此結(jié)論與國(guó)內(nèi)外相關(guān)研究相吻合[15-16]，然而判斷這2個(gè)回歸模型好壞的標(biāo)準(zhǔn)是多方面的。就數(shù)理統(tǒng)計(jì)角度而言，根據(jù)回歸模型來(lái)計(jì)算相關(guān)系數(shù)(R)及擬合值與實(shí)測(cè)值之間的決定系數(shù)(R2)，R和R2值越大模型的相關(guān)程度越高、擬合效果也就越好。從表3可以看出公犢牛四元回歸模型的R和R2值大于二元回歸模型的R和R2值，因此估測(cè)公犢牛體重時(shí)候以選擇四元回歸方程更為適宜，即：Y=-323.10+1.40X1+0.83X2+1.45X3-0.06X4+7.02X5，這與付雪峰[17]得出的結(jié)果相似，但如果同時(shí)考慮體尺指標(biāo)測(cè)量的繁瑣程度時(shí)可以考慮應(yīng)用二元回歸方程，即：Y=-290.26+1.83X3+2.45X1；對(duì)母犢牛四元回歸模型的R和R2值大于二元回歸模型的R和R2值，因此估測(cè)母犢牛體重時(shí)候以選擇四元回歸方程更為適宜，即：Y=-282.70+1.29X1+0.03X2+1.89X3+0.14X4+5.69X5，但如果同時(shí)考慮體尺指標(biāo)測(cè)量的繁瑣程度時(shí)，可以考慮應(yīng)用二元回歸方程，即：Y=-247.94+2.01X3+1.73X1來(lái)估測(cè)新疆地區(qū)西門塔爾犢牛體重。

4 結(jié) 論

新疆地區(qū)西門塔爾犢牛體重與體高、體長(zhǎng)、胸圍、腿圍、管圍五個(gè)指標(biāo)之間的相關(guān)系數(shù)分別為0.924、0.759、0.945、0.853、0.763。構(gòu)建的估測(cè)新疆地區(qū)西門塔爾犢牛體重的四元與二元回歸模型，決定系數(shù)均大于0.9。