邵學(xué)芬，耿文誠(chéng)，向必勇，朱勝全，趙安杰

(云南省種羊繁育推廣中心，云南 尋甸 655204)

云上黑山羊是云南省自主培育的第二個(gè)國(guó)家審定的羊新品種，是中國(guó)第一個(gè)肉用黑山羊新品種、第三個(gè)肉用山羊新品種，2019年4月28日農(nóng)業(yè)農(nóng)村部第168號(hào)公告該新品種通過(guò)國(guó)家畜禽遺傳資源委員會(huì)審定，新品種證書號(hào)為“農(nóng)03新品種證字第18號(hào)”。云上黑山羊由云南省畜牧獸醫(yī)科學(xué)院聯(lián)合云南省種羊繁育推廣中心等單位以努比山羊?yàn)楦副尽⒃茙X黑山羊?yàn)槟副?，?jīng)過(guò)22年5個(gè)世代系統(tǒng)選育而成，具有被毛全黑、生長(zhǎng)發(fā)育快、常年發(fā)情、繁殖力高、產(chǎn)肉性能好、適應(yīng)性強(qiáng)和耐粗飼等優(yōu)良特性。

在種羊繁育生產(chǎn)過(guò)程中，為了便于生產(chǎn)管理，減少人力物力的消耗，降低生產(chǎn)成本，提高種羊生產(chǎn)經(jīng)濟(jì)效益，滿足社會(huì)對(duì)良種的需求，種羊生產(chǎn)單位通常采用批量生產(chǎn)的方式，主要表現(xiàn)在集中批量投入配種、批量產(chǎn)羔、批量斷奶、批量培育以及疫苗批量免疫注射、體重體尺等生產(chǎn)技術(shù)指標(biāo)集中批量監(jiān)測(cè)等。然而，我們既不能保證每批次投產(chǎn)的母羊同時(shí)發(fā)情，更無(wú)法讓它們同時(shí)產(chǎn)羔，致使同一批次生產(chǎn)的種羊存在一定的年齡差異，而體重體尺等生產(chǎn)技術(shù)指標(biāo)測(cè)定卻通常在同一日內(nèi)完成。生產(chǎn)技術(shù)指標(biāo)是判斷種羊優(yōu)劣的主要依據(jù)，而不同日齡的種羊卻用同一日測(cè)定的指標(biāo)來(lái)判斷其優(yōu)劣顯然缺乏說(shuō)服力。因此，將不同日齡種羊的生產(chǎn)技術(shù)指標(biāo)校正為同一日齡，可更為準(zhǔn)確地評(píng)判一只種羊的優(yōu)劣，使種羊生產(chǎn)部門的選種選育工作更加科學(xué)，也能讓養(yǎng)殖戶(場(chǎng))更準(zhǔn)確地了解其所購(gòu)買的種羊質(zhì)量。

斷奶重是種羊質(zhì)量評(píng)判的重要指標(biāo)之一，尤其是在種羊供不應(yīng)求的年份，養(yǎng)殖戶(場(chǎng))會(huì)在羔羊斷奶后迅速將種羊買走。此時(shí)，斷奶重幾乎成了評(píng)判種羊優(yōu)劣的唯一指標(biāo)。作為種羊生產(chǎn)單位，為了維持繁殖群體的數(shù)量及合理的年齡結(jié)構(gòu)，每年均需按比例選留更新羊。在羔羊斷奶時(shí)預(yù)選更新羊是云南省種羊繁育推廣中心近年來(lái)的常規(guī)做法，預(yù)選的更新羊在終選前不向客戶銷售。將不同斷奶日齡羔羊的斷奶重校正為標(biāo)準(zhǔn)斷奶日齡的斷奶重，可將最優(yōu)秀的羔羊選拔到預(yù)留更新羊群中進(jìn)行培育，讓最優(yōu)秀的種羊留在種羊生產(chǎn)單位，為社會(huì)培育出更多的良種。

羔羊斷奶重受多種因素的影響，如性別、胎次、同胎羔羊數(shù)、斷奶日齡、疾病、飼養(yǎng)管理水平等[1]，在此僅研究斷奶日齡與斷奶重的關(guān)系。在選種過(guò)程中，根據(jù)選育目標(biāo)可加入同胎羔羊數(shù)等指標(biāo)與斷奶重一起進(jìn)行綜合評(píng)判，來(lái)彌補(bǔ)單一用斷奶重來(lái)評(píng)判種羊優(yōu)劣的缺陷。

