沙烏頭豬育肥階段生長曲線擬合的研究

陸雪林，吳昊旻，雷勝輝，李何君，張新生，沈 東

（1.上海市動(dòng)物疫病預(yù)防控制中心，上海 201103；2.上海市松江區(qū)動(dòng)物疫病預(yù)防控制中心，上海 201600；3.上海市崇明區(qū)種畜場，上海 202150）

沙烏頭豬又名沙河頭豬、沙胡頭豬，其頭部大、耳下垂，形似黃砂茶壺，鬃毛烏黑色，其中心產(chǎn)區(qū)啟東市在1928 年建縣之前稱為外沙（屬上海市崇明縣管轄），因而得名。沙烏頭豬主要分布于江蘇省海門市及上海市的崇明區(qū)[1]。沙烏頭豬具有早熟、繁殖力高、泌乳力強(qiáng)、使用年限長和肉質(zhì)鮮美、嫩而多汁等特點(diǎn)，屬于江海型豬種，是我國豬種中具有繁殖力強(qiáng)、產(chǎn)仔數(shù)多的著名地方品種?！吨袊笄葸z傳資源志·豬志》（2011）將沙烏頭豬從太湖豬類群里面劃出成為一個(gè)獨(dú)立品種[2]。2014 年2 月，沙烏頭豬被列入《國家級畜禽遺傳資源保護(hù)名錄》。

畜禽生長曲線的擬合和分析是畜禽育種和生產(chǎn)的基礎(chǔ)工作，是研究畜禽生長發(fā)育規(guī)律主要方法之一。理想的生長曲線模型在指導(dǎo)飼養(yǎng)管理和育種實(shí)踐及控制畜禽生長發(fā)育過程中具有重要作用。上海市的4 個(gè)地方畜禽品種被列入《國家級畜禽遺傳資源保護(hù)名錄》，其中梅山豬、浦東白豬和浦東雞已通過Logistic、Gompertz、Von Bertalanffy 模型擬合分析了生長曲線，摸清了其生長發(fā)育規(guī)律，取得了良好的效果，而沙烏頭豬至今還沒有這方面的研究和分析。

本研究運(yùn)用2 種典型的非線性生長曲線模型Logistic、Von Bertalanffy 對沙烏頭豬進(jìn)行生長曲線擬合，通過分析其生長速度和性別間的差異，揭示沙烏頭豬的生長發(fā)育規(guī)律，為該品種保護(hù)、選育和開發(fā)利用提供參考依據(jù)。

1 材料與方法

1.1 實(shí)驗(yàn)動(dòng)物 實(shí)驗(yàn)豬由上海市崇明區(qū)種畜場提供，2016—2018 年，每年從純繁沙烏頭豬選取3～4 窩，每窩隨機(jī)選擇體重適中的仔豬2～3 頭作為待測豬，佩戴電子識(shí)別耳標(biāo)。

1.2 飼養(yǎng)管理和數(shù)據(jù)采集 仔豬在場內(nèi)完成哺乳和保育階段后，達(dá)約20 kg 體重送至上海市種豬測定中心，經(jīng)過隔離觀察后轉(zhuǎn)入測定舍集中測定。采用奧斯本全自動(dòng)種豬生產(chǎn)性能測定系統(tǒng)每天測定豬個(gè)體重、采食量等數(shù)據(jù)。測定期間的飼料采用指定公司生產(chǎn)的商品飼料，歷年使用的飼料配方未變且營養(yǎng)標(biāo)準(zhǔn)一致，3 種飼料的主要營養(yǎng)成分見表1。采用自由采食和飲水，常規(guī)免疫接種，集中測定時(shí)間安排在每年的9—12 月，累計(jì)測定數(shù)量沙烏頭公豬（閹割）42 頭、母豬39 頭。

1.3 生長曲線模型 采用Logistic、Von Bertalanffy 2種生長曲線模型，模型中待定參數(shù)A 為極限生長量（終極生長量或成熟體重），k 為瞬時(shí)相對生長率，B 為常數(shù)尺度。2 種生長曲線模型表達(dá)式、拐點(diǎn)體重、拐點(diǎn)日齡、最大日增重及相對生長率見表2。統(tǒng)計(jì)分析采用SPSS 19.0 軟件，利用不同日齡體重?cái)?shù)據(jù)擬合算出模型參數(shù)的最優(yōu)估計(jì)值A(chǔ)、B、k，建立生長模型，算出拐點(diǎn)體重、拐點(diǎn)日齡、最大日增重和擬合度（R2）。

表1 商品豬飼料主要營養(yǎng)成分

表2 本研究所采用的2 種生長曲線模型

2 結(jié)果

2.1 沙烏頭公豬和母豬的實(shí)際體重測量值 由表3 可知，112～238 日齡母豬體重略大于公豬，252 日齡后公豬體重超過母豬。

表3 沙烏頭豬各日齡的體重 kg

2.2 沙烏頭豬生長模型的擬合分析 由表4 可知，2 種模型均能很好地?cái)M合沙烏頭公豬和母豬的生長曲線，其擬合度均達(dá)到0.99 以上。2 種模型擬合的生長曲線與實(shí)際生長曲線基本吻合。就2 種模型而言，Von Bertalanffany 模型的擬合效果更好，公、母豬的擬合度分別達(dá)到0.995 和0.999，高于Logistic 模型的擬合度。由圖1 和圖2 可以看出，用2 種模型擬合的沙烏頭豬公豬和母豬各日齡體重值與實(shí)際體重測量值（實(shí)測值）均非常接近。

