楊岸奇，吳買生，向擁軍，陳 斌，*，李朝輝，張善文，劉 偉，吳攀峰

(1.湖南農(nóng)業(yè)大學(xué)動(dòng)物科學(xué)技術(shù)學(xué)院， 長(zhǎng)沙 410128； 2.湘潭市畜牧水產(chǎn)局， 湖南 湘潭 411104；3.湘潭市家畜育種站， 湖南 湘潭 411104； 4.畜禽遺傳改良湖南省重點(diǎn)實(shí)驗(yàn)室， 長(zhǎng)沙 410128)

沙子嶺豬生長(zhǎng)性狀校正公式的制定

楊岸奇1，4，吳買生2，3，向擁軍3，陳 斌1，4，*，李朝輝3，張善文3，劉 偉3，吳攀峰3

(1.湖南農(nóng)業(yè)大學(xué)動(dòng)物科學(xué)技術(shù)學(xué)院， 長(zhǎng)沙 410128； 2.湘潭市畜牧水產(chǎn)局， 湖南 湘潭 411104；3.湘潭市家畜育種站， 湖南 湘潭 411104； 4.畜禽遺傳改良湖南省重點(diǎn)實(shí)驗(yàn)室， 長(zhǎng)沙 410128)

生長(zhǎng)性狀校正公式對(duì)于動(dòng)物遺傳評(píng)估具有重要的意義。標(biāo)準(zhǔn)的種豬測(cè)定多見(jiàn)于瘦肉型豬種，如大約克夏、杜洛克、長(zhǎng)白等，且在此基礎(chǔ)上都已經(jīng)有了深入地研究和較為成熟的理論體系。但就目前狀況而言，部分地方豬種在這方面的研究還未有完善的性狀校正公式，后續(xù)的遺傳評(píng)估等工作也無(wú)法開展。

沙子嶺豬屬于國(guó)家級(jí)畜禽遺傳資源保護(hù)對(duì)象，具有較快的生長(zhǎng)速度，性成熟早，5月齡左右即可進(jìn)行配種，沙子嶺豬適配體重在30～50 kg左右。相對(duì)于瘦肉型豬而言，沙子嶺豬配種時(shí)體重較輕，在此階段進(jìn)行性能測(cè)定保證了沙子嶺豬生長(zhǎng)性狀資料的連續(xù)性而不至于中途數(shù)據(jù)缺失。根據(jù)育種原理和沙子嶺豬實(shí)際生長(zhǎng)環(huán)境，越早對(duì)沙子嶺豬生長(zhǎng)性狀進(jìn)行測(cè)定可以獲得更豐富的性能數(shù)據(jù)，不僅能保證資料的完整性與連續(xù)性，還可以降低疾病風(fēng)險(xiǎn)和育種成本。

本文通過(guò)對(duì)沙子嶺豬進(jìn)行性能測(cè)定，結(jié)合沙子嶺豬實(shí)際生長(zhǎng)過(guò)程，對(duì)動(dòng)物生長(zhǎng)模型進(jìn)行擬合而選擇較為合理模型方程，從而制定目標(biāo)性狀校正公式。

1??材料與方法

1．1 試驗(yàn)材料

本試驗(yàn)研究對(duì)象為湘潭市沙子嶺豬資源場(chǎng)2013年11月至2014年2月出生的70頭后備母豬，試驗(yàn)豬要求有完整的系譜，且有著相同的飼養(yǎng)管理?xiàng)l件和統(tǒng)一的測(cè)定標(biāo)準(zhǔn)與儀器（由于公豬數(shù)目少，因而僅對(duì)母豬生長(zhǎng)發(fā)育性能進(jìn)行了連續(xù)測(cè)定）。

圖1 日齡與體重散點(diǎn)圖

1．2 試驗(yàn)儀器電子秤，籠稱，A超（測(cè)定活體背膘厚，測(cè)定三點(diǎn)取平均值），計(jì)算機(jī)軟件MATLAB（矩陣實(shí)驗(yàn)室）。

1．3 研究方法

1.3.1 動(dòng)物生長(zhǎng)模型的最優(yōu)選擇

其中，Y為因變量觀測(cè)值向量，矩陣X為自變量觀測(cè)值矩陣，β為模型未知參數(shù)的向量。本文研究的生長(zhǎng)日齡與對(duì)應(yīng)體重屬于一元非線性回歸，是上述回歸模型的一般形式。