1 材料與方法

1.1 羊群概況

2019年云南省種羊繁育推廣中心肉牛隊(duì)飼養(yǎng)兩群云上黑山羊，繁殖母羊年初存欄數(shù)406只、年末存欄數(shù)351只，年均飼養(yǎng)數(shù)384只；分批次投入繁殖生產(chǎn)，全年共生產(chǎn)10批次；全年新生羔羊707只，斷奶羔羊612只；年平均產(chǎn)羔率214.89%，每胎次單羔率20.67%、雙羔率61.40%、三羔率16.72%、四羔率1.22%。為避免飼養(yǎng)員飼養(yǎng)管理水平差異，在此只以其中一個(gè)飼養(yǎng)員所管理的羔羊進(jìn)行分析。

1.2 樣本羊只

樣本羊只為肉牛隊(duì)繁育生產(chǎn)并在2019年度斷奶的全部云上黑山羊羔羊，共計(jì)327只(其中：♂163只、♀164只)。327只樣本羊只分為5批次生產(chǎn)，分別為2019年2月13日斷奶的104只、5月20日33只、8月16日35只、11月13日80只、12月25日75只。

1.3 分析方法

1.3.1公羔與母羔分別進(jìn)行分析

眾所周知，公羔與母羔生長(zhǎng)發(fā)育差異較大，通常分開(kāi)進(jìn)行分析。

1.3.2制作樣本散點(diǎn)趨勢(shì)圖

為了直觀地看出兩個(gè)相關(guān)變量間的變化趨勢(shì)，將每一對(duì)觀測(cè)值在平面直角坐標(biāo)系上描點(diǎn)，作出散點(diǎn)圖[2]。從斷奶日齡(x)與斷奶重(y)散點(diǎn)趨勢(shì)圖(圖1、圖2)可初步判斷出斷奶日齡與斷奶重呈直線回歸關(guān)系。

圖1 公羔斷奶日齡與斷奶重散點(diǎn)趨勢(shì)圖

圖2 母羔斷奶日齡與斷奶重散點(diǎn)趨勢(shì)圖

1.3.3樣本處理

由于羔羊斷奶重個(gè)體間差異較大，為了更準(zhǔn)確地?cái)M合出斷奶日齡與斷奶重間的直線回歸方程，使其更具有實(shí)用價(jià)值，本研究將樣本分為4種處理。

全樣本：原始樣本，不作處理。

等日齡平均：將原始樣本中斷奶日齡相同羊只歸類為同一個(gè)樣本，計(jì)算出其平均斷奶重。

期間5日平均：在327個(gè)原始樣本中，最小斷奶日齡為74 d、最大斷奶日齡為113 d，以間隔5日齡為一個(gè)梯度劃分樣本，分別計(jì)算每個(gè)梯度內(nèi)樣本的平均斷奶日齡和平均斷奶重。樣本梯度期間分別為：71～75 d、76～80 d、……、111～115 d。

升序10羊平均：將原始樣本按主要關(guān)鍵字“斷奶日齡”、次要關(guān)鍵字“斷奶重”升序排序后，按順序以每10個(gè)原始樣本的平均斷奶日齡和平均斷奶重重新生成一組新樣本。

1.3.4線性擬合及分析

以斷奶日齡為自變量x、斷奶重為依變量y，按直線回歸擬合方法建立數(shù)學(xué)模型y=a+bx方程式，計(jì)算相關(guān)系數(shù)并進(jìn)行顯著性檢驗(yàn)。比較分析樣本不同處理方法的檢驗(yàn)效果，找出最具實(shí)用價(jià)值的數(shù)學(xué)模型。

1.3.5斷奶重校正系數(shù)的計(jì)算

利用數(shù)學(xué)模型y=a+bx方程式計(jì)算出羔羊不同日齡的體重，稱為理論體重；根據(jù)種羊生產(chǎn)單位生產(chǎn)管理的需求核定標(biāo)準(zhǔn)斷奶日齡(云南省種羊繁育推廣中心云上黑山羊標(biāo)準(zhǔn)斷奶日齡為90日齡)，各日齡與標(biāo)準(zhǔn)斷奶日齡理論體重比值即為該日齡斷奶羔羊的校正系數(shù)。

2 結(jié)果與分析

2.1 各處理樣本的線性擬合及相關(guān)分析