2.3 沙烏頭豬相對生長率 2 種模型對沙烏頭豬擬合得到的相對生長率均是隨日齡增加而降低（表5），公、母豬基本一致，Logistic 模型112～126 日齡母豬比公豬略大，隨后母豬的相對生長率較公豬略??；Von Bertalanffany 模型112～168 日齡母豬比公豬略大，隨后母豬的相對生長率較公豬略小，Von Bertalanffany 模型的相對生長率比Logistic 模型大。

2.4 沙烏頭豬生長曲線擬合優(yōu)度的檢驗(yàn) 根據(jù)崔黨群等[3]報(bào)道，可采用適合性χ2檢驗(yàn)進(jìn)行曲線回歸方程擬合度的檢驗(yàn)。由表6 可知，2 種模型的χ2檢驗(yàn)值均明顯小于χ2（0.05）和χ2（0.01）（df=11），這說明2 種模型均能較好地?cái)M合沙烏頭豬公、母豬育肥階段生長曲線。綜合R2和χ2檢驗(yàn)得出，Von Bertalanffany 模型比Logistic 模型擬合效果更優(yōu)。

3 討 論

豬的生長發(fā)育主要受遺傳和環(huán)境影響，但生長發(fā)育的基本特征具有相對穩(wěn)定性，這種穩(wěn)定性可以采用生長曲線進(jìn)行擬合分析，為豬的生長發(fā)育、精準(zhǔn)飼養(yǎng)和選育提供基礎(chǔ)數(shù)據(jù)。目前，最常用的非線性生長曲線模型有Logistic、Von Bertalanffy 和Gompertz，不同品種的生長曲線模型存在差異。

本研究主要是針對沙烏頭豬的育肥階段（20～85 kg），采用奧斯本全自動(dòng)種豬生產(chǎn)性能測定系統(tǒng)測定，數(shù)據(jù)精確可靠，采用Logistic、Von Bertalanffy 2 種生長曲線模型擬合分析均取得很好的效果，通過綜合R2和χ2檢驗(yàn)得出，Von Bertalanffany 模型比Logistic 模型擬合效果更優(yōu)。Von Bertalanffy 模型不僅擬合度最高，公、母豬的擬合度分別達(dá)到0.995 和0.999，而且拐點(diǎn)日齡、拐點(diǎn)體重和體重最大值均符合沙烏頭豬育肥階段的生長發(fā)育規(guī)律，因此，用于預(yù)測沙烏頭豬育肥階段體重最佳。這與陳明飛等[4]研究的貴州小型豬生長曲線和張順華等[5]研究的涼山豬的生長發(fā)育規(guī)律及其生長曲線擬合分析的結(jié)果一致。但在對不同品種豬的研究中3 種生長曲線模型的擬合效果也有不一致報(bào)道。鐘土木等[6]對金華豬I 系生長曲線分析與擬合研究和李慶崗等[7]對圩豬體重非線性生長曲線擬合分析，Gompertz 模型在擬合度和預(yù)測體重效果方面相對較好。肖煒等[8]對長白豬生長肥育期生長規(guī)律的研究和高鵬飛等[9]對五指山豬近交系生長曲線進(jìn)行的分析和擬合研究中Logistic 模型擬合效果較好。這可能與豬的品種、飼料配方、飼養(yǎng)環(huán)境、飼養(yǎng)周期、研究方法等不同有關(guān)，也可能與測定的生長階段不同有關(guān)。本研究盡管在沙烏頭豬育肥階段的測定非常精確，取得了令人滿意的效果，但豬完整的生長發(fā)育測定時(shí)間應(yīng)該從初生開始更加理想，這樣才能更好地反映沙烏頭豬整個(gè)生長發(fā)育階段的規(guī)律，因此有待于在今后的研究中進(jìn)一步完善。

表4 2 種生長曲線擬合參數(shù)估計(jì)值和擬合度

表5 沙烏頭豬的相對生長率

表6 體重模型估計(jì)值與實(shí)測值適合性χ2 檢驗(yàn)

4 結(jié) 論

本研究可得出，在沙烏頭豬112～266 日齡，Logistic、Von Bertalanffany 2 種生長曲線擬合度均超過0.99，均能較好地對公、母豬生長曲線進(jìn)行擬合，Von Bertalanffy 模型擬合度高于Logistic 模型。在Von Bertalanffy 模型擬合生長曲線中，沙烏頭豬公豬最大理論體重為124.227 kg，生長拐點(diǎn)體重是36.808 kg，拐點(diǎn)日齡是146 日齡，最大日增重為454.579 g，沙烏頭豬母豬最大理論體重為111.979 kg，生長拐點(diǎn)體重是33.179 kg，拐點(diǎn)日齡是133 日齡，最大日增重為467.824 g。通過擬合沙烏頭豬的生長曲線，可為動(dòng)態(tài)分析沙烏頭豬的生長發(fā)育、指導(dǎo)飼養(yǎng)管理和選育提供可靠依據(jù)。