高原型藏羊早期生長(zhǎng)曲線擬合分析

張軍霞，方德成，賀 娜

（青海大學(xué)農(nóng)牧學(xué)院，青海 西寧 810016）

藏系綿羊簡(jiǎn)稱藏羊，主要分布在青海、西藏、甘肅、四川等海拔3 000 m 及以上的高原地區(qū)[1-2]。藏羊長(zhǎng)期在高寒缺氧的環(huán)境條件下生存，逐漸形成了體質(zhì)結(jié)實(shí)、肉質(zhì)良好、適應(yīng)性和抗病力強(qiáng)等特征[3-4]。是廣大牧民賴以生存的生活資料[5]。家畜的早期生長(zhǎng)發(fā)育情況對(duì)于整個(gè)生長(zhǎng)發(fā)育過程有著至關(guān)重要的作用，因此，了解家畜早期的生長(zhǎng)發(fā)育狀況是非常重要的。生長(zhǎng)曲線的擬合分析是研究家畜在不同生長(zhǎng)發(fā)育階段特點(diǎn)的一種常用方法，是家畜育種和生產(chǎn)的基礎(chǔ)，理想的生長(zhǎng)曲線模型對(duì)家畜飼養(yǎng)管理、育種實(shí)踐及控制生長(zhǎng)發(fā)育具有重要指導(dǎo)意義[6-8]。為此，通過建立羔羊早期生長(zhǎng)曲線模型，動(dòng)態(tài)分析高原型藏羊的生長(zhǎng)發(fā)育強(qiáng)度和規(guī)律，以期為高原型藏羊的生產(chǎn)潛力、體格、體重增長(zhǎng)的研究提供數(shù)據(jù)積累，為高原型藏羊飼養(yǎng)管理和選育提供理論依據(jù)。

1 材料與方法

1.1 試驗(yàn)羔羊的選擇

試驗(yàn)羊只來自青海省海北藏族自治州剛察縣吉泉鄉(xiāng)寧夏村養(yǎng)殖合作社，隨機(jī)選取年齡、體重和胎次相近的高原型藏羊妊娠母羊100 只，待母羊分娩后隨機(jī)選擇公、母羔羊各40 只，測(cè)定正常生長(zhǎng)發(fā)育羔羊的體重。

1.2 飼養(yǎng)管理

采用自然放牧加補(bǔ)飼的飼養(yǎng)方式，早、晚補(bǔ)飼精料，白天放牧。放牧地為天然草地。初生羔羊以母乳為主，從15 日齡開始每只羔羊每日補(bǔ)飼精料50 g，2 月齡前根據(jù)羔羊采食情況每周調(diào)整一次精料補(bǔ)飼量，2 月齡斷奶后每天補(bǔ)飼精料200 g，每2 周調(diào)整一次，每天清掃羊舍和定時(shí)清洗舍內(nèi)水槽，保證飲水清潔，自由飲水。羔羊2 月齡斷奶后采取分群飼養(yǎng)，日糧配方及營(yíng)養(yǎng)水平如表1 所列。

1.3 測(cè)定方法

參照綿、山羊生產(chǎn)性能測(cè)定技術(shù)規(guī)范(NY/T 1236–2006)[9]測(cè)定高原型藏羊體重指標(biāo)、初生重及1 ? 6 月齡的體重指標(biāo)，所有測(cè)定工作均在08: 00 之前羔羊空腹?fàn)顟B(tài)下進(jìn)行。

1.4 數(shù)據(jù)處理

日增重(G)、相對(duì)生長(zhǎng)率(R)和擬合度的計(jì)算公式如下：

式中：Wt為結(jié)束測(cè)定時(shí)的體重， W0為開始測(cè)定時(shí)體重，t 測(cè)定天數(shù)[5]。

式中：R2為曲線擬合度，Wi為體重實(shí)測(cè)值，W 為擬合曲線體重估計(jì)值，為體重實(shí)測(cè)平均值。

采 用Logistic、Gompertz 和Von Bertalanffy 3 種非線性擬合模型進(jìn)行生長(zhǎng)曲線擬合。利用Excel 和SPSS 19.0 對(duì)數(shù)據(jù)進(jìn)行初步整理和統(tǒng)計(jì)分析，應(yīng)用SPSS 19.0 軟件非線性回歸方程進(jìn)行分析(表2)。

表1 試驗(yàn)日糧組成及營(yíng)養(yǎng)水平Table 1 Ingredients and nutrient levels of experimental diets

2 結(jié)果與分析

2.1 高原型藏羊早期生長(zhǎng)發(fā)育規(guī)律

