云上黑山羊斷奶日齡與斷奶重相關(guān)性分析及斷奶重的校正

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

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

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

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

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

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

1 材料與方法

1.1 羊群概況

2019年云南省種羊繁育推廣中心肉牛隊飼養(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)管理水平差異，在此只以其中一個飼養(yǎng)員所管理的羔羊進行分析。

1.2 樣本羊只

樣本羊只為肉牛隊繁育生產(chǎn)并在2019年度斷奶的全部云上黑山羊羔羊，共計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公羔與母羔分別進行分析

眾所周知，公羔與母羔生長發(fā)育差異較大，通常分開進行分析。

1.3.2制作樣本散點趨勢圖

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

圖1 公羔斷奶日齡與斷奶重散點趨勢圖

圖2 母羔斷奶日齡與斷奶重散點趨勢圖

1.3.3樣本處理

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

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

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

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

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

1.3.4線性擬合及分析

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

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

利用數(shù)學(xué)模型y=a+bx方程式計算出羔羊不同日齡的體重，稱為理論體重；根據(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)分析

對4種處理的樣本進行直線回歸擬合并進行相關(guān)系數(shù)計算及顯著性檢驗。結(jié)果：4種樣本處理方法羔羊斷奶日齡與斷奶重的線性擬合度均很高，t檢驗均為極顯著(P<0.01)。見表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é)模型擬合越好，實用價值越高。在4種樣本處理中，相關(guān)系數(shù)r值為期間5日平均>升序10羊平均>等日齡平均>全樣本。僅從r值判斷，“期間5日平均”處理擬合度最高、“全樣本”處理擬合度最低。然而，4種處理自由度ν差異較大；自由度ν越大，出現(xiàn)大幅偏離趨勢線的樣本數(shù)也就會越多，相關(guān)系數(shù)r值就會變小。為了減小自由度ν值對擬合度評判的影響，此處采用了“t檢驗值”與“tν,0.01值”比值作為輔助評判指標(biāo)。t/tν,0.01值越大，說明相關(guān)系數(shù)r值t檢驗越顯著。t/tν,0.01值為全樣本>升序10羊平均>等日齡平均>期間5日平均。

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

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

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

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

2.2.3比選結(jié)論

通過綜合比較分析，結(jié)論為“升序10羊平均”處理所擬合的數(shù)學(xué)模型最佳，可應(yīng)用于實際生產(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é)模型圖示見圖3。

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

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

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

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

3 討論

3.1 原始樣本及其處理方法

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

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

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

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

3.3 斷奶時間的確定

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

3.4 公羔與母羔的增重變化

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

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

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