崔登帥，熊三亞，鄭浩，李龍?jiān)?，余乃?biāo)，黃智勇，肖石軍，郭源梅

里岔黑豬和巴里雜交豬母豬100kg體重日齡校正方法的比較

崔登帥1，熊三亞1，鄭浩2，李龍?jiān)?，余乃標(biāo)1，黃智勇2，肖石軍1，郭源梅1

1江西農(nóng)業(yè)大學(xué)省部共建種豬遺傳改良與養(yǎng)殖技術(shù)國(guó)家重點(diǎn)實(shí)驗(yàn)室，南昌 330045；2江西山下華系種豬養(yǎng)殖有限公司，江西贛州 341900

【目的】比較3種達(dá)100 kg體重日齡（D100kg）的校正系數(shù)的校正效果，為種豬場(chǎng)根據(jù)自身實(shí)際情況選擇合適的D100kg校正系數(shù)提供參考?！痉椒ā恳岳锊砗谪i群體及其與巴克夏雜交1代（巴里1代）和2代（巴里2代）共1 344頭豬只為研究對(duì)象，利用其生長(zhǎng)肥育期體重和日齡的數(shù)據(jù)，以始測(cè)體重≥25 kg、終測(cè)體重在80—125 kg、過(guò)始測(cè)點(diǎn)和終測(cè)點(diǎn)直線的斜率在群體平均斜率±2個(gè)標(biāo)準(zhǔn)差之內(nèi)等標(biāo)準(zhǔn)進(jìn)行質(zhì)控，然后按照下面描述的三種方法估計(jì)里岔黑豬、巴里1代、巴里2代以及它們混合群體達(dá)100 kg校正系數(shù)，并評(píng)估它們的校正效果。一是回歸校正系數(shù)（reg）：利用測(cè)定期間個(gè)體的體重和日齡數(shù)據(jù)，擬合體重對(duì)日齡的一元回歸方程，把其斜率記1i；然后利用它終測(cè)點(diǎn)和原點(diǎn)（日齡和體重均為0）作一條直線，把其斜率記為2i；1i/2i的平均數(shù)就是reg。二是兩點(diǎn)校正系數(shù)（2point）：利用個(gè)體初測(cè)點(diǎn)和終測(cè)點(diǎn)作一條直線，得到斜率1i，斜率2i及2point計(jì)算方法與reg計(jì)算方法相同。三是單點(diǎn)校正系數(shù)（1point）：擬合所有個(gè)體的終測(cè)體重對(duì)終測(cè)日齡的一元回歸方程，得到斜率1；第二條直線斜率與reg相同；1/2i的平均數(shù)就是1point。分別將CF、2point和1point代入校正公式，計(jì)算出相應(yīng)的D100kg.CF.reg、D100kg.2point和D100kg.1point。根據(jù)前面擬合的日齡對(duì)體重的一元回歸，計(jì)算個(gè)體達(dá)100kg的回歸日齡（D100kg.reg），并計(jì)算它與其他D100kg的相關(guān)系數(shù)和離差。根據(jù)相關(guān)系數(shù)和離差的大小，評(píng)估3種方法的優(yōu)劣。為比較3種校正系數(shù)的群體適用性，用某個(gè)群體的去校正其他3個(gè)群體的D100kg，并根據(jù)其校正日齡與該群體D100kg.reg相關(guān)系數(shù)的高低來(lái)評(píng)估其群體適用性的優(yōu)劣。通過(guò)分析終測(cè)體重和校正誤差之間的關(guān)系，探究終測(cè)體重對(duì)D100kg的影響?？紤]單點(diǎn)校正誤差與終測(cè)體重之間存在很強(qiáng)的線性關(guān)系，利用終測(cè)體重對(duì)D100kg.1point進(jìn)行二次校正，以提高單點(diǎn)校正的準(zhǔn)確性?！窘Y(jié)果】質(zhì)控后共1 181頭母豬用于后續(xù)分析。在同一群體中，reg和2point的值相近，它們極顯著大于1point。D100kg.CF.reg和D100kg.2point與D100kg.reg很接近，相關(guān)系數(shù)均大于0.98，說(shuō)明這兩種方法的校正效果很好。D100kg.1point與D100kg.reg的差異較大，相關(guān)系數(shù)僅為0.89，校正效果較差。reg和2point的群體適應(yīng)性很好，均可用于其他群體D100kg的校正。1point的群體適應(yīng)性較差，不宜用于其他群體D100kg的校正。4個(gè)群體中，終測(cè)體重離目標(biāo)體重（100 kg）越遠(yuǎn)，D100kg.CF.reg校正誤差越大，D100kg.2point校正誤差沒(méi)有明顯的變化。D100kg.1point的誤差與終測(cè)體重均存在極顯著的負(fù)相關(guān)，分別為-0.9812、-0.9627、-0.9786和-0.9352。經(jīng)過(guò)二次校正后，1point校正效果有顯著提升，但與reg校正效果仍有較大差距，優(yōu)先選擇reg?！窘Y(jié)論】reg和2point的校正效果和群體適應(yīng)性都很好，可用于生產(chǎn)實(shí)踐。單點(diǎn)校正系數(shù)（1point）的校正效果不理想，不宜用于生產(chǎn)實(shí)踐。在校正D100kg時(shí)，可以根據(jù)豬只的數(shù)據(jù)量，采用合適的方法，以提高估計(jì)準(zhǔn)確性，降低估計(jì)誤差。

里岔黑豬；達(dá)100 kg體重日齡；校正系數(shù)；回歸日齡

0 引言

