元江本地黄牛在不同饲养方式下生长曲线拟合及体重与体尺指标线性回归分析研究

孙 峰，龚国安，龙云发，温昌饶，金福况，杨菊琴，白来章，杨灼雄

(1.元江县畜牧与饲草饲料站，云南 元江 653300；2.云南省元江县那诺乡农业农村综合服务中心，云南 元江 653305；3.云南省元江县咪哩乡农业农村综合服务中心，云南 元江 653303；4.曼来镇农业农村综合服务中心，云南 元江 6533001；5.云南省元江县甘庄街道农业农村综合服务中心，云南 元江 653399；6.云南省玉溪市畜牧站，云南 玉溪 653100)

1 材料与方法

1.1 试验设计

随机选取规模化牛场舍饲饲养的本地黄牛以及牛场周边农户放牧饲养的本地黄牛犊牛各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 多元线性回归方程

对本地放牧黄牛、本地舍饲黄牛的体重与体尺外貌指标、月龄进行回归分析，获得五元线性回归方程。回归方程决定系数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模型最优；两组牛体重、体尺与月龄之间均呈极显著相关，五元线性回归方程中月龄的回归系数都大于体尺性状的回归系数，其次体尺性状中管围也有很高的回归系数。