魯中肉羊體重與體尺指標(biāo)的相關(guān)性及回歸分析

楊存明，張夢(mèng)華，何軍敏，張曉雪，劉桂芬，魏晨，任一帆，毛靜藝，李雪，張文靜，田可川，黃錫霞

(1. 新疆農(nóng)業(yè)大學(xué)動(dòng)物科學(xué)學(xué)院，新疆 烏魯木齊 830052；2. 山東省農(nóng)業(yè)科學(xué)院畜牧獸醫(yī)研究所，山東 濟(jì)南 250100)

綿羊育種以提高其重要經(jīng)濟(jì)性狀的生產(chǎn)水平為主要目的[1]。 魯中肉羊是以白頭杜泊羊?yàn)楦副?，湖羊?yàn)槟副?，雜交繁育十余年育成的肉用綿羊新品種，具有成長速度快、繁殖率高、產(chǎn)肉性能好、適應(yīng)性強(qiáng)等優(yōu)良特性[2-3]。 生長性狀是肉羊育種的重要考察因素，體尺和體重都是衡量生長發(fā)育的重要指標(biāo)，體尺數(shù)據(jù)反映畜體的生長發(fā)育狀況及各部位間的相對(duì)發(fā)育關(guān)系，體現(xiàn)家畜個(gè)體或群體的生產(chǎn)性能，而體重則是檢驗(yàn)飼養(yǎng)管理水平的主要依據(jù)[4]。 對(duì)體重與體尺指標(biāo)進(jìn)行相關(guān)性分析、通徑系數(shù)分析及逐步回歸分析，有助于提高相關(guān)性狀的選擇精度，為保障綿羊生產(chǎn)系統(tǒng)的整體效率和盈利能力提供依據(jù)[5-6]。 前人采用這些方法對(duì)不同品種羊，如杜泊羊、薩?？四秆?、陶賽特母羊[7]、青海無角歐拉羊羔羊[8]、大尾寒羊[9]、科羅基羊[10]、烏達(dá)羊[11]、杜湖雜種羔羊[12]、云上黑山羊[13]、澳洲白羊[14]、公灘羊[15]等進(jìn)行了研究。

本研究以周歲魯中肉羊?yàn)閷?duì)象，測(cè)定體重和體尺指標(biāo)參數(shù)，對(duì)相關(guān)數(shù)據(jù)進(jìn)行了初步統(tǒng)計(jì)，并對(duì)體重和體尺指標(biāo)進(jìn)行了相關(guān)性分析、通徑分析以及多元回歸分析，建立了最優(yōu)回歸方程，以期為優(yōu)化魯中肉羊的飼養(yǎng)管理結(jié)構(gòu)奠定基礎(chǔ)，并為后期生產(chǎn)實(shí)踐以及品種選育工作提供參考。

1 材料與方法

1.1 試驗(yàn)動(dòng)物

試驗(yàn)于山東嬴泰農(nóng)牧科技有限公司進(jìn)行。 隨機(jī)選取周歲魯中肉羊共計(jì)1 289 只，其中公羊330只(平均體重83.49 kg)，母羊959 只(平均體重60.73 kg)，所有羊只均發(fā)育良好，健康無病，實(shí)行常規(guī)免疫，舍飼圈養(yǎng)。

1.2 測(cè)定指標(biāo)及方法

參照《家畜育種學(xué)》[1]和《綿、山羊生產(chǎn)性能測(cè)定技術(shù)規(guī)范》[16]對(duì)周歲魯中肉羊公羊和母羊的體重(Y)、體高(X1)、體長(X2)、胸圍(X3)、胸寬(X4)、胸深(X5)、管圍(X6)進(jìn)行測(cè)定，每個(gè)指標(biāo)2次重復(fù)。

1.3 數(shù)據(jù)處理與分析

通過Microsoft Excel 2019 對(duì)試驗(yàn)數(shù)據(jù)進(jìn)行初步整理，采用拉依達(dá)準(zhǔn)則(平均值±3 倍標(biāo)準(zhǔn)差)將指標(biāo)異常值進(jìn)行剔除。 采用SPSS 25.0 軟件對(duì)數(shù)據(jù)進(jìn)行描述性統(tǒng)計(jì)，以各體尺指標(biāo)為自變量，體重為因變量進(jìn)行相關(guān)性分析，并采用多元回歸分析中的逐步回歸法建立周歲魯中肉羊體重和體尺指標(biāo)的最優(yōu)回歸方程，最后通徑分析各體尺性狀對(duì)體重的直接和間接作用。

2 結(jié)果與分析

2.1 魯中肉羊體重和體尺指標(biāo)的測(cè)定

周歲魯中肉羊公羊和母羊的體重和體尺指標(biāo)測(cè)定結(jié)果見表1、表2。 結(jié)果顯示，周歲魯中肉羊公羊體重的變異系數(shù)較大，為9.17%，母羊體重變異系數(shù)為2.67%，公羊和母羊體尺性狀中變異系數(shù)最大的分別是胸寬和體高，分別為5.01%和6.91%。

