哺乳期放养云岭山羊羔羊生长曲线拟合分析

和贵文，和云凤，杨亚晋，郭爱伟

（西南林业大学生命科学学院，云南 昆明 650224）

云岭山羊是肉皮兼用型山羊品种， 主要分布于云南省境内云岭山系及其余脉的哀牢山、 无量山和乌蒙山延伸地区，是云南省数量最多、分布最广的山羊地方良种，具有耐寒、耐热、耐粗饲、抗病力和适应性强等特点［1］。 云岭山羊羔羊在哺乳期间的生长发育会直接影响断奶后的生产潜能，继而影响后期的生产性能， 羔羊在哺乳期发育良好可使养殖者获得更好的生产和经济效益［2］。 羔羊断奶体重受多种因素的影响， 初生体重是关键因素之一，此外母羊在妊娠期营养状况、产羔年龄以及饲养环境、 饲粮营养水平等均会影响羔羊断奶体重。研究畜禽生长发育规律，通常需要测定体重和体尺，分析不同时期的生长曲线，进行体重和体尺间的相关性分析。通过对畜禽生长发育分析，掌握畜禽生长发育规律， 可以提高选育效果和养殖效益，对家畜的合理饲养和补饲具有重要意义［3-5］。目前对云岭山羊的研究主要集中在整个生产周期以及生产性能和生理生化指标等方面， 对哺乳期羔羊生长情况的研究缺少基础数据。 生长曲线可以直观描述畜禽在生长过程中随年龄增长生长发育数据的累积变化轨迹。 常用的非线性模型有Logistic、Gompertz 和Von Bertalanffy 模型等，不同畜禽的生长轨迹符合不同的非线性模型。 目前对放牧加补饲后山羊生长发育规律的研究较多，或对哺乳期不同补饲效果进行了探究， 但对单一放牧条件下哺乳期云岭山羊羔羊生长发育的研究缺乏报道。 研究放养条件下哺乳期云岭山羊羔羊生长发育规律可以找出哺乳期放养云岭山羊羔羊最佳的补饲时间， 充实放养条件下哺乳期云岭山羊羔羊生长发育的基础数据。 该研究分析了哺乳期放养云岭山羊羔羊体重和体尺增长规律， 采用Logistic、Gompertz 和Von Bertalanffy 3 种非线性模型［6］拟合哺乳期放养云岭山羊羔羊生长曲线，并计算了拐点体重以及拐点周龄， 分析了初生体重与试验末体重的相关性， 为今后哺乳期放养云岭山羊羔羊进行合理补饲提供依据。

1 材料与方法

1.1 试验羊及饲养管理

试验羊饲养于云南省丽江市玉龙纳西族自治县宝山乡。选胎次、体重接近的母羊24 只，共产羔羊36 只，公羔17 只，母羔19 只。 试验期间每天8：00—19：00 母羊和羔羊在外自由采食， 晚上在半开放式圈舍休息，期间无补饲。 试验期为56 d。

1.2 测定指标

主要测定羔羊0～8 周龄的体重、体尺指标。每周测定空腹体重，隔1 周测定体高、体斜长、胸围、胸深、胸宽等，具体测定方法参考《绵、山羊生产性能测定技术规范》［7］（NY/T 1236—2006）。

体高：肩胛最高点到地面垂直距离，单位为cm。

体斜长：胸突到坐骨结节后端直线距离，单位为cm。

胸围：肩胛后端绕胸一周长度，单位为cm。

胸深：肩胛最高到胸突直线距离，单位为cm。

胸宽：肩胛最宽处左右两侧直线距离，单位为cm。

1.3 生长曲线拟合与分析

采用SPSS 21.0 软件中Logistic、Gompertz 和Von Bertalanffy 3 种非线性模型分别对公羔、 母羔生长曲线进行拟合，用拟合度（R2）对3 条曲线进行验证，通过Excel 绘制出生长曲线。 通过拟合公式生成的A、B、K 值计算拐点体重、拐点日龄以及预测体重等值， 经线性相关分析对羔羊初生体重与羔羊断奶体重密切程度与相关方向进行分析得到Pearson（皮尔逊）相关系数。 拟合曲线公式见表1。

相关性分析公式如下：

2 结果与分析

2.1 哺乳期放养云岭山羊羔羊生长发育规律分析

哺乳期放养云岭山羊羔羊体重增长曲线如图1 所示， 在哺乳期间放养云岭山羊羔羊体重随周龄的增加而增长。 由表2 可知，8 周龄时公羔平均体重为（10.16±2.26）kg，母羔平均体重为（9.54±1.68）kg，公羔平均体重高于母羔；周增重在第1 周最大（公羔1.96 kg，母羔1.15 kg），之后增重速率较第1 周减缓。 由表3 可知，随着周龄的增加，体高、体斜长、胸围、胸宽、胸深有增加的趋势，但胸宽和胸深增加较小。 断奶前羔羊体尺各项指标在哺乳期间均表现出较为稳定的增长 （见图2），体重和体尺各项指标间均表现出良好的相关性 （相关系数均大于0.990）（见表4）。

