南江黄羊生长曲线拟合分析

李桂贤

福建省清流县嵩口畜牧兽医站，福建清流 365307

南江黄羊原产于四川省南江县，是我国培育的著名的肉用山羊品种，具有体格高大、生长快、肉质好、繁殖力高、适应性强等优点，适宜我国大部分农区饲养，福建省从1985 年起陆续到南江县引进种羊[1]并进行扩繁推广，现已成为福建省肉羊养殖主要品种。近年来，由于改变养殖方式，在自然放牧基础上进行精料补饲，同时福建省与四川省在气候、地理、牧草植被等方面存在较大的差异，生长性能与原产地可能会有差别。梁学武等[2]、李文杨等[3]、刘远等[4-5]在福建省相继开展波尔山羊、戴云山羊、福清山羊生长拟合方面研究，但鲜见南江黄羊在福建省生长发育的报道。为研究南江黄羊在福建省内的生长发育规律，对南江黄羊开展不同阶段的体重测定，并利用Von Bertalalanffy、Gom pertz、Logistic 3种数学模型建立生长曲线模型，以期为在福建省养殖南江黄羊制定合理的饲养方案提供依据。

1 材料与方法

1.1 试验羊选择

试验在福建省清流县某南江黄羊养殖场进行，选择出生日龄相差不超过10 d 的南江黄羊羔羊60 只，其中公母各半。待试验羊达到2 月龄时断奶。

1.2 饲养管理

试验羔羊断奶后，将公母统一编号并分栏饲养。每天放牧4～5 h，放牧时间从13:00 开始，到傍晚归牧。归牧后进行补饲精料，精料采购福建友容牧业公司生产的羔羊精料和育肥精料。补充量按前期（2～4 月龄）补充羔羊精料100 g/d，中期（5～9 月龄）、后期（10～12 月龄）分别补充育肥精料150、200 g/d。日常饲养管理及免疫、驱虫按常规进行。羔羊精料与育肥精料营养见表1。

表1 羔羊精料与育肥期精料营养水平 %

1.3 体重测定

每只试验羊在初生至12 月每月进行称重，除初生重外，各月龄体重均于月末07:00 放牧前空腹测定。

1.4 生长模型选择

应 用Von Bertalalanffy、Gom pertz、Logistic 3种非线性生长模型对各月龄体重进行拟合与分析，模型公式见表2。

表2 3 种非线性模型公式

1.5 数据处理及分析方法

原始数据用Excel 表格进行整理，并绘制南江黄羊的累积生长、绝对生长折线图。利用SPSS 24 非线性回归分析Von Bertalalanffy、Gom pertz、Logistic 3 种模型的生长曲线，求出函数公式的参数A、B、K的参数估算值和拟合度（R2），建立生长模型，同时计算出拐点体重、拐点月龄。

2 结果与分析

2.1 南江黄羊累积生长

南江黄羊累积生长见图1 所示。南江黄羊公、母羊累积生长随月龄增大而持续增长；0～2 月龄，公、母羊体重差异不大，公、母羊初生重分别为2.57、2.30 kg；3 月龄后公羊体重逐渐比母羊大，12月龄时公、母羊体重分别为36.62、31.23 kg，相差5.39 kg。

图1 南江黄羊累积生长

2.2 南江黄羊绝对生长

在绝对生长方面，南江黄羊公、母羊生长趋势基本相同，但公羊生长比母羊更快。1～3 月龄绝对生长下降，特别在2～3 月龄下降较大，可能与这时期羔羊受断奶应激、由母乳喂养向食草过渡影响有关；4 月龄后生长速度加快，6 月龄达到生长最高峰，公、母羊生长速度分别达到4.15、3.69 kg/月，7～12 月龄生长速度呈下降趋势（见图2）。

图2 南江黄羊绝对生长

2.3 南江黄羊生长曲线的拟合分析