【研究意義】我國(guó)是全球養(yǎng)豬第一大國(guó)，生豬存欄量和豬肉產(chǎn)量均占世界的半壁江山[1-5]。目前，國(guó)內(nèi)肉豬以“外三元”為主[6-7]，雖然生長(zhǎng)速度快，瘦肉率高，但是肉質(zhì)較差，難以滿足人們對(duì)于優(yōu)質(zhì)豬肉的需求[8-9]。我國(guó)地方豬種肉質(zhì)好，能夠滿足人們對(duì)于優(yōu)質(zhì)豬肉的需求，但生長(zhǎng)速度有待提高[10-12]。達(dá)100 kg體重日齡（D100kg）是度量生長(zhǎng)速度最常用的指標(biāo)[13-14]，也是種豬選育的目標(biāo)性狀之一[15]。由于在飼養(yǎng)過(guò)程中很難對(duì)豬只的體重進(jìn)行實(shí)時(shí)監(jiān)測(cè)，且豬只的體重受采食、飲水、排泄等行為的影響，因此不能準(zhǔn)確測(cè)定豬只D100kg的實(shí)際日齡[16]。在生產(chǎn)實(shí)踐中，對(duì)75—115 kg豬只進(jìn)行稱重，然后利用校正系數(shù)通過(guò)校正公式來(lái)計(jì)算D100kg[17]。精準(zhǔn)校正D100kg有利于加快我國(guó)地方豬種生長(zhǎng)速度的遺傳改良，克服我國(guó)地方豬種生長(zhǎng)緩慢的不足?！厩叭搜芯窟M(jìn)展】校正系數(shù)主要有兩種：一種是2016年以前國(guó)家生豬產(chǎn)業(yè)體系使用的校正系數(shù)[18]，另一種是目前國(guó)家生豬產(chǎn)業(yè)體系使用的校正系數(shù)[19]。筆者前期的研究結(jié)果顯示這兩種校正系數(shù)的校正效果沒(méi)有顯著差異，因此本研究只選用第一種校正系數(shù)進(jìn)行后續(xù)的分析[16]。目前，只有少數(shù)地方品種對(duì)D100kg校正系數(shù)進(jìn)行了估計(jì)[16, 20]，絕大部分地方豬種D100kg校正系數(shù)未見(jiàn)報(bào)道?！颈狙芯壳腥朦c(diǎn)】目前估計(jì)D100kg校正系數(shù)需要用到豬只至少2個(gè)時(shí)間點(diǎn)體重和日齡的數(shù)據(jù)，如杜洛克、長(zhǎng)白、大白豬、里岔黑豬等[21-30]。時(shí)間點(diǎn)越多，估計(jì)的D100kg校正系數(shù)就越準(zhǔn)確，但耗費(fèi)人力物力也就越多。如果可以在不影響準(zhǔn)確性的前提下只用一個(gè)時(shí)間點(diǎn)的數(shù)據(jù)來(lái)估計(jì)，則可以大幅度節(jié)約成本?！緮M解決的關(guān)鍵問(wèn)題】本研究以里岔黑豬及其與巴克夏公豬雜交的F1和F2代為研究對(duì)象，利用生長(zhǎng)肥育期測(cè)定的體重和日齡數(shù)據(jù)，探究數(shù)據(jù)量與的校正效果之間關(guān)系。希望探索出一種能夠簡(jiǎn)化測(cè)定，節(jié)約成本，并且適用于推廣應(yīng)用的D100kg校正方法。

1 材料與方法

1.1 試驗(yàn)豬只

研究所用的試驗(yàn)豬只均來(lái)自于江西山下華系種豬養(yǎng)殖有限公司定南陽(yáng)林山下原種豬場(chǎng)（江西，定南）[31]，實(shí)驗(yàn)數(shù)據(jù)采集于2014—2019年。飼養(yǎng)管理和營(yíng)養(yǎng)水平參見(jiàn)文獻(xiàn)[16]，豬只自由采食和飲水。通過(guò)潤(rùn)農(nóng)性能測(cè)定系統(tǒng)（深圳，潤(rùn)農(nóng)科技），測(cè)定試驗(yàn)豬只生長(zhǎng)肥育期的體重，并記錄測(cè)定日期。本研究利用前期報(bào)道的1 344頭豬只生長(zhǎng)肥育期的體重和日期[16]，其中里岔黑豬204頭（閹公豬23頭，母豬181頭）、巴克夏和里岔黑豬（巴里）雜交1代644頭（閹公豬31頭，母豬613頭）和2代496頭（閹公豬64頭，母豬432頭）。

1.2 試驗(yàn)方法

1.2.1 數(shù)據(jù)質(zhì)控 這1 344頭豬已經(jīng)通過(guò)了以下質(zhì)控標(biāo)準(zhǔn)：測(cè)定期≥30 d；始測(cè)體重≥25 kg，終測(cè)體重在80—125 kg；始測(cè)日齡>60 d和終測(cè)日齡<280 d；優(yōu)化后體重對(duì)日齡一元回歸方程擬合度≥0.95的個(gè)體[16]。本研究新增一個(gè)質(zhì)控標(biāo)準(zhǔn)，即經(jīng)過(guò)始測(cè)點(diǎn)和終測(cè)點(diǎn)直線的斜率在群體平均斜率2個(gè)標(biāo)準(zhǔn)差之內(nèi)。

1.2.2 D100kg校正系數(shù)的估計(jì) 由于群體對(duì)于校正系數(shù)有顯著的影響[16]，因此分群體估計(jì)D100kg校正系數(shù)。根據(jù)之前的報(bào)道，巴里雜交1代和2代母豬的非常接近，分別為1.389和1.380，且與里岔黑豬的（1.346）差異也很小[16]。另外，由于里岔黑豬是巴里雜交豬的一個(gè)親本，與它們有遺傳關(guān)系。因此，將3個(gè)群體合并一個(gè)群體，即混合群體，加大樣本量，提高估計(jì)的準(zhǔn)確性。根據(jù)使用數(shù)據(jù)量的不同，采用3種方法計(jì)算豬只D100kg的校正系數(shù)：

一是回歸校正系數(shù)（reg）：利用自動(dòng)飼喂系統(tǒng)生長(zhǎng)肥育期記錄的個(gè)體所有體重和日齡數(shù)據(jù)，擬合體重對(duì)日齡的一元回歸方程，把其斜率記為1i。然后利用它的終測(cè)點(diǎn)和原點(diǎn)（日齡和體重均為0）作一條直線，把其斜率記為2i，reg就等于1i/2i的平均數(shù)。

二是兩點(diǎn)校正系數(shù)（2point）：利用個(gè)體始測(cè)點(diǎn)（始測(cè)體重，始測(cè)日齡）和終測(cè)點(diǎn)（終測(cè)體重，終測(cè)日齡）的數(shù)據(jù)，作一條直線，把其斜率記為1i。剩下的步驟和回歸校正系數(shù)的步驟完全相同。

三是單點(diǎn)校正系數(shù)（1point）：使用所有個(gè)體的終測(cè)數(shù)據(jù)，擬合終測(cè)體重對(duì)終測(cè)日齡的一元回歸方程，把其斜率記為1。然后利用個(gè)體終測(cè)點(diǎn)和原點(diǎn)（日齡和體重均為0）作為一條直線，把其斜率記為2i，1point為1/2i的平均數(shù)。

1.2.3 D100kg的計(jì)算 利用生長(zhǎng)期的體重和日齡，對(duì)每個(gè)個(gè)體擬合日齡對(duì)體重的一元回歸方程。把目標(biāo)體重100 kg代入回歸方程即可計(jì)算出回歸D100kg，記為D100kg.reg，并將其看作為真實(shí)的D100kg。

把估計(jì)的校正系數(shù)（reg、2point和1point）、終測(cè)體重和終測(cè)日齡代入以下公式[18]，就得到個(gè)體的D100kg（D100kg.CF.reg、D100kg.2point和D100kg.1point）：

利用R語(yǔ)言的lm和aov進(jìn)行方差分析，通過(guò)agricolae包中的duncan.test函數(shù)進(jìn)行多重比較的校正，比較不同群體D100kg之間的差異。

1.2.4 校正日齡之間的相關(guān)系數(shù) 利用R語(yǔ)言的pairs函數(shù)，作不同方法估計(jì)的D100kg之間的散點(diǎn)圖，然后用cor.test函數(shù)，計(jì)算D100kg之間的相關(guān)系數(shù)，并進(jìn)行顯著性檢驗(yàn)。將D100kg.reg看作為真實(shí)的D100kg，其他校正日齡與它的相關(guān)系數(shù)越高，說(shuō)明該方法估計(jì)的校正日齡越準(zhǔn)確；反之，越不準(zhǔn)確。

