黃俻華，尹華東

（1.賀州市水產(chǎn)畜牧獸醫(yī)局，廣西賀州 542899；2.四川農(nóng)業(yè)大學動物科技學院，成都 611130）

彭縣黃雞產(chǎn)蛋率和累計產(chǎn)蛋數(shù)曲線擬合分析

黃俻華1，尹華東2*

（1.賀州市水產(chǎn)畜牧獸醫(yī)局，廣西賀州 542899；2.四川農(nóng)業(yè)大學動物科技學院，成都 611130）

【目的】探討彭縣黃雞產(chǎn)蛋率和產(chǎn)蛋數(shù)變化規(guī)律，為提高其產(chǎn)蛋性能提供理論依據(jù)。【方法】統(tǒng)計彭縣黃雞選育群體21～43周齡的每周產(chǎn)蛋率和累計產(chǎn)蛋數(shù)，采用伍德、分室和楊寧模型對產(chǎn)蛋率曲線進行擬合分析，以及采用Logistic、Gompertz和Von Bertalanffy模型對累計產(chǎn)蛋數(shù)曲線進行擬合分析，并與實際生產(chǎn)值進行比較?！窘Y(jié)果】楊寧模型的擬合度最高（R2=0.968），分室模型的次之（R2=0.908），伍德模型的效果最差（R2=0.787），楊寧模型的預測值與實際觀察值最接近；Von Bertalanffy模型對累積產(chǎn)蛋數(shù)曲線的擬合度最高（R2=0.998），預測的拐點產(chǎn)蛋周齡為13.386周，拐點產(chǎn)蛋數(shù)為48.534枚，預計累積產(chǎn)蛋數(shù)85.13枚，比實際低1.31枚?！窘Y(jié)論】楊寧模型和Von Bertalanffy模型分別對彭縣黃雞的產(chǎn)蛋率和累計產(chǎn)蛋數(shù)擬合效果最佳，可用于彭縣黃雞產(chǎn)蛋規(guī)律的預測和分析，可為彭縣黃雞的選育和飼養(yǎng)管理提供參考。

彭縣黃雞；數(shù)學模型；產(chǎn)蛋率；累計產(chǎn)蛋數(shù)；曲線擬合

彭縣黃雞為我國優(yōu)秀的地方雞種之一，入選了1989年版的《中國家禽品種志》，為廣泛分布于成都平原和丘陵地區(qū)黃雞類群中的代表型。20世紀80—90年代，由于國外高產(chǎn)雞品種的大量引入，彭縣黃雞養(yǎng)殖數(shù)量驟減，種質(zhì)資源快速萎縮，瀕臨滅絕。隨著國家對畜禽遺傳種質(zhì)資源的不斷重視，成都市農(nóng)委于2006年在全市范圍內(nèi)開展了彭縣黃雞種質(zhì)資源調(diào)查，并收集活體交由專門公司進行保種。經(jīng)過近年來的努力，彭縣黃雞群里數(shù)量也得到較大程度的改善，保種成效顯著，并逐步進行開發(fā)利用，但產(chǎn)蛋較低，實際生產(chǎn)效果不佳[1-2]。