對(duì)4種處理的樣本進(jìn)行直線回歸擬合并進(jìn)行相關(guān)系數(shù)計(jì)算及顯著性檢驗(yàn)。結(jié)果：4種樣本處理方法羔羊斷奶日齡與斷奶重的線性擬合度均很高，t檢驗(yàn)均為極顯著(P<0.01)。見(jiàn)表1、表2。

表1 公羔斷奶日齡與斷奶重線性擬合及相關(guān)分析

表2 母羔斷奶日齡與斷奶重線性擬合及相關(guān)分析

2.2 樣本處理方式的比選

2.2.1數(shù)學(xué)模型相關(guān)系數(shù)r值比較

相關(guān)系數(shù)r值越接近1，數(shù)學(xué)模型擬合越好，實(shí)用價(jià)值越高。在4種樣本處理中，相關(guān)系數(shù)r值為期間5日平均>升序10羊平均>等日齡平均>全樣本。僅從r值判斷，“期間5日平均”處理擬合度最高、“全樣本”處理擬合度最低。然而，4種處理自由度ν差異較大；自由度ν越大，出現(xiàn)大幅偏離趨勢(shì)線的樣本數(shù)也就會(huì)越多，相關(guān)系數(shù)r值就會(huì)變小。為了減小自由度ν值對(duì)擬合度評(píng)判的影響，此處采用了“t檢驗(yàn)值”與“tν,0.01值”比值作為輔助評(píng)判指標(biāo)。t/tν,0.01值越大，說(shuō)明相關(guān)系數(shù)r值t檢驗(yàn)越顯著。t/tν,0.01值為全樣本>升序10羊平均>等日齡平均>期間5日平均。

2.2.2用數(shù)學(xué)模型校正樣本羊只斷奶重比較

根據(jù)對(duì)原始樣本的統(tǒng)計(jì)，公羔平均斷奶日齡為96.93 d，母羔平均斷奶日齡為97.66 d，用各處理擬合的數(shù)學(xué)模型將原始樣本羊只校正為平均斷奶日齡的斷奶重，并與原始樣本平均斷奶重進(jìn)行對(duì)比分析(表3、表4)。結(jié)果：4個(gè)數(shù)學(xué)模型校正的羔羊平均斷奶重與原始樣本羔羊平均斷奶重間的差異很小，公羔差異最大的為“期間5日平均”處理，平均斷奶重較原始樣本增加34.61 g、0.18%；母羔差異最大的為“等日齡平均”處理，平均斷奶重較原始樣本增加15.62 g、0.10%。差異絕對(duì)值排序,公羔為全樣本<升序10羊平均<等日齡平均<期間5日平均;母羔為全樣本<升序10羊平均<期間5日平均<等日齡平均。

表3 公羔校正為平均斷奶日齡的斷奶重統(tǒng)計(jì)比較表

表4 母羔校正為平均斷奶日齡的斷奶重統(tǒng)計(jì)比較表

2.2.3比選結(jié)論

通過(guò)綜合比較分析，結(jié)論為“升序10羊平均”處理所擬合的數(shù)學(xué)模型最佳，可應(yīng)用于實(shí)際生產(chǎn)中，回歸方程為公羔y=-6.73+0.27x、r=0.8945(P<0.01)，母羔y=3.49+0.13x、r=0.8515(P<0.01)。數(shù)學(xué)模型圖示見(jiàn)圖3。

圖3 “升序10羊平均”處理羔羊斷奶日齡與斷奶重直線回歸數(shù)學(xué)模型圖

2.3 90日齡斷奶重校正系數(shù)參照表

根據(jù)回歸方程y=-6.73+0.27x(公羔)及y=3.49+0.13x(母羔)計(jì)算出云上黑山羊標(biāo)準(zhǔn)斷奶日齡為90 d的斷奶重校正系數(shù)，見(jiàn)表5。

表5 云上黑山羊90日齡斷奶重校正系數(shù)參照表

3 討論

3.1 原始樣本及其處理方法

在統(tǒng)計(jì)學(xué)上，樣本數(shù)量越大，所分析出來(lái)的結(jié)論越可靠。然而，在大樣本中，往往會(huì)存在一些極端的數(shù)據(jù)，影響到統(tǒng)計(jì)分析者的判斷力(如對(duì)一組數(shù)據(jù)是直線相關(guān)或是曲線相關(guān)的判斷)或是影響到分析結(jié)果(如相關(guān)系數(shù)、顯著性檢驗(yàn))。因此，多數(shù)情況下需要對(duì)原始樣本進(jìn)行分類處理來(lái)形成新的統(tǒng)計(jì)分析數(shù)據(jù)源。