1.2.5 校正日齡的無(wú)偏性 利用校正誤差，即校正日齡與回歸日齡的差值，來(lái)評(píng)估校正方法的無(wú)偏性。校正誤差的平均數(shù)離原點(diǎn)越近越無(wú)偏，越遠(yuǎn)偏差越大。如果平均數(shù)為負(fù)，即偏??；反之，則偏大。校正誤差與終測(cè)體重之間的散點(diǎn)圖能進(jìn)一步顯示校正日齡的無(wú)偏性：如果校正日齡關(guān)于=0對(duì)稱，則說(shuō)明是無(wú)偏的，反之則有偏。

1.2.6 校正系數(shù)的適用性 以某群體的校正系數(shù)去估計(jì)其他群體的D100kg，利用R語(yǔ)言cor.test函數(shù)，計(jì)算它與該群體D100kg.reg之間的相關(guān)系數(shù)，以此來(lái)評(píng)估校正系數(shù)的適用性。相關(guān)系數(shù)越大，則說(shuō)明校正系數(shù)的適用性越好；反之，校正系數(shù)的適用性越差。

1.2.7 終測(cè)體重對(duì)D100kg的影響 據(jù)前人的研究結(jié)果，終測(cè)體重對(duì)D100kg準(zhǔn)確性有顯著的影響[32]。終測(cè)體重與目標(biāo)體重的差異越大，校正D100kg的誤差也越大。通過(guò)終測(cè)體重和校正誤差（校正D100kg與D100kg.reg的離差）的散點(diǎn)圖，可以看出D100kg.1point的誤差與終測(cè)體重存在線性關(guān)系。利用R語(yǔ)言的cor.test函數(shù)，計(jì)算D100kg.1point的誤差與終測(cè)體重的相關(guān)系數(shù)，并用lm函數(shù)擬合它們之間的一元回歸公式。

1.2.8 D100kg.1point的二次校正 由于單點(diǎn)校正誤差與終測(cè)體重（BWend）之間存在很強(qiáng)的線性關(guān)系，且單點(diǎn)校正效果并不理想，因此建立單點(diǎn)校正誤差對(duì)終測(cè)體重的回歸方程，把體重和回歸方程的斜率（）代入下列公式，對(duì)D100kg.1point進(jìn)行二次校正（D100kg.1point.cor），以提高單點(diǎn)校正的準(zhǔn)確性：

D100kg.1point.cor=D100kg.1point+(BWend-100)。

2 結(jié)果

2.1 數(shù)據(jù)質(zhì)控結(jié)果

在這1 344個(gè)個(gè)體中，一共1 290個(gè)個(gè)體通過(guò)了質(zhì)控，其中里岔黑豬197頭（母豬175頭，閹公豬22頭），巴里1代621頭（母豬593頭，閹公豬28頭），巴里2代472頭（母豬413頭，閹公豬59頭）。由于性別和群體對(duì)校正系數(shù)影響顯著，且3個(gè)群體中閹公豬數(shù)量都較少，影響閹公豬校正系數(shù)估計(jì)的準(zhǔn)確性，因此，僅1 181頭母豬用于后續(xù)分析。通過(guò)質(zhì)控個(gè)體始測(cè)點(diǎn)和終測(cè)點(diǎn)連線的趨勢(shì)基本相似（圖1），即豬只的生長(zhǎng)規(guī)律相似，說(shuō)明質(zhì)控后的數(shù)據(jù)質(zhì)量較好。

2.2 校正系數(shù)CF的估計(jì)

在里岔黑豬、巴里1代和2代以及混合群體中，3種方法估計(jì)的校正系數(shù)的簡(jiǎn)單統(tǒng)計(jì)量如表1所示。在同一個(gè)群體中，reg和2point的值相近，但它們極顯著（＜2.2×10-16）大于1point。

2.3 校正日齡D100kg的估計(jì)

回歸日齡（D100kg.reg）和三種方法估計(jì)的D100kg見(jiàn)表2。在這4個(gè)群體中，D100kg.reg、reg校正的D100kg（D100kg.CF.reg）和2point校正的D100kg（D100kg.2point）比較接近，特別是D100kg.reg和D100kg.CF.reg（≥0.0613）。它們與1point校正的D100kg（D100kg.1point）差異均比較大，均達(dá)到了極顯著水平（≤2.56×10-13）。在里岔黑豬中，D100kg.reg和D100kg.CF.reg與D100kg.2point差異不顯著（≥0.4500），但在另外2個(gè)群體中，它們之間的差異均達(dá)到極顯著水平（≤2.94×10-6）。

里岔黑豬4種方法估計(jì)的D100kg均極顯著（＜2.2×10-16）大于巴里1代和巴里2代相應(yīng)方法估計(jì)的D100kg。除D100kg.1point（=0.0568）外，剩余3種方法校正的巴里2代D100kg均極顯著（≤3.76×10-8）大于巴里1代的D100kg。

圖1 質(zhì)控后個(gè)體始測(cè)點(diǎn)和終測(cè)點(diǎn)的連線圖

表1 校正系數(shù)CF的簡(jiǎn)單統(tǒng)計(jì)量

表2 達(dá)100kg日齡的簡(jiǎn)單統(tǒng)計(jì)量

2.4 校正日齡之間的相關(guān)性

4種方法估計(jì)的D100kg散點(diǎn)圖見(jiàn)圖2，它們之間的簡(jiǎn)單相關(guān)系數(shù)見(jiàn)表3。除D100kg.1point外，D100kg之間均存在極強(qiáng)的正相關(guān)，≥0.9796，＜2.2×10-16。D100kg.1point與另外3種D100kg之間也存在極強(qiáng)的正相關(guān)，但相關(guān)系數(shù)沒(méi)有后三者之間那么強(qiáng)，0.8873≤≤ 0.8999，＜2.2×10-16。

D100kg與D100kg.reg相關(guān)系數(shù)可以用來(lái)度量其估計(jì)的準(zhǔn)確性。D100kg.CF.reg和D100kg.2point與D100kg.reg相關(guān)系數(shù)均在0.95以上，分別為0.9925和0.9796。因此，這兩種方法估計(jì)的D100kg很準(zhǔn)確。D100kg.1point與D100kg.reg相關(guān)系數(shù)為0.8873，顯著的小于D100kg.CF.reg和D100kg.2point與D100kg.reg相關(guān)系數(shù)，說(shuō)明1point校正的D100kg的準(zhǔn)確性要差一些。

表3 達(dá)100kg日齡之間的簡(jiǎn)單相關(guān)系數(shù)

上三角為相關(guān)系數(shù)，對(duì)角線為觀察數(shù)，下三角為相關(guān)系數(shù)對(duì)應(yīng)的值

