不同品種的后備種豬體尺預(yù)測體重校正公式的研究

孟慶利，張 軍，王志云，周海深

（北京養(yǎng)豬育種中心，北京 海淀 100194）

畜禽各表型性狀是遺傳因素和環(huán)境因素共同作用的結(jié)果，各性狀的遺傳力也差距很大。體尺是一種重要的性狀，通過體尺的測定，可以有效地預(yù)測種豬的體重，目前國內(nèi)外常用到的利用體尺預(yù)測豬體重的預(yù)測公式如下：蘭巨生（1963年）提出的預(yù)測公式，李劍秋（2002）[1]提出的預(yù)測公式。本文通過大量的后備豬體尺測定，通過典型相關(guān)分析可以發(fā)現(xiàn)，胸圍、體長這兩個體尺性狀對體重的預(yù)測能力最高，很多研究者也多用胸圍×體長或胸圍2×體長來預(yù)估體重。本研究利用原始數(shù)據(jù)分別計算了這兩個指標(biāo)，然后利用SAS的GLM過程得到其對體重性狀的回歸方程。

1 材料與方法

1.1 試驗(yàn)時間與地點(diǎn)

試驗(yàn)時間：2018年1月至2019年3月；試驗(yàn)地點(diǎn)：北京養(yǎng)豬育種中心種豬場。

1.2 試驗(yàn)動物與測定指標(biāo)

選擇200日齡左右的大白、長白和杜洛克后備豬1 258、372和236頭進(jìn)行種豬體尺指標(biāo)體長、體高、胸圍、管圍和體重的測定工作。

體長：在種豬正常站立時，從兩耳根中點(diǎn)聯(lián)線的中部起，用卷尺沿背脊量到尾根的第一自然輪紋為止。

體高：在種豬正常站立時，肩部最高到地面的垂直距離。

胸圍：在種豬正常站立時，在肩胛骨后緣用皮尺測量胸部的垂直周徑，松緊度以皮尺自然貼緊毛皮為宜。

管圍：在種豬正常站立時，左前肢管部最細(xì)處的周徑。

2 結(jié)果與分析

2.1 數(shù)據(jù)質(zhì)控

2.1.1 表型數(shù)據(jù)奇異點(diǎn)篩選

（1）同一個耳號，同一測定日期，有兩次記錄，個別性狀兩次記錄有差異，刪除16條記錄，編號分別是：227、237、662、673、764、814、836、850、1401、1405、1416、1429、1553、1563、1844、1869。

（2）編號為104、1491個體的胸圍分別是11.2 cm和14 cm，刪除此2條記錄。

（3）編號為 147、959、1037個體的體重分別是193 kg、1 226.5 kg和13.5 kg，刪除此3條記錄。

（4）編號為870個體的管圍為9.6 cm，刪除此條記錄。

（5）編號為833、1081個體的體高分別是8 cm和7 cm，刪除此2條記錄。

（6）編號為280、1932個體的體長分別是1 026.2 cm和15 cm，刪除此2條記錄。

2.1.2 體重與測定日齡的關(guān)系 根據(jù)圖1測定日齡和體重的散點(diǎn)圖，發(fā)現(xiàn)編號為1835、1828、1830、1826、1823、1822、1836、1829、1825、1831、1818、1821、1824、1820、1819、1827、1834、1833、1832 的 個 體 為獨(dú)立一群，刪除19條記錄。

2.1.3 體重與體長關(guān)系 根據(jù)圖2體長和體重的散點(diǎn)圖，發(fā)現(xiàn)編號為 1676、1680、1679、1670、1671、1681、1675、1658、1677、1640、1645、1655、1672、1668、1639、1648、1654、1685、1688、1644、1656、1665、1666、1674、1643、1651、1682、1684、1660、1661、1667、1662、1649、1669、1641、1646、1647、1652、1673、1657、1683、1689、1686、1659、1687、1642、1663、1653、1678、1650、1664的個體為獨(dú)立一群，刪除51條記錄。體長的預(yù)測能力最高。從該結(jié)果發(fā)現(xiàn)，可以用胸圍和體長兩個體尺性狀來預(yù)估體重性狀。

2.3 體重預(yù)估模型構(gòu)建

2.3.1 所有個體

