閩清毛腳雞生長曲線擬合與分析

林秋娟 羅曉瑩 王文凱 曾彥欽 林福忠 陳 震 福州市農業(yè)科學研究所 福州 350018

閩清毛腳雞是福建省福州市稀有的肉蛋兼用型地方品種，因腿羽、脛羽、趾羽發(fā)達，看似“腳穿褲子”而得名。 該品種耐粗飼、肉質細嫩鮮美，具有較強的覓食力、適應性和抗病力[1]。2014 年閩清毛腳雞作為優(yōu)秀的地方原始品種被批準列入《國家級畜禽遺傳資源保護名錄》，具有較大的開發(fā)、保護和利用價值，但目前對該品種的深入研究尚處于起步階段， 國內外還沒有閩清毛腳雞早期發(fā)育規(guī)律的相關研究。

畜禽育種和生產中常運用Logistic、Gompertz 和Von Bertalanffy 三種模型研究不同環(huán)境下或不同品種畜禽生長發(fā)育的基本模式， 通過結合分析3 種模型曲線的參數估計值與畜禽的飼養(yǎng)相關數據， 預測畜禽各生長期的體重生長率、飼料消耗量、飼料轉化率，研究營養(yǎng)、溫度等外界因素對生長的影響，從而指導畜禽的飼養(yǎng)管理和選種選育工作[2]。 因此，本研究通過對閩清毛腳雞0～22 周齡時的體重情況進行曲線擬合， 以期為該種質資源的保護、 育種、 飼養(yǎng)管理及資源的合理開發(fā)與利用等方面的工作奠定基礎。

1 材料與方法

1.1 試驗材料 在福州市農業(yè)科學研究所孵化場入孵400 枚閩清毛腳雞種蛋， 從孵化出殼的雞群中選取健康、初生重均勻的200 羽雛雞（公母各半），出雛時戴上腳號進行個體編號。

1.2 飼養(yǎng)管理 6 周齡前網上平養(yǎng) （公母混合），6 周齡后地面、網上混合平養(yǎng)（公母分群），在同一階段內，全部雞群使用同樣的飼料和飼喂方法，整個過程中自由采食和飲水，并遵循常規(guī)的免疫流程，進行科學的免疫。

1.3 體重測定 用電子秤測量閩清毛腳雞初生重，在2、4、6、8、10、12、14、16、18、20、22 周的周末對試驗雞禁食12 h 后進行空腹體重測量并記錄，計算群體每兩周的平均體重。

1.4 統計方法

1.4.1 擬合曲線模型 運用Logistic、Gompertz 和Von Bertalanffy 三種非線性生長模型分析閩清毛腳雞的體重及生長發(fā)育特征。 具體模型表達式及相關參數見表1， 其中A 、B、k、t 分別代表極限生長量、常數尺度、瞬時相對生長率和周齡。

表1 用于擬合的3 種非線性數學模型

1.4.2 數據統計分析 利用EXCEL 軟件對雞群不同周齡的體重相應進行計算，并制作相應的圖表。使用SPSS27.0 軟件進行非線性模型擬合，使用逐次迭代的方式， 計算3 個模型參數A、 B、 k 的最優(yōu)估算值和擬合度指數R2，構建出生長模型，并對推算出極限體重、拐點周齡、拐點體重進行計算。當R2值趨于1 時，該模型的擬合效果較好。

2 結果與分析

2.1 閩清毛腳雞實際觀測生長發(fā)育情況分析 結合圖1、表2 可以看出，不同性別的閩清毛腳雞累積生長基本一致，均呈近“S”型，基本符合家禽正常的生長發(fā)育規(guī)律，0～4 周期間公雞和母雞體重相差不大，6 周齡以后公雞體重大于母雞。

圖1 0～22 周齡閩清毛腳雞體重累積生長曲線

表2 3 種生長曲線模型的擬合度和參數估計值

從圖2 分析可知， 公雞在0～14 周期間絕對生長呈現快速上升趨勢，14 周齡達到最高峰， 周平均增重為173.05 g，14 周以后絕對生長呈現逐漸下降趨勢； 母雞0～10 周期間絕對生長呈現快速上升趨勢，10 周齡達到最高峰， 周平均增重為149.97 g，10周以后絕對生長呈現逐漸下降趨勢。0～4 周期間，不同性別的閩清毛腳雞絕對生長差異不大， 公雞與母雞周平均增重都小于80 g，4 周齡后周平均增重大于80 g，6 周以后公雞的周平均增重明顯大于母雞。