高原型藏羔羊體重隨年齡的增長(zhǎng)而增加(表3)，2 月齡公、母羔羊日增重和1 月齡相對(duì)生長(zhǎng)率達(dá)到峰值，隨后有明顯降低趨勢(shì)，表明羔羊初生時(shí)生長(zhǎng)強(qiáng)度較大，隨著年齡的增加而逐漸減緩。

2.2 高原型藏羊早期生長(zhǎng)曲線擬合

利用 Gompertz、Logistic 和 Von Bertalanffy 模型擬合高原型藏羔羊生長(zhǎng)曲線結(jié)果如表4 所列。3 種模型的擬合度均較高，均能較好地?cái)M合高原型藏羔羊的生長(zhǎng)曲線，公、母羔羊擬合模型均以Von Bertalanffy 數(shù)學(xué)模型擬合度最高，擬合效果較另外兩個(gè)模型更佳，在Von Bertalanffy 模型擬合生長(zhǎng)曲線中，公羔拐點(diǎn)月齡、拐點(diǎn)體重和最大日增重分別為1.35、12 kg 和246.67 g·d?1，母羔拐點(diǎn)月齡、拐點(diǎn)體重和最大日增重分別為1.36、11.96 kg 和239.33 g·d?1。

表2 3 種曲線模型表達(dá)式及相關(guān)參數(shù)Table 2 The formulae and related parameters of three growth curve models

表3 高原型藏羔羊各月齡體重Table 3 Monthly body weight for plateau-type Tibetan sheep

表4 3 種生長(zhǎng)曲線擬合參數(shù)估計(jì)值及擬合度Table 4 Estimated parameter values and fitting degree of the three growth curve models

2.3 高原型藏羊早期生長(zhǎng)曲線繪制

根據(jù)高原型藏羊公、母羔早期生長(zhǎng)階段的體重、日增重和相對(duì)生長(zhǎng)率的Gompertz、Logistic 和Von Bertalanffy 模型擬合理論值和實(shí)際測(cè)定分別繪制出累積生長(zhǎng)曲線(圖1)、絕對(duì)生長(zhǎng)曲線(圖2)和相對(duì)生長(zhǎng)曲線(圖3)，可以看出，高原型藏羊公、母羔實(shí)際體重值與3 種模型擬合的理論值接近。 公羔在2 月齡斷奶前實(shí)際日增重高于理論值，斷奶后低于理論值，在5 ? 6 月齡間開始逐漸高于理論值。母羔在斷奶前日增重高于理論值，斷奶后逐漸低于理論值，在4 ? 5 月齡間開始趨于一致。

圖1 公、母羔羊?qū)崪y(cè)體重與模型估計(jì)體重累積生長(zhǎng)曲線Figure 1 Growth curves based on the measured and estimated values of the body weights of ram and ewe lambs

圖2 公、母羔羊?qū)崪y(cè)體重與模型估計(jì)體重絕對(duì)生長(zhǎng)曲線Figure 2 Growth curves based on the measured and estimated values of gains in ram and ewe lambs

圖3 公、母羔羊?qū)崪y(cè)體重與模型估計(jì)體重相對(duì)生長(zhǎng)曲線Figure 3 Growth curves based on the measured and estimated values of the relative growth of ram and ewe lambs

3 討論

