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

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

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

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

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

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

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

1 材料與方法

1.1 試驗(yàn)豬只

研究所用的試驗(yàn)豬只均來自于江西山下華系種豬養(yǎng)殖有限公司定南陽林山下原種豬場（江西，定南）[31]，實(shí)驗(yàn)數(shù)據(jù)采集于2014—2019年。飼養(yǎng)管理和營養(yǎng)水平參見文獻(xiàn)[16]，豬只自由采食和飲水。通過潤農(nóng)性能測定系統(tǒng)（深圳，潤農(nóng)科技），測定試驗(yàn)豬只生長肥育期的體重，并記錄測定日期。本研究利用前期報(bào)道的1 344頭豬只生長肥育期的體重和日期[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)通過了以下質(zhì)控標(biāo)準(zhǔn)：測定期≥30 d；始測體重≥25 kg，終測體重在80—125 kg；始測日齡>60 d和終測日齡<280 d；優(yōu)化后體重對日齡一元回歸方程擬合度≥0.95的個(gè)體[16]。本研究新增一個(gè)質(zhì)控標(biāo)準(zhǔn)，即經(jīng)過始測點(diǎn)和終測點(diǎn)直線的斜率在群體平均斜率2個(gè)標(biāo)準(zhǔn)差之內(nèi)。

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

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

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

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

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

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

1.2.4 校正日齡之間的相關(guān)系數(shù) 利用R語言的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ù)越高，說明該方法估計(jì)的校正日齡越準(zhǔn)確；反之，越不準(zhǔn)確。

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

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

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

1.2.8 D100kg.1point的二次校正 由于單點(diǎn)校正誤差與終測體重（BWend）之間存在很強(qiáng)的線性關(guān)系，且單點(diǎn)校正效果并不理想，因此建立單點(diǎn)校正誤差對終測體重的回歸方程，把體重和回歸方程的斜率（）代入下列公式，對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è)體通過了質(zhì)控，其中里岔黑豬197頭（母豬175頭，閹公豬22頭），巴里1代621頭（母豬593頭，閹公豬28頭），巴里2代472頭（母豬413頭，閹公豬59頭）。由于性別和群體對校正系數(shù)影響顯著，且3個(gè)群體中閹公豬數(shù)量都較少，影響閹公豬校正系數(shù)估計(jì)的準(zhǔn)確性，因此，僅1 181頭母豬用于后續(xù)分析。通過質(zhì)控個(gè)體始測點(diǎn)和終測點(diǎn)連線的趨勢基本相似（圖1），即豬只的生長規(guī)律相似，說明質(zhì)控后的數(shù)據(jù)質(zhì)量較好。

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

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

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

回歸日齡（D100kg.reg）和三種方法估計(jì)的D100kg見表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è)體始測點(diǎn)和終測點(diǎn)的連線圖

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

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

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

4種方法估計(jì)的D100kg散點(diǎn)圖見圖2，它們之間的簡單相關(guān)系數(shù)見表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ù)沒有后三者之間那么強(qiáng)，0.8873≤≤ 0.8999，＜2.2×10-16。

D100kg與D100kg.reg相關(guān)系數(shù)可以用來度量其估計(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ù)，說明1point校正的D100kg的準(zhǔn)確性要差一些。

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

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

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

2.5 校正日齡無偏性評估

在里岔黑豬、巴里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沒有顯著差異。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），這說明三種校正系數(shù)均有很好的適用性。

2.7 終測體重對D100kg的影響

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

在這3個(gè)群體中，D100kg.1point的誤差與終測體重均存在極顯著的負(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)過二次校正后，單點(diǎn)校正準(zhǔn)確性得到非常明顯的改善（圖2，表2和3），誤差與終測體重之間的線性關(guān)系也消失了（圖3）。但是，誤差與終測體重仍然一定曲線關(guān)系，尤其在巴里1世代和巴里2世代：隨著終測體重與目標(biāo)體重差異的絕對值越來越大，校正誤差的值就負(fù)得越多，即單點(diǎn)二次校正日齡比回歸日齡偏小越來越多。

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

3 討論

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

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

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

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

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

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

3.3 終測體重對CF的影響

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

3.4 D100kg的無偏性

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

3.5 CF適用范圍討論

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

根據(jù)豬只可用數(shù)據(jù)的多少，采用最佳的方法來校正D100kg可以提高校正準(zhǔn)確性，降低估計(jì)誤差。對于有肥育期體重?cái)?shù)據(jù)的個(gè)體，計(jì)算D100kg.reg比較合適。對于只有始測和終測數(shù)據(jù)的豬只，如果reg可以估計(jì)，則先估計(jì)reg，然后計(jì)算D100kg.CF.reg；否則先估計(jì)2point，然后獲得D100kg.2point。對于只有終測數(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ù)國家重點(diǎn)實(shí)驗(yàn)室主任黃路生教授的大力支持。

[1] 熊信威. 豬12號染色體上增加肌內(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ì)效益. 國外畜牧學(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ī)模化養(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ā)展趨勢[J]. 農(nóng)家致富顧問, (12): 208.

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

[5] 周道雷, 李保明, 施正香. 2003. 我國養(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)金良, 曹建國. 中國土種黑豬的開發(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] 陳鵬, 楊在賓, 張慶, 姜淑貞. 里岔黑豬與 “杜長大” 三元豬生長性能和免疫性能的對比研究. 豬業(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ì)性狀比較分析. 中國畜牧獸醫(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)? 賴昭勝, 黃黎斌, 黃智勇, 涂金敏, 邱恒清, 黃路生, 肖石軍. 中國地方豬種利用現(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] 劉國信. 加強(qiáng)地方豬種保護(hù)利用提高生豬產(chǎn)業(yè)競爭優(yōu)勢. 浙江畜牧獸醫(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] 于福清. 我國地方豬種資源最新狀況與保護(hù)利用建議. 中國豬業(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] 劉艷花, 郭建華, 李書勛, 紀(jì)志和. 關(guān)于加快里岔黑豬新品種研究開發(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), 李新云, 付亮亮, 趙書紅. 豬主要經(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ù)研究. 中國農(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] 全國畜牧獸醫(yī)總站. 全國種豬遺傳評估方案(試行) [Z]. 2000: 60.

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

[18] 822-2004" N T. 中華人民共和國農(nóng)業(yè)行業(yè)標(biāo)準(zhǔn)-種豬生產(chǎn)性能測定規(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] 全國生豬遺傳改良計(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] 李加琪. 中國瘦肉型豬性能測定校正公式. 中國豬業(yè)科技大會暨中國畜牧獸醫(yī)學(xué)會2015年學(xué)術(shù)年會論文集. 廈門. 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)著. 瘦肉型種豬早期生長性狀校正公式的制定及評估. 華中農(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] 陳斌. 瘦肉型豬的場內(nèi)遺傳評估及遺傳分析研究[D]. 長沙: 湖南農(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] 張哲, 張豪, 陳贊謀, 李加琪. 種豬育種性能測定校正公式研究. 中國畜牧雜志, 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] 楊岸奇, 吳買生, 向擁軍, 陳斌, 李朝輝, 張善文, 劉偉, 吳攀峰. 沙子嶺豬生長性狀校正公式的制定. 豬業(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] 胡閃耀, 陳歡, 閆德超, 張子君, 毛元良, 毛慧敏, 王利通, 王國水, 左波. 大白、長白達(dá)100 kg體重時(shí)體長校正公式的研究. 豬業(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)性能測定與應(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] 賀婕妤, 王斌虎, 廖柱, 謝紅濤, 易國強(qiá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

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

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

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