本文所采用的參數(shù)擬合方法為麥夸特法，又稱為阻尼因子法，是基于牛頓—高斯法得到的矩陣加上阻尼因子而形成。若P=0，則麥夸特法與牛頓-高斯法的計(jì)算結(jié)果一致；若P＞0，則麥夸特法比牛頓-高斯法具有更好的收斂性。

麥夸特法可以從任意初始值開始計(jì)算全局最優(yōu)解，并能有效地減小擬合結(jié)果誤差，在大規(guī)模數(shù)學(xué)優(yōu)化問(wèn)題上，能更簡(jiǎn)便且快速地完成全局計(jì)算任務(wù)[2]。麥夸特法在參數(shù)估計(jì)領(lǐng)域應(yīng)用廣泛，是較為公認(rèn)的最優(yōu)參數(shù)估計(jì)方法。

本文所選擇的動(dòng)物生長(zhǎng)模型為常用的擬合動(dòng)物生長(zhǎng)模型：Logistic模型、Compertz曲線、Saturation模型（飽和生長(zhǎng)模型），指數(shù)增長(zhǎng)模型等等。

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

沙子嶺豬體重與其活體背膘厚關(guān)系模型在此前未有研究，基于圖2所顯示的體重與活體背膘厚的散點(diǎn)圖和畜牧研究者對(duì)其他豬種的研究經(jīng)驗(yàn)，沙子嶺豬與其活體背膘厚模型方程可進(jìn)行線性擬合，其模型方程為：Bf=α+βw，其中，Bf為沙子嶺豬活體背膘厚，w為對(duì)應(yīng)體重， α，β為待定參數(shù)。

圖2 體重與背膘厚散點(diǎn)圖

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

1.3.3 顯著性檢驗(yàn)與擬合優(yōu)度比較

本文只需要對(duì)動(dòng)物生長(zhǎng)模型進(jìn)行擬合優(yōu)度比較，而動(dòng)物生長(zhǎng)模型屬于非線性回歸，因此可以通過(guò)構(gòu)建檢驗(yàn)統(tǒng)計(jì)量R2(Coefficient of Determination，可決系數(shù))來(lái)比較各模型的擬合優(yōu)越性。

1.3.4 性狀校正公式的制定

根據(jù)1.3.1和1.3.2對(duì)生長(zhǎng)模型進(jìn)行最優(yōu)選擇和體重與活體背膘厚線性擬合，再進(jìn)行校正公式的制定。其中，校正公式：目標(biāo)性狀校正值=目標(biāo)性狀實(shí)測(cè)值×C；C為校正系數(shù)，C=性狀的目標(biāo)預(yù)測(cè)值/性狀的實(shí)際預(yù)測(cè)值[3]。

2??結(jié)果與分析

2．1 動(dòng)物生長(zhǎng)模型擬合結(jié)果與分析

根據(jù)1.3.1 所提供的動(dòng)物生長(zhǎng)模型方程，利用數(shù)學(xué)軟件MATLAB進(jìn)行擬合，結(jié)果如表2和圖3所示。

在對(duì)數(shù)據(jù)進(jìn)行多種擬合的情形中，通過(guò)比較RSS、Adjusted R2、RMSE可以判定擬合結(jié)果的優(yōu)劣。RSS越小，RMSE也越小，Adjusted R2則會(huì)越大，說(shuō)明擬合的效果越好。根據(jù)表2結(jié)果顯示，RSS從小到大依次為：Logistic模型，1 892；Saturation模 型，8 237；Compertz模型，9 544；Brody模型，21 870。RMSE從小到大順序依次為：Logistic模 型，2.017；Saturation模型，4.204；Compertz模 型，4.53；Brody模型為6.859。Adjusted R2從大到小依次為：Logistic模型，0.988 9；Compertz模 型，0.954 4；Saturation模型，0.943 9；Brody模型，0.871 5。Logistic模型3個(gè)指標(biāo)優(yōu)于其他3種模型的擬合結(jié)果，而Saturation模型雖在RSS和RMSE兩方面優(yōu)于Compertz模型，Compertz的Adjusted R2卻略大于Saturation模型，但是區(qū)別不大。Brody模型的參數(shù)擬合結(jié)果與其置信區(qū)間上下限距離較遠(yuǎn)，根據(jù)切比雪夫不等式可知參數(shù)的估計(jì)方差較大，說(shuō)明其離散程度大，而RSS、Adjusted R2、RMSE3項(xiàng)指標(biāo)均遜于其他3種模型方程，表明Brody模型的模擬結(jié)果不理想。