圖2 0～22 周齡閩清毛腳雞平均周增重曲線

2.2 閩清毛腳雞生長曲線擬合分析 結合表2-表3 及圖3-圖4 分析，3 種生長模型的擬合度都非常高，均達到0.99 以上，基本符合閩清毛腳雞的實際生長。 以擬合度做比較，Logistic 曲線模型的擬合度＜Von Bertalanffy 曲線模型的擬合度＜Gompertz 曲線模型的擬合度，其中Logistic 和Von Bertalanffy 曲線模型0 周齡的擬合值與公母雞的實測值相差較大。3 種曲線模型中，Gompertz 曲線模型公雞、母雞極限生長體重即A 值分別為2 650.863 g、1 892.757 g，與閩清毛腳雞成年公雞 （2.2～2.6 kg）、 成年母雞（1.7～2.2 kg）體重接近。 綜合考慮擬合度、3 個模擬參數、拐點體重及拐點周齡等指標，Gompertz 曲線模型與實際曲線最接近，為最佳擬合模型。

圖3 三種模型擬合曲線和實際曲線（公雞）

圖4 三種模型擬合曲線和實際曲線（母雞）

表3 閩清毛腳雞實測值與3 種生長曲線模型擬合值

3 結論與討論

經對閩清毛腳雞的體重累積及平均周增重曲線進行研究， 結果表明： 閩清毛腳雞的公雞與母雞在4 周前體重無顯著性差別，因此，可以采取混合飼養(yǎng)方式。4 周后公雞的平均周增重顯著高于母雞，可采取分群養(yǎng)殖的方式， 并在此基礎上配備不同的飼料營養(yǎng)水平，以達到一定的飼養(yǎng)目標要求。

從生長曲線擬合結果可知，Logistic、Gompertz和Von Bertalanffy 三個生長模型的擬合度都達到0.99 以上， 說明均能很好地擬合閩清毛腳雞生長曲線，描述其早期體重的生長過程和規(guī)律。對比3 種生長模型，Gompertz 模型擬合度大于Logistic 和Von Bertalanffy 模型， 且模型的擬合值和實測值較為接近，Gompertz 模型能更好地模擬閩清毛腳雞生長發(fā)育規(guī)律。 與周玉剛等[3]、吳艷等[4]、柳儉強等[5]、王克文等[6]、林伯全等[7]分別研究得出的不同品種雞群最佳生長曲線結果相一致。 但與李桂賢[10]、馬猛等[11]、萬火福等[12]、譚玉文等[13]、姜自琴等[14]以Von Bertalanffy生長曲線模型最佳的研究結果不一致。與曹贊等[15]、趙月平等[16]、李健[17]以Logistic 擬合生長曲線最優(yōu)的結論不同。這可能與種雞品種、性別、遺傳背景、飼養(yǎng)管理、生長環(huán)境及其交互效應等的差異所造成，說明不同生長模型增長模型對各種家禽的適應性存在有差異。

將A、B、k 值代入方程，閩清毛腳雞公雞、母雞Gompertz 曲線方程分別為y=2650.863e-4.472e-0.157t、y=1892.757e-4.551e-0.175t，模型曲線擬合結果表明：閩清毛腳雞公雞、 母雞極限體重分別為2 650.863 g、1 892.757 g，拐點體重分別為975.20 g、693.31 g，拐點周齡分別為9.54 周、8.66 周， 即閩清毛腳雞在飼養(yǎng)過程中，9.54 周、8.66 周分別是公雞、 母雞體重增長速率從快速增長到緩慢增長的臨界值。 公雞的拐點周齡晚于母雞。 因此，在實際生產中，應加強9 周齡前的飼養(yǎng)管理工作，需提供高能、高蛋白日糧，滿足閩清毛腳雞生長需要，同時保證鈣、磷飼料供給，保證閩清毛腳雞骨骼生長發(fā)育良好， 最大限度地發(fā)揮其生長潛能。