柏云心 ，譚 紅 ，李玉蓮 ，吳買(mǎi)生 ，楚隨安 ，陳 斌

（1.湖南農(nóng)業(yè)大學(xué)動(dòng)物科技學(xué)院，湖南 長(zhǎng)沙 410128；2.湘潭家畜育種站，湖南 湘潭 411100；3.湘潭縣農(nóng)業(yè)農(nóng)村局，湖南 湘潭 411228）

隨著人們生活水平的提高，人們喜食瘦肉的同時(shí)也注重肉的品質(zhì)與風(fēng)味，高品質(zhì)的豬肉需求急劇增大，按照如今市場(chǎng)需求導(dǎo)向進(jìn)行豬種培育和豬肉供應(yīng)已經(jīng)成為當(dāng)今養(yǎng)豬業(yè)發(fā)展的重要部分。

為了培育出可滿足現(xiàn)在市場(chǎng)對(duì)高品質(zhì)豬肉需求的優(yōu)秀豬種，打造具有持久市場(chǎng)競(jìng)爭(zhēng)力的優(yōu)質(zhì)豬品牌，湘沙豬配套系選育課題組經(jīng)歷10余年，經(jīng)過(guò)5個(gè)世代的持續(xù)選育，成功培育了湘沙豬配套系。該課題組研究發(fā)現(xiàn)湘沙豬配套系的父母代豬肉質(zhì)好，于是計(jì)劃利用父母代豬進(jìn)行橫交固定，繼代選育出新品種—湘嶺豬，同時(shí)繼承了沙子嶺豬的肉質(zhì)風(fēng)味和巴克夏豬的生長(zhǎng)優(yōu)勢(shì)，為規(guī)?；B(yǎng)殖提供了可能。

湘嶺豬作為一個(gè)尚在培育的新品種，目前已經(jīng)培育到Ⅱ世代，許多資料和相關(guān)數(shù)據(jù)都十分欠缺。在生產(chǎn)中，體重是反映豬的生長(zhǎng)、健康、生產(chǎn)性能情況以及品種性能的重要指標(biāo)，科學(xué)測(cè)定與記錄體重變化，可以實(shí)現(xiàn)對(duì)生長(zhǎng)曲線擬合的研究，繼而了解豬的生長(zhǎng)規(guī)律。目前有許多關(guān)于動(dòng)物生長(zhǎng)曲線擬合研究的文獻(xiàn)，基于這些研究基礎(chǔ)，文章采用了Logistic模型、Gompertz模型和Saturation模型來(lái)對(duì)湘嶺豬體重與背膘生長(zhǎng)規(guī)律進(jìn)行擬合。

如今對(duì)于瘦肉型豬種性狀的校正公式已有較為可靠的制定方法與參考數(shù)據(jù)，而地方豬種的校正公式制定方法主要參照瘦肉型豬種，相關(guān)研究十分匱乏。本文對(duì)湘嶺豬的校正公式制定具有重要意義，不僅填補(bǔ)了湘嶺豬校正公式的空白，同時(shí)也為地方豬種和雜交豬種的校正公式制定提供了參考。

1 材料與方法

1.1 試驗(yàn)材料

本試驗(yàn)以湘潭市沙子嶺豬資源場(chǎng)2021年7月出生的Ⅱ世代湘嶺豬為試驗(yàn)對(duì)象，具有完整的系譜，飼養(yǎng)管理?xiàng)l件相同，出生日期相近。公豬14頭、母豬11頭，測(cè)定每頭豬的初生、21日齡、35日齡、60日齡、90日齡、120日齡、150日齡、180日齡時(shí)的體重以及120日齡、150日齡、180日齡的活體背膘厚。

1.2 試驗(yàn)工具

電子籠秤、獸用B超儀（百勝公司）、統(tǒng)計(jì)分析軟件SAS 9.4。

1.3研究方法

1.3.1 生長(zhǎng)曲線研究方法