The numbers in the upper triangle, diagonal and lower triangle are the correlation coefficients, observation numbers andvalues, respectively

2.5 校正日齡無(wú)偏性評(píng)估

在里岔黑豬、巴里1代、巴里2代和混合群體，reg校正誤差的均值±標(biāo)準(zhǔn)差（值）分別為-0.522± 3.6655（=0.0613）、0.107±2.1040（=0.2177）、-0.033± 1.8764（=0.7219）和-0.035±2.3402（=0.6044），均與0沒(méi)有顯著差異。2point校正誤差的均值±標(biāo)準(zhǔn)差（值）分別為-0.208±4.2646（=0.5199）、-0.629± 3.2440（=2.94×10-6）、-1.097±2.9674（=3.61×10-13）和-0.730±3.3360（=1.06×10-13）。除里岔黑豬群體外，其他3個(gè)群體均與0存在極顯著差異。1point校正誤差的均值±標(biāo)準(zhǔn)差（值）分別為19.626±26.6659（＜2.2×10-16）、-5.175±15.1430（=6.0×10-16）、-7.810±18.4939（＜2.2×10-16）和-2.422±20.6247（=5.80×10-5），均與0有極顯著差異。1point二次校正誤差的均值±標(biāo)準(zhǔn)差（值）分別為-1.472± 5.1416（=0.0002）、-1.745±4.0992（＜2.2×10-16）、-2.421±3.8033（＜2.2×10-16）和-0.043±7.3030（= 0.8384）。除混合群體外，其他三個(gè)群體均與0存在極顯著差異。

2.6 校正系數(shù)的適用性

用其他群體的校正系數(shù)和用本群體的校正系數(shù)估計(jì)的D100kg與該群體D100kg.reg間的相關(guān)系數(shù)非常接近（表4），這說(shuō)明三種校正系數(shù)均有很好的適用性。

2.7 終測(cè)體重對(duì)D100kg的影響

校正誤差與終測(cè)體重之間的散點(diǎn)圖見(jiàn)圖3。在這3個(gè)群體中，終測(cè)體重離目標(biāo)體重（100 kg）越遠(yuǎn)，D100kg.CF.reg校正誤差越大，D100kg.2point校正誤差沒(méi)有明顯的變化，D100kg.1point校正誤差呈負(fù)線性變化。終測(cè)體重越小，D100kg.1point校正日齡就比回歸日齡大的越多；反之，校正日齡就比回歸日齡小的越多。

在這3個(gè)群體中，D100kg.1point的誤差與終測(cè)體重均存在極顯著的負(fù)相關(guān)，分別為-0.9812、-0.9627和-0.9786，值均＜2.2×10-16。擬合一元回歸模型的截距分別為299.1974、135.6703和161.0781，斜率分別為-3.0067、-1.3742和-1.6350，擬合度分別為0.9628、0.9267和0.9577。

表4 校正系數(shù)通用性結(jié)果

所有相關(guān)系數(shù)均達(dá)到了極顯著水平（＜2.2×10-16） All correlation coefficients are highly significant (＜2.2×10-16)

2.8 D100kg.1point的二次校正

經(jīng)過(guò)二次校正后，單點(diǎn)校正準(zhǔn)確性得到非常明顯的改善（圖2，表2和3），誤差與終測(cè)體重之間的線性關(guān)系也消失了（圖3）。但是，誤差與終測(cè)體重仍然一定曲線關(guān)系，尤其在巴里1世代和巴里2世代：隨著終測(cè)體重與目標(biāo)體重差異的絕對(duì)值越來(lái)越大，校正誤差的值就負(fù)得越多，即單點(diǎn)二次校正日齡比回歸日齡偏小越來(lái)越多。

第1-4列分別為里岔黑豬、巴里1世代、巴里2世代和混合群體；第1 - 4行分別為回歸、兩點(diǎn)、單點(diǎn)和單點(diǎn)二次校正日齡

3 討論

本研究以里岔黑豬及其與巴克夏雜交（巴里）群體為研究對(duì)象，采用3種方法估計(jì)豬只達(dá)100 kg體重日齡的校正系數(shù)，并以回歸日齡為標(biāo)準(zhǔn)，評(píng)估其他3種方法校正的準(zhǔn)確性。除單點(diǎn)（終測(cè)點(diǎn)）估計(jì)的校正效果不理想外，其余兩種方法的準(zhǔn)確性都很高，均可以用于生產(chǎn)實(shí)踐。本研究結(jié)果為種豬場(chǎng)選擇達(dá)100 kg體重日齡校正系數(shù)的估計(jì)方法提供參考。

3.1 選擇校正系數(shù)CF的理由

根據(jù)前人的報(bào)道，雖然校正系數(shù)和的校正效果沒(méi)有顯著差異，但是校正的D100kg與D100kg.reg的相關(guān)系數(shù)（0.9918）比校正的D100kg與D100kg.reg的相關(guān)系數(shù)（0.9917）略高，且它與D100kg.reg差異的絕對(duì)值（-0.1385d）比校正的D100kg與D100kg.reg差異的絕對(duì)值（-0.2741d）小[16]，所以本研究選用校正效果相對(duì)較好的。

目前，我國(guó)對(duì)于豬只達(dá)100 kg體重日齡的校正系數(shù)方面研究較少，尤其是對(duì)符合我國(guó)地方豬種生長(zhǎng)規(guī)律的研究更少。在我國(guó)《全國(guó)種豬遺傳評(píng)估方案（試行）》中推薦使用來(lái)自于加拿大種豬遺傳改良中心的校正公式，但是加拿大瘦肉型豬目標(biāo)體重校正公式只分性別，不分群體，且我國(guó)地方豬種與西方豬種在生長(zhǎng)速度上也存在一定的差異，因此導(dǎo)致對(duì)于我國(guó)地方豬種D100kg的估計(jì)不夠準(zhǔn)確，難以提供參考[23]。在之前的研究中，鄭浩等人通過(guò)研究發(fā)現(xiàn)里岔黑豬的小于加拿大種豬遺傳改良中心杜洛克、長(zhǎng)白和大白D100kg的[16]。筆者認(rèn)為借鑒這些豬種的校正系數(shù)，會(huì)降低對(duì)于D100kg估計(jì)的準(zhǔn)確性，從而影響生長(zhǎng)速度的選擇效果。

3.2 數(shù)據(jù)量對(duì)CF的影響

從估計(jì)校正系數(shù)使用的數(shù)據(jù)量來(lái)說(shuō)，reg使用的數(shù)據(jù)量最大。在計(jì)算reg時(shí)，使用了豬只肥育期每天的體重，受異常體重的影響小。此外，肥育期豬只體重變化還可以監(jiān)測(cè)個(gè)體的健康狀況，有助于數(shù)據(jù)的質(zhì)控。因此，筆者認(rèn)為這3種方法中reg最好，D100kg.CF.reg與D100kg.reg的相關(guān)系數(shù)最大也印證了這一點(diǎn)。2point只利用了始測(cè)和終測(cè)的體重，容易受始測(cè)和終測(cè)體重的影響。在始測(cè)和終測(cè)的體重測(cè)量比較準(zhǔn)確的情況下，2point和reg沒(méi)有本質(zhì)差異，D100kg.2point與D100kg.CF.reg之間的相關(guān)系數(shù)高達(dá)0.9796。1point只利用了終測(cè)數(shù)據(jù)，使用數(shù)據(jù)量最少，校正效果也是最差的，D100kg.1point與D100kg.reg的相關(guān)系數(shù)僅為0.8873。但經(jīng)過(guò)二次校正后，與D100kg.reg的相關(guān)系數(shù)為0.9349。