在家畜批量繁殖生產(chǎn)中，由于不可能做到所有母畜同日發(fā)情、同日配種、同日受孕、同日分娩，而斷奶卻在同一天，導(dǎo)致仔幼畜斷奶日齡參差不齊?！暗热正g平均”“期間5日平均”處理樣本的最大缺陷是小樣本數(shù)量不一致，即“試驗(yàn)設(shè)計(jì)中重復(fù)數(shù)不一致”，這也導(dǎo)致了其分析結(jié)果遜于“升序10羊平均”。當(dāng)然，“10羊平均”并不是絕對(duì)的，可根據(jù)原始樣本的數(shù)量及其分布采用“5羊平均”“20羊平均”或是其他重復(fù)數(shù)。

3.2 校正系數(shù)的使用區(qū)間

斷奶日齡越接近標(biāo)準(zhǔn)斷奶日齡，所校正出來(lái)的斷奶重越準(zhǔn)確。本研究使用的樣本最小斷奶日齡為74 d、最大斷奶日齡為113 d，理論上只對(duì)74～113日齡的斷奶羔羊有效。但是，在實(shí)際生產(chǎn)中，經(jīng)常會(huì)出現(xiàn)一些主觀原因或客觀原因?qū)е赂嵫蛭茨馨雌跀嗄?，致使個(gè)別羔羊斷奶時(shí)不足74日齡或是大于113日齡，建議可將校正系數(shù)的使用區(qū)間定為61～120日齡，超越此范圍的不建議使用該方法進(jìn)行斷奶重校正。

3.3 斷奶時(shí)間的確定

在生產(chǎn)中，除了一些主觀原因或客觀原因?qū)е赂嵫蛭茨馨雌跀嗄掏?，還有一個(gè)原因是沒(méi)有計(jì)算羔羊平均出生日期。正確的斷奶時(shí)間為“平均出生日期+標(biāo)準(zhǔn)斷奶日齡”，如某批次羔羊平均出生日期為2020年1月10日，標(biāo)準(zhǔn)斷奶日齡90 d，則應(yīng)該在2020年4月9日這天斷奶。

3.4 公羔與母羔的增重變化

在各處理擬合的回歸方程中，母羔的斜率(回歸系數(shù))約為公羔的50%，說(shuō)明此期間公羔的相對(duì)增重比率遠(yuǎn)大于母羔。從圖3可以看出，2條回歸線在73日齡處交匯，也就是說(shuō)，73日齡前母羔的體重大于公羔，73日齡后公羔的體重大于母羔。也可以理解為，哺乳前期母羔生長(zhǎng)速度大于公羔，哺乳后期公羔生長(zhǎng)速度大于母羔。是否有這種現(xiàn)象，有待今后進(jìn)一步研究。

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

對(duì)于斷奶重的校正，不少科技人員會(huì)使用“日均生長(zhǎng)速度”，因?yàn)檫@種方法容易被絕大多數(shù)科技人員理解并熟練使用。但是，哺乳期的仔幼畜顯然不是每天都有相同的增重，這樣校正出來(lái)的斷奶重誤差較大。使用數(shù)學(xué)模型來(lái)校正斷奶重，只要數(shù)學(xué)模型擬合度高、相關(guān)系數(shù)較大、t檢驗(yàn)顯著，誤差是可控的。然而，建立數(shù)學(xué)模型需要科技人員熟練掌握數(shù)理統(tǒng)計(jì)知識(shí)，目前多數(shù)養(yǎng)殖場(chǎng)(戶)的科技人員實(shí)際操作能力很強(qiáng)，但數(shù)理統(tǒng)計(jì)分析能力相對(duì)較弱，不易建立數(shù)學(xué)模型。使用日均生長(zhǎng)速度校正斷奶重也許是基層科技人員的無(wú)奈之舉，建議科研單位、高校的專家學(xué)者能夠在這方面給予更多的技術(shù)支持，讓最基層的養(yǎng)殖場(chǎng)(戶)早日實(shí)現(xiàn)科學(xué)化、精準(zhǔn)化的生產(chǎn)管理。