利用Von Bertalalanffy、Gom pertz、Logist ic 3种生长曲线模型对南江黄羊公母羊的生长曲线拟合，计算出拟合参数估计值及拟合度（R2），结果见表3。从表3 可以看出，3 种生长曲线模型的拟合度均较高（均达0.99 以上），拟合度（R2）是评价非线性生长模型与实际生长发育拟合的重要指标，拟合度值越接近1，说明与实际越相符。通过体重实测值与3 种模型的体重估计值比较发现（表4），3 种模型对体重的拟合曲线估计与实际观测值基本一致，但初生实际值均低于3 种模型估计值，差异较大。

表3 南江黄羊生长3 种模型拟合参数

表4 南江黄羊公母羊体重实测值和各拟合曲线估计值 kg

3 种模型中，Von Bertalalanffy 模型不仅拟合度最高，而且在成熟体重、拐点月龄参数与实际更相符，表示该模型的拟合更与实际吻合。应用该模型建立南江黄羊公、母体重的生长曲线方程分别为W=46.77（1-0.57e-0.17t）3、W=39.08（1-0.56e-0.18t）3，拐点月龄分别为3.15、2.92 月；拐点体重分别为13.86、11.58 kg。

3 讨 论

3.1 养殖方式对南江黄羊生长的影响

品种、营养水平、饲养管理是影响畜禽生长的主要因素。本次试验南江黄羊公、母羊6 月龄体重分别为23.27、20.34 kg，徐永锋等[6]在原产地测定放牧的同月龄公、母羊体重分别为17.01、15.84 kg，分别增加6.26、4.50 kg；试验公、母羊12 月龄体重分别为36.62、31.23 kg。陈瑜等[7]在四川测定放牧的12 月龄公、母羊体重分别为35.67、29.63 kg，分别增加2.18、1.15 kg。南江黄羊在福建省生长性能得到提高，主要是进行补饲精料，同时福建省比四川省的年平均气温更高、无霜期更长，牧草更茂盛，特别在冬天枯草季节能采食更多的牧草，更能满足肉羊生长营养需求。

3.2 生长曲线对生产的指导作用

生长曲线的拟合对分析畜禽生长、制订合理饲养管理、选育研究等有重要的作用。通过3 种曲线模型对南江黄羊生长拟合，羔羊初生重实际值均低于估算值，可能是在母羊妊娠期没有补饲精料有关。南江黄羊个体较大，对营养水平要求较高，单纯放牧无法满足妊娠期营养需求，导致羔羊初生重偏小；羔羊在初生到3 月龄的绝对生长折线几乎呈垂直下降，可能与母羊泌乳不足、断奶应激有关。在生产实践中，应加强母羊妊娠期和哺乳期的营养，提高羔羊初生重和哺乳期的泌乳能力；羔羊断奶前应饲喂优质代乳料，尽早促进瘤胃发育，保证羔羊平稳断奶。本研究表明，南江黄羊公、母羊的生长拐点在3 月龄左右，在养殖中应加强拐点前的饲养管理与营养供给，最大限度满足其生长需求。

3.3 生长曲线的选择

不同畜禽品种所选最佳生长曲线有所差别。如姜勋平等[8]研究林细杂交羊生长曲线拟合，结果以Gom pertz 模型拟合度（R2）最大；梁学武等[2]研究波尔山羊的曲线拟合，结果以Logistic 模型拟合最好。本次在福建省内开展南江黄羊生长曲线研究，结果以Von Bertalalanffy 模型拟合最好，表现在拟合度最大，成熟体重、拐点月龄与实际更相符，这与李文杨等[3]对福清山羊的研究一致。为了模型公式各参数更准确，在测量时应尽可能提供大量翔实的数据分析；在最佳模型的选择上，除了比较拟合度大小外，还要结合估算出成熟体重、拐点日龄、拐点体重与实际情况是否相符，来选择最佳的模型。

4 结 论

采 用Von Bertalalanffy、Gom pertz、Logistic 3种常用的生长曲线模型拟合南江黄羊生长，各模型拟合度均≥0.99，其中以Von Bertalalanffy 模型拟合效果最好。建立南江黄羊公、母体重的生长曲线方程分别为W=46.77（1-0.57e-0.17t）3、W=39.08（1-0.56e-0.18t）3，拐点月龄分别为3.15、2.92 月；拐点体重为分别13.86、11.58 kg。南江黄羊在福建省养殖采用放牧+精料补饲的方式，生长性能得到提高。