孫 峰，龔國安，龍云發(fā)，溫昌饒，金福況，楊菊琴，白來章，楊灼雄

(1.元江縣畜牧與飼草飼料站，云南 元江 653300；2.云南省元江縣那諾鄉(xiāng)農(nóng)業(yè)農(nóng)村綜合服務中心，云南 元江 653305；3.云南省元江縣咪哩鄉(xiāng)農(nóng)業(yè)農(nóng)村綜合服務中心，云南 元江 653303；4.曼來鎮(zhèn)農(nóng)業(yè)農(nóng)村綜合服務中心，云南 元江 6533001；5.云南省元江縣甘莊街道農(nóng)業(yè)農(nóng)村綜合服務中心，云南 元江 653399；6.云南省玉溪市畜牧站，云南 玉溪 653100)

1 材料與方法

1.1 試驗設計

隨機選取規(guī)?；錾犸曪曫B(yǎng)的本地黃牛以及牛場周邊農(nóng)戶放牧飼養(yǎng)的本地黃牛犢牛各20頭、公母各10頭，從初生到12月齡逐月測定體重體尺，并擬合適宜的生長發(fā)育模型、體重與體尺指標的回歸模型。

1.2 試驗牛

試驗放牧的元江本地犢牛為隨機選取牛場周邊1個自然村若干肉牛養(yǎng)殖戶出生月齡相近的犢牛，采取放牧的方式；舍飼的本地黃牛為牛場牛只，飼養(yǎng)管理同牛場其他雜交牛一致。母牛均為世代在元江當?shù)仫曫B(yǎng)，犢牛通過本交繁育，無凍精改良或引種等引進外來血緣的情況。

1.3 飼養(yǎng)管理

1.3.1 本地黃牛農(nóng)戶放牧飼養(yǎng) 本地放牧黃牛白天放牧，牧歸補充作物秸稈等粗飼料，犢牛隨母哺乳、天然放牧、自然斷奶。粗飼料均為本地當季的新鮮飼草或作物秸稈，母牛在妊娠期及產(chǎn)犢后補充數(shù)量玉米面。

1.3.2 本地黃牛舍飼 母牛：母牛產(chǎn)后提供營養(yǎng)全面、質(zhì)量較好的精、粗飼料，精料以配合飼料為主，粗料為當季所儲存的干草料、青貯、新鮮飼草或作物秸稈等，早晚各投喂1次，每頭母牛精料投喂量每天不少于2 kg，根據(jù)不同母牛差異可做調(diào)整，自由飲水，自由運動，早期斷奶，早期配種。

犢牛：犢牛出生后保證及時吃到初乳，7 d后用優(yōu)質(zhì)青干草誘飼，15 d左右飼喂少量顆粒飼料，1個月后自由采食青干草及顆粒飼料，當犢牛日平均采食量達自身體重的1%以上，2～3月齡時，母犢分離，完全斷奶。斷奶后自由采食粗飼料，自由飲水。精料：犢牛開食顆粒料，選用云南龍谷生物科技有限公司生產(chǎn)的犢牛專用顆粒料，約6月齡后逐漸轉(zhuǎn)為正常的飼糧(日糧)飼喂。

1.4 測定指標及方法

所有試驗牛只從初生到12月齡每月定期在早飼前，空腹狀態(tài)下進行性能測定數(shù)據(jù)為試驗犢牛的體重、體高、體斜長、胸圍及管圍。

1.5 數(shù)據(jù)處理與分析

根據(jù)元江本地黃牛的體重、體尺隨月齡的增加情況，運用Gompertz模型、Logistic模型、Von Bertalanffy模型、Brody模型，選擇出最契合的生長曲線模型方程。用相關(guān)指數(shù)R2(擬合度)的大小評價擬合的結(jié)果。

根據(jù)元江本地黃牛的體重體尺隨月齡的增加情況，運用體尺、月齡與體重進行逐步回歸分析，獲得五元線性回歸方程。并用相關(guān)系數(shù)R2及估計標準誤差的大小評價擬合的結(jié)果。相關(guān)系數(shù)R2越大或估計標準誤差越小，說明該多元回歸方程越能反應該體重與體尺、月齡的規(guī)律。

2 結(jié)果與分析

2.1 體重生長曲線

由表1可知，Gompertz、Logistic、Von Bertalanffy、Brody 4種模型擬合肉牛體重的擬合度(R2)均較高，均達0.9以上，但是均低于。其中以Logistic模型擬合度(R2)最高，本地舍飼黃牛均高于本地放牧黃牛。本地放牧黃牛和舍飼黃牛公牛拐點月齡分別為4.75月、5.23月，拐點體重分別為70.23 kg、101.28 kg。體重生長模型為：本地放牧黃牛，y=140.456/(1+5.012×e-0.339×月齡)，R2=0.911；本地舍飼黃牛，y=202.55/(1+7.064×e-0.374×月齡)，R2=0.962；母牛拐點月齡分別為5.54月、5.30月，拐點體重分別為68.15 kg、93.34 kg。體重生長模型為：本地放牧黃牛，y=136.291/(1+5.6×e-0.311×月齡)，R2=0.904；本地舍飼黃牛，y=186.681/(1+7.271×e-0.374×月齡)，R2=0.969；本地放牧和舍飼黃牛體重生長曲線如圖1所示，各組生長曲線均呈現(xiàn)出一定程度的“S”形。

圖1 不同飼養(yǎng)方式對本地黃牛體重生長曲線