表2 各動(dòng)物生長(zhǎng)模型方程擬合結(jié)果

圖3 動(dòng)物生長(zhǎng)模型方程擬合對(duì)比圖

從圖3可以更為直觀地看出，Brody模型在起始點(diǎn)位置與終點(diǎn)要高出其他3種模型，而其他3種模型在起點(diǎn)幾乎位于相同位置。終點(diǎn)處，Saturation模型也稍微偏離點(diǎn)的密集處，只有Logistic與Compertz模型集中在點(diǎn)較多的位置。整體而言Logistic模型、Compertz模型呈較為明顯的“S”型曲線，Saturation模型為類似于直線的曲線，但都很好地穿過(guò)了點(diǎn)較為集中的地方。因此，從擬合結(jié)果來(lái)看，除Brody模型以外，其他3種模型理論上均可以進(jìn)行生長(zhǎng)性狀的校正公式制定。

2．2 體重與活體背膘厚擬合結(jié)果與分析

根據(jù)1.3.2內(nèi)容，通過(guò)利用MATLAB統(tǒng)計(jì)工具箱中自帶的regstats函數(shù)，并調(diào)用reglm函數(shù)作一元線性回歸[4]，分析結(jié)果如表3、表4所示。

方差分析表中的P=0.000 1，說(shuō)明整個(gè)回歸方程是極顯著的；參數(shù)估計(jì)表中P值均小于0.01，同時(shí)說(shuō)明回歸方程中2個(gè)參數(shù)的估計(jì)值均為極顯著。但后備母豬群體所呈現(xiàn)的性能離散程度太大，引起這種情況的可能原因?yàn)椋孩偌膊?dǎo)致部分后備母豬生長(zhǎng)發(fā)育情況不正常；②測(cè)定人員的測(cè)定水平有待提高；③后備母豬群體不夠大?？傮w而言，后備母豬活體背膘厚與質(zhì)量的回歸方程是合理的。

2．3 校正公式制定結(jié)果與分析

根據(jù)湘潭市沙子嶺豬資源場(chǎng)實(shí)際情況，要考慮其產(chǎn)仔日期的集中度，很多母豬初配月齡要超過(guò)5個(gè)月。根據(jù)圖1所呈現(xiàn)的日齡與體重的關(guān)系，5月齡左右，沙子嶺后備母豬體重在30～50 kg，同時(shí)兼顧沙子嶺豬資源場(chǎng)的實(shí)際情況，可將50 kg作為校正公式的基準(zhǔn)。

表3 方差分析表

表4 參數(shù)估計(jì)表

2.3.1 達(dá)50 kg體重日齡的校正公式

根據(jù)2.1的擬合結(jié)果可知Logistic模型，Compertz模型及Saturation模型均可以作為校正公式制定的待定曲線。根據(jù)1.3.4內(nèi)容可知，目標(biāo)性狀校正值=目標(biāo)性狀實(shí)測(cè)值×C，則×C，校正系數(shù)

根據(jù)上述公式推導(dǎo)過(guò)程可發(fā)現(xiàn)如果實(shí)測(cè)體重大于k值，則校正公式會(huì)無(wú)意義，因此此校正公式的限制條件為：w實(shí)＜52.81 。利用相同的方法可以推導(dǎo)出另外2種模型的校正公式（表5）。

2.3.2 達(dá)50 kg的活體背膘厚校正公式

根據(jù)上述內(nèi)容可知沙子嶺后備母豬活體背膘厚與體重的線性回歸方程為：

表5 各模型制定的達(dá)50 kg日齡校正公式

Bf=15.0786+0.1837w

參照2.3.1的推導(dǎo)過(guò)程，利用基于Saturation方程制定的達(dá)50 kg體重日齡的校正公式結(jié)合活體背膘厚與體重的線性回歸方程可以得到活體背膘厚的校正公式為：

3??討論

3．1 模型選用的探討

Logistic模型、Compertz模型與Saturation模型在數(shù)據(jù)擬合過(guò)程中有著較高的擬合度，與杜洛克、長(zhǎng)白豬及金華豬等品種的研究有著類似結(jié)果，但Saturation模型不能很好地刻畫動(dòng)物生長(zhǎng)發(fā)育規(guī)律，其他兩者在這方面有優(yōu)勢(shì)，也并不說(shuō)明這幾種模型是最優(yōu)的，可能還存在其他模型能夠得到更好地?cái)M合。體重與活體背膘厚的線性模型的擬合度不高，但是線性模型方程各項(xiàng)分析均為極顯著，國(guó)內(nèi)其他研究表明在瘦肉型豬上，二者關(guān)系也可以用線性模型來(lái)描述。

