太行鸡麻羽系产蛋率与累计产蛋量曲线拟合分析

刘玉芳，范业锴，郭思武，韩海银，杨俊琦，张伟峰

(河北工程大学生命科学与食品工程学院，河北 邯郸 056038)

太行鸡又称河北柴鸡，是河北省重要的地方品种，具有适应性强、体型小、肉质鲜美、蛋重小而营养价值高等特点[1]。太行鸡产于太行山脉，饲喂及驯化时间距今已有8 000多年，1960年开始与外来品种进行杂交培育新品系，直至2015年11月通过了国家畜禽遗传资源鉴定，被正式列入国家地方畜禽遗传资源名录[2]。在太行鸡的高产群体中500日龄平均产蛋数量在187枚左右，仍然处于高产品系培育与选育的阶段[3]。目前，太行鸡种鸡存栏量7万套，河北省饲养量在1 200万只左右，有17家1万只以上规模的饲养场[4]，结合不同的饲养技术正在进行品系和技术推广。由于太行鸡麻羽系养殖尚处于探索阶段，因此，开展太行鸡麻羽系产蛋曲线的研究，有利于了解太行鸡麻羽系的生产性能，对合理利用种质资源具有重要意义。

家禽产蛋曲线拟合分析常用的非线性模型包括：Wood、McMillan、Yang和Ali-Schaeffer 4种模型[5]；累计产蛋数曲线模拟一般使用生长曲线模型进行拟合，包括冈珀茨(Gompertz)、逻辑模型(Logistic)和贝塔朗菲(von Bertalanffy)3种模型[6]。本研究采用上述模型对太行鸡麻羽系第2世代蛋鸡进行产蛋率和累计产蛋曲线拟合分析，为太行鸡麻羽系选育及饲养技术推广工作提供理论基础。

1 材料与方法

1.1 数据来源

太行鸡麻羽系第2世代母鸡于河北省邢台沙河市金凯养殖公司在相同的笼养环境条件下进行饲养。育雏期、育成期笼内平养，140日龄见蛋后进行单笼饲养。0～4周龄、5～15周龄和16周龄后日粮蛋白质水平分别为19%、15%～16.5%和14%～17.5%。人工喂料，自由饮水，每天喷淋10～15 min。系谱记录完整、留有血液样本的太行鸡麻羽系母鸡3 500只，记录个体开产到500日龄产蛋量(21～73周龄)，统计21～73周龄产蛋率和累计产蛋量用于曲线拟合分析。

1.2 产蛋率曲线拟合模型

使用伍德(Wood)、分室模型(McMillan)、杨(Yang)和阿里-谢弗(Ali-Schaeffer)4种模型进行产蛋率曲线拟合，4种模型的公式见表1。

表1 产蛋率曲线的4种拟合模型

1.3 累计产蛋曲线拟合模型

使用Gompertz、Logistic和von Bertalanffy 3种模型对累计产蛋量曲线进行拟合。3种模型的公式见表2。比较不同的产蛋曲线，主要看初始产蛋数、拐点产蛋数、拐点周龄这3个参数，并进行产蛋率曲线差别分析。

表2 累计产蛋量曲线3种拟合模型

1.4 统计与分析

数据使用Excel软件进行初步处理，使用SPSS Statistics 19.0中的非线性回归程序进行曲线拟合，各模型拟合效果用拟合度R2评价。

2 结果与分析

2.1 太行鸡麻羽系产蛋率曲线拟合分析

太行鸡麻羽系产蛋率曲线拟合结果见表3、表4和图1。太行鸡麻羽系在19周(27～45周龄)的产蛋期间保持着较高的产蛋率，维持在60%以上；其中7周(29～35周龄)的产蛋率维持在70%以上的高水平。对4种模型的拟合度(R2)进行计算，有3种模型的拟合度均在0.9以上，其中Yang模型拟合效果最好，拟合度为0.989。Yang模型基本能够体现太行鸡麻羽系的产蛋率的变化趋势，可以用于太行鸡麻羽系产蛋率的曲线拟合，为笼养条件下太行鸡麻羽系的生产管理和产蛋性能的遗传改良提供理论基础。

表3 太行鸡麻羽系产蛋率曲线拟合结果

表4 产蛋率实测值与拟合曲线估计值

图1 太行鸡麻羽系产蛋率曲线拟合比较

2.2 太行鸡麻羽系累计产蛋量曲线拟合分析