產(chǎn)蛋量是家禽生產(chǎn)的重要數(shù)量性狀，受飼養(yǎng)環(huán)境、營養(yǎng)狀況、管理水平等多種因素的影響[3]。產(chǎn)蛋曲線不僅能顯示產(chǎn)蛋率的高低，而且可以反映不同雞生產(chǎn)時期的飼養(yǎng)管理水平，用于指導生產(chǎn)，還可以為準確評估產(chǎn)蛋量性狀的育種值提供參考依據(jù)[4]。產(chǎn)蛋曲線數(shù)學模型表示產(chǎn)蛋與時間的函數(shù)，是從動態(tài)的角度研究產(chǎn)蛋率隨時間而變化的規(guī)律[5]。通過模型能夠較好地描述雞群產(chǎn)蛋特點，以早期產(chǎn)蛋數(shù)據(jù)預測全期產(chǎn)蛋量，提高種雞早期選擇的準確性，同時對雞場管理和生產(chǎn)狀況的系統(tǒng)分析也具有重要的指導意義[6]。產(chǎn)蛋率曲線數(shù)學模型受到了國外內(nèi)學者的廣泛關(guān)注，先后提出了分室模型[7]、伍德模型[8]、楊寧模型[9]、三次或四次方程[10]等。然而，產(chǎn)蛋率曲線雖可以描繪不同階段的產(chǎn)蛋特點，但不能預測整個生產(chǎn)期的總產(chǎn)蛋數(shù)以及產(chǎn)蛋拐點。為了解決這一問題，早期學者采用Gompertz、Logistic以及Von Bertalanffy等非線性模型對產(chǎn)蛋總數(shù)進行擬合，并取得了良好的效果[11-12]。

本研究將分別采用分室模型，伍德模型，楊寧模型對彭縣黃雞的產(chǎn)蛋率曲線以及采用Gompertz，Logistic和Von Bertalanffy模型對累計產(chǎn)蛋數(shù)進行擬合分析，并與實際統(tǒng)計值進行比較，以期掌握彭縣黃雞的產(chǎn)蛋變化規(guī)律，為提高其產(chǎn)蛋性能提供理論依據(jù)，也為彭縣黃雞的飼養(yǎng)管理和生產(chǎn)狀況分析提供指導。

1 材料和方法

1.1 試驗動物與飼養(yǎng)管理

以四川彭州天華禽業(yè)有限公司彭縣黃雞保種場飼養(yǎng)的彭縣黃雞為材料，選取2000只個體，同一棟雞舍單籠飼養(yǎng)，自動喂料，自動飲水，人工撿蛋，并按照《彭縣黃雞種雞飼養(yǎng)管理規(guī)程》進行日常管理。

1.2 產(chǎn)蛋數(shù)據(jù)記錄

每只雞從產(chǎn)第一個蛋開始，按個體記錄每天產(chǎn)蛋情況至300日齡，統(tǒng)計群體的300日齡累積產(chǎn)蛋數(shù)，并計算21至43周的周產(chǎn)蛋率。

1.3 統(tǒng)計與分析

1.3.1 產(chǎn)蛋率曲線擬合模型

伍德模型：y（t）=atbe-ct

分室模型：y（t）=a[1-e-c（t-d）]e-bt

楊寧模型：y（t）=ae-bt/[1+e-c（t-d）]

以上模型中，各式中a、b、c和d為待定參數(shù)，e為自然對數(shù)底，t為產(chǎn)蛋周齡，y（t）為第 t周的產(chǎn)蛋率。利用SPSS 20.0計算參數(shù)a、b、c、d的最優(yōu)值以及擬合度R2。

1.3.2 累計產(chǎn)蛋數(shù)擬合模型

Logistic：y（t）=A/(1+Be-kt)

Gompertz：y（t）=Ae-Be（-kt）

Von Bertalanffy：y（t）=A(1-Be-kt)3

各式中A為極限產(chǎn)蛋數(shù)、B為參數(shù)、K為接近極限速度、t為產(chǎn)蛋周齡（t=1表示雞為第21周齡）、y（t）為第 t周累計產(chǎn)蛋數(shù)。拐點周齡分別為：（lnB）/k、（lnB）/k和（ln3*B）/k；拐點產(chǎn)蛋數(shù)分別為：A/2、A/e和8A/27。用SPSS 20.0計算模型參數(shù)A、B、k最優(yōu)值以及擬合度R2。

2 結(jié)果與分析

2.1 彭縣黃雞產(chǎn)蛋率曲線擬合

