李佳瑩，梁 龍，賀 鑫，巴音花，何鵬飛，張 寧，賀三剛*

（1.新疆農(nóng)業(yè)大學(xué)動(dòng)物科學(xué)學(xué)院，烏魯木齊 830052；2.新疆畜牧科學(xué)院生物技術(shù)研究所/農(nóng)業(yè)農(nóng)村部草食家畜遺傳育種與繁殖重點(diǎn)實(shí)驗(yàn)室/新疆動(dòng)物生物技術(shù)重點(diǎn)實(shí)驗(yàn)室，烏魯木齊 830011）

在綿羊的選種選育中，為了便于生產(chǎn)管理，減少人力資源的經(jīng)濟(jì)成本，管理者通常采取統(tǒng)一管理的方式進(jìn)行統(tǒng)一配種、統(tǒng)一注射疫苗、統(tǒng)一飼喂、統(tǒng)一培養(yǎng)、統(tǒng)一接羔，統(tǒng)一斷奶，統(tǒng)一測量體重和體尺體況。這些指標(biāo)檢測都要在同一天進(jìn)行，但是我們不能保證母羊同時(shí)發(fā)情，也不能保證羔羊統(tǒng)一出生，更不能計(jì)算羔羊的斷奶天數(shù)讓它們剛好滿100 d進(jìn)行斷奶，蘇保全[1]指出：斷乳羔羊育肥是羊肉生產(chǎn)的主要形式。斷奶重作為生產(chǎn)發(fā)育和評定羔羊重要的指標(biāo)，同一天斷奶的數(shù)據(jù)并不能作為評價(jià)羔羊斷奶重的指標(biāo)。根據(jù)出生類型和性別進(jìn)行劃分羔羊類型進(jìn)行斷奶重的校正，使斷奶體重校正到同一水平，可以為接下來的研究提供理論基礎(chǔ)，使羔羊的選種選育更加科學(xué)。

綿羊不僅是人類需要的主要肉類食品和輔助畜產(chǎn)品的來源，也是人類農(nóng)業(yè)文明傳播的主要組成部分[2]。斷奶重是反映羔羊早期生長發(fā)育健康狀況、羔羊生產(chǎn)性能的重要性狀，受到遺傳、環(huán)境和母體效應(yīng)等多種因素的影響[3]。羔羊的斷奶鑒定是綿羊育種工作中早期選種的重要環(huán)節(jié)，尤其是在種公羊的選種時(shí)，科研人員和農(nóng)場養(yǎng)殖戶會(huì)根據(jù)斷奶時(shí)期羔羊的體重作為評判羔羊好壞的唯一指標(biāo)。斷奶體重作為一項(xiàng)重要的中間經(jīng)濟(jì)性狀，是影響家畜生產(chǎn)能力的決定指標(biāo)，倍受畜牧工作者的重視[4]。此時(shí)便突出了羔羊斷奶體重的重要性，成年母羊配種周期長，通常為45 d左右，在日后的產(chǎn)羔過程中，難免會(huì)造成羔羊出生日期的差異，而牧場和農(nóng)戶為了節(jié)約時(shí)間、經(jīng)濟(jì)、人力的成本通常會(huì)選擇統(tǒng)一斷奶，彼時(shí)測定的斷奶體重嚴(yán)重缺乏說服力，不能作為評判羔羊優(yōu)劣的標(biāo)準(zhǔn)。將不同日齡斷奶羔羊體重校正到標(biāo)準(zhǔn)化的同一水平上，可將優(yōu)秀的羔羊集中在優(yōu)秀的羊群中飼喂培養(yǎng)，為社會(huì)挑選出優(yōu)良的羊種。

1 材料與方法

1.1 試驗(yàn)材料

1.1.1 數(shù)據(jù)來源

采用新疆維吾爾自治區(qū)伊犁哈薩克自治州昭蘇馬場薩?？朔N羊場2021年7月中旬?dāng)嗄惕b定數(shù)據(jù)共1 030條，羔羊品種為特哈雜交F2代羔羊，即以特克塞爾羊?yàn)楦副?，以特克塞爾×哈薩克F1代羊?yàn)槟副?，雜交繁育獲得級進(jìn)F2代羊。

1.1.2 試驗(yàn)數(shù)據(jù)