3.3 終測(cè)體重對(duì)CF的影響

根據(jù)D100kg的校正公式，體重離差（終測(cè)體重與目標(biāo)體重的離差）越小，D100kg的誤差也越小。由于真實(shí)的D100kg不可獲取，因此把D100kg.reg看作為真實(shí)的D100kg，3種方法估計(jì)的D100kg與它的離差看作為它們的誤差。從圖3可以看出D100kg.2point的誤差與體重離差之間沒(méi)有明顯的關(guān)系。通過(guò)分析體重離差小、而誤差大的數(shù)據(jù)后，發(fā)現(xiàn)造成誤差大是由于個(gè)體在終測(cè)時(shí)終測(cè)體重和回歸體重的差異較大造成的。由于終測(cè)體重和目標(biāo)體重很接近，因此D100kg.2point和終測(cè)日齡也很接近。但是回歸體重與終測(cè)體重差異較大，因此D100kg.reg與終測(cè)日齡差異較大，即與D100kg.2point的差異較大，從而導(dǎo)致D100kg.2point誤差也較大。

3.4 D100kg的無(wú)偏性

采用線性回歸的方法來(lái)估計(jì)D100kg的假設(shè)是體重和日齡在測(cè)定期間是線性的。從體重對(duì)日齡的簡(jiǎn)單回歸的擬合度非常接近1（質(zhì)控后2≥0.95）來(lái)看[16]，這個(gè)假設(shè)即使不完全滿足，也應(yīng)該近似滿足。在滿足這個(gè)假設(shè)的前提下，D100k.reg是無(wú)偏的。D100kg.CF.reg與D100k.reg之間沒(méi)有顯著差異，因此D100kg.CF.reg也是無(wú)偏的。D100kg.2point與D100k.reg之間存在極顯著差異（里岔群體除外），但這種差異很可能是起始重和終測(cè)重的測(cè)定誤差造成的（見(jiàn)3.3）。從圖3可以看出D100kg.2point的誤差關(guān)于=0對(duì)稱，則說(shuō)明D100kg.2point是無(wú)偏的。D100kg.1point是有偏的：當(dāng)終測(cè)體重大于100 kg時(shí)偏小；當(dāng)終測(cè)體重小于100 kg時(shí)偏大。即使經(jīng)過(guò)二次校正，單點(diǎn)校正的D100kg仍然是有偏的：總體上來(lái)看是偏小，尤其是終測(cè)體重與目標(biāo)體重差異的絕對(duì)值較大時(shí)。

3.5 CF適用范圍討論

如果有自動(dòng)飼喂系統(tǒng)，能夠自動(dòng)記錄肥育期豬只的體重，建議優(yōu)先選用reg。如果沒(méi)有自動(dòng)飼喂系統(tǒng)，很難準(zhǔn)確獲取肥育期豬只所有的體重，建議測(cè)定始測(cè)和終測(cè)體重，用2point來(lái)計(jì)算D100kg。在任何情況下均不建議使用1point，因?yàn)槭褂闷渌后w的reg和2point校正效果均優(yōu)于本群體的1point校正效果（表4）。雖然經(jīng)過(guò)二次校正后單點(diǎn)校正效果顯著提升，但是進(jìn)行二次校正需要校正誤差，即校正日齡減回歸日齡。如果能夠得到回歸日齡，就能計(jì)算出reg，因此可以選擇使用reg進(jìn)行校正，因?yàn)閞eg的校正效果最佳。

根據(jù)豬只可用數(shù)據(jù)的多少，采用最佳的方法來(lái)校正D100kg可以提高校正準(zhǔn)確性，降低估計(jì)誤差。對(duì)于有肥育期體重?cái)?shù)據(jù)的個(gè)體，計(jì)算D100kg.reg比較合適。對(duì)于只有始測(cè)和終測(cè)數(shù)據(jù)的豬只，如果reg可以估計(jì)，則先估計(jì)reg，然后計(jì)算D100kg.CF.reg；否則先估計(jì)2point，然后獲得D100kg.2point。對(duì)于只有終測(cè)數(shù)據(jù)的豬只，如果reg或2point可以估計(jì)，則先估計(jì)reg或2point，然后計(jì)算D100kg.CF.reg或D100kg.2point；否則借鑒其他群體的reg或2point，然后計(jì)算D100kg.CF.reg或D100kg.2point。

4 結(jié)論

回歸校正系數(shù)（reg）和兩點(diǎn)校正系數(shù)（2point）的校正效果都很好，均可用于生產(chǎn)實(shí)踐。單點(diǎn)校正系數(shù)（1point）的校正效果不理想，不適合用于生產(chǎn)實(shí)踐。在估計(jì)D100kg時(shí)，可以根據(jù)豬只的數(shù)據(jù)量，采用合適的方法，可以提高估計(jì)準(zhǔn)確性，降低估計(jì)誤差。

致謝：本研究得到江西農(nóng)業(yè)大學(xué)省部共建種豬遺傳改良與養(yǎng)殖技術(shù)國(guó)家重點(diǎn)實(shí)驗(yàn)室主任黃路生教授的大力支持。

[1] 熊信威. 豬12號(hào)染色體上增加肌內(nèi)脂肪含量的因果基因及其因果突變的鑒別與作用機(jī)理研究[D]. 南昌: 江西農(nóng)業(yè)大學(xué), 2017.

XIONG X W. Identification of the causative gene and mutation for increasingcontent at SSC12 and its mechanism study[D]. Nanchang: Jiangxi Agricultural University, 2017. (in Chinese)

[2] 張?zhí)m勇. 養(yǎng)豬業(yè)怎樣提高經(jīng)濟(jì)效益. 國(guó)外畜牧學(xué)(豬與禽), 2009, 29(3): 93-94. doi:10.3969/j.issn.1001-0769.2009.03.040.

ZHANG L Y. How to improve the economic efficiency of swine industry. Pigs and Poultry, 2009, 29(3): 93-94. doi:10.3969/j.issn. 1001-0769.2009.03.040. (in Chinese)

[3] 博亞. 規(guī)?；B(yǎng)豬業(yè)發(fā)展探析. 農(nóng)業(yè)技術(shù)與裝備, 2011(19): 4-6.

BO Y. Study on the development of swine industry in large-scale. Agricultural Technology & Equipment, 2011(19): 4-6. (in Chinese)

