连城白鸭青壳新品系生长曲线拟合及分析

缪中纬 梁阿政 庄晓东 李 丽 朱志明 辛清武 章琳俐 郑嫩珠＊

(1.福建省农业科学院畜牧兽医研究所 福州 350013；2.福建省畜禽遗传育种重点实验室 福州 350013；3.漳州昌龙农牧有限公司 福建漳州 360000)

福建省拥有4 个国家级蛋鸭品种资源， 其中连城白鸭因其特殊的外观特征而全国知名，拥有“唯一的药用鸭”的美誉[1]。 原种连城白鸭体型偏小，对环境敏感，受应激产蛋量陡降，蛋壳颜色为纯白色，蛋重59 g 左右，不利于蛋品加工。 福建省农业科学院畜牧兽医研究所通过多年的努力， 以连城白鸭为选育素材，培育出了连城白鸭青壳新品系。该鸭不仅省料，产蛋持久性强，而且蛋壳颜色为青色，蛋重适中（62～68 g），受到广大养殖户和消费者的青睐。 生长拟合曲线模型可以拟合生长时间与体重变化之间的规律，研究两者之间的关系，有助于了解连城白鸭青壳新品系前期的生长性能和养殖成本的控制。 禽类生长规律的研究是禽类育种的基础， 目前国内外学者经常采用生长拟合曲线模型来研究禽类的生长规律[2-3]。 2000 年至今，国内外学者已经初步建立了一些禽类的产蛋性能、 生长性能、 肌肉发育等的拟合曲线模型[4-9]，目前主要有Logistic、Gompertz 和Von Bertalanffy[10]三种生长拟合曲线模型在国内外禽类育种研究方面采用较多，3 种曲线的参数估计值可以反映不同禽类品种在生长率、 拐点体重和拐点周龄的差异[11]。 这些参数的差异性可直观地比较不同禽类品种之间的差异， 更有利于了解不同禽类品种之间的生长特点和遗传规律[12]。 连城白鸭青壳新品系是一个新培育的品系， 相关的生长发育数据和生长拟合曲线模型还未见报道，研究其生长性能，对于连城白鸭青壳新品系的培育成本、 实际生产过程中饲料的调控和蛋鸭新品系的推广具有重要作用。 本试验以连城白鸭青壳新品系母鸭出雏（0 周龄）称重后，每隔7 d 为一个称重周期，称取0～17 周的实际体重， 通过SPSS 22 软件建立Logistic、Gompertz 和Von Bertalanffy 三种不同生长拟合曲线模型对连城白鸭青壳新品系进行分析和研究。

1 材料与方法

1.1 试验动物 试验用连城白鸭青壳新品系种蛋由福建省农业科学院畜牧兽医研究所家禽育种研究中心提供，孵化并饲养于漳州昌龙农牧有限公司，试验采用同批次出雏的连城白鸭青壳新品系母鸭1 200 羽。

1.2 饲养管理和体重测定 试验鸭出雏（0 周龄）后通过人工性别判定， 对母雏进行脚号标记并称重，0～17 周试验期间，每隔7 d，清晨空腹条件下随机称重150 羽试验鸭，计算其平均体重。试验鸭采取统一的饲养管理条件，网上平养，自由采食，常规免疫程序。饲料由福建省漳州市华龙饲料有限公司提供，不同饲养阶段的饲料营养成分见表1。

表1 不同时期饲料的营养成分

1.3 生长拟合曲线模型表达式 连城白鸭青壳新品系母鸭出雏称重后，每7 d 为一个称重周期，称取0～17 周鸭的实际体重， 通过SPSS 22 软件建立Logistic、Gompertz 和Von Bertalanffy 三种不同生长拟合曲线模型对连城白鸭青壳新品系的生长曲线进行拟合。 各模型表达式[13-15]见表2。

表2 用于拟合的三种非线性模型

1.4 生长拟合曲线模型处理及分析方法 参考戴国俊等[16]介绍的方法，利用Excel2016 进行实际观测值的数据整理和绘制曲线图， 采用统计学软件SPSS22 对试验数据进行分析，分析结果以平均值表示。通过SPSS22 软件，算出3 个模型各自的A 值、B值和K 值，通过各表达式的A 值、B 值和K 值，计算拐点体重和拐点体重参数， 并通过计算每周龄的生长估值，建立生长曲线模型，根据拟合度（R2）高低来评价生长曲线模型拟合度。

2 结果与分析