采集的數(shù)據(jù)包括羔羊個(gè)體號、性別、出生日期（5月出生羔羊數(shù)量較少，出生日期太晚，不在人工配種預(yù)期的產(chǎn)羔周期內(nèi)，故將其剔除）、母羊年齡、單雙羔（由于3胎數(shù)量極少<10只，故將3胎數(shù)據(jù)劃歸到雙胎數(shù)據(jù)之中）、初生重、斷奶重、死亡記錄等。將數(shù)據(jù)錄入Excel表格中進(jìn)行整理，為保證數(shù)據(jù)的準(zhǔn)確性，剔除死亡、缺失數(shù)據(jù)。羔羊斷奶重?cái)?shù)據(jù)記為平均數(shù)加減3倍標(biāo)準(zhǔn)差來確定，篩選刪除異常值，其中包括羔羊死胎數(shù)據(jù)101條；出生于2月份羔羊數(shù)據(jù)3條；出生于5月份羔羊數(shù)據(jù)36條；不明羔羊性別的數(shù)據(jù)1條；不明羔羊母親年齡數(shù)據(jù)17條；共刪除異常值158條，最終有效數(shù)據(jù)872條。

1.2 試驗(yàn)方法

1.2.1 水平劃分

將試驗(yàn)類別按照出生類型和性別進(jìn)行劃分為四個(gè)類型，分別為：單胎公羊、單胎母羊、雙胎公羊、雙胎母羊。根據(jù)其各個(gè)類型的日齡和斷奶重進(jìn)行校正。

1.2.2 統(tǒng)計(jì)分析

運(yùn)用SPSS26.0統(tǒng)計(jì)軟件進(jìn)行描點(diǎn)繪制散點(diǎn)圖，運(yùn)用雙變量相關(guān)程序求得斷奶日齡與斷奶重的相關(guān)系數(shù)，并進(jìn)行線性回歸分析，最終得到理想校正值。

2 結(jié)果與分析

2.1 斷奶體重校正

2.1.1 斷奶體重與斷奶日齡的散點(diǎn)圖

為了簡潔直觀的觀察斷奶日齡與斷奶重之間的變化趨勢，特制定平面散點(diǎn)圖將每一對觀測值描點(diǎn)，形成散點(diǎn)圖，見圖1，橫軸x軸代表斷奶日齡（d），斷奶重（kg）在y軸上表示。初步判斷斷奶日齡與斷奶重

圖1 斷奶日齡與斷奶重的散點(diǎn)圖

[4]呈直線回歸關(guān)系。

2.1.2 相關(guān)系數(shù)與回歸方程

為了將影響斷奶校正的因素降至最低，故將出生類型與性別分類校正。在斷奶日齡中顯示，四組類別的斷奶日齡差距跨度較大，范圍較廣，增加了不同日齡斷奶羔羊同一時(shí)間斷奶造成羔羊斷奶重的誤差，體現(xiàn)了校正斷奶重的必要性。

表1、2、3、4的結(jié)果顯示，單胎公羊、單胎母羊、雙胎公羊、雙胎母羊的斷奶日齡與斷奶重的回歸方程分別為y=0.205x+12.404；y=0.195x+11.472；y=0.212x+3.218；y=0.216x+1.203。四組類別的相關(guān)系數(shù)分別為0.468、0.571、0.503、0.516，均有極顯著的影響（P＜0.01），并大于總體相關(guān)系數(shù)。按照系數(shù)顯著，回歸方程一定顯著；系數(shù)不顯著，回歸方程一定不顯著的原則，表中的結(jié)果數(shù)據(jù)可信度高，為接下來的分析提供了理論基礎(chǔ)。

表1 羔羊斷奶日齡和斷奶體重相關(guān)系數(shù)和回歸方程

2.1.3 斷奶體重理論值和校正系數(shù)

斷奶日齡每間隔5 d劃分為一組，取中位數(shù)和目標(biāo)斷奶日齡100 d，將其代入單胎公羊表2、單胎母羊表3、雙胎公羊表4、雙胎母羊表5的回歸方程，分別得出斷奶重的理論值和校正日齡體重值。校正系數(shù)為校正日齡體重值除以斷奶重的理論值。最后，用校正系數(shù)乘以實(shí)際斷奶重得出校正體重值。

表2 單胎公羊不同日齡斷奶重理論值和校正系數(shù)

表3 單胎母羊不同日齡斷奶重理論值和校正系數(shù)

表4 雙胎公羊不同日齡斷奶重理論值和校正系數(shù)

表5 雙胎母羊不同日齡斷奶重理論值和校正系數(shù)

3 討 論