（1）多變量模型。多變量模型是構(gòu)建多個自變量對1個依變量回歸關(guān)系的模型。本研究使用SAS的REG過程，首先構(gòu)建多變量回歸方程，之后基于Cp標(biāo)準(zhǔn)（概念預(yù)測標(biāo)準(zhǔn)）和AIC準(zhǔn)則篩選出最佳模型，具體模型見表2，表中的R2值為擬合度，體重的單位為千克（kg），4個體尺性狀的單位為厘米（cm）。

（2）單變量分析。通過典型相關(guān)分析可以發(fā)現(xiàn)，胸圍、體長這兩個體尺性狀對體重的預(yù)測能力最高，很多研究者也多用胸圍×體長或胸圍2×體長來預(yù)估體重。本研究利用原始數(shù)據(jù)分別計算了這兩個指標(biāo)，然后利用SAS的GLM過程得到其對體重性狀的回歸方程，見表2。

2.2 體重與體尺的相關(guān)分析

2.2.1 簡單相關(guān)系數(shù)與偏相關(guān)系數(shù) 表1為所有豬只體重與各體尺性狀的相關(guān)系數(shù)，其中簡單相關(guān)系數(shù)一列為簡單相關(guān)（Pearson相關(guān)），是表示兩變量間的直線關(guān)系。當(dāng)對多個變量的相關(guān)性進(jìn)行分析時，僅僅根據(jù)兩個變量間的簡單相關(guān)系數(shù)往往不能很好地說明彼此間的相關(guān)關(guān)系，只有去除其他變量影響的情況下，計算出它們之間的相關(guān)系數(shù)，才能更好地解釋變量間的相關(guān)情況，故而計算了體重與各體尺性狀的偏相關(guān)系數(shù)（剔除其他所有變量）。表1中所有的相關(guān)系數(shù)和偏相關(guān)系數(shù)經(jīng)統(tǒng)計檢驗(yàn)后，其P值均小于10-4。從結(jié)果中可以看到體重與體長、胸圍的相關(guān)程度高于體高和管圍性狀。

表2 不同體尺性狀預(yù)估體重性狀的回歸方程

2.3.2 不同品系的體重預(yù)估模型構(gòu)建 根據(jù)不同品種分別估計體重性狀的回歸方程，見表3。

表1 體重與各體尺性狀的相關(guān)系數(shù)

2.2.2 體重與各體尺性狀的典型相關(guān)分析 典型相關(guān)分析是度量兩組變量間的整體相關(guān)性。利用SAS的CANCORR過程，構(gòu)建體尺性狀組的典型變量V1=0.407 2×體長+0.117 1×體高+0.583 4×胸圍-0.091 1×管圍，經(jīng)過Wilks'Lambda檢驗(yàn)，P為2.78×10-6，達(dá)到統(tǒng)計顯著。體重性狀的變異可以被V1變量解釋64.57%。豬胸圍、體長、體高和管圍的判定系數(shù)分別是 0.582 1、0.521 7、0.369 9 和 0.005 6，說明胸圍、體長和體高對體重有一定的預(yù)測能力，且胸圍和

表3 不同品種利用體尺性狀預(yù)估體重的回歸方程

3 分析與討論

本研究結(jié)果表明，胸圍、體長和體高對體重有一定的預(yù)測能力，且胸圍和體長的預(yù)測能力最高。通過SAS的GLM過程得到不同體尺性狀預(yù)估體重性狀的回歸方程和不同品種利用體尺性狀預(yù)估體重的回歸方程。結(jié)果為：杜洛克豬通過體尺指標(biāo)預(yù)測體重R2值最高，都在0.7以上；長白豬通過體尺指標(biāo)預(yù)測體重R2值居中，都在0.6以上、0.7以下；大白豬通過體尺指標(biāo)預(yù)測體重R2值最低都在0.5～0.6之間?？傮w來說，通過細(xì)分品種，通過體尺預(yù)測體重R2值都在0.5以上，校正公式都具有實(shí)用性。將所有數(shù)據(jù)不分品種進(jìn)行分析得到不同體尺性狀預(yù)估體重性狀的回歸方程R2值在0.6～0.7之間，同樣具有較強(qiáng)的實(shí)用性，但沒有細(xì)分品種所得出的R2值更精細(xì)，說明不同品種之間通過體尺預(yù)測體重具有一定的差異性。在實(shí)際應(yīng)用過程中，最好根據(jù)品種制定校正公式，準(zhǔn)確性會更高。