根據(jù)彭縣黃雞不同周齡的產(chǎn)蛋率，采用伍德模型、分室模式以及楊寧模型對其產(chǎn)蛋率曲線進行了擬合，結(jié)果見表1。將估計的參數(shù)帶入模型，分別得到三個模型的擬合方程，其中伍德模型：y（t）=2.2×10-6t7.129e-0.321t；分室模型：y（t）=87.129[1-e-0.705（t-18.417）]e-0.065t；楊寧模型：y（t）=1.051e-0.016t/[1+e-0.968（t-23.169）]。三種模型的F檢驗結(jié)果都達到差異極顯著水平（P＜0.01），表明這三種模型均可用于彭縣黃雞產(chǎn)蛋率曲線的擬合。然而，三種模型的擬合效果存在差異，其中楊寧模型的擬合度最高（R2=0.968），分室模型的次之（R2=0.908），伍德模型的效果最差（R2=0.787）。

表1 不同模型對產(chǎn)蛋率曲線擬合的參數(shù)估計Table 1 Parameters estimation of different models for the fitting curves of laying rate

2.2 產(chǎn)蛋率觀察值與實際值比較

由表2可知，彭縣黃雞產(chǎn)蛋率較低，27至31周齡為產(chǎn)蛋高峰期，產(chǎn)蛋率均超過65%，但最高峰期產(chǎn)蛋率也未超過70%，其中第28周齡產(chǎn)蛋率最高，達到68.59%。從三個模型的擬合值來看，伍德模型的擬合效果最差，最高產(chǎn)蛋率僅為59.12%；分室模型擬合產(chǎn)蛋率最高值為67.34%，但出現(xiàn)在第30周齡；楊寧模型的擬合效果最佳，最高峰產(chǎn)蛋率為68.97%，比實際觀察值略高，也出現(xiàn)在第28周齡。由圖1的曲線也可以看出，楊寧模型擬合出來的曲線與實際的產(chǎn)蛋率曲線吻合度最高；分室模型擬合曲線前期較好，但產(chǎn)蛋后期比實際曲線略低；伍德模型擬合的曲線產(chǎn)蛋率比實際提前，但總體產(chǎn)蛋率比實際低。

圖1 不同模型擬合產(chǎn)蛋率曲線與實際產(chǎn)蛋率曲線比較Figure 1 Comparison of different fitting curves and actual curve of laying rate

2.3 彭縣黃雞累計產(chǎn)蛋數(shù)曲線擬合

根據(jù)彭縣黃雞群體各周齡的累計產(chǎn)蛋情況，采用Logistic、Gompertz以及Von Bertalanffy三種模型對累積產(chǎn)蛋數(shù)曲線進行了擬合。由表3可知，三種非線性模型均可以對彭縣黃雞的累積產(chǎn)蛋數(shù)曲線進行擬合，且擬合度都超過0.98，其中以Von Bertalanffy模型的擬合度最高，達到了0.998。結(jié)合表4和圖2來看，不同模型擬合出的拐點周齡以及相對應(yīng)的拐點產(chǎn)蛋數(shù)與實際值相差不大；雖然在開產(chǎn)前期，不同模型的擬合值均大于實際，但擬合的總產(chǎn)蛋數(shù)均比實際數(shù)低，其中Von Bertalanffy模型的擬合數(shù)據(jù)與實際值最為接近，僅比實際低1.31枚。

表2 產(chǎn)蛋率觀察值和擬合值的比較Table 2 Comparison of observed values and fitted values of laying rate

3 討論與結(jié)論

3.1 討論