遺傳因素、營(yíng)養(yǎng)水平和飼養(yǎng)管理等均能影響動(dòng)物的生長(zhǎng)發(fā)育過程[9]。高原型藏羊公羔羊在生長(zhǎng)早期體重增長(zhǎng)要高于母羔羊，表明高原型藏羊公羔羊體格較母羔羊體格大，說明公羔羊在0 ? 6 月齡的生長(zhǎng)發(fā)育速度高于母羔羊。高原型藏羔羊在初生至2 月齡階段是體重增長(zhǎng)的旺盛期，是高原型藏羊生長(zhǎng)發(fā)育過程中的重要環(huán)節(jié)，改善這一階段的飼養(yǎng)管理可大大提高原型藏羊后期生長(zhǎng)發(fā)育速率[10]。利用Gompertz、Logistic 和Von Bertalanffy 3 種非線性生長(zhǎng)曲線擬合模型均能對(duì)高原型藏羊公、母羔羊生長(zhǎng)發(fā)育進(jìn)行生長(zhǎng)曲線擬合分析，R2均在0.990 以上，能較好地模擬其體重增長(zhǎng)的趨勢(shì)和特點(diǎn)[11]。3 個(gè)數(shù)學(xué)模型擬合結(jié)果表明，對(duì)高原型藏羊公、母羔的成熟體重、拐點(diǎn)月齡、拐點(diǎn)體重的擬合，Von Bertalanffy模型較Gompertz 和Logistic 模型更為準(zhǔn)確，擬合度更高。陳玲等[12]運(yùn)用 Logistic、Gompertz 和 Von Bertalanffy 3 種模型對(duì)湖羊的早期生長(zhǎng)發(fā)育規(guī)律進(jìn)行了擬合研究，結(jié)果表明Logistic 和Gompertz 模型對(duì)湖羊的生長(zhǎng)發(fā)育擬合效果最好，并根據(jù)此模型建立了湖羊的生長(zhǎng)發(fā)育曲線。劉遠(yuǎn)等[13]利用上述3 種模型對(duì)戴云山羊進(jìn)行分析，結(jié)果發(fā)現(xiàn)Von Bertalanffy模型的擬合結(jié)果更接近實(shí)際測(cè)量值，擬合效果最佳。索效軍等[14]對(duì)湖北黑頭羊早期生長(zhǎng)發(fā)育規(guī)律進(jìn)行了擬合研究，其中公羊最佳擬合模型為 Gompertz，母羊最佳擬合模型為 Von Bertalanffy。本研究中，從3 種數(shù)學(xué)模型擬合估計(jì)值來看，Von Bertalanffy 模型對(duì)高原型藏羊羔羊體重生長(zhǎng)擬合結(jié)果更接近于實(shí)測(cè)值；Logistic、Gompertz 兩個(gè)模型與實(shí)測(cè)值之間有一定的波動(dòng)。各模型擬合生長(zhǎng)過程相比，絕對(duì)生長(zhǎng)和相對(duì)生長(zhǎng)的實(shí)測(cè)值在生長(zhǎng)早期呈現(xiàn)較大波動(dòng)。梁新亮等[15]測(cè)定了在高寒地區(qū)舍飼條件下杜泊羊♂ ×湖羊♀(杜湖)、澳洲白羊♂ × 湖羊♀(澳湖)和特克塞爾羊♂ × 湖羊♀(特湖)3種不同雜交組合F1代羔羊早期生長(zhǎng)性能，比較分析了3組雜交羔羊的生長(zhǎng)觀測(cè)值和模型估計(jì)值，結(jié)果發(fā)現(xiàn)，兩種模型擬合度(R2)都在0.980 以上，Gopmpertz 模型對(duì)生長(zhǎng)擬合效果優(yōu)于 Logistic 模型。袁飛等[16]對(duì)南江黃羊生長(zhǎng)發(fā)育規(guī)律的研究結(jié)果與本研究結(jié)果類似。本研究中絕對(duì)生長(zhǎng)曲線圖顯示，公、母羔羊在2 月齡體重增長(zhǎng)速度均出現(xiàn)了明顯的變化，羔羊在2 月齡斷奶前營(yíng)養(yǎng)主要來源于母羊的乳汁，母羊產(chǎn)后28 d 泌乳量達(dá)到高峰，之后逐漸下降，而2 月齡后羔羊瘤胃并未發(fā)育完全，尚不能完全消化吸收其所采食飼料中的營(yíng)養(yǎng)成分，致使羔羊2 月齡后日增重增長(zhǎng)速度減緩，加之?dāng)嗄虘?yīng)激的影響，最終導(dǎo)致2 月齡羔羊增重改變。3 種數(shù)學(xué)模型擬合曲線中公、母羔羊均是以Von Bertalanffy 模型擬合曲線最為切合，在2 月齡后實(shí)測(cè)生長(zhǎng)曲線與Von Bertalanffy 模型擬合曲線逐漸趨于重合，這表明隨著對(duì)環(huán)境條件的逐漸適應(yīng)以及飼養(yǎng)管理水平的改善，加快了高原型藏羊的生長(zhǎng)。在生產(chǎn)中應(yīng)綜合考慮飼養(yǎng)管理的各個(gè)環(huán)節(jié)，保證羔羊早期平穩(wěn)、快速生長(zhǎng)，充分發(fā)揮高原型藏羊生產(chǎn)潛力。

4 結(jié)論

在 高 原 型 藏 羊0 ? 6 月 齡 階 段，Logistic、Von Bertalanffany 3 種生長(zhǎng)曲線擬合度均超過 0.99，能較好地對(duì)高原型藏羊公、母羊生長(zhǎng)曲線進(jìn)行擬合，其中，Von Bertalanffy 模型擬合度高于 Logistic 和 Gompertz模型。在 Von Bertalanffy 模型擬合生長(zhǎng)曲線中，公羊拐點(diǎn)月齡、拐點(diǎn)體重和最大日增重分別為1.35、12 kg和246.67 g·d?1，母羊拐點(diǎn)月齡、拐點(diǎn)體重和最大日增重分別為1.36、11.96 kg 和239.33 g·d?1。通過擬合高原型藏羊的生長(zhǎng)曲線，可為動(dòng)態(tài)分析高原型藏羊的生長(zhǎng)發(fā)育、指導(dǎo)飼養(yǎng)管理和為羔羊選育提供可靠依據(jù)。