[4] 陳日. 2021. 生態(tài)豬養(yǎng)殖技術(shù)及發(fā)展趨勢(shì)[J]. 農(nóng)家致富顧問(wèn), (12): 208.

CHEN R. The Technology of Raising Pigs Ecologically and Its Developing Trend., 2021 (12): 208. (in Chinese)

[5] 周道雷, 李保明, 施正香. 2003. 我國(guó)養(yǎng)豬產(chǎn)業(yè)生產(chǎn)現(xiàn)狀及其發(fā)展方向. 農(nóng)業(yè)工程學(xué)報(bào), 19(z1): 101-105.

ZHOU D L, LI B M, SHI Z X. Current Status and Development Strategy of Swine Industry in China. Transactions of the Chinese Societyof Agricultural Engineering, 2003, 19(z1): 101-105. (in Chinese)

[6] 孫泉云, 衛(wèi)金良, 曹建國(guó). 中國(guó)土種黑豬的開(kāi)發(fā)和利用現(xiàn)狀. 上海畜牧獸醫(yī)通訊, 2013(3): 45. doi:10.3969/j.issn.1000-7725.2013.03. 026.

SUN Q Y, WEI J L, CAO J G. Present Status of Exploitation and Utilization of Indigenous Black Pig Breeds in China. Shanghai Journal of Animal Husbandry and Veterinary Medicine, 2013(3): 45. doi:10.3969/j.issn.1000-7725.2013.03.026. (in Chinese)

[7] 陳鵬, 楊在賓, 張慶, 姜淑貞. 里岔黑豬與 “杜長(zhǎng)大” 三元豬生長(zhǎng)性能和免疫性能的對(duì)比研究. 豬業(yè)科學(xué), 2017, 34(1): 131-133. doi:10.3969/j.issn.1673-5358.2017.01.051.

CHEN P, YANG Z B, ZHANG Q, JIANG S Z. Comparative the growth and immunity performances of Licha black pig to Duroc× Landrace×Yorkshire pigs. Swine Industry Science, 2017, 34(1): 131-133. doi:10.3969/j.issn.1673-5358.2017.01.051. (in Chinese)

[8] 曾檢華, 羅艷鳳, 胡杰, 陳瑤生, 劉小紅. 不同黑豬雜交組合生產(chǎn)性能研究與分析. 養(yǎng)豬, 2015(1): 61-64. doi:10.13257/j.cnki.21- 1104/s.2015.01.024.

ZENG J H, LUO Y F, HU J, CHEN Y S, LIU X H. The research and analysis of production performance in different hybrid black pigs. Swine Production, 2015(1): 61-64. doi:10.13257/j.cnki.21-1104/s. 2015.01.024. (in Chinese)

[9] 李文通, 劉穎, 周榮, 沈自泉, 董勝華, 宋孝濤, 王月江, 唐輝. 江泉黑豬與杜長(zhǎng)大三元雜交豬屠宰性能及肉質(zhì)性狀比較分析. 中國(guó)畜牧獸醫(yī), 2017, 44(7): 2057-2064. doi:10.16431/j.cnki.1671-7236. 2017.07.022.

LI W T, LIU Y, ZHOU R, SHEN Z Q, DONG S H, SONG X T, WANG Y J, TANG H. A comparative study on the carcass and meat quality traits of Jiangquan black pig and DLY pig. China Animal Husbandry & Veterinary Medicine, 2017, 44(7): 2057-2064. doi:10. 16431/j.cnki.1671-7236.2017.07.022. (in Chinese)

[10] 郭源梅, 李龍?jiān)? 賴昭勝, 黃黎斌, 黃智勇, 涂金敏, 邱恒清, 黃路生, 肖石軍. 中國(guó)地方豬種利用現(xiàn)狀與展望. 江西農(nóng)業(yè)大學(xué)學(xué)報(bào), 2017, 39(3): 427-435. doi:10.13836/j.jjau.2017055.

GUO Y M, LI L Y, LAI Z S, HUANG L B, HUANG Z Y, TU J M, QIU H Q, HUANG L S, XIAO S J. Present status and prospect of utilization of indigenous pig breeds in China. Acta Agriculturae Universitatis Jiangxiensis, 2017, 39(3): 427-435. doi:10.13836/j. jjau.2017055. (in Chinese)

[11] 劉國(guó)信. 加強(qiáng)地方豬種保護(hù)利用提高生豬產(chǎn)業(yè)競(jìng)爭(zhēng)優(yōu)勢(shì). 浙江畜牧獸醫(yī), 2021, 46(1): 14-15.

LIU G X. Strengthening the protection and utilization of native pig breeds and improving the competitive advantage of pig industry. Zhejiang Journal Animal Science and Veterinary Medicine, 2021, 46(1): 14-15. (in Chinese)

[12] 于福清. 我國(guó)地方豬種資源最新?tīng)顩r與保護(hù)利用建議. 中國(guó)豬業(yè), 2012, 7(3): 21-23. doi:10.16174/j.cnki.115435.2012.03.012.

YU F Q. The latest situation of Chinese native pig breeds and suggestions on protection and utilization of them. China Swine Industry, 2012, 7(3): 21-23. doi:10.16174/j.cnki.115435.2012.03.012. (in Chinese)

[13] 劉艷花, 郭建華, 李書(shū)勛, 紀(jì)志和. 關(guān)于加快里岔黑豬新品種研究開(kāi)發(fā)的思考. 山東畜牧獸醫(yī), 2017, 38(6): 20-21.

LIU Y H, GUO J H, LI S X, JI Z H. Thoughts on accelerating the research and development of new breed using Licha black pig. Shandong Journal of Animal Science and Veterinary Medicine, 2017, 38(6): 20-21. (in Chinese)

[14] 彭瀟, 尹立林, 梅全順, 王海燕, 劉小磊, 朱猛進(jìn), 李新云, 付亮亮, 趙書(shū)紅. 豬主要經(jīng)濟(jì)性狀的基因組選擇研究. 畜牧獸醫(yī)學(xué)報(bào), 2019, 50(2): 439-445. doi:10.11843/j.issn.0366-6964.2019.02.023.

PENG X, YIN L L, MEI Q S, WANG H Y, LIU X L, ZHU M J, LI X Y, FU L L, ZHAO S H. A study of genome selection based on the porcine major economic traits. Acta Veterinaria et Zootechnica Sinica, 2019, 50(2): 439-445. doi:10.11843/j.issn.0366-6964.2019.02.023. (in Chinese)

[15] WANG K J, LIU D W, HERNANDEZ-SANCHEZ J, CHEN J, LIU C K, WU Z F, FANG M Y, LI N. Genome wide association analysis reveals new production trait genes in a male duroc population. PLoS ONE, 2015, 10(9): e0139207. doi:10.1371/journal.pone.0139207.

[16] 鄭浩, 涂金敏, 熊秀萍, 張則凱, 李龍?jiān)? 黃黎斌, 黃智勇, 肖石軍, 郭源梅. 里岔黑豬和巴里雜交豬達(dá)100 kg體重日齡校正系數(shù)研究. 中國(guó)農(nóng)業(yè)科學(xué), 2020, 53(12): 2493-2501. doi:10.3864/j.issn.0578- 1752.2020.12.015.