綜合以前的文獻(xiàn)資料可知湘嶺豬的日齡與相應(yīng)的體重之間是非線性關(guān)系，不能用線性模型去統(tǒng)計(jì)分析，本研究選擇了三種常用的模型即Logistic模型、Gompertz模型和Saturation模型對(duì)湘嶺豬體重的生長(zhǎng)規(guī)律進(jìn)行了擬合，如表1。

表1 待定的生長(zhǎng)模型

1.3.2 體重與活體背膘厚模型擬合

本文采用了線性模型分別對(duì)湘嶺豬公豬、母豬的體重與活體背膘厚進(jìn)行了擬合。以下是線性模型的表達(dá)式:

其中，BF為活體背膘厚；a為待定參數(shù)；b為斜率。

1.3.3 模型的方差分析與擬合度比較

本文采用了方差分析法對(duì)模型的擬合效果進(jìn)行了顯著性檢驗(yàn)，其原理為：

令SSR為回歸平方和，SSe為殘差平方和，則，

若要回歸方程有意義，就需要證明因變量的變化主要是由自變量的變化引起的，所以SSR越大并且SSe越小，則F值越大即差異越顯著，表示自變量對(duì)因變量的影響越大。

Logistic模型和Gompertz模型的擬合優(yōu)越性也常用可決系數(shù)R2（Coefficient of Determination）來(lái)表示，其值越接近1，代表擬合度越高。

1.3.4 校正公式的制定

目前對(duì)性狀校正公式制定普遍采用的是國(guó)外的計(jì)算方法，也主要是對(duì)瘦肉型豬進(jìn)行了校正公式的制定，其中校正因子的數(shù)值有較為準(zhǔn)確的參考，但是對(duì)于地方豬以及其他新培育的豬種的校正公式制定的研究較少，缺少合適的計(jì)算公式以及參考數(shù)值。綜合對(duì)瘦肉型豬的校正公式的制定方法，以及公式的實(shí)際應(yīng)用情況分析，本試驗(yàn)采用Saturation模型擬合結(jié)果來(lái)制定湘嶺豬的達(dá)50 kg體重日齡的校正公式，利用線性模型對(duì)體重與活體背膘厚的擬合結(jié)果來(lái)制定50 kg體重背膘厚的校正公式。

根據(jù)試驗(yàn)對(duì)生長(zhǎng)模型的擬合結(jié)果，得到相關(guān)參數(shù)，然后制定達(dá)50 kg體重日齡與50 kg體重背膘厚的校正公式。

其中，C 為校正系數(shù)，C=性狀的目標(biāo)預(yù)測(cè)值÷性狀的實(shí)測(cè)值。

2 結(jié)果與分析

2.1 湘嶺豬體重生長(zhǎng)曲線模型擬合

根據(jù)動(dòng)物生長(zhǎng)模型，利用統(tǒng)計(jì)軟件SAS 9.4分別對(duì)湘嶺豬公豬、母豬的體重與日齡進(jìn)行擬合，各結(jié)果如表2、表3，結(jié)果均以“估計(jì)值±標(biāo)準(zhǔn)誤”表示。

表2 湘嶺豬體重生長(zhǎng)曲線模型擬合結(jié)果

表3 生長(zhǎng)曲線模型參數(shù)的95%置信區(qū)間

在公豬體重與日齡的模型擬合中，Logistic模型、Gompertz模型和Saturation模型可決系數(shù)R2的值分別為0.994 5、0.994 4、0.917 1，可以看出Logistic模型、Gompertz模型擬合度非常接近，都比Saturation模型要高。母豬的生長(zhǎng)模型擬合結(jié)果中，Logistic模型、Gompertz模型和Saturation模型可決系數(shù)R2的值分別為0.989 8、0.990 1、0.918 3，都大于0.9，對(duì)湘嶺豬的體重生長(zhǎng)規(guī)律擬合程度極高，Logistic模型、Gompertz模型擬合度非常接近，均優(yōu)于Saturation模型。