表1 周歲魯中肉羊公羊體重、體尺指標(biāo)測(cè)定結(jié)果(n＝330)

表2 周歲魯中肉羊母羊體重、體尺指標(biāo)測(cè)定結(jié)果(n＝959)

2.2 魯中肉羊體重與體尺指標(biāo)的相關(guān)性分析

從表3 可知，魯中肉羊公羊體重與胸深呈顯著負(fù)相關(guān)(P＜0.05)，胸深與胸圍呈顯著正相關(guān)(P＜0.05)，體高與胸深呈極顯著負(fù)相關(guān)(P＜0.01)，體高與管圍間無顯著相關(guān)性(P＞0.05)，其余各性狀間均呈極顯著正相關(guān)(P＜0.01)。 各體尺指標(biāo)與體重間的相關(guān)系數(shù)依次為體高＞胸圍＞胸寬＞體長＞管圍＞胸深。 母羊體長與各體尺指標(biāo)及體重之間均呈極顯著負(fù)相關(guān)(P＜0.01)，其余各體尺指標(biāo)之間以及各體尺指標(biāo)與體重之間均呈現(xiàn)極顯著正相關(guān)(P＜0.01)，各體尺指標(biāo)與體重的相關(guān)系數(shù)依次為胸寬＞胸圍＞體高＞胸深＞管圍＞體長。

表3 魯中肉羊各性狀間相關(guān)系數(shù)分析

2.3 魯中肉羊體重與體尺指標(biāo)的通徑分析

由表4 可知，在魯中肉羊公羊各體尺指標(biāo)中，體高對(duì)體重的直接作用最大，通徑系數(shù)為0.325，其他指標(biāo)對(duì)體重的直接作用由高到低依次為胸圍、胸寬、體長和胸深，通徑系數(shù)分別為0.289、0.228、0.218 和-0.157。 胸圍對(duì)體重的間接作用最大，通徑系數(shù)為0.615，其他指標(biāo)對(duì)體重的間接作用由高到低依次為體高、胸寬、體長、胸深，通徑系數(shù)分別為0.552、0.494、0.375 和-0.011。 由此可見，體高、胸圍和胸寬對(duì)魯中肉羊公羊體重的影響無論是直接作用還是間接作用均較大，說明其體高、胸圍和胸寬與體重密切相關(guān)。

表4 魯中肉羊公羊體尺對(duì)體重的通徑分析

由表5 可知，在魯中肉羊母羊各體尺指標(biāo)中，胸寬對(duì)體重的直接作用最大，通徑系數(shù)為0.689，其次為胸圍、管圍、體長和體高，通徑系數(shù)分別為0.342、-0.151、-0.279 和-0.294。 胸寬對(duì)體重的間接作用最大，通徑系數(shù)為0.776，其他指標(biāo)由高到低依次為體長、胸圍、管圍和體高。 說明魯中肉羊母羊胸圍和胸寬同樣與體重密切相關(guān)。

表5 魯中肉羊母羊體尺對(duì)體重的通徑分析

2.4 魯中肉羊體尺指標(biāo)與體重的回歸分析

以周歲魯中肉羊體重為因變量，各體尺指標(biāo)為自變量，采用多元回歸分析中的逐步回歸分析法建立多元線性回歸方程。 由表6 可知，模型1～5 中魯中肉羊公羊體重分別與體高、胸圍、體長、胸深、胸寬呈極顯著線性回歸關(guān)系(P＜0.01)。 由此建立以下5 個(gè)線性回歸方程:模型1:Y ＝2.166X1-67.778，校正R2＝0.693；模型2:Y ＝1.381X1+1.181X3-132.499，校正R2＝0.776；模型3:Y ＝1.273X1+0.482X2+0.990X3-144.109，校正R2＝0.801；模型4:Y ＝1.073X1+0.577X2+1.125X3-0.640X5-129.685，校正R2＝0.810；模型5:Y ＝0.845X1+0.554X2+0.818X3+1.245X4-0.868X5-107.903，校正R2＝0.826。 模型5 的校正R2值最高，因此其建立的方程為最優(yōu)回歸方程。

表6 魯中肉羊公羊體重與體尺的回歸分析

由表7 可知，模型1～5 中魯中肉羊母羊體重分別與胸寬、管圍、胸圍、體長、體高呈極顯著線性回歸關(guān)系(P＜0.01)。 由此建立以下5 個(gè)線性回歸方程:模型1:Y ＝1.537X4+22.071，校正R2＝0.800；模型2:Y＝1.616X4-0.998X6+28.625，校正R2＝0.819；模型3:Y ＝0.345X3+1.299X4-1.792X6+11.829，校正R2＝0.836；模型4:Y ＝-0.032X2+0.342X3+1.220X4-1.770X6+15.911，校正R2＝0.839；模型5:Y ＝-0.104X1-0.125X2+0.420X3+1.184X4-1.033X6+16.088，校正R2＝0.844。 模型5 的校正R2值最高，為最優(yōu)回歸方程。