一個成熟的家禽品種，其產(chǎn)蛋具有規(guī)律性，反映在整個產(chǎn)蛋期內(nèi)產(chǎn)蛋率的變化具有一定的模式，雞群開產(chǎn)后，在最初5～6周內(nèi)產(chǎn)蛋率迅速上升，高峰期過后平穩(wěn)地下降至產(chǎn)蛋末期[13]。將每周的產(chǎn)蛋率連接起來形成曲線，則可直觀反映出該品種的產(chǎn)蛋性能，通過與標準產(chǎn)蛋曲線的比對，找出實際生產(chǎn)中存在的問題，可有效地指導生產(chǎn)實踐。通過數(shù)學模型來擬合產(chǎn)蛋曲線最早是由McMillan等在1970年提出的分室模型[7]，隨后McNally提出了伍德模型及其修正式[8]，1994年我國的楊寧教授在對產(chǎn)蛋規(guī)律性變化進行深入分析后提出了楊寧模型[9]，這三種模型成為目前應(yīng)用范圍最廣的產(chǎn)蛋曲線擬合模型。

表3 累積產(chǎn)蛋數(shù)曲線模型擬合參數(shù)估計值Table 3 Estimation parameters of fitted curve of total egg number

表4 實際累積產(chǎn)蛋數(shù)與不同模型擬合值的比較Table 4 Comparison of observed values and fitted values of total egg numbers

圖2 不同模型擬合累計產(chǎn)蛋數(shù)與實際產(chǎn)蛋數(shù)實測曲線比較Figure 2 Comparison of different fitting curves and actual curve of total egg number

本研究采用分室模型、伍德模型以及楊寧模型對彭縣黃雞群選育群體的產(chǎn)蛋率曲線進行了擬合分析，結(jié)果發(fā)現(xiàn)楊寧模型的擬合效果最佳，其產(chǎn)蛋率擬合曲線與實際產(chǎn)蛋率曲線的吻合度最高，這與大多數(shù)的研究一致，如：顧玉萍等研究京海黃雞[14]，蔣松等研究二郎山山地雞[15]，李俊營等研究淮南麻黃雞[16]，董晶等研究信宜懷鄉(xiāng)雞[17]，但馮敏山等[18]和李鴻飛[19]的研究卻發(fā)現(xiàn)分室模型對AA肉雞和豫粉1號商品代蛋雞的產(chǎn)蛋曲線模擬效果最佳，因此說明了盡管楊寧模型為綜合前人研究提出的改進型模型，也并不是適用于所有的品種，但也不排除可能由于飼養(yǎng)管理過程不當造成的產(chǎn)蛋異常，從而導致模型差異。

在家禽飼養(yǎng)過程中除了要關(guān)注不同周齡產(chǎn)蛋率以外，累計產(chǎn)蛋數(shù)以及產(chǎn)蛋拐點也必須掌握。眾多的研究已經(jīng)表明，Logistic、Gompertz以及 Von Bertallanffy模型是對家禽生長發(fā)育規(guī)律進行曲線擬合與分析的理想曲線模型，可動態(tài)地預測家禽的生長隨時間而變化的規(guī)律[20-22]。與此同時，研究者發(fā)現(xiàn)，這三種模型也適用于對家禽累計產(chǎn)蛋數(shù)曲線的擬合。為此，本研究采用Gompertz，Logistic以及Von Bertalanffy三種非線性模型對彭縣黃雞的累積產(chǎn)蛋數(shù)進行了擬合，并發(fā)現(xiàn)三種模型均適用于累積產(chǎn)蛋數(shù)的擬合，Von Bertalanffy模型的擬合數(shù)據(jù)與實際生產(chǎn)數(shù)據(jù)更為接近，這與王曉峰等[23]，李國輝等[24]，駱磊等[4]的研究結(jié)果一致。

從實際生產(chǎn)數(shù)據(jù)來看，彭縣黃雞的產(chǎn)蛋性能較低，除了通過持續(xù)不斷的品種選育，精細化的飼養(yǎng)管理也是提高產(chǎn)蛋性能的有效手段。在飼養(yǎng)過程中，可以采用楊寧模型和Von Bertalanffy模型對群體的產(chǎn)蛋率曲線和產(chǎn)蛋數(shù)曲線進行擬合，并與實際曲線進行比較，及時發(fā)現(xiàn)飼養(yǎng)中存在的問題，通過改變光照時間和調(diào)整喂料量等方法來避免外界環(huán)境因素導致的生產(chǎn)性能降低，影響選育效果。

