吳有林 陽文攀,3 何清輝 張文勁 李文靜 劉鳴慧 羅 誠 丁能水,3＊

（1.福建傲農(nóng)生物科技集團(tuán)股份有限公司 福建 漳州 363000；2.上杭傲農(nóng)槐豬產(chǎn)業(yè)發(fā)展有限公司 福建 龍巖 364200；3.江西農(nóng)業(yè)大學(xué)豬遺傳改良與養(yǎng)殖技術(shù)國家重點(diǎn)實(shí)驗(yàn)室 南昌 330045；4.漳浦縣農(nóng)業(yè)農(nóng)村局 福建 漳州 363200）

上杭槐豬是福建地方優(yōu)良品種豬，其性情溫順，母性好，能適應(yīng)粗放的飼養(yǎng)管理[1]，屬國家級(jí)畜禽品種保護(hù)資源。與瘦肉型豬（外來種豬）相比，槐豬肌纖維細(xì)、羥脯氨酸含量高，肌內(nèi)脂肪含量高、膽固醇含量低[2]。但上杭槐豬生長速度較慢，個(gè)體間變異較大，有充足的選育空間[3-4]。由于上杭槐豬初生重較?。s為0.6 kg），個(gè)體間特征差異不明顯，初生選留效果較差；而在斷乳時(shí)對(duì)上杭槐豬進(jìn)行精選，可增加優(yōu)秀上杭槐豬使用比例。研究表明，斷乳日齡對(duì)仔豬生長性能、養(yǎng)分消化率及消化酶活性有顯著的影響，斷乳越早仔豬生長及消化方面的機(jī)能越不完善[5]。斷乳日齡可顯著影響仔豬抗氧化功能，且斷乳日齡越早，機(jī)體抗氧化功能降低的程度越大，恢復(fù)時(shí)間越長[6]。斷乳仔豬個(gè)體重越大，生產(chǎn)性能越高，且仔豬斷乳重與腹瀉率和死亡率呈顯著負(fù)相關(guān)關(guān)系[7]。故選擇合適的斷乳日齡十分必要，在現(xiàn)有的本地豬研究中，發(fā)現(xiàn)臨高仔豬在35日齡或42日齡時(shí)斷乳較好[8]，五指山仔豬35日齡或42日齡斷乳可改善五指山仔豬生長性能[9]。與其他土豬斷乳日齡類似，上杭槐豬合適斷乳日齡為42日齡。然而對(duì)于上杭槐豬斷乳前生長性能研究較少，因此，了解上杭槐豬斷乳前生長趨勢，有助于加強(qiáng)上杭槐豬斷乳前飼養(yǎng)管理。通過對(duì)上杭槐豬斷乳時(shí)生長性能的衡量，可制定適宜的斷乳選留策略以替代出生選留，篩選性能優(yōu)越的仔豬。本研究對(duì)上杭槐豬斷乳前不同日齡體重進(jìn)行統(tǒng)計(jì)，以確定上杭槐豬斷乳前生長曲線。同時(shí)創(chuàng)建達(dá)8 kg體重日齡校正公式，制定斷乳選留策略，以此挑選斷乳時(shí)生長性能優(yōu)越的仔豬。

1 材料與方法

1.1 試驗(yàn)動(dòng)物、稱重與記錄 選擇上杭傲農(nóng)槐豬產(chǎn)業(yè)發(fā)展有限公司2021年3月31日至2021年5月1日出生的103頭仔豬，依據(jù)《上杭傲農(nóng)槐豬分娩舍飼養(yǎng)管理標(biāo)準(zhǔn)》進(jìn)行飼喂管理。仔豬出生第2 d進(jìn)行出生登記，15 d后進(jìn)行體重稱量，記錄仔豬耳號(hào)、性別、出生日期、記錄日期和體重。初次稱重至斷乳期間進(jìn)行多次體重稱量、記錄，共獲得624條數(shù)據(jù)。

1.2 統(tǒng)計(jì)分析

1.2.1 數(shù)據(jù)處理 從稱重記錄中篩選高質(zhì)量的稱重記錄，個(gè)體耳號(hào)、性別、出生日期、記錄日期和體重記錄完整的為一條合格記錄，使用合格數(shù)據(jù)進(jìn)行早期生長性能描述。

用于達(dá)8 kg體重日齡校正公式研究的仔豬合格稱重記錄至少為5次，體重在8 kg附近的稱重記錄間隔最多為7 d（4月1日至4月2日間隔1 d）。

1.2.2 達(dá)8 kg體重日齡校正系數(shù)計(jì)算 根據(jù)加拿大達(dá)30 kg體重日齡校正公式，建立槐豬達(dá)8 kg體重日齡校正公式，校正日齡=實(shí)際日齡-（實(shí)際體重-8）×CF。其中，CF為達(dá)8 kg體重日齡校正系數(shù)。

首先使用R（v4.1.0）通過線性插值方式從距8 kg最近的2條數(shù)據(jù)中獲得8 kg體重日齡真值，通過R語言lm函數(shù)使用每個(gè)個(gè)體的多條實(shí)測值對(duì)校正公式參數(shù)進(jìn)行求解，從而獲取槐豬達(dá)8 kg體重日齡校正系數(shù)。