湘嶺豬公豬的Logistic模型估計(jì)的拐點(diǎn)體重和拐點(diǎn)日齡分別為49.3 kg和124.1 d，Gompertz模型估計(jì)的拐點(diǎn)體重和拐點(diǎn)日齡分別為51.7 kg和127.9 d；湘嶺豬母豬的Logistic模型估計(jì)的拐點(diǎn)體重和拐點(diǎn)日齡分別為46.7 kg、121.7 d，Gompertz模型估計(jì)的拐點(diǎn)體重和拐點(diǎn)日齡分別47.8 kg 和122.9 d。兩種模型估計(jì)的拐點(diǎn)體重與拐點(diǎn)日齡均非常接近，公豬相對(duì)于母豬來(lái)說(shuō)，拐點(diǎn)體重與拐點(diǎn)日齡數(shù)值較大。

Logistic模型、Gompertz模型和Saturation模型的方差分析結(jié)果如表4所示。

從表4中可以看出，公豬與母豬的Logistic模型、Gompertz模型和Saturation模型的方差分析結(jié)果Pr值都小于0.01，說(shuō)明三種模型都與試驗(yàn)數(shù)據(jù)極顯著擬合。

表4 方差分析結(jié)果

2.2 體重與背膘厚回歸方程擬合

利用線性回歸模型對(duì)湘嶺豬體重與背膘厚兩者生長(zhǎng)關(guān)系進(jìn)行了擬合，并同時(shí)對(duì)模型進(jìn)行了方差分析，結(jié)果如表5、表6所示。

表5 生長(zhǎng)模型參數(shù)估計(jì)值

表6 方差分析結(jié)果

根據(jù)上述結(jié)果可知，公豬與母豬的線性模型參數(shù)b的Pr值均小于0.000 1，表示體重與背膘厚的線性關(guān)系極顯著成立。

公豬的體重與活體背膘厚線性模型參數(shù)估計(jì)值為： a = -11.890 3；b = 0.585 0，由此可知，公豬的體重與活體背膘厚線性回歸方程為：

母豬的體重與活體背膘厚線性模型參數(shù)估計(jì)值為： a = -10.993 2；b = 0.594 9，由此可知，母豬的體重與活體背膘厚線性回歸方程為：

對(duì)公豬、母豬體重與背膘厚的線性模型擬合的結(jié)果進(jìn)行方差分析，F(xiàn)值分別為61.18和32.72，公豬，母豬的Pr 都小于0.000 1，表示公豬、母豬體重與背膘厚的線性關(guān)系極顯著成立。

2.3 校正公式的制定

2.3.1 達(dá)50 kg體重日齡的校正公式制定結(jié)果

因?yàn)?，校正系?shù)C=性狀的預(yù)測(cè)值÷性狀的實(shí)測(cè)值。

那么t50=t實(shí)×C，校正系數(shù)C=t50/t實(shí)，其中t50為達(dá)50 kg體重日齡模型估計(jì)值，t實(shí)為實(shí)測(cè)日齡。因?yàn)镾aturation 模型表達(dá)式為Wt=t/（a+b·t），可知，，其中W實(shí)為實(shí)測(cè)體重，校正因子=1-50b，帶入公豬Saturation 模型參數(shù)值，整理得：

利用同樣原理可推出母豬Saturation 模型下的校正公式。

由公式可看出當(dāng)實(shí)測(cè)體重大于50 kg時(shí)，達(dá)50 kg校正日齡只為會(huì)比實(shí)測(cè)日齡更大；反之，校正日齡比實(shí)測(cè)日齡更小，這與事實(shí)相悖。根據(jù)實(shí)際情況，本文將校正因子（1-50b）取其絕對(duì)值，最終得到湘嶺豬公、母豬的校正公式為：

在代入50 kg體重附近的實(shí)際值校正后，發(fā)現(xiàn)與實(shí)際情況接近，本試驗(yàn)制定的校正公式對(duì)于湘嶺豬在一定體重范圍內(nèi)是具有可信度的。

2.3.2 50 kg體重活體背膘厚校正公式制定結(jié)果