3.2 結(jié)論

楊寧模型對彭縣黃雞產(chǎn)蛋率曲線的擬合效果最好，擬合度為0.968；Von Bertalanffy模擬累計產(chǎn)蛋數(shù)曲線的效果最佳，這兩種模型綜合對彭縣黃雞的產(chǎn)蛋規(guī)律進行預測和評估，可為彭縣黃雞的選育和生產(chǎn)實踐提供參考。

［1］韓順順，陸登宇，鄒云峰，等.彭縣黃雞體尺性狀與屠宰性能的典型相關(guān)分析［J］.中國家禽，2016，38（10）：7-10.

［2］劉子豪，韓順順，鄒云峰，等.彭縣黃雞不同產(chǎn)蛋日齡蛋品質(zhì)比較分析［J］.中國家禽，2016，38（14）：13-16.

［3］邱祥聘.家禽學［M］.成都：四川科學技術(shù)出版社，1993.

［4］駱磊，張龍，黃欽柯，等.地方蛋雞新品系（HS）基礎(chǔ)群產(chǎn)蛋率與累計產(chǎn)蛋數(shù)曲線擬合［J］.四川農(nóng)業(yè)大學學報，2013，31（4）：438-443.

［5］劉天飛，瞿浩，馬杰，等.商業(yè)化模式下的優(yōu)質(zhì)肉雞產(chǎn)蛋曲線數(shù)學模型研究［J］.四川農(nóng)業(yè)大學學報，2012，30（1）：110-113.

［6］賀興龍，王克華，杜房房，等.蘇禽黑羽青殼蛋雞產(chǎn)蛋性能及蛋品質(zhì)分析［J］.中國家禽，2013，35（9）：11-14.

［7］MCMILLAN I，F(xiàn)ITZ-EARLE M，ROBSON D S.Quantitative genetics of fertility.I.Lifetime egg production of Drosophila melanogaster—theoretica［lJ］.Genetics，1970，65（2）：349-353.

［8］MCNALLY D H.Mathematical model for poultry egg production［J］.Biometrics，1971，27（3）：735-738.

［9］YANG N，WU CHANGXIN，MCMILLAN I.New mathematical model of poultry egg production［J］.Poultry Science，1989，68（4）：476-481.

［10］杜文興.蛋雞產(chǎn)蛋率分布規(guī)律的研究［J］.畜牧與獸醫(yī)，1994，26（6）：262-263.

［11］GANG F Y，ZHEN Y S.A study on the growth curve of and maximum profit from layer-type cockerel chicks［J］.British Poultry Science，1997（38）：445-446.

［12］MIGNON-GRASTEAU S，BEAUMONT C，LE E，et al.Genetic parameters of growth curve parameters in male and female chickens［J］.British Poultry Science，1999（40）：44-51.

［13］鄧炳真，溫麗華，麥劍威，等.南海黃麻雞不同階段蛋品質(zhì)研究［J］.中國家禽，2016，38（5）：10-13.

［14］顧玉萍，侯啟瑞，王金玉，等.京海黃雞生長和產(chǎn)蛋率曲線的擬合［J］.中國家禽，2010，32（20）：26-29.

［15］蔣松，肖禮華，張龍，等.二郎山山地雞SD02和SD03品系一世代產(chǎn)蛋率曲線比較分析［J］.黑龍江畜牧獸醫(yī)，2012（7）：1-3.

［16］李俊營，詹凱，劉偉，等.籠養(yǎng)淮南麻黃雞產(chǎn)蛋曲線的數(shù)學模型研究［J］.福建農(nóng)業(yè)學報，2013，28（12）：1204-2107.