3 討論

3.1 魯中肉羊體重與體尺性狀之間的相關(guān)性分析

本試驗(yàn)分析結(jié)果表明，魯中肉羊的體尺指標(biāo)與體重有極顯著的相關(guān)性(P＜0.01)，其中公羊體高與體重的相關(guān)性最強(qiáng)，其次為胸圍；母羊胸寬與體重的相關(guān)性最強(qiáng)，其次為胸圍。 可能由于本研究所選羊只均為周歲魯中肉羊，體型基本穩(wěn)定，所以其體高、胸寬和胸圍等指標(biāo)與羊只的肌肉、骨骼、內(nèi)臟的生長發(fā)育以及膘情等情況密切相關(guān)。在青海無角歐拉羊羔羊[8]中，體高對(duì)體重的影響最大，其次為胸圍；大尾寒羊[9]體高和胸圍不僅對(duì)體重有較大的直接作用，還輔助其他性狀有較大的間接作用，這均與本研究結(jié)果類似。 而對(duì)科羅基羊[10]、烏達(dá)羊[11]、杜湖雜種羔羊[12]、云上黑山羊[13]、澳洲白羊[14]、灘羊[15]、杜泊羊、特克塞爾羊[17]等研究中，均顯示胸圍與體重的相關(guān)系數(shù)最大，其差異可能是由于品種不同和羊只年齡不一致所導(dǎo)致的。 因此，在魯中肉羊的品種選育過程中，應(yīng)以體重為主，兼顧體高、胸圍、胸寬等其他體尺性狀。

3.2 魯中肉羊體重與體尺指標(biāo)的通徑分析

本研究通過對(duì)周歲魯中肉羊體尺體重間的通徑分析發(fā)現(xiàn)，公羊和母羊體尺指標(biāo)中對(duì)體重的直接作用最大的分別為體高和胸寬，其次均為胸圍，間接作用最大的均為胸寬，其次分別為胸圍和體高。 對(duì)成年歐拉羊[18]、云上黑山羊[13]、酉州烏羊[19]、海南黑山羊[20]的分析均表明胸圍是影響體重的主要體尺指標(biāo)，與本研究中對(duì)魯中肉羊公羊的分析結(jié)果類似，但存在一定差異，可能是由所選品種和飼養(yǎng)水平不同造成的。 在青海無角歐拉羊羔羊[8]的研究中，體高對(duì)體重的直接影響最大，與本研究中對(duì)魯中肉羊公羊的分析結(jié)果一致。本研究中，通徑分析結(jié)果和相關(guān)性分析結(jié)果均表明魯中肉羊公羊主要通過體高、胸圍，母羊主要通過胸寬、胸圍直接或間接影響體重。 因此，在進(jìn)行魯中肉羊體重性狀選育時(shí)，應(yīng)區(qū)分性別，公羊、母羊分別以體高、胸圍和胸寬、胸圍指標(biāo)為主要參考。

3.3 魯中肉羊體尺指標(biāo)與體重的回歸分析

本研究通過多元回歸分析建立了周歲魯中肉羊公羊與母羊體尺指標(biāo)與體重的最優(yōu)回歸方程，與對(duì)青海無角歐拉羊[8]、杜湖雜種羔羊[12]、巴特肯羊[21]、巴美肉羊[22]、布拖黑綿羊育成母羊[23]等品種羊的研究結(jié)果相似。 剩余因子(e ＝1-R2)分別為0.174 和0.156，說明入選回歸方程的體尺性狀對(duì)體重有較大的決定作用。 但還有其他因素未考慮，因此下一步需要對(duì)影響魯中肉羊體重的其他因素進(jìn)行進(jìn)一步分析。

4 結(jié)論

周歲魯中肉羊公羊體高與體重的相關(guān)性最強(qiáng)，母羊胸寬與體重的相關(guān)性最強(qiáng)，其次均為胸圍。 體高和胸寬主要通過直接作用影響體重，而胸圍、體長、胸深和管圍主要通過間接作用影響體重。 因此，在魯中肉羊品種選育過程中，公羊應(yīng)重點(diǎn)關(guān)注體高、胸圍，母羊應(yīng)重點(diǎn)關(guān)注胸寬，從而提高選育速度。 魯中肉羊公羊體尺指標(biāo)與體重的最優(yōu)方程為Y ＝0.845X1+0.554X2+0.818X3+1.245X4-0.868X5-107.903(P＜0.01，R2＝0.826)，母羊?yàn)閅＝-0.104X1- 0. 125X2+ 0. 420X3+ 1. 184X4-1.033X6+16.088(P＜0.01，R2＝0.844)，可作為利用體尺指標(biāo)評(píng)估體重的參考依據(jù)。