本文所采用“日齡校正方法”結(jié)果顯示回歸方程達(dá)到極顯著的水平（P<0.01），說明用本文的四組回歸方程校正系數(shù)是可靠的。而散點(diǎn)圖顯示斷奶日齡與斷奶體重呈正相關(guān)，說明了校正斷奶體重的必要性，能夠進(jìn)行更好的選種選育。王瑞金[5]鄭浩[6]等人同樣運(yùn)用了求得相關(guān)系數(shù)與回歸方程最后通過日齡來校正體重。

3.1 校正斷奶重的研究進(jìn)展

在早些年前，當(dāng)時(shí)的計(jì)算機(jī)系統(tǒng)還不算發(fā)達(dá)的時(shí)候，李蒙召[7]張延齡[8]通過影響羔羊斷奶重的因素如單雙羔、胎次、哺乳期等進(jìn)行羔羊分類，提出校正羔羊斷奶重的方法，運(yùn)用日增重法、回歸法、多項(xiàng)式的方法校正羔羊的體重。齊侃虎等[9]運(yùn)用最小二乘法制定校正初生重與斷奶重的系數(shù)。為后續(xù)新興的計(jì)算機(jī)統(tǒng)計(jì)軟件算法奠定了基礎(chǔ)。劉望宏[10]將三種不同品種的豬進(jìn)行了早期和全期生長性狀的體重測定，運(yùn)用線性回歸及非線性回歸的模型，根據(jù)數(shù)據(jù)的特點(diǎn)進(jìn)行了抽樣、分階段、連續(xù)測定的方法，選取了最適合的模型。結(jié)果顯示：Logistic、Gompertz、Saturation等模型適合體重性狀與日齡的關(guān)系，相關(guān)擬合平均指數(shù)達(dá)到0.9，達(dá)到極顯著（P<0.01）其中Gompertz模型更具有生物學(xué)意義，Saturation的模型回歸更接近一條直線，最適合做校正公式。呂鳳林[11]對內(nèi)江仔豬運(yùn)用不同的校正方法并進(jìn)行比較分析，結(jié)果以平均日增重和Tayl or公式對內(nèi)江仔豬體重校正較理想。邵學(xué)芬等[12]采用“全樣本”、“等日齡平均”、“期間5日平均”、“升序10羊平均”四種模型斷對奶日齡的分割方法進(jìn)行斷奶重的校正，結(jié)果顯示，四種方法的擬合度均很高，t檢驗(yàn)的結(jié)果均為極顯著（P<0.01），通過相關(guān)系數(shù)的比較和與原始樣本斷奶重的差異分析，“升序10羊平均”的模型最適合。此種方法是近年來國內(nèi)鮮見的僅用斷奶日齡劃分模型的方法，具有創(chuàng)新性，對斷奶重校正提供了諸多靈感。

3.2 斷奶重校正在實(shí)際生產(chǎn)應(yīng)用中存在的問題

對于斷奶重校正的方法，大多數(shù)的科研技術(shù)人員會(huì)采用最簡便的“日增重校正法”，通過羔羊斷奶前的平均日增重乘以校正天數(shù)。但此種方法的弊端在于不能保證羔羊每天的日增重相同，誤差較大。而建立數(shù)學(xué)模型進(jìn)行斷奶重校正，保證其相關(guān)系數(shù)顯著、擬合值接近1，t檢驗(yàn)顯著，產(chǎn)生的誤差遠(yuǎn)遠(yuǎn)小于日增重法。而建立數(shù)學(xué)模型需要在養(yǎng)殖場和牧民中的相關(guān)科技人員掌握一定的理論知識(shí)和學(xué)科基礎(chǔ)。但目前看來大部分人員實(shí)際動(dòng)手操作能力強(qiáng)，但理論知識(shí)都達(dá)不到相關(guān)要求，很難建立起數(shù)學(xué)模型。

4 結(jié)論

通過劃分生長類型和性別不同分別建立四種斷奶天數(shù)與斷奶重雙變量相關(guān)系數(shù)，得出斷奶天數(shù)和斷奶重均有極顯著相關(guān)（P<0.01）的結(jié)果，建立散點(diǎn)回歸圖，通過趨勢線初步判斷斷奶天數(shù)與斷奶重呈線性回歸，再建立斷奶重與斷奶天數(shù)的線性回歸方程，按照5日間隔將斷奶天數(shù)劃分14組，計(jì)算各組中數(shù)，分別代入目標(biāo)校正天數(shù)100 d求出校正系數(shù)，與實(shí)際斷奶重相乘得出校正斷奶重。為后續(xù)的選種選育建立特哈雜交F2代種羊遺傳評估體系提供了一定的參考標(biāo)準(zhǔn)。