根據(jù)2.2節(jié)可知，湘嶺豬體重與活體背膘厚線性回歸模型表達(dá)式為：

母豬豬50 kg體重活體背膘厚校正公式：

3 討論

3.1 關(guān)于體重生長(zhǎng)曲線擬合結(jié)果

通過(guò)前面的統(tǒng)計(jì)結(jié)果可知，Logistic模型和Gompertz模型擬合度優(yōu)于Saturation模型，三者對(duì)湘嶺豬均有極高的擬合度，與杜洛克、梅山豬公豬及從江香豬等品種的研究有著類似結(jié)果。徐如海等對(duì)14頭杜洛克公豬連續(xù)稱重后，對(duì)其生長(zhǎng)曲線進(jìn)行了Logistic模型和Gompertz模型擬合，最終選定了Gompertz模型。許棟等通過(guò)Logistic模型和Gompertz模型擬合梅山豬公豬的生長(zhǎng)發(fā)育數(shù)據(jù)，結(jié)果以Logistic模型的擬合優(yōu)度最高，擬合結(jié)果最接近真實(shí)情況。徐永健等采用Logistic模型和Gompertz模型擬合了從江香豬公豬的生長(zhǎng)曲線，得出Logistic模型為最適模型。通過(guò)查閱對(duì)豬的生長(zhǎng)曲線擬合研究結(jié)果來(lái)看，對(duì)公豬的生長(zhǎng)曲線模型更多采用了Logistic模型，與本文的研究結(jié)果一致。為了更好地掌握湘嶺豬生長(zhǎng)規(guī)律，本試驗(yàn)還研究了母豬的生長(zhǎng)曲線模型，通過(guò)對(duì)比三種模型的決定系數(shù)R2發(fā)現(xiàn)，母豬和公豬擬合結(jié)果類似，R2數(shù)值都很接近1，方差分析三種模型的擬合度都極為顯著，最后通過(guò)與實(shí)際生長(zhǎng)發(fā)育數(shù)據(jù)比較分析，母豬最適生長(zhǎng)模型采用了Gompertz模型。

由2.1節(jié)可以看出，Logistic模型統(tǒng)計(jì)出的湘嶺豬生長(zhǎng)速度略快于Gompertz模型，母豬要比公豬更早進(jìn)入體重快速增長(zhǎng)階段，即母豬的發(fā)育略早于公豬。利用拐點(diǎn)日齡，可以指導(dǎo)湘嶺豬的飼養(yǎng)管理，達(dá)拐點(diǎn)日齡后，可通過(guò)提高豬只的營(yíng)養(yǎng)水平和飼養(yǎng)供給量滿足豬只在拐點(diǎn)日齡后體重快速增長(zhǎng)的營(yíng)養(yǎng)需要。

3.2 關(guān)于背膘厚與體重生長(zhǎng)曲線擬合結(jié)果

經(jīng)查閱文獻(xiàn)，目前并沒(méi)有被普遍接受的對(duì)體重與活體背膘厚兩者生長(zhǎng)規(guī)律擬合程度較高的生長(zhǎng)模型，國(guó)外對(duì)于豬的體重與背膘厚生長(zhǎng)規(guī)律也研究得較少。褚青坡等首次參考 Schinckel等的模型對(duì)蘇淮豬進(jìn)行了體重與背膘的生長(zhǎng)曲線擬合，擬合程度并不高。劉望宏等對(duì)長(zhǎng)白豬、杜洛克豬、大約克豬的體重與背膘厚進(jìn)行了Logistic模型、Gompertz模型和直線模型擬合，結(jié)果直線回歸模型擬合度最高，并且Logistic模型、Gompertz模型擬合曲線圖接近直線，同時(shí)直線模型具有制定校正公式的優(yōu)勢(shì)。本試驗(yàn)預(yù)計(jì)對(duì)湘嶺豬的體重與背膘厚的生長(zhǎng)曲線進(jìn)行Logistic模型、Gompertz模型擬合，但利用SAS軟件初步統(tǒng)計(jì)分析后，發(fā)現(xiàn)無(wú)法迭代出結(jié)果，不滿足收斂準(zhǔn)則，因此只采用了直線回歸模型。在母豬的直線回歸模型擬合結(jié)果中，參數(shù)a的估計(jì)值t檢驗(yàn)Pr=0.119 8>0.05，不顯著，代表a的估計(jì)值誤差較大，在實(shí)際生產(chǎn)中，參數(shù)a的意義并不大。但回歸系數(shù)b的估計(jì)值經(jīng)t檢驗(yàn)后，公豬、母豬都為極顯著，代表了該系數(shù)具有很大可信度。