ZHENG H, TU J M, XIONG X P, ZHANG Z K, LI L Y, HUANG L B, HUANG Z Y, XIAO S J, GUO Y M. Estimating the correction coefficient of days to 100 kg inblack pig and its intercross with Berkshire. Scientia Agricultura Sinica, 2020, 53(12): 2493-2501. doi:10.3864/j.issn.0578-1752.2020.12.015. (in Chinese)

[17] 全國(guó)畜牧獸醫(yī)總站. 全國(guó)種豬遺傳評(píng)估方案(試行) [Z]. 2000: 60.

National Animal Husbandry & Veterinary Service.(). Beijing, 2000: 60. (in Chinese)

[18] 822-2004" N T. 中華人民共和國(guó)農(nóng)業(yè)行業(yè)標(biāo)準(zhǔn)-種豬生產(chǎn)性能測(cè)定規(guī)程. 北京: 中華人民共和國(guó)農(nóng)業(yè)部, 2004.

NY/T 822-2004. People's Republic of China agricultural industry standards – Rules for performance testing of breeding pig. Beijing: Ministry of Agriculture of the People's Republic of China, 2004. (in Chinese)

[19] 全國(guó)生豬遺傳改良計(jì)劃工作小組領(lǐng)導(dǎo)辦公室. 關(guān)于調(diào)整種豬達(dá)100Kg體重日齡和達(dá)100Kg體重活體背膘厚校正公式的通知. 2016: 1-2.

LEADING GROUP OF WORK OFFICE OF CHINA SWINE GENETIC IMPROVEMENT PROGRAM. Notice on adjusting the calibration formula of days to 100 kg and backfat thickness at 100 kg in vivo for breeding pigs. 2016: 1-2. (in Chinese)

[20] 李加琪. 中國(guó)瘦肉型豬性能測(cè)定校正公式. 中國(guó)豬業(yè)科技大會(huì)暨中國(guó)畜牧獸醫(yī)學(xué)會(huì)2015年學(xué)術(shù)年會(huì)論文集. 廈門(mén). 2015: 5-5.

LI J Q. Correction Formula for Production Performance Testing of Chinese Lean Pigs. Proceedings of China Pig Science and Technology Conference and 2015 Academic Annual Meeting of China Animal Husbandry and Veterinary Society, Xiamen. 2015: 5-5. (in Chinese)

[21] 劉望宏, 張勤, 胡軍勇, 倪德斌, 熊遠(yuǎn)著. 瘦肉型種豬早期生長(zhǎng)性狀校正公式的制定及評(píng)估. 華中農(nóng)業(yè)大學(xué)學(xué)報(bào), 2010, 29(4): 469-474. doi:10.13300/j.cnki.hnlkxb.2010.04.024.

LIU W H, ZHANG Q, HU J Y, NI D B, XIONG Y Z. Design and evaluation of adjusted equations for early growth performance of commercial breeding swines. Journal of Huazhong Agricultural University, 2010, 29(4): 469-474. doi:10.13300/j.cnki.hnlkxb.2010. 04.024. (in Chinese)

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

CHEN B. Studies on the on-farm genetic evaluation and genetic analysis of lean type swine[D]. Changsha: Hunan Agricultural University, 2005. (in Chinese)

[23] 張哲, 張豪, 陳贊謀, 李加琪. 種豬育種性能測(cè)定校正公式研究. 中國(guó)畜牧雜志, 2015, 51(16): 49-54. doi:10.3969/j.issn.0258-7033. 2015.16.011.

ZHANG Z, ZHANG H, CHEN Z M, LI J Q. Development of correction formula for production traits in swine breeding. Chinese Journal of Animal Science, 2015, 51(16): 49-54. doi:10.3969/j.issn. 0258-7033.2015.16.011. (in Chinese)

[24] 王林云. 豬育種工作中某些數(shù)據(jù)的校正方法. 畜牧與獸醫(yī), 2000, 32(1): 20-21. doi:10.3969/j.issn.0529-5130.2000.01.010.

WANG L Y. Correction of some indexes in pig breeding. Animal Husbandry & Veterinary Medicine, 2000, 32(1): 20-21. doi:10.3969/ j.issn.0529-5130.2000.01.010. (in Chinese)

[25] 呂鳳林, 黎繼華. 內(nèi)江仔豬體重校正方法的研究. 四川畜牧獸醫(yī), 1993, 20(2): 1-3.

Lü F L, LI J H. A study of the method to adjust bodyweights of Neijiang pig. Sichuan Animal & Veterinary Sciences, 1993, 20(2): 1-3. (in Chinese)

[26] 楊岸奇, 吳買(mǎi)生, 向擁軍, 陳斌, 李朝輝, 張善文, 劉偉, 吳攀峰. 沙子嶺豬生長(zhǎng)性狀校正公式的制定. 豬業(yè)科學(xué), 2015, 32(2): 132-134. doi:10.3969/j.issn.1673-5358.2015.02.048.

YANG A Q, WU M S, XIANG Y J, CHEN B, LI Z H, ZHANG S W, LIU W, WU P F. Formulating the correction formula for growth traits of Shaziling pig. Swine Industry Science, 2015, 32(2): 132-134. doi:10.3969/j.issn.1673-5358.2015.02.048. (in Chinese)

[27] 胡閃耀, 陳歡, 閆德超, 張子君, 毛元良, 毛慧敏, 王利通, 王國(guó)水, 左波. 大白、長(zhǎng)白達(dá)100 kg體重時(shí)體長(zhǎng)校正公式的研究. 豬業(yè)科學(xué), 2016, 33(1): 110-111. doi:10.3969/j.issn.1673-5358.2016.01.036.

HU S Y, CHEN H, YAN D C, ZHANG Z J, MAO Y L, MAO H M, WANG L T, WANG G S, ZUO B. Study on the correction formula of body length for body weight at 100 kg in Large White and Landrace. Swine Industry Science, 2016, 33(1): 110-111. doi:10.3969/j.issn.1673- 5358.2016.01.036. (in Chinese)

[28] 孫賢,陳衛(wèi)東, 梁鳳蓮, 左濤, 居瑞海, 魏傳德. 豐海種豬生產(chǎn)性能測(cè)定與應(yīng)用報(bào)告. 上海畜牧獸醫(yī)通訊, 2010(3): 16-17. doi:10.3969/j. issn.1000-7725.2010.03.009.

SUN X, CHEN W D, LIANG F L, ZUO T, JU R H, WEI C D. Reports on production performance testing and its application of fenghai breeding pig. Shanghai Journal of Animal Husbandry and Veterinary Medicine, 2010(3): 16-17. doi:10.3969/j.issn.1000-7725.2010.03.009. (in Chinese)