本文模型方程的選擇是在前人對(duì)其他品種進(jìn)行研究經(jīng)驗(yàn)上進(jìn)行的，且在最初數(shù)據(jù)選用上，并不為了刻意呈現(xiàn)某種狀態(tài)而人為地?zé)o根據(jù)地去掉部分?jǐn)?shù)據(jù)，由于湘潭市沙子嶺豬場(chǎng)沙子嶺后備母豬群體數(shù)目相對(duì)有限，對(duì)沙子嶺豬其他分布區(qū)的資料收集較為困難，在散點(diǎn)圖上呈現(xiàn)了沙子嶺豬各階段體重離散程度較大的狀態(tài)。擬合過(guò)程中，基于殘差及其置信區(qū)間的分析，通過(guò)MATLAB軟件將異常點(diǎn)剔除后再重新擬合。因此，所有的模型參數(shù)值是“分析—剔除—再分析”的結(jié)果。

3．2 校正公式制定結(jié)果討論

基于Saturation模型制定沙子嶺豬達(dá)50 kg體重日齡校正公式，在考慮綜合因素下，優(yōu)于其他2種結(jié)果，而達(dá)50 kg體重的活體背膘厚校正公式，也是通過(guò)對(duì)以前研究經(jīng)驗(yàn)加以總結(jié)的結(jié)果，2個(gè)公式變換后與美國(guó)和加拿大的早期性狀校正公式結(jié)構(gòu)類似。對(duì)上述公式，利用假設(shè)檢驗(yàn)對(duì)真實(shí)值與校正值進(jìn)行比較，P＞0.05，說(shuō)明2組數(shù)據(jù)差異不顯著，公式在現(xiàn)有的條件下有一定的可信度。楊飛來(lái)[5]的研究分析了大白豬以50 kg作為標(biāo)準(zhǔn)體重代替以100 kg作為標(biāo)準(zhǔn)體重可行性，但用了線性模型來(lái)描述日齡與體重的關(guān)系。陳斌[6]利用非線性模型將豬早期的體重直接校正到達(dá)100 kg體重日齡，研究方法值得借鑒，但通用性不強(qiáng)。

本研究制定的沙子嶺豬性能測(cè)定資料校正公式為沙子嶺豬性狀的遺傳評(píng)估做了鋪墊，保證了沙子嶺豬待選目標(biāo)性狀的統(tǒng)一性，便于計(jì)算相應(yīng)性狀的遺傳參數(shù)，對(duì)于后續(xù)育種方案與優(yōu)化等工作也有著重要作用。

[1] 劉望宏， 胡軍勇， 倪德斌，等. 瘦肉型豬早期生長(zhǎng)發(fā)育規(guī)律的擬合模型研究[J]. 華中農(nóng)業(yè)大學(xué)學(xué)報(bào)， 2010， 29(3)：335-340.

[2] 曾蒙秀，宋友桂.麥夸特算法在X射線衍射物相定量分析中的應(yīng)用[J].地球科學(xué)：中國(guó)地質(zhì)大學(xué)學(xué)報(bào)，2013，38(2)：431-440.

[3] 張勤， 張啟能. 生物統(tǒng)計(jì)學(xué)[M]. 北京：中國(guó)農(nóng)業(yè)大學(xué)出版社，2003：139-140.

[4] 謝中華， 李國(guó)棟， 劉煥進(jìn)， 等.MATLAB從零到進(jìn)階[M]. 北京：北京航空航天大學(xué)出版社， 2012：375-387.

[5] 楊飛來(lái).應(yīng)用多性狀動(dòng)物模型BLUP法進(jìn)行豬的遺傳評(píng)定研究[D].長(zhǎng)沙：湖南農(nóng)業(yè)大學(xué)， 2002.

[6] 陳斌.瘦肉型豬的場(chǎng)內(nèi)遺傳評(píng)估及遺傳分析研究[D].長(zhǎng)沙：湖南農(nóng)業(yè)大學(xué)，2005.

2014-09-25 ）

國(guó)家現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術(shù)體系建設(shè)專項(xiàng)資金（CARS-36）