［17］董晶，杜炳旺，劉宜林，等.信宜懷鄉(xiāng)雞C品系第3世代產(chǎn)蛋率和累計產(chǎn)蛋數(shù)曲線擬合分析［J］.南方農(nóng)業(yè)學報，2017，48（1）：158-162.

［18］馮敏山，高山林，楊增芳，等.AA父母代肉雞產(chǎn)蛋曲線的數(shù)學模型分析［J］.中國家禽，2003，7（1）：20-21.

［19］李鴻飛.“豫粉1號”商品代蛋雞主要生產(chǎn)性能研究［D］.鄭州：河南農(nóng)業(yè)大學，2015.

［20］鄺智祥，王俐智，朱慶，等.朗德鵝早期生長發(fā)育模型研究［J］.四川農(nóng)業(yè)大學學報，2008，26（2）：191-196.

［21］付睿琦，趙桂蘋，劉冉冉，等.北京油雞生長曲線擬合與分析［J］.中國家禽，2014，36（3）：6-9.

［22］申志超，劉麗仙，榮華，等.鹽津烏骨雞生長曲線分析與擬合研究［J］.中國家禽，2016，38（15）：63-65.

［23］王曉峰，王成麗，施會強.“京海黃雞1號”配套系父母代產(chǎn)蛋率與累計產(chǎn)蛋數(shù)曲線擬合研究［J］.中國家禽，2015，37（18）：9-13.

［24］李國輝，韓威，張會永，等.蘆花雞H系生長及累計產(chǎn)蛋數(shù)曲線擬合分析［J］.中國家禽，2016，38（20）：60-62.

Curve Fitting Analysis for Laying Rate and Total Egg Number of Pengxian Yellow Chicken

HUANG Bei-hua1，YIN Hua-dong2*
（1.Hezhou Husbandry-Veterinary and Fishery Bureau，Hezhou 542899，Guangxi，China；2.College of Animal Science and Technology，Sichuan Agricultural University，Chengdu 611130，China)

【Objective】In order to provide theoretical basis for improving laying performance,the present study investigated the change rule of laying rate and egg production of Pengxian Yellow Chicken.【Method】Laying rate and total egg number of Pengxian Yellow Chicken from 21 to 43 weeks were collected.Three models，namely Wood，Yang-Ning and McMillan were adopted to fitting analysis on laying rate.Logistic，Von Bertalanffy and Compertz were employed in fitting analysis on total egg number.Additionally，the fitted values of laying rate and total egg number were compared to the actual observed values.【Results】The degree of fitting of Yang-Ning was the highest(R2=0.968)followed by McMillan(R2=0.908)and Wood(R2=0.787)，and the fitted value of Yang-Ning was the closest to the actual observed values.The degree of fitting of Von Bertalanffy for total eggproductionwasthehighest(R2=0.998)，and the inflection point of laying week and egg production were 13.386 and 48.534，respectively.The fitted total egg production was 85.13，and 1.31 less than the actual observed number.【Conclusion】Yang-Ning and Von Bertalanffy have the optimum fitting effects on laying rate and total egg production for Pengxian Yellow Chicken，which can be used for prediction and analysis of egg production rules，and can provide scientific references for breeding and feeding and management of Pengxian Yellow Chicken.

PengxianYellowChicken；mathematicalmodel；layingrate；totaleggproduction；curvefitting

S831

A

1000-2650（2017）03-0420-05

10.16036/j.issn.1000-2650.2017.03.021

2017-04-05

成都市農(nóng)業(yè)技術(shù)成果應(yīng)用示范項目（2015-NY01-00036-NC）。

黃俻華，畜牧師，主要從事畜牧生產(chǎn)工作。*責任作者：尹華東，副研究員，主要從事家禽遺傳育種工作，E-mail：yinhuadongsicau@163.com。

（本文審稿：劉益平；責任編輯：秦碧雯；英文編輯：劉益平）