表1 不同飼養(yǎng)方式對本地黃牛體重生長曲線模型及擬合參數(shù)

2.2 體重、體尺相關(guān)性分析

由表2可見，元江本地放牧及舍飼黃牛公母牛體重、體尺與月齡之間均呈極顯著相關(guān)(P<0.01)。

表2 不同飼養(yǎng)方式與本地黃牛體重，月齡，體尺相關(guān)性

2.3 多元線性回歸方程

對本地放牧黃牛、本地舍飼黃牛的體重與體尺外貌指標、月齡進行回歸分析，獲得五元線性回歸方程。回歸方程決定系數(shù)R2表明回歸平方和占總變異量的百分比，數(shù)據(jù)越大，表明擬合度越好，表明體高、體斜長、胸圍、管圍、月齡均對體重有較大的決定作用，具有較高的參考價值。由表3可知，五元線性回歸方程中，月齡的回歸系數(shù)都大于其余性狀的回歸系數(shù)，其五元線性回歸方程如下。

表3 不同飼養(yǎng)方式與本地黃牛五元線性回歸分析

公牛：本地放牧黃牛，體重=-78.988+3.015×月齡-0.123×體高+0.291×體斜長+1.223×胸圍+0.778×管圍，R2=0.975；本地舍飼黃牛，體重=-99.127+5.308×月齡+0.408×體高-0.214×體斜長+1.146×胸圍+3.031×管圍，R2=0.987。

母牛：本地放牧黃牛，體重=-76.739+3.713×月齡+0.632×體高+0.214×體斜長-0.057×胸圍+6.013×管圍，R2=0.975；本地舍飼黃牛，體重=-67.661+7.566×月齡+0.107×體高+0.385×體斜長+0.227×胸圍+5.088×管圍，R2=0.980。

3 討論

3.1 生長曲線模型

動物機體的生長曲線擬合是分析和研究其生長發(fā)育規(guī)律的重要方法，目前最常見的生長曲線模型主要有Gompertz、Logistic、Von Bertallanffy、Brody 4種模型，其特點是擬合效果較好，參數(shù)估計值等能夠反映實際生產(chǎn)中的一些問題[1]。生長曲線拐點是反映畜禽機體生長率最大的時間點，根據(jù)曲線方程可推算出拐點體重及達到拐點體重時的月齡等重要指標[2]。擬合度較高的生長曲線模型在確定種畜禽適配年齡，指導畜禽改良及選育工作，和確定出欄、屠宰年齡及控制生長發(fā)育上均具有主要的參考價值[3]，目前國內(nèi)外許多肉牛品種的生長曲線都已建立。但也有相關(guān)研究表明對于不同的品種、性狀等，最佳的生長曲線擬合模型不盡相同[4]，故應該根據(jù)實際情況具體分析。通過對元江本地黃牛生長曲線的擬合，求出最佳生長曲線并得出各自的拐點月齡，而拐點是肉牛生長率的轉(zhuǎn)折點，我們應該在實際的生產(chǎn)活動中更加注意拐點階段的飼糧營養(yǎng)水平[5]，以達到犢牛個體生長速度減速時間后移及拐點體重最大化的目的，為元江本地黃牛的育肥及出欄、選種選育等提供理論依據(jù)。

本研究結(jié)果顯示拐點體重舍飼牛高于放牧牛，拐點體重的大小可以作為比較一個牛牛群飼養(yǎng)水平的重要依據(jù)，我們可以在今后的肉牛養(yǎng)殖中可以通過拐點月齡階段飼糧營養(yǎng)水平的調(diào)節(jié)、因地制宜選擇適合的養(yǎng)殖模式實現(xiàn)對元江肉牛養(yǎng)殖的提質(zhì)增效。

3.2 多元回歸方程

多元線性回歸在生物統(tǒng)計中具有舉足輕重的地位。畜禽體尺數(shù)據(jù)一方面反映了其體軀大小和結(jié)構(gòu)、發(fā)育狀況等情況，另一方面也可以推算出機體組織、器官的發(fā)育程度及體重，它與畜禽的生產(chǎn)性能、生理機能、抗病性能以及對適應能力等息息相關(guān)[6]。我們可以通過回歸關(guān)系顯著的回歸方程進行估算畜禽體重，指實際生產(chǎn)活動，提高經(jīng)濟效益。

R2作為判斷回歸模型好壞及選擇的標準，其值越大模型的相關(guān)程度越高、擬合效果也越好，回歸方程中各項指標回歸系數(shù)越高，證明該指標對體重的影響越大。本研究中，五元線性回歸方程均有較高的擬合度(R2>0.97)，其中月齡的回歸系數(shù)都大于體尺性狀的回歸系數(shù)，體尺性狀中管圍除本地放牧黃牛公牛組外也有很高的回歸系數(shù)，故建議對元江不同品種肉牛建立多元線性回歸方程時，除月齡外還應該充分考慮管圍的影響，從而使得其多元線性回歸方程更加完善，這與何鴻源[7]研究香格里拉牦牛的結(jié)果類似。

4 結(jié)論

元江本地黃牛體重生長曲線的擬合以Logistic模型最優(yōu)；兩組牛體重、體尺與月齡之間均呈極顯著相關(guān)，五元線性回歸方程中月齡的回歸系數(shù)都大于體尺性狀的回歸系數(shù)，其次體尺性狀中管圍也有很高的回歸系數(shù)。