元江本地黃牛在不同飼養方式下生長曲線擬合及體重與體尺指標線性回歸分析研究

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

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

1 材料與方法

1.1 試驗設計

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

1.2 試驗牛

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

1.3 飼養管理

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

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

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

1.4 測定指標及方法

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

1.5 數據處理與分析

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

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

2 結果與分析

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所示，各組生長曲線均呈現出一定程度的“S”形。

圖1 不同飼養方式對本地黃牛體重生長曲線

表1 不同飼養方式對本地黃牛體重生長曲線模型及擬合參數

2.2 體重、體尺相關性分析

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

表2 不同飼養方式與本地黃牛體重，月齡，體尺相關性

2.3 多元線性回歸方程

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

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

公牛：本地放牧黃牛，體重=-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 生長曲線模型

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

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

3.2 多元回歸方程

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

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

4 結論

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