2.1 连城白鸭青壳新品系母鸭的实际观测生长曲线分析 从图1-图2 可知，连城白鸭青壳新品系母鸭3 周龄前平均体重从35.1 g 增长到348.7 g，平均周增重104.53 g；3 周龄至7 周龄， 平均体重从348.7 g 增长到1 158.9 g，平均周增重202.55 g；7 周龄至10 周龄， 平均体重从1158.9 g 增长到1 459.4 g，平均周增重100.17 g；10 周龄至17 周龄，平均体重从1 459.4 g 增长到1 520.7 g，平均周增重8.76 g。

图1 连城白鸭青壳新品系体重累积生长曲线

图2 连城白鸭青壳新品系周增重曲线

2.2 连城白鸭青壳新品系母鸭生长拟合曲线模型分析 由表3 可知，Gompertz、Logistic 和Von Bertalanffy 三种拟合曲线模型的极限体重 （A） 分别为1 542.032 g、1 510.456 g 和599.036 g， 拐点周龄分别3.97 周、4.99 周和8.26 周， 拐点体重为别为567.28 g、755.23 g 和177.49 g，拟合度（R2）值分别为0.997、0.998 和0.954。

表3 3 种模型连城白鸭青壳新品系生长曲线参数估计值和拟合度

从图3 来看，Gompertz 和Logistic 两个生长拟合曲线模型均能很好地拟合连城白鸭青壳新品系母鸭的生长发育过程， 而Von Bertalanffy 生长拟合曲线模型并不能很好地拟合青壳新品系母鸭的生长发育过程，与实际观测生长曲线出入较大。表4 结果表明：除0 周龄和1 周龄外，Gompertz 模型拟合连城白鸭青壳新品系母鸭其余周龄与实际观测值基本符合。 除0 周龄外，Logistic 模型拟合连城白鸭青壳新品系母鸭其余周龄与实际观测值基本符合。 Logistic模型拟合度（R2=0.998）是三者中最高的，也说明Logistic 的拟合曲线模型效果最好。

图3 连城白鸭青壳新品系生长曲线与三种生长拟合曲线比较

表4 体重实际观测值与拟合曲线估计值比较

3 讨 论

本试验发现，Gompertz 和Logistic 两个生长拟合曲线模型均能很好地拟合连城白鸭青壳新品系母鸭的生长发育过程， 而Von Bertalanffy 生长拟合曲线模型与实际观测生长曲线偏差较大，综合比较后，Logistic 模型拟合度（R2=0.998）拟合效果更好。 相关文献表明：唐雪峰等[17]、陈清等[18]、张超[19]、王晋[4]在研究淮南麻鸭、籽鹅、金陵白鸭、连城白鸭、北京鸭时认为最优生长拟合曲线模型为Logistic；张超[19]、肖慎华等[20]在研究莆田黑鸭、樱桃谷鸭时认为最优生长拟合曲线模型为Gompertz；黄骏腾等[21]、苗小猛等[22]、王聪等[23]、吉文林等[24]、丁贤群等[25]在研究雷州黑鸭、三穗鸭、高邮鸭、金定鸭和绍兴鸭时认为最优生长拟合曲线模型为Von Bertalanffy；吴艳等[26]在研究荆江蛋鸭青壳系时认为最优生长拟合曲线模型为Gompertz 和Von Bertalanffy。 说明不同的水禽品种，在不同季节、不同营养水平和不同饲养管理条件下，其生长曲线的最优拟合模型不尽相同， 也说明我国不同地域水禽资源的多样性。

本试验Logistic 生长拟合曲线模型结果表明：连城白鸭青壳系母鸭极限体重为1 510.456 g，拐点体重为755.23 g，拐点周龄为4.99 周。 张超[19]、黄骏腾等[21]、吴艳等[26]研究结果表明：金陵白鸭、连城白鸭、莆田黑鸭、雷州黑鸭、荆江蛋鸭青壳系的拐点周龄分别为4.1～5.3 周、3.7～5 周、3.8～5.1 周、4.66 周、4.97～5.47 周，与本试验结果接近，说明虽然3 种生长拟合曲线模型不同， 但是最优生长拟合曲线模型得到的参数具有可比性。 本试验3 个生长拟合曲线模型拟合度（R2）的高低，均是根据实际观测生长值估算出3 个生长拟合曲线模型的A 值、B 值和K值，然后根据方程式算出估算值，该方法使得模型拟合度（R2）不同程度偏高，在试验过程中由于采用常规的免疫程序， 部分疫苗对于鸭的生长性能有一定的应激性，影响其生长体重。

通过生长拟合曲线模型的建立， 可以更好地了解连城白鸭青壳新品系母鸭各个阶段的生长发育情况，合理控制培育成本，为提高连城白鸭青壳新品系母鸭市场竞争力奠定基础。