[29] WU P X, WANG K, YANG Q, ZHOU J, CHEN D J, LIU Y H, MA J D, TANG Q Z, JIN L, XIAO W H, LOU P E, JIANG A N, JIANG Y Z, ZHU L, LI M Z, LI X W, TANG G Q. Whole-genome re-sequencing association study for direct genetic effects and social genetic effects of six growth traits in Large White pigs. Scientific Reports, 2019, 9: 9667. doi:10.1038/s41598-019-45919-0.

[30] LI F G, LEI M G, ZHENG R, ZUO B, JIANG S W, DENG C Y, XIONG Y Z. The effects of estrogen receptor locus on reproductive tracts components and performance traits in large white × Meishan F2offspring. Asian-Australasian Journal of Animal Sciences, 2004, 17(9): 1223-1226. doi:10.5713/ajas.2004.1223.

[31] 李龍?jiān)? 肖石軍, 黃黎斌, 黃智勇, 涂金敏, 熊秀萍, 郭源梅, 黃路生. 巴克夏和里岔黑豬不同雜交組合生產(chǎn)性能的比較研究. 畜牧獸醫(yī)學(xué)報(bào), 2018, 49(1): 26-35. doi:10.11843/j.issn.0366-6964.2018.01. 004.

LI L Y, XIAO S J, HUANG L B, HUANG Z Y, TU J M, XIONG X P, GUO Y M, HUANG L S. Comparative studies on productive performances among different cross combinations of Berkshire andblack pig. Acta Veterinaria et Zootechnica Sinica, 2018, 49(1): 26-35. doi:10.11843/j.issn.0366-6964.2018.01.004. (in Chinese)

[32] 賀婕妤, 王斌虎, 廖柱, 謝紅濤, 易國(guó)強(qiáng), 劉毓文, 敖紅, 唐中林. 長(zhǎng)白和大白豬主要生長(zhǎng)性狀的遺傳參數(shù)估計(jì). 畜牧獸醫(yī)學(xué)報(bào), 2021, 52(8): 2115-2123. doi:10.11843/j.issn.0366-6964.2021.08.005.

HE J Y, WANG B H, LIAO Z, XIE H T, YI G Q, LIU Y W, AO H, TANG Z L. Estimation of genetic parameters of main growth traits in Landrace and large white pigs. Acta Veterinaria et Zootechnica Sinica, 2021, 52(8): 2115-2123. doi:10.11843/j.issn.0366-6964.2021.08.005. (in Chinese)

Comparing Methods for Correcting Days to 100 kg of Sows in Licha Black Pig and Its Intercross with Berkshire

CUI DengShuai1, XIONG SanYa1, ZHENG Hao2, LI LongYun1, YU NaiBiao1, HUANG ZhiYong2, XIAO ShiJun1, GUO YuanMei1

1State Key Laboratory of Pig Genetic Improvement and Production Technology, Jiangxi Agricultural University, Nanchang 330045;2Jiangxi Shanxia Huaxi Pig Breeding Company Limited, Ganzhou 341900, Jiangxi

【Objective】 The aims of this study were to compare the powers of three correction coefficients for days to 100 kg (D100kg), so as to provide the reference for pig breeding farms to select the appropriate correction coefficient for D100kgaccording to the data they have. 【Method】A total of 1 344 pigs were used in this study, including Licha black pig, F1(BL1) and F2(BL2) from its intercross with Berkshire boars. Body weight and age in days for each pig were recorded everyday during its growth and fattening period. The data were qualified according to following criteria: the start body weight ≥25kg, the final body weight among 80 – 125 kg, the slope of the regression equation of body weight to age in days within±2 standard deviations of its average and so on. Three methods described as following were used to estimate the correction coefficient of the days to 100 kg () into Licha Black pig, BL1, BL2and their mixed population. First one was the regressive correction coefficient (reg):1iwas the slope of the simple regression equation of body weight to age in days for theth individual, and2iwas the slop of the straight line through the testing end point of theth individual and the origin (age and weight were both 0), then the average of1i/2iin the population wasreg. Second one is the two-point correction coefficient (2point):1iwas the slope of the straight line through the testing start and end points of theth individual, and the methods to calculate2iand2pointwere as same as those ofreg. The last one was the single-point correction coefficient (1point):1was the slope of the simple regression equation of body weight to age in days at the testing end point for all pigs, and the methods to calculate2iand1pointwere as same as those ofreg. Putreg,2pointand1pointinto the correction equation, respectively, and got D100kg.CF.reg, D100kg.2pointand D100kg.1point. The simple regression of age in days to body weight was used to calculate the regressive age in days to 100 kg (D100kg.reg), and the correlation coefficients and deviations between D100kg.regand other D100kgwere calculated. The merits and demerits of the three methods were evaluated according to their correlation coefficients and deviations to D100kg.reg. D100kgof a population was estimated by using the correction coefficients from another population, and its correlation coefficient with D100kg.regof the population was calculated to evaluate whether the correction coefficient estimated from one population could be used in other populations. The relationship between body weight at the end of performance test (BWend) and correction error was also investigated to explore whether BWendhad an effect on D100kg. Taking consideration of the strong linear relationship between BWendand the correction error of1point, BWendwas used to correct D100kg.1pointto improve the correction accuracy of1point.【Result】 A total of 1 181 sows were used for following analyses after quality control. In the same population, the values ofregand2pointwere almost the same, and both of which were highly significantly greater thanCF. D100kg.CF.regand D100kg.2pointof each pig were very close to its D100kg.reg, and all of the correlation coefficients among D100kg.CF.reg, D100kg.2pointand D100kg.regwere greater than 0.98, indicating that the correction powers ofregand2pointwere very good. The difference between D100kg.1pointand D100kg.regwas large, and their correlation coefficient was 0.89, so the correction power was not very good.regand2pointhad good population transplantations, and could be used to correct D100kgof other populations, while2pointhad poor population transplantations, and it was unsuitable to correct D100kgof other populations. With the difference of BWendto the target body weight increasing, the correction error of D100kg.CF.regwas also increasing in the three populations, and the correction error of D100kg.2pointwas no significant changing. D100kg.1pointwas significant negative correlation with BWend, the correlation coefficients were -0.9812, -0.9627, -0.9786 and -0.9352 in the four populations, respectively. The correction power of1pointhas been significantly improved after adjusting with BWend, butregwas still better than1point, and it was preferred.【Conclusion】The correction powers and transplantations of correction coefficientsregand2pointwere very good, soregand2pointwere suitable to be used to correct Din practice. The correction power of1pointwas not very good, therefore,1pointwas unsuitable to be used in practice. According to the datum amount of a pig, its Dwas estimated with the appropriate correction coefficient in order to improve the estimation accuracy and reduce the estimation error.

Licha Black pig; days to 100 kg; correction coefficient; regressive days to 100 kg

10.3864/j.issn.0578-1752.2023.06.013

2021-11-29；

2022-05-17

國(guó)家自然科學(xué)基金（31972542和31660303）

崔登帥，E-mail：1305568175@qq.com。通信作者郭源梅，Tel/Fax：0791-83813080；E-mail：gyuanmei@hotmail.com

（責(zé)任編輯 林鑒非）