太行鸡麻羽系累计产蛋量曲线拟合结果见表5、表6和图2。3种曲线模型拟合太行鸡麻羽系累计产蛋数效果较好，拟合度均在0.98以上，其中von Bertalanffy模型拟合度最高，为0.998。von Bertalanffy模型与实际曲线最为接近，累计产蛋量为185.73枚，更接近实际值。因此，von Bertalanffy模型为太行鸡麻羽系累计产蛋量最佳拟合模型，其预测的累计产蛋量的拐点周龄为38.37，拐点产蛋数为66.95枚。

表5 太行鸡麻羽系累计产蛋量曲线拟合结果

表6 累计产蛋量实测值与拟合曲线估计值

图2 太行鸡麻羽系累计产蛋量曲线拟合比较

3 讨论

在产蛋家禽的育种过程中，统计产蛋家禽产蛋率和累计产蛋量数据是一项重要的工作，对所得数据进行拟合和分析是家禽育种和生产的工作基础，而选择合适的统计模型动态的描述产蛋量和产蛋率随时间变化的规律对指导产蛋家禽饲养管理和育种实践具有重要作用。产蛋率曲线的研究中常用的模型有4种，分别为Wood模型、McMillan模型、Yang模型和Ali-Schaeffer模型，判断模型优劣的标准是拟合度，拟合度越好越接近实际观测值，R2值就越接近于1，反之则拟合效果差。本研究采用了上述4种常用产蛋率曲线数据模型，结果表明Wood模型估计结果与实际观测值差距较大，拟合度仅为0.662，与前人结果相似[7-9]。而孙艳发等[5]在龙岩山麻鸭高产系产蛋率曲线的研究中发现：Ali-Schaeffer模型的拟合效果最好，适合龙岩山麻鸭高产系产蛋曲线的拟合。本研究结果表明，采用上述4种模型对太行鸡麻羽系产蛋率曲线进行非线性回归拟合，根据模型参数、拟合度及预测数据与实际观测值的相近程度等指标综合分析，Yang模型的拟合效果最好，拟合度在0.989，是太行鸡麻羽系产蛋率曲线的最佳模型，模型预测曲线与实际产蛋率曲线拟合度较高，该结果与李岩等[10]和蒋松等[11]研究结果一致。产蛋性状受到多种因素的调控影响，尽管均为产蛋家禽，但其产蛋性能差异较大，造成这一现象的原因可能与产蛋家禽的物种(比如鸡和鸭，尽管均为产蛋家禽，其产蛋率曲线差异较大)、品种品系及饲养模式等相关。对太行鸡麻羽系产蛋率曲线的拟合映射出该系育种与饲养管理的特点，该系产蛋率与引进品种相比产蛋率较低，达到产蛋高峰产蛋率为73.77%，且维持时间较短，因此仍需在选种选育中进行高产品系的选择，以提高太行鸡的生产性能。

常见用于累计产蛋量拟合的模型有3种，包括Gompertz模型、Logistic模型和von Bertalanffy模型。累计产蛋量在产蛋家禽的研究中是一个非常重要的数量性状，反应了产蛋家禽的生产能力，遵循一定的规律和曲线分布。近年来，对地方品种家禽的产蛋量研究比较广泛，国内育种人员先后对京海黄鸡、二郎山山地鸡、丝羽乌骨鸡和坝上长尾鸡等蛋鸡品种的产蛋量情况进行了拟合分析，对提高蛋鸡产蛋性能有一定的指导意义[8,11-13]。对于太行鸡麻羽系的累计产蛋量拟合曲线研究中发现von Bertalanffy模型的拟合程度由于其他2种模型，拟合度为0.998，根据累计产蛋量拐点周龄数、拐点产蛋量、拟合度和与实际观测值之间的相近程度等指标综合分析，von Bertalanffy模型更适合太行鸡麻羽系预测整个产蛋周期的平均个体累计产蛋量，更好地预测整个产蛋过程，预测的拐点周龄为38.37，拐点产蛋量为66.95枚，与实际产蛋规律最为接近。研究表明，von Bertalanffy模型是当前使用较多的用以描述产蛋家禽累计产蛋量规律的曲线模型[14-15]，本研究结果与前人结果一致。然而该结果也反应出太行鸡产蛋周期内平均累计产蛋量较低，与高产品系相比仍存在一定的差距，需要进一步的选种选配，以提高产蛋周期的产蛋数量，改善太行鸡产蛋性能。

4 结论

Yang模型在太行鸡麻羽系的产蛋率曲线拟合效果最好，拟合度为0.989；而von Bertalanffy模型在太行鸡麻羽系的累计产蛋量曲线拟合效果最佳，拟合度为0.998，预测的累计产蛋量拐点周龄为38.37，拐点产蛋量为66.95枚。