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

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

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

1 材料與方法

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

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

1.2 試驗(yàn)動(dòng)物與測(cè)定指標(biāo)

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

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

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

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

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

2 結(jié)果與分析

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

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

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

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

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

（4）編號(hào)為870個(gè)體的管圍為9.6 cm，刪除此條記錄。

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

（6）編號(hào)為280、1932個(gè)體的體長(zhǎng)分別是1 026.2 cm和15 cm，刪除此2條記錄。

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

2.1.3 體重與體長(zhǎng)關(guān)系 根據(jù)圖2體長(zhǎng)和體重的散點(diǎn)圖，發(fā)現(xiàn)編號(hào)為 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的個(gè)體為獨(dú)立一群，刪除51條記錄。體長(zhǎng)的預(yù)測(cè)能力最高。從該結(jié)果發(fā)現(xiàn)，可以用胸圍和體長(zhǎng)兩個(gè)體尺性狀來(lái)預(yù)估體重性狀。

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

2.3.1 所有個(gè)體

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

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

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

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

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

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

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

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

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

3 分析與討論

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