宁都黄公鸡3 种生长曲线拟合的比较分析

谭玉文，朱学农，章 逸，张仁颠，刘慧琳，周玉强，周 敏*

（1.南昌师范学院生物技术研究所，江西南昌 330032；2.南昌师范学院数学与计算机科学系，江西南昌 330032）

宁都黄鸡是江西省著名的优良肉用黄鸡品种，也是我国宝贵的家禽品种资源，不仅肉质鲜嫩、营养丰富，还具有耐粗饲、抗应激能力强等特点，深受广大消费者和养殖户的青睐。近年来，江西省对其保种育种投入了大量的人力物力，也对相关课题进行了深入研究。生长曲线能描述畜禽体重随年龄的增长表现出S 型变化这一规律，它不仅能动态地反映畜禽的生长过程和生长规律[1]，还能指导畜禽的饲养管理，提高畜禽的生长效率，也能为加快畜禽品种的选育提供可靠的理论依据。本研究采用Logistic、Gompertz 和Bertalanffy 3 种非线性模型对宁都黄公鸡进行生长曲线拟合，为宁都黄公鸡的生产实践、品种选育和种质资源的保护提供理论依据。

1 材料与方法

1.1 实验动物 实验鸡由江西南师科技有限公司提供，选取同批次出雏的未经选育的宁都黄鸡，选择503 只1日龄健康公鸡笼养，1 周龄时随机对鸡只进行翅号编号，为防止翅号掉落，在5 周龄时进行脚号编号。

1.2 饲养管理 实验于2018 年4 月中旬至8 月上旬在相同的饲养条件下进行，鸡只自由采食和饮水，饲粮参照美国NRC 饲养标准并根据我国《鸡饲养标准》进行配制。按照常规的免疫接种，实施科学的规范管理，定时观察鸡群健康状况并做好日常记录。

1.3 实验方法 分别在0～16 周各周进行空腹称重，统计数据，计算平均体重，并参照《NY/T 83-2004 家禽生产性能名词术语和度量统计方法》计算相对生长率和绝对生长率。

1.4 生长曲线的拟合 采用Logistic、Gompertz 和Bertalanffy 3 种模型对实验鸡只生长情况进行拟合，模型表达式见表1。

1.5 统计分析 采用Excel 2013 进行数据整理计算，并绘制曲线图，运用SPSS 22.0 进行生长曲线的拟合，算出模型参数的最优估计值 A、B、K，建立生长模型，算出拐点体重、拐点周龄和最大周增重，根据拟合度（R2）和残差平方和（E）预测初生重与实测体重的差距评价生长模型。

2 结果与分析

2.1 宁都黄公鸡生长情况 由图1 可以看出，宁都黄公鸡体重随年龄增加而增加，在 1～3 周龄前生长相对缓慢，4 周龄后生长速度逐渐加快，6～7 周龄生长速度最快、周增重也最高；另外11～12 周龄体重上升明显，14 周龄后体重增加缓慢，生长强度减缓，体重趋于稳定。

表1 3 种常见非线性生长曲线模型

图1 宁都黄公鸡实际生长曲线和周增重曲线

2.2 宁都黄公鸡3 种生长曲线的比较 由图2 可知，宁都黄公鸡3 种曲线模型估计的最大绝对生长率均出现在第8 周，8 周以后都出现不同程度的下降。

图2 宁都黄公鸡3 种曲线模型的绝对生长率

由图3 可知，相对生长率的3 种模型曲线基本相似：前6 周，Bertalanffy 模型的相对生长率高于其他2 种模型，6 周以后3 种模型曲线基本趋于一致。这一结果表明6 周龄之前应当采取合理的配方满足宁都黄公鸡的营养需要。

2.3 宁都黄公鸡不同生长曲线拟合与分析 由表2 可知，3 种模型拟合曲线的拟合度R2均大于0.985，说明3 种模型均适合拟合宁都黄公鸡的生长情况。其中，Bertalanffy 模型拟合度最高且残差平方和最小，且Bertalanffy 模型初生重和各周拟合值与实测值最为接近（图4、表3）。综合以上数据，宁都黄公鸡生长曲线的拟合以Bertalanffy 模型最佳。

图3 宁都黄公鸡3 种曲线模型的相对生长率

图4 宁都黄公鸡实际生长曲线与3 种拟合曲线的比较

表2 0～16 周宁都黄公鸡3 种生长拟合曲线参数估计值与拟合度

表3 3 种模型预测值与实测值之间的比较

3 讨 论

拐点周龄是家禽生长曲线拟合的重要参数，动物的生长环境、营养、品种、性别及其交互效应都会影响模型估计参数及拐点[2]。宁都黄公鸡拐点周龄较早，本研究中Bertalanffy 模型宁都黄公鸡拟合的拐点周龄是7.32 周，3 种模型预测的拐点周龄平均在7.93 周龄。结合绝对生长率和相对生长率的结果，在饲养过程为满足宁都黄公鸡生长发育的需求，应特别注意8 周前宁都黄公鸡饲料营养的全价性和平衡性。在用模型进行曲线拟合时发现，宁都黄公鸡在8～10 周龄的周增重明显下降，这主要是因为鸡群正处于夏季高温和降低饲养密度换笼造成鸡群产生应激所致。

本研究中，3 种非线性模型均能很好地模拟宁都黄公鸡的生长曲线，其各方程的拟合度都达到了0.985 以上，但综合R2、残差平方和以及模拟的体重与实际观察体重这3 个因素，Bertalanffy 模型拟合效果要优于Logistic 模型和Gompertz 模型，这与蒋世远等[3-4]、何兰花等[5]、牛建芹等[6]的研究结果相一致。由于本实验采用的是未经选育的宁都黄公鸡，生长发育较为迟缓，这也与杨海明等[7]提到的Bertalanffy 模型更适用于增重相对较迟缓的品种这一结论相符合。在众多学者研究中，Logistic、Gompertz 和Bertalanffy 这3 种模型对家禽生长曲线的拟合效果也不尽相同。刘海斌等[8]、张乐超等[9]、白优等[10]的研究结果表明，Logistic 模型优于Gompertz、Bertalanffy 模型；而刘双艳等[11]、张华杰等[12]、余春林等[13]研究则认为Gompertz 模型更好，这可能是由于品种、性别、饲养管理、亲缘关系等原因造成。

4 结 论

本实验条件下，Logistic、Gompertz 和Bertalanffy模型都能较好地拟合宁都黄公鸡的生长曲线，其中Bertalanffy 模型的拟合效果最佳，宁都黄公鸡的模型方程为Yt=1336.053（1-0.714e-0.104t）3，拟合度为0.994，拐点周龄为7.32 周龄，拐点体重为395.57 g。