蔣世遠(yuǎn) ,陳小萍 ,陳福祥 ,張跟喜
(1.廣西鳳翔集團(tuán)股份有限公司,廣西 北海 536100;2.揚(yáng)州大學(xué)動(dòng)物科學(xué)與技術(shù)學(xué)院,江蘇 揚(yáng)州 225009)
動(dòng)物生長(zhǎng)過程中,體重的增長(zhǎng)一般呈現(xiàn)非線性的規(guī)律,因此,動(dòng)物生長(zhǎng)過程中體重的增長(zhǎng)規(guī)律可以用非線性的數(shù)學(xué)模型進(jìn)行擬合[1]。目前,常用于動(dòng)物生長(zhǎng)發(fā)育規(guī)律研究的非線性模型一般有3種:Logistic、Gompertz和 Bertalanffy,這 3種模型擬合的結(jié)果普遍較好,準(zhǔn)確度也較高,對(duì)動(dòng)物的生長(zhǎng)發(fā)育描述較為準(zhǔn)確[2]。 該研究采用 Logistic、Gompertz、Bertalanffy 3種動(dòng)物生長(zhǎng)曲線模型,分別對(duì)桂香雞母雞的生長(zhǎng)曲線進(jìn)行擬合分析,通過該研究能對(duì)桂香雞的生長(zhǎng)過程有一個(gè)動(dòng)態(tài)的了解,從而為桂香雞制定合理的飼養(yǎng)和管理方案,使其生長(zhǎng)能力充分發(fā)揮,進(jìn)而降低養(yǎng)殖成本,加速桂香雞的開發(fā)與推廣。
以2016年飼養(yǎng)在廣西鳳翔集團(tuán)股份有限公司育種場(chǎng)的桂香雞母雞為實(shí)驗(yàn)動(dòng)物。根據(jù)公司飼養(yǎng)管理技術(shù)規(guī)程,按照飼養(yǎng)日齡將飼養(yǎng)過程分為1~5周齡(雛雞料)、6~10周齡(育成 1號(hào)料)、11~16周齡(育成2號(hào)料)3個(gè)不同階段,飼喂不同營(yíng)養(yǎng)水平的飼糧。桂香雞母雞3個(gè)階段飼糧營(yíng)養(yǎng)指標(biāo)見表1。飼養(yǎng)過程自由采食、飲水,執(zhí)行常規(guī)免疫程序,實(shí)施規(guī)范的飼養(yǎng)管理。在飼養(yǎng)過程中,每周隨機(jī)選取120只母雞,稱量并記錄體重(空腹稱重),數(shù)據(jù)統(tǒng)計(jì)至16周齡。
表1 桂香雞母雞生長(zhǎng)各階段飼糧營(yíng)養(yǎng)指標(biāo)
該 研 究 采 用 Logistic、Gompertz、Bertalanffy 3種動(dòng)物生長(zhǎng)曲線模型,運(yùn)用SPSS 19.0軟件,對(duì)桂香雞母雞的生長(zhǎng)曲線進(jìn)行擬合分析。3種畜禽常用生長(zhǎng)曲線模型公式等參數(shù)見表2。
表2 3種畜禽常用生長(zhǎng)曲線模型
表3 桂香雞母雞生長(zhǎng)曲線模型參數(shù)
采用Excel 2013計(jì)算桂香雞母雞每周的平均體重,利用SPSS 19.0軟件計(jì)算每種模型的擬合度(R2)和參數(shù)估計(jì)值(A,B,k)。 根據(jù)參數(shù)估計(jì)值 A、B、k,以及表2中的公式,計(jì)算各模型的拐點(diǎn)體重、拐點(diǎn)周齡和最大周增重。
桂香雞母雞累積生長(zhǎng)曲線(圖1)顯示,桂香雞母雞在1~3周齡時(shí)生長(zhǎng)較慢,在4周齡時(shí)生長(zhǎng)速度加快,在5~6周齡時(shí)生長(zhǎng)速度最快,在6周齡時(shí),周增重達(dá)到 99.15 g。
圖1 桂香雞母雞累積生長(zhǎng)曲線
表3為桂香雞母雞3種生長(zhǎng)曲線擬合的模型參數(shù)(A、B、k)、擬合度(R2)、拐點(diǎn)體重、拐點(diǎn)周齡、最大周增重。由表3可知,3種模型對(duì)桂香雞母雞體重?cái)M合度:Bertalanffy>Gompertz>Logistic, 且均達(dá)到0.981以上,其中,Bertalanffy模型的擬合度R2最高(0.993),說明在桂香雞母雞中該模型的擬合效果最好。Bertalanffy模型擬合的桂香雞母雞的拐點(diǎn)體重為333.10 g,拐點(diǎn)周齡4.55周,最大周增重為 78.94 g。
表4為桂香雞母雞體重實(shí)際觀察值與3種生長(zhǎng)模型預(yù)測(cè)值的數(shù)據(jù)對(duì)比。圖2為桂香雞母雞實(shí)際體重增長(zhǎng)曲線與3種生長(zhǎng)模型擬合曲線的對(duì)比。從整體上看,Logistic預(yù)測(cè)效果最差,在1、6和7周齡,預(yù)測(cè)值和實(shí)際值相差都在50 g以上;Bertalanffy預(yù)測(cè)值與實(shí)際值更加接近,更適合用于預(yù)測(cè)桂香雞母雞的體重。但從表3和圖2可以看出,3種模型預(yù)測(cè)的 6~7周齡體重以及 10~11周齡體重與實(shí)際觀察值都存在一定的差距,Bertalanffy預(yù)測(cè)值與實(shí)際值相差 30~41 g,Logistic的相差 44~60 g,Gompertz的相差 35~46 g。
生長(zhǎng)曲線直接將動(dòng)物體的生長(zhǎng)過程用線性的數(shù)據(jù)直觀地呈現(xiàn)在人們面前,利用數(shù)學(xué)模型對(duì)所獲得的生長(zhǎng)曲線進(jìn)行擬合分析,進(jìn)而從中找出畜禽的生長(zhǎng)規(guī)律,為日后畜禽生產(chǎn)發(fā)展提供寶貴的參考資料。目前,生長(zhǎng)曲線的擬合在家禽、家畜等方面已有較多報(bào)道[3-6]。 張燕等[7]、李青青等[8]、劉雙艷等[9]、伍維高等[10]運(yùn)用了 3 種模型對(duì)不同雞品種的生長(zhǎng)進(jìn)行了曲線擬合,結(jié)果一致發(fā)現(xiàn)Gompertz 模型的擬合效果最佳;而李?。?1]、牛建芹等[2]對(duì)不同雞品種生長(zhǎng)進(jìn)行擬合發(fā)現(xiàn)Bertalanffy模型的擬合效果最佳。該研究發(fā)現(xiàn),3種模型在桂香雞母雞中的擬合度均達(dá)到了0.981以上,但Bertalanffy模型的擬合效果最佳,預(yù)測(cè)值與實(shí)際值更加接近。綜合以上研究結(jié)果可知,生長(zhǎng)曲線擬合模型本身具有不同的特點(diǎn),不同的雞品種有著不同的生長(zhǎng)規(guī)律,所以選擇的最優(yōu)曲線擬合模型會(huì)有所不同,針對(duì)不同的品種選擇合適的模型,有助于了解其生長(zhǎng)發(fā)育規(guī)律。
表4 桂香雞母雞3種模型體重預(yù)測(cè)值
圖2 桂香雞母雞實(shí)際體重增長(zhǎng)曲線與3種擬合曲線的比較
該研究中,雖然3種模型在桂香雞中的擬合度均達(dá)到了0.981以上,但通過與實(shí)際值進(jìn)行比較,發(fā)現(xiàn)6~7周齡體重以及10~11周齡體重與實(shí)際觀察值都存在一定的差距,Bertalanffy預(yù)測(cè)值與實(shí)際值相差最小,在 30~41 g,Logistic和Gompertz預(yù)測(cè)值和實(shí)際值相差均在35~60 g。這可能是在桂香雞3階段法飼養(yǎng)的過程中,6和11周齡時(shí)日糧的替換造成桂香雞母雞的生長(zhǎng)受到一定程度的影響,因此,有必要對(duì)桂香雞母雞不同階段的營(yíng)養(yǎng)配方進(jìn)行研究,分析目前配方的合理性。
在生長(zhǎng)發(fā)育過程中,生長(zhǎng)速率最快的時(shí)間被稱為生長(zhǎng)拐點(diǎn)。李莉等[12]對(duì)桂香雞母雞1~17周齡的生長(zhǎng)曲線擬合表明,Logistic、Gompertz和 Bertalanffy模型擬合的拐點(diǎn)周齡分別為10.85、10.34和10.75周。在該研究中,3種模型擬合的生長(zhǎng)拐點(diǎn)分別為 6.78、5.36和 4.55周。 在李莉等[12]的研究中,桂香雞母雞16周齡的實(shí)際體重是1 501.40 g,而在該研究中桂香雞母雞16周齡的實(shí)際體重是977.39 g,這主要是由于在桂香雞母雞實(shí)際飼養(yǎng)過程中,為防止種雞長(zhǎng)速過快,改變了原有的飼養(yǎng)方式,對(duì)母雞限制飼喂造成的。
3種模型對(duì)桂香雞母雞生長(zhǎng)規(guī)律的擬合度均達(dá)到了0.981以上,但Bertalanffy模型擬合效果最好。在Bertalanffy模型下,預(yù)測(cè)的母雞拐點(diǎn)周齡為4.55周,拐點(diǎn)體重為333.10 g。該研究結(jié)果揭示了桂香雞的生長(zhǎng)發(fā)育規(guī)律,同時(shí)為桂香雞的生產(chǎn)管理和育種提供了參考依據(jù)。