由2.2節(jié)中公豬、母豬的體重與背膘線性回歸方程可計(jì)算出，初生重為1 kg左右的小豬活體背膘厚為負(fù)值，這明顯與實(shí)際相悖。造成這一現(xiàn)象的原因，筆者認(rèn)為由于本試驗(yàn)測(cè)定活體背膘厚的體重范圍在45 kg到95 kg這一階段，背膘生長(zhǎng)速度遠(yuǎn)大于初生重到45 kg這一階段，所以導(dǎo)致模型與豬只前面階段的生長(zhǎng)規(guī)律不擬合這一現(xiàn)象，也是造成縱截距a的估計(jì)值不顯著原因所在。由實(shí)際情況可知，體重與背膘厚兩者的生長(zhǎng)關(guān)系并不是直線關(guān)系，只是前期背膘生長(zhǎng)相較于中后期過(guò)于緩慢，而中后期二者生長(zhǎng)關(guān)系更像直線關(guān)系，而本文只研究了湘嶺豬120日齡到180日齡的生長(zhǎng)規(guī)律，也就出現(xiàn)縱截距估計(jì)值誤差過(guò)大而回歸系數(shù)b在本試驗(yàn)數(shù)據(jù)下可信度極顯著的現(xiàn)象。關(guān)于體重與背膘厚生長(zhǎng)關(guān)系還有待更多的研究，以期尋找到與實(shí)際生長(zhǎng)規(guī)律擬合度更高的模型。

3.3 關(guān)于校正公式制定

經(jīng)查閱相關(guān)文獻(xiàn)發(fā)現(xiàn)，Saturation模型在制定校正公式時(shí)更具優(yōu)勢(shì)，利用其制定的校正公式更符合實(shí)際，限制條件也更少。根據(jù)Saturation 模型對(duì)湘嶺豬體重與日齡的擬合參數(shù)估計(jì)值來(lái)制定的達(dá)50 kg體重日齡校正公式，在實(shí)際的數(shù)據(jù)帶入計(jì)算后，發(fā)現(xiàn)在一定范圍內(nèi)是具有實(shí)際參考意義的。50 kg體重背膘厚校正公式則根據(jù)線性模型擬合結(jié)果制定而成，兩個(gè)校正公式變換形式后與美國(guó)早期性狀校正公式結(jié)構(gòu)相似。本試驗(yàn)校正公式制定方法與結(jié)果和劉望宏等對(duì)長(zhǎng)白豬、杜洛克豬、大約克豬制定的性狀校正公式結(jié)果相似，公式在現(xiàn)有的條件下有一定的可信度。湘嶺豬性狀校正公式，能將在一定范圍內(nèi)波動(dòng)的數(shù)據(jù)校正至同一標(biāo)準(zhǔn)下，有利于客觀的評(píng)價(jià)以及計(jì)算相應(yīng)性狀的遺傳參數(shù)，對(duì)于育種工作有著重要實(shí)際意義。

4 結(jié)論

通過(guò)對(duì)湘嶺豬公豬、母豬的體重與日齡生長(zhǎng)曲線擬合研究發(fā)現(xiàn)，公豬最適生長(zhǎng)模型采用Logistic模型，母豬最適生長(zhǎng)模型采用Gompertz模型。在制定達(dá)50 kg體重日齡的校正公式時(shí)，運(yùn)用Saturation模型制定的校正公式更具有實(shí)際意義。