1.2.3 槐豬達(dá)8 kg體重日齡校正公式擬合度分析通過槐豬達(dá)8 kg體重日齡校正公式，計(jì)算不同日齡、體重下槐豬達(dá)8 kg體重校正日齡，通過8 kg體重日齡真值減去校正日齡從而獲取殘差。使用R語言ggplot2包作箱線圖，比較不同體重箱線圖中位線與殘差=0的標(biāo)準(zhǔn)線。

2 結(jié)果與分析

2.1 早期生長性能統(tǒng)計(jì) 使用103頭上杭槐豬624次合格稱量記錄進(jìn)行不同日齡體重箱線圖描繪，使用局部多項(xiàng)式回歸進(jìn)行生長曲線擬合。由圖1可知，槐豬在斷乳前體重增長擬合線近似直線，槐豬在42 d時(shí)體重約為8 kg。

圖1 槐豬不同日齡體重箱線圖

2.2 數(shù)據(jù)質(zhì)控 去除稱重記錄只包含大于8 kg或小于6 kg的槐豬后，剩余76頭槐豬仔豬503條合格記錄。為保證進(jìn)行達(dá)8 kg體重日齡校正公式的豬只稱量質(zhì)量，去除測量記錄少于5條或在8 kg附近的稱重記錄間隔大于7 d的豬只，剩余49頭槐豬仔豬383條合格記錄。質(zhì)控后槐豬日齡體重箱線圖見圖2，可知槐豬不同日齡體重箱線圖中線近似處于一條直線上。通過線性插值方式從距8 kg最近的2條數(shù)據(jù)中獲得8 kg體重日齡真值，準(zhǔn)確性較高?；必i達(dá)8 kg體重日齡分布見圖3，無極值豬只存在。

圖2 質(zhì)控后槐豬日齡體重箱線圖

2.3 達(dá)8 kg體重日齡校正系數(shù) 使用R語言l m函數(shù)基于線性擬合的方式，用每個(gè)個(gè)體的多條實(shí)測值對(duì)校正公式參數(shù)進(jìn)行求解。獲取槐豬達(dá)8 kg體重日齡校正系數(shù)為4.643，P值為2×10-16，判定系數(shù)為0.96。由此可知，通過線性擬合獲取的達(dá)8 kg體重日齡校正系數(shù)校正效果較好。

2.4 達(dá)8 kg體重日齡校正公式擬合度分析 以不同日齡、體重下槐豬達(dá)8 kg體重校正日齡與通過線性插值方式獲取的8 kg體重日齡真值殘差作箱線圖。由圖3可知，體重在6～8 kg時(shí)，殘差中值與0最接近，而當(dāng)體重小于6 kg或大于8 kg時(shí)殘差中值與0偏差較大。

圖3 槐豬達(dá)8 kg體重日齡分布

圖3 達(dá)8 kg體重日齡殘差箱線圖

3 討 論

在規(guī)?；?、集約化養(yǎng)豬生產(chǎn)中，仔豬早期斷乳能夠增加產(chǎn)房欄舍使用率，提高養(yǎng)殖場經(jīng)濟(jì)效益[10]，外來種豬實(shí)行早期斷乳，如21日齡、28日齡斷乳[11-12]。斷乳日齡對(duì)仔豬體重、日增重、相對(duì)生長速率和日采食量均有顯著影響，且斷乳越早，這些指標(biāo)于斷乳后下降的幅度越大[13-14]。上杭槐豬屬于地方優(yōu)良品種豬，由于其出生重較小，早期斷乳死亡率較高，斷乳日齡較晚為42日齡。上杭槐豬斷乳前生長曲線近似直線，在42日齡時(shí)體重約為8 kg。

為了有效利用上杭槐豬種質(zhì)資源，斷乳時(shí)對(duì)上杭槐豬生長性能進(jìn)行初步評(píng)估，可擴(kuò)大優(yōu)秀上杭槐豬的利用率?，F(xiàn)有對(duì)外來種豬體重校正公式的研究較多[15-16]，對(duì)土豬品種體重校正公式的研究較少[17]，且不能直接進(jìn)行套用。為了對(duì)斷乳時(shí)上杭槐豬生長性能進(jìn)行初步評(píng)估，設(shè)立上杭槐豬達(dá)8 kg體重日齡校正公式：校正日齡=實(shí)際日齡-（實(shí)際體重-8）×CF，其中CF=4.643。由達(dá)8 kg體重日齡校正公式擬合度分析可知，該公式在6～8 kg時(shí)有良好的校正效果。公式在9 kg時(shí)校正效果較差，可能是由于收集的數(shù)據(jù)量較少的原因。

通過對(duì)上杭槐豬斷乳前不同日齡體重統(tǒng)計(jì)，繪制上杭槐豬斷乳前生長曲線，確定42日齡斷乳重約為8 kg。通過上杭槐豬達(dá)8 kg體重日齡校正公式的建立，可為上杭槐豬斷乳選留策略制定提供參考。

4 結(jié) 論

通過對(duì)上杭槐豬斷乳前生長性能的描述，確定了上杭槐豬42日齡斷乳時(shí)目標(biāo)體重。進(jìn)而完成對(duì)上杭槐豬校正8 kg體重日齡公式的建立，有助于衡量斷乳時(shí)上杭槐豬生長性能，并在斷乳時(shí)對(duì)上杭槐豬進(jìn)行初次精選，增強(qiáng)槐豬群體性能及保種利用。