2.2 生长曲线拟合结果

表1 3 种生长曲线模型

图1 哺乳期羔羊总体实测平均体重及周增重

表2 哺乳期羔羊生长情况 单位：kg

表3 哺乳期羔羊体尺变化 单位：cm

表4 体重体尺相关系数（Pearson）

由表5 可知， 对哺乳期放养云岭山羊羔羊体重进行拟合，3 种模型R2均大于0.99，说明用3 种曲线均能较好地拟合放养云岭山羊羔羊从初生到断奶期间的生长发育规律；Logistic 模型拐点周龄和拐点体重在3 种模型中均是最大的，公、母羔拐点体重均超过5 kg，体重的拐点出现在20 日龄以后。在拐点周龄后母羔生长发育速度比公羔缓慢 （见图3）。公羔和母羔拐点体重的差异表现为Logistic模型最小，Gompertz 模型次之，Von Bertalanffy 模型最大。在拐点周龄方面公、母羔的差异则表现为Von Bertalanffy 模型最大，Gompertz 模型次之，且拐点周龄相差1 周以上，Logistic 模型在拐点周龄方面受性别影响最小，差异低于1 周。初生体重对羔羊试验末体重有极显著（P<0.01）影响（见表6），通过相关性系数皮尔逊值可知， 初生体重对公羔和母羔试验末体重影响有一定差异， 对公羔的影响比母羔大。

图2 体尺变化曲线

表5 拟合曲线参数

3 讨论

在研究羊的生长发育规律时用非线性模型拟合生长曲线的报道较多，主要有Gompertz、Brody、Richards、Logistic、Polynomial 和Von Bertalanffy 等非线性模型［8-9］。 而Von Bertalanffy、Gompertz 和Logistic 模型精确性好、拟合度高，是目前采用较多的生长曲线非线性模型， 主要用于拟合畜禽生长发育规律［10］。 该研究运用Von Bertalanffy、Gompertz 和Logistic 3 种非线性模型均能较好地模拟放养云岭山羊哺乳期羔羊的生长曲线， 拟合度均在0.99 以上；Gompertz 和Von Bertalanffy 模型对体重拟合度均达到0.999，拟合效果优于Logistic，说明Gompertz 和Von Bertalanffy 模型拟合结果更接近实际测量值。研究结果表明，Gompertz 和Von Bertalanffy 模型能很好地模拟云岭山羊哺乳期羔羊生长发育规律。 该研究结果为在云岭山羊不同发育阶段采取相应的饲养措施， 最大限度地发挥放养云岭山羊生长潜力与提高经济效益提供帮助。Gompertz 和Von Bertalanffy 模型曲线拟合结果表明，公羔的生长发育曲线拐点时间分别为3.212周和2.862 周，拐点体重分别为5.691、5.352 kg，母羔生长发育曲线拐点时间分别为1.927 周和1.325 周，拐点体重分别为4.532、3.960 kg，公、母羔拐点周龄和拐点体重存在差异，公羔平均3.037周龄体重增长速度达到最大，而母羔平均在1.626周龄体重增长速度达到最大。 公羔拐点体重比母羔高，而公羔与母羔相比拐点日龄也延迟，表明公羔达到最大生长速度的时间要晚于母羔。 拐点体重与拐点日龄的变化与袁飞等［10］研究发现的早期南江黄羊生长发育结果相似。 在出生后进入快速生长期会导致母羊泌乳量无法满足羔羊营养需要，这对于出生后生长快速期羔羊极为不利，应该给予及时补饲， 这样才能稳定羔羊后期的生产成绩。 在整个哺乳期，体重增长曲线接近线性，且羔羊生长速率随着时间的推移降低， 在出生第1 周内公羔、母羔均出现了最大周增重，第2 周到第3周以及第5 周到第6 周增重下降最快， 这与陕北绒山羊哺乳期体重生长变化相似［11］。 在8 周龄时公羔体重高于母羔， 与马鸣山种山羊场羔羊的生长结果一致，但在初生体重上与叶瑞卿等［12］人的试验结果有较大差异， 这可能是该次试验未对母羊进行补饲和不同地区环境间的差异导致的结果。此外，该次试验发现羔羊的初生体重与断奶体重呈正相关，吴庶青等［13］对妊娠后期苏尼特羊进行限制饲养研究发现， 低能量水平饲养会降低羔羊初生体重， 结合该次试验结果说明在母羊妊娠期间应适当补饲， 通过增加母羊能量摄入提高羔羊的初生体重， 进而提高羔羊后期生长性能和经济效益。

图3 羔羊生长曲线拟合

表6 羔羊初生体重与试验末体重相关性分析 单位：kg

4 结论

Logistic、Gompertz 和Von Bertalanffy 3 种非线性模型对放养云岭山羊羔羊生长曲线具有良好的拟合效果，Gompertz 和Von Bertalanffy 模型对体重的拟合度最高，拟合结果更接近实际测量值。该研究结果为今后对哺乳期放养云岭山羊羔羊进行科学补饲提供依据。