肉用西門塔爾牛群體生長曲線擬合及體重與體尺相關(guān)性分析的研究

梁永虎，朱　波，金生云，寶金山，徐凌洋，陳　燕，高　雪，張路培，高會江，李俊雅*

(1.中國農(nóng)業(yè)科學(xué)院北京畜牧獸醫(yī)研究所，北京 100193； 2.錫林郭勒盟烏拉蓋管理區(qū)獸醫(yī)局，烏拉蓋 026321)

西門塔爾牛早在20世紀(jì)初作為兼用品種引入我國[1]，主要用于改良我國黃牛，其雜交后代大約占我國雜交改良牛的50%[2]。由于其良好的生產(chǎn)性能[3]，育種工作者對其進行了大量的基礎(chǔ)性研究工作。H.Z.Fan等[4-8]利用全基因組關(guān)聯(lián)分析方法對肉用西門塔爾牛群體的肉質(zhì)性狀和生長性狀進行了初步研究。B.Zhu等[8-9]首次系統(tǒng)地評估了肉用西門塔爾牛群體的生長發(fā)育、胴體和肉質(zhì)性狀的基因組育種值估計的準(zhǔn)確性，為我國肉牛全基因組選擇研究提供理論基礎(chǔ)。牛紅等[10]采用非求導(dǎo)約束最大似然法估計遺傳力、遺傳相關(guān)和方差組分對肉用西門塔爾牛群體的重要經(jīng)濟性狀進行了遺傳參數(shù)和方差組分估計，為將來制定育種方案和遺傳評估奠定基礎(chǔ)。

畜禽生長曲線擬合是研究和分析畜禽生長發(fā)育規(guī)律的重要方法之一，鑒于動物生長的非線性特性，可以通過幾個非線性生物學(xué)參數(shù)描述動物的生長發(fā)育規(guī)律，即生長曲線模型。Logistic、Brody、Gompertz和Bertallanffy等4種模型是常用的用于描述畜禽生長規(guī)律的非線性數(shù)學(xué)模型，不僅對畜禽生長曲線擬合效果好，而且其中的參數(shù)估計值都能夠反映不同種群在生長率和最大體重方面的差異[11-12]。該模型在指導(dǎo)牛的育種、確定適配年齡和屠宰年齡及控制生長發(fā)育上均起到了重要作用。國外許多肉牛品種，如海福特牛和安格斯牛的生長曲線均已建立[13-14]，國內(nèi)秦川牛、延邊牛和安西牛等品種的生長曲線方程也已建立[15-18]。馬光輝等[19]應(yīng)用Wood不完全伽瑪函數(shù)模型、Nelder逆多項式模型和Wilmink模型等3個數(shù)學(xué)模型,分別對新疆呼圖壁種公牛場肉用西門塔爾牛第一胎、第二胎、第三胎、

第四胎、四胎及以上所有胎次的泌乳曲線進行了擬合,通過SAS統(tǒng)計軟件進行分析，發(fā)現(xiàn)Wood模型的擬合度最好,擬合效果最優(yōu)。而關(guān)于肉用西門塔爾牛群體的生長曲線尚未見報道。

本研究通過對內(nèi)蒙古烏拉蓋地區(qū)肉用西門塔爾牛群體的體重、體高、體斜長、十字部高、胸圍和腹圍等指標(biāo)數(shù)據(jù)進行分析，研究其生長發(fā)育規(guī)律并揭示各性狀間的相互關(guān)系，同時肉用西門塔爾牛體重的預(yù)測方程可為實際的體重測量提供參考依據(jù)。

1　材料與方法

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

本研究所采用數(shù)據(jù)來自于中國農(nóng)業(yè)科學(xué)院北京畜牧獸醫(yī)研究所牛遺傳育種創(chuàng)新團隊。該團隊從2008年開始構(gòu)建肉用西門塔爾牛資源群體，該資源群體位于內(nèi)蒙古錫林郭勒盟烏拉蓋管理區(qū)，經(jīng)過逐年擴群，至2016年肉用西門塔爾牛群體基礎(chǔ)母牛已超過2 000頭，資源群體的父親來自于39頭肉用西門塔爾牛種公牛。每年7～8月份，測定該群體犢牛的生長發(fā)育性狀，包括體重和體尺性狀。為了便于數(shù)據(jù)測量，每年10月份將烏拉蓋地區(qū)5～9月齡的肉用西門塔爾牛運送到北京金維福仁清真食品有限公司，并按照統(tǒng)一的飼養(yǎng)管理方法進行集中育肥，在集中育肥過程中，每3個月需測量一次體重和體尺數(shù)據(jù)。

1.2　擬合生長曲線模型

分別選擇Logistic、Brody、Gompertz和Bertallanffy 4種生長曲線模型擬合其生長曲線(表1)，根據(jù)曲線方程可反映出生長發(fā)育規(guī)律的曲線拐點體重及達到拐點體重的拐點月齡等指標(biāo)[20]。

表1生長曲線模型及特征

Table1Growthcurvemodelsandcharacteristics

生長曲線模型Growthcurvemodel公式Formula拐點體重/kgInflectionpointweight拐點月齡Inflectionpointmonth相對生長率Therelativegrowthrate最大日增重/kgThebiggestdailygainLogisticY=A/(1+Be(-kt))A/2(LnB)/kk(1-wi/A)kw/2GompertzY=A×e-B×exp(-kt)A/e(LnB)/kk(lnA-lnwi)kwBertallanffyY=A×(1-B×e-kt)38A/27(Ln3B)/k3k((A/wi)1/3-1)3kw/2BrodyY=A×(1-B×e-kt)__k((A/wi)-1)kw/2

方程中，A為成熟體重或體尺，B為常數(shù)尺度，k為順勢相對生長率，wi為觀測值體重，t為生長時間,e為自然對數(shù)

Ais mature body weight or body measurements,Bis constant scale,kis the relative growth rate,wiis observed weight,tis the growth time，eis the natural logarithm in equations

R2是用來評價擬合曲線方程優(yōu)劣的指標(biāo)，即擬合度：

1.3　統(tǒng)計分析

使用SAS 9.4軟件中的NLIN過程對整理后的數(shù)據(jù)進行曲線模型方程參數(shù)的最優(yōu)估計[21]，采用Gauss-Newton法進行環(huán)迭代來估計參數(shù)使殘差平方和最小，收斂標(biāo)準(zhǔn)為10-8。同時對肉用西門塔爾牛體尺及體重進行相關(guān)性分析，以及利用體尺性狀對體重進行回歸分析，獲得對肉用西門塔爾牛體重的回歸預(yù)測方程。采用于志等[22]對中國荷斯坦牛育成牛生長發(fā)育規(guī)律的相關(guān)分析方法，計算公式：

Y=ba+b1X1+b2X2+b3X3+b4X4+b5X5

公式中，Y為體重，X1為體高，X2為十字部高，X3為體斜長，X4為胸圍，X5為腹圍，ba為常數(shù)，b1、b2、b3、b4、b5為回歸系數(shù)。

相關(guān)系數(shù)r根據(jù)Pearson公式進行分析，計算模型：

上式中，i=1，2，3，4，5，N為試驗牛頭數(shù)，X1、X2、X3、X4、X5同上，Y為體重。

2　結(jié)　果

2.1　肉用西門塔爾牛不同月齡體尺體重統(tǒng)計結(jié)果

如表2 所示，肉用西門塔爾牛從出生到20月齡體重、體尺平均值呈逐增態(tài)勢，但各階段增幅不一，體重自出生到12月齡增幅較大，13～20月齡增幅緩慢。對于體高、體斜長、十字部高、胸圍和腹圍均是在出生到6月齡增速最快，之后增速緩慢，15～20月齡增幅最小。從標(biāo)準(zhǔn)誤可以看出，隨著月齡的增加，肉用西門塔爾牛個體間的體尺波動較小。而體重是由出生最小的標(biāo)準(zhǔn)誤6.99逐漸增加到18月齡的64.09，個體間的差異較大。

2.2　肉用西門塔爾牛4種數(shù)學(xué)模型各參數(shù)估計值

通過SAS 9.4運行結(jié)果匯總4種模型的參數(shù)估計值如表3所示，可以看出4種模型都能很好擬合肉用西門塔爾牛的各項指標(biāo)，除體重外其他各項擬合度均在0.99以上。Bertallanffy模型在20月齡時最大體重可達831.40 kg, 而Logistic模型在20月齡是僅為617.20 kg, Gompertzhe模型和 Brody模型則分別為764.00 kg和651.20 kg，處于前兩者之間。但順勢相對生長率最大的Logistic模型較最小的Bertallanffy模型高0.15。

生長曲線拐點是反映動物生長率最大的時間點，Gompertz和Bertallanff生長曲線模型的拐點是固定的，如表4所示，其中Logistic模型在9.46月齡時到達拐點，體重為308.55 kg，最大日增重為1.140 kg,對于模型Gompertz，生長至8.47月齡時，最大日增重是1.051 kg，隨后逐漸下降，此時體重為281.81 kg。Bertallanff模型的最大日增重最小1.007 kg,在7.15月齡時，體重達到268.62 kg。因為Brody模型的拐點月齡和體重是不穩(wěn)定的，所以不存在相關(guān)參數(shù)值。

表2肉用西門塔爾牛體尺、體重統(tǒng)計結(jié)果

Table2ThestatisticsofbodymeasurementsandbodyweightofSimmentalbeefcattle

月齡Month個數(shù)/頭Number體重/kgBodyweight體高/cmWithersheight十字部高/cmHipheight體斜長/cmHodylength胸圍/cmChestcircumference腹圍/cmAbdominalgirth054438.33±6.9972.57±6.5178.92±7.1067.81±6.4276.35±6.4775.83±7.86166369.75±25.4883.00±9.9388.10±9.3384.73±9.5399.77±11.80109.14±22.172487105.57±11.8590.79±6.3297.00±8.5795.88±9.04114.56±9.26126.39±12.073820135.48±23.0798.00±7.83101.79±5.30103.97±7.62122.17±11.48136.34±11.974318166.86±24.63101.08±6.94105.11±4.75108.18±6.70129.15±9.26147.81±10.355842184.67±26.37102.36±6.72107.40±6.46112.57±8.48134.08±9.26155.73±10.696684210.49±28.86104.98±5.13110.29±4.60115.89±7.39140.27±8.76164.10±12.247336221.23±34.70105.29±5.58111.18±4.67117.55±5.56141.68±8.13167.77±13.178587261.82±44.16112.11±4.68117.44±5.00122.31±6.01154.34±13.80181.00±14.629589266.45±38.93113.63±4.66117.86±4.88120.04±6.09156.81±10.74185.10±11.510128312.60±47.75116.79±4.97122.79±6.36124.28±6.58163.01±10.38193.44±12.3411349356.16±44.44118.24±4.68124.62±5.68129.79±8.15168.73±9.08198.63±11.3712775401.78±56.33121.36±4.63128.12±4.92133.39±6.99174.31±9.24204.48±12.3513692441.88±53.36123.39±4.29130.20±4.59136.19±8.14180.67±7.75212.82±9.4614545475.70±47.24124.97±4.40131.63±4.88141.69±7.45184.44±7.84216.24±8.8915657485.56±54.62126.51±4.09132.89±4.99144.42±8.67187.38±8.91219.09±9.9616664494.86±54.10127.07±4.77133.55±5.68144.79±8.02189.30±7.08219.60±17.5017857500.91±56.22128.27±4.39135.10±4.14145.28±7.11190.98±7.44219.65±8.6018384518.36±64.09128.97±4.30135.70±2.91146.86±9.22193.44±10.93221.23±9.0319223536.90±41.77131.25±4.07135.71±2.98147.92±5.60195.83±8.31223.42±8.1820189557.43±41.36132.20±2.17140.20±3.49147.97±6.79199.00±5.62228.67±3.78

表3Gompertz、Bertallanff、Logistic和Brody模型各指標(biāo)的相關(guān)估計值

Table3EstimatedvaluesofGompertz,Bertallanff,LogisticandBrodymodelforbodymassindexes

模型Model參數(shù)Parameter體重/kgBodyweight體高/cmWithersheight體斜長/cmBodylength十字部高/cmHipheight胸圍/cmChestcircumference腹圍/cmAbdominalgirthGompertzA764.00131.20147.60139.80202.20227.40B2.720.17420.72139.800.921.04k0.120.16520.180.160.160.20R20.98760.99510.99580.99830.99690.9969LogisticA617.20130.50145.40138.00196.60222.00B8.760.56740.99580.731.391.68k0.230.17750.230.190.210.27R20.98740.99440.99550.99820.99670.9966BertallanffyA831.40129.00148.60140.60205.000.27B0.630.73350.220.170.270.30k0.080.21450.170.150.140.18R20.98760.99180.99590.99830.99700.9970BrodyA651.20132.70151.20142.50213.10238.30B0.530.44540.530.450.630.67k0.140.14070.140.130.110.13R20.98740.99660.99600.99830.99720.9972

表4Gompertz、Bertallanff和Logistic模型對體重指標(biāo)的相關(guān)估計值

Table4EstimatedvaluesofGompertz,BertallanffandLogisticmodelsforbodyweightindexes

生長模型Growthmodel拐點月齡Inflectionpointmonth拐點體重/kgInflectionpointbodyweight最大日增重/kgMaximumdailygain相對生長率RelativegrowthrateGompertz8.47281.811.0510.1251Bertallanff7.15268.621.0070.1328Logistic9.46308.551.1400.1235

2.3　利用4種模型繪制西門塔爾肉牛體尺和體重的生長曲線

4種數(shù)學(xué)模型對肉用西門塔爾牛體重、體尺的擬合結(jié)果如圖1所示，出生到20月齡的肉用西門塔爾牛生長階段4種模型理論值呈逐漸穩(wěn)步上升的趨勢。犢牛從出生到12月齡上升較快之后則逐漸平緩，且4種模型對不同體尺性狀的擬合效果較一致，導(dǎo)致圖中出現(xiàn)了4種曲線重疊的現(xiàn)象。而Gompertz和Bertallanffy兩種模型對體重的擬合優(yōu)于其他兩種模型，這說明，對于該地區(qū)的肉用西門塔爾牛來說，這兩個模型更能反映出本群體的生長發(fā)育規(guī)律。而且，Gompertz模型得到的最大體重可達到831.40 kg，是4種模型中最佳值，因此，對烏拉蓋地區(qū)肉用西門塔爾牛體重擬合的最佳模型是Gompertz。Brody模型對體高、體斜長、胸圍、腹圍和十字部高的擬合效果均最優(yōu)，結(jié)果分別為132.7、151.20、213.10、238.30、217.65和142.50 cm，較其他模型高出0.3～7.1 cm。

圖中GOMP_M、Logi_M、Bert_M、Brody_M 分別代表模型Gompertz、Logistic、Bertallanff和Brody, BW、BL、CC、WH、HH、AG分別代表體重、體斜長、胸圍、體高、十字部高和腹圍GOMP_M, Logi_M, Bert_M, Brody_M represent models Gompertz, Logistic, Bertallanff and Brody, respectively, the BW,BL,CC,WH,HH,AG represent body weight, body length, chest circumference, withers height, hip height, abdominal grith, respectively圖1　利用4種模型繪制西門塔爾肉牛體尺和體重的生長曲線Fig.1　Using the 4 models to draw the growth curves of the body measurements and body weight of Simmental beef cattle

2.4　肉用西門塔爾牛體重與體尺的相關(guān)性分析

由表5可知，右上側(cè)為各性狀間的相關(guān)系數(shù)，左下側(cè)是各性狀相關(guān)關(guān)系顯著性的檢驗結(jié)果。體重與體高、體斜長、十字部高、胸圍和腹圍的相關(guān)性較高，其中體重與胸圍的相關(guān)性最高(r=0.958 55)，腹圍次之；體高與體斜長相關(guān)性最高(r=0.983 87)，胸圍次之。因此，肉用西門塔爾牛6個性狀指標(biāo)間呈極顯著相關(guān)，其中胸圍和腹圍是可直接用來反映體重大小的最重要的兩個指標(biāo)。

2.5　肉用西門塔爾牛體尺與體重的線性回歸分析

2.5.1肉用西門塔爾牛體尺與體重的線性回歸分析通過SAS軟件分析表明，剔除與體重回歸關(guān)系不顯著的體高和十字部高，體斜長、胸圍和腹圍與體重的回歸關(guān)系極顯著(P<0.000 1)，檢驗結(jié)果見表6和表7。建立回歸方程：Y=-366.485 70+1.307 37X3+3.896 12X4-0.417 50X5。

表5肉用西門塔爾牛各性狀指標(biāo)間的相關(guān)性

Table5ThecorrelationanalysisamongdifferenttraitsofSimmentalbeefcattle

性狀Trait體重Bodyweight體高Withersheight體斜長Bodylength十字部高Hipheight胸圍Chestcircumference腹圍Abdominalgirth體重Bodyweight1.0000.939180.940950.939120.958550.94539體高Withersheight<0.00011.0000.983870.960730.972030.96604體斜長Bodylength<0.0001<0.00011.0000.961150.969010.96373十字部高Hipheight<0.0001<0.0001<0.00011.0000.963680.96021胸圍Chestcircumference<0.0001<0.0001<0.0001<0.00011.0000.98538腹圍Abdominalgirth<0.0001<0.0001<0.0001<0.0001<0.00011.000

表6肉用西門塔爾牛體斜長、胸圍、腹圍與體重回歸關(guān)系的顯著性檢驗

Table6Thesignificancetestoftheregressionrelationshipbetweenthebodylength,chestcircumference,abdominalgirthandbodyweightofSimmentalbeefcattle

源Source自由度Degreeoffreedom(df)平方和Sumofsquares均方MeansquareF值F-valuePr>F模型Model4462920911573021608.76<0.0001誤差Error541389182719.37515校正合計Total5455018391

表7肉用西門塔爾牛體斜長、胸圍、腹圍與體重的回歸系數(shù)

Table7Theregressioncoefficientsbetweenthebodylength,chestcircumference,abdominalgirthandbodyweightofSimmentalbeefcattle

變量Variable參數(shù)估計Parameterestimation標(biāo)準(zhǔn)誤差StandarderrorII型SSIIstyleSSF值F-valuePr>FIntercept-366.485706.8980960142292822.65<0.0001X31.307370.1772411592954.41<0.0001X43.896120.21357709072332.79<0.0001X5-0.417500.15892147066.09<0.0001

2.5.2肉用西門塔爾牛體重回歸預(yù)測方程的驗證為了能夠進一步驗證由以上分析所得出的肉用西門塔爾牛體重回歸預(yù)測方程的實用性，選取該地區(qū)12月齡肉用西門塔爾牛共計125頭個體數(shù)據(jù)代入所得回歸預(yù)測方程計算出回歸預(yù)測體重。通過SAS 9.4軟件的配對樣本t檢驗當(dāng)μ=0 時用t分布函數(shù)95%區(qū)間內(nèi)來判斷肉用西門塔爾牛實測體重與回歸預(yù)測體重之間的差異性，兩組體重數(shù)據(jù)間的相關(guān)系數(shù)為0.927 4(P<0.001),檢驗結(jié)果如表8所示，實測體重與回歸預(yù)測體重之間無顯著性差異(P=0.095),證實了烏拉蓋地區(qū)肉用西門塔爾牛體重回歸預(yù)測方程Y=-366.48 570+1.307 37X3+3.896 12X4-0.417 50X5具有統(tǒng)計學(xué)意義。

表8肉用西門塔爾牛實測體重與回歸方程預(yù)測體重樣本配對t檢驗

Table8SamplepairingttestofthemeasuredbodyweightandpredictedbodyweightbyregressionequationofSimmentalbeefcattle

性狀Trait平均值Mean標(biāo)準(zhǔn)差SD標(biāo)準(zhǔn)誤SE上限Lower下限Uppert自由度Df雙側(cè)P值Pvaluesonbothsides實際體重Actualbodyweight401.2547.792.303.44-7.52-0.741250.095預(yù)測體重Predictedbodyweight402.2130.282.77

3　討　論

3.1　肉用西門塔爾牛體重生長曲線模型的討論

前人通過研究比較不同模型來分析母牛的生長發(fā)育[23],一致認(rèn)為用Brody模型描述母牛生長發(fā)育過程效果較差，這些牛品種包括安西母牛、延邊母牛、秦川母牛、徽成母牛、安格斯母牛、海福特母牛和娟姍母牛。而本研究選取烏拉蓋地區(qū)肉用西門塔爾牛從出生到20月齡的體重和體尺測量數(shù)據(jù)進行分析發(fā)現(xiàn)，4種生長曲線模型對體重和體尺性狀都能夠得到很好的擬合。動物生長曲線模型實質(zhì)上是把動物生長過程中各個階段的體重和體尺的數(shù)據(jù)資料提供的信息綜合成少數(shù)幾個可表現(xiàn)具體生物學(xué)意義的參數(shù)[24]。動物生長曲線的分析和擬合也是研究畜禽生長發(fā)育規(guī)律的主要方法之一，尤其是國外許多知名牛品種已通過多種數(shù)學(xué)模型建立了適合該品種的生長曲線方程[13,18]，為個體選育提供數(shù)據(jù)支持。熊飛等[25]通過對檳榔江水牛的體重和體尺性狀進行累計生長曲線的預(yù)測方程擬合發(fā)現(xiàn)，檳榔江水牛在犢牛時期的生長發(fā)育較快，其新陳代謝能力相對旺盛，而成年后則逐步趨于穩(wěn)定。本研究中所應(yīng)用的4種非線性模型中，Brody模型僅僅是描述漸近生長過程；Logistic、Gompertz是拐點固定的S形生長曲線模型[26]；Bertalantfy是拐點可變的S形生長曲線模型。Logistic模型拐點在極限體重的1/2處拐點體重為0.5A，該模型適用于早期生長慢、生長拐點較遲的牛品種。Gompertz模型拐點在極限體重的1/e處，拐點體重為0.368A，所以Gompertz模型適合描述早期生長迅速的生長過程，也是當(dāng)前使用較多的生長曲線模型。對于本試驗群體，由于早期生長迅速，Gompertz模型擬合結(jié)果達到0.987 6，是4種模型中最好的，能很好的反映該品種的生長特性從而指導(dǎo)選育。

本研究所選試驗群體數(shù)據(jù)是來自內(nèi)蒙古烏拉蓋地區(qū)肉用西門塔爾牛不同年份各生長階段的體重體尺記錄，烏拉蓋地區(qū)屬森林草原向典型草原過渡地帶，以典型草原為主，氣候?qū)侔霛駶櫋敫珊荡箨懶詺夂颍瑲鉁氐停浼据^長，年均氣溫為-0.9 ℃，水資源豐富。本地肉用西門塔爾牛一般以放養(yǎng)為主，自由采食，生長速度較快，體型大，是優(yōu)良的牛品種。通過對體重體尺生長曲線模型的分析發(fā)現(xiàn)，犢牛的體尺和體重在出生至1周歲上升速度較快。烏拉蓋地區(qū)肉用西門塔爾牛的成熟體重為617.2～831.4 kg，這與陳江凌等[27]對荷斯坦、西門塔爾、夏洛萊、利木贊和紅安格斯牛用Bert模型得到的成熟體重793.7～1 081.2 kg存在差距，主要是因為所測牛只進行3月齡早期斷奶及6月齡分群飼養(yǎng)，由于受飼養(yǎng)、管理、天氣、飼喂條件等因素制約，致使生長過程出現(xiàn)較大偏差。根據(jù)張麗[18]和左福元[15]等分別對南陽母牛和金黃阿奎頓雜交后代牛與黑安格斯雜交牛體重生長拐點的研究發(fā)現(xiàn)，南陽母牛(4.41月齡，113.77 kg)、金黃阿奎頓雜交牛后代(4.16月齡，71.83 kg)和黑安格斯雜交牛(4.97月齡，78.38 kg)生長拐點都出現(xiàn)在4月齡，而德國對西門塔爾雜交后代牛體重生長拐點的研究發(fā)現(xiàn)，拐點月齡為9.58月齡時體重是115.23 kg,本研究的肉用西門塔爾牛群體利用Gompertz 模型得到拐點月齡為8.47，此時的拐點體重為281.81 kg, Logistic模型則分別為9.46月齡和308.55 kg，是前人所報道的拐點體重的2倍。品種早熟性決定牛體重拐點時間，拐點時間越晚則表示該品種成熟越晚，體型更大。

3.2　肉用西門塔爾牛體尺與體重相關(guān)性和回歸性的分析

畜禽體型外貌是遺傳因素與環(huán)境因素相互作用的結(jié)果。與其他牛品種相比，新疆褐牛的體型外貌、尻部、后軀對品種特征、生產(chǎn)性能方面的影響非常大[28-29]。而肉用西門塔爾牛后備牛以達到理想的生長指標(biāo)且14～15月齡可以配種(體高127 cm、體重360 kg)為主要目標(biāo)[30-31]。肉用西門塔爾牛的體重和體尺是反映該群體身體結(jié)構(gòu)及發(fā)育程度的重要指標(biāo)，是用來衡量其經(jīng)濟用途最直接的依據(jù)之一。體尺測量是采用專門的工具對家畜各部位進行準(zhǔn)確度量，研究本品種外貌特征、估計體重及生產(chǎn)性能的重要手段之一，不同月齡肉用西門塔爾牛體尺生長發(fā)育速度不同。通過對肉用西門塔爾牛進行相關(guān)性分析可知，體尺性狀與體重之間均呈極顯著正相關(guān)，各個體尺性狀通過直接或間接的作用對體重產(chǎn)生影響。體尺與體重的相關(guān)性系數(shù)大小依次為胸圍>腹圍>體斜長>體高>十字部高。由此可以發(fā)現(xiàn)，胸圍和腹圍與體重的相關(guān)性系數(shù)分別達到0.958 55和0.945 39，是對體重產(chǎn)生影響最大的指標(biāo)。在生產(chǎn)實踐中應(yīng)加以重視與體重相關(guān)性高的指標(biāo)，提高經(jīng)濟效益。

多元線性分析是統(tǒng)計分析及建立預(yù)測方程的重要手段。通過對烏拉蓋地區(qū)肉用西門塔爾牛體高、十字部高、體斜長、胸圍、腹圍與體重的線性回歸分析,利用向前回歸剔除顯著性不明顯的X1、X2和X3后得到體斜長、胸圍、腹圍與體重新的預(yù)測方程：Y=-366.485 70+1.307 37X3+3.896 12X4-0.417 50X5，并一步驗證了該方程的實用性，預(yù)測體重與實際體重間無顯著差異，該擬合方程有較好的估測能力。通過該研究可以看出，不同年齡的肉用西門塔爾牛在所研究階段體重與體尺性狀均呈一定程度的波動狀態(tài)。

4　結(jié)　論

本研究通過對內(nèi)蒙古烏拉蓋地區(qū)肉用西門塔爾牛群體的體重、體高、體斜長、十字部高、胸圍和腹圍等指標(biāo)數(shù)據(jù)的生長發(fā)育規(guī)律進行研究和分析發(fā)現(xiàn)，Gompertz和Bertallanffy兩種模型對體重的擬合優(yōu)于其他兩種模型，Brody模型對體高、體斜長、胸圍、腹圍和十字部高的擬合效果均最優(yōu)。通過SAS 9.4軟件的配對樣本t檢驗，當(dāng)μ=0的t分布函數(shù)用95%區(qū)間內(nèi)來判斷肉用西門塔爾牛實測體重與回歸預(yù)測體重之間的差異性，兩組體重數(shù)據(jù)之間的相關(guān)系數(shù)為0.927(P<0.001)，實測體重與回歸預(yù)測體重之間無顯著性差異(P<0.001)，預(yù)測體重回歸方程：Y=-366.485 70+1.307 37X3+3.896 12X4-0.417 50X5，可應(yīng)用于估測肉用西門塔爾牛體重的生產(chǎn)實踐中，為烏拉蓋地區(qū)肉用西門塔爾牛的選育提供相應(yīng)的數(shù)據(jù)支持。

參考文獻(References)：

[1]王雅春, 陳幼春. 世界西門塔爾牛育種現(xiàn)狀及其在中國的應(yīng)用方向[J].中國牛業(yè)科學(xué), 2008, 34(6): 1-4.

WANG Y C, CHEN Y C. Breeding progress of simmental World-wide and the Principle of its application in China[J].ChinaIndustryScience, 2008, 34(6): 1-4. (in Chinese)

[2]許尚忠, 李俊雅, 任紅艷, 等. 中國西門塔爾牛選育及其進展[J].中國畜禽種業(yè), 2008, 4(5): 13-15.

XU S Z, LI J Y, REN H Y, et al. Chinese simmental breeding and its progress[J].TheChineseLivestockandPoultryBreeding, 2008, 4(5): 13-15. (in Chinese)

[3]方曉敏, 許尚忠, 張英漢, 等. 我國新的牛種資源——中國西門塔爾牛[J].黃牛雜志, 2002, 28(5): 67-69.

FANG X M, XU S Z, ZHANG Y H, et al. The new kind of resources-Chinese simmental cattle[J].JournalofYellowCattleScience, 2002, 28(5): 67-69. (in Chinese)

[4]FAN H Z, WU Y, ZHOU X J, et al. Pathway-based genome-wide association studies for two meat production traits in simmental cattle[J].SciRep, 2015, 5: 18389.

[5]WU Y, FAN H Z, WANG Y H, et al. Genome-wide association studies using haplotypes and individual SNPs in Simmental cattle[J].PLoSOne, 2014, 9(10): e109330.

[6]NIU H, ZHU B, GUO P, et al. Estimation of linkage disequilibrium levels and haplotype block structure in Chinese Simmental and Wagyu beef cattle using high-density genotypes[J].LivestSci, 2016, 190: 1-9.

[7]XIA J W, QI X, WU Y, et al. Genome-wide association study identifies loci and candidate genes for meat quality traits in Simmental beef cattle[J].MammGenome, 2016, 27(5-6): 246-255.

[8]ZHU B, NIU H, ZHANG W G, et al. Genome wide association study and genomic prediction for fatty acid composition in Chinese Simmental beef cattle using high density SNP array[J].BMCGenomics,2017, 18: 464.

[9]ZHU B, ZHU M, JIANG J C, et al. The impact of variable degrees of freedom and scale parameters in bayesian methods for genomic prediction in Chinese simmental beef cattle[J].PLoSOne, 2016, 11(5): e0154118.

[10]牛紅, 寶金山, 吳洋, 等. 中國西門塔爾牛肉用群體重要經(jīng)濟性狀遺傳參數(shù)估計[J].畜牧獸醫(yī)學(xué)報, 2016, 47(9): 1817-1823.

NIU H, BAO J S, WU Y, et al. Estimation of genetic parameters for economic important traits in Chinese Simmental beef cattle[J].ActaVeterinariaetZootechnicaSinica, 2016, 47(9): 1817-1823. (in Chinese)

[11]楊運清, 繆堯源. 動物生長曲線擬合方法的探討[J].畜牧獸醫(yī)學(xué)報, 1992, 23(3): 219-224.

YANG Y Q, MIAO Y Y. Animal growth curve fitting method discussion[J].ActaVeterinariaetZootechnicaSinica,1992, 23(3): 219-224. (in Chinese)

[12]楊運清. 動物生長模型的優(yōu)化擬合方法[J].數(shù)理統(tǒng)計與管理, 1991, 10(5): 43-48.

YANG Y Q. Animal growth model of the optimum fitting method[J].ApplicationofStatisticsandManagement,1991, 10(5): 43-48. (in Chinese)

[14]BULLOCK K D, BERTRAND J K, BENYSHEK L L. Genetic and environmental parameters for mature weight and other growth measures in Polled Hereford cattle[J].JAnimSci, 1993, 71(7): 1737-1741.

[15]左福元, 王玲, 周沛, 等. 三峽庫區(qū)川南山地黃牛及雜種后代生長曲線擬合分析[J].中國畜牧雜志, 2010, 46(9): 1-4.

ZUO F Y, WANG L, ZHOU P, et al. Analysis of growth curve of Chuannan Mountainous cattles and their hybrids in Three Gorges Reservoir Area[J].ChineseJournalofAnimalScience,2010, 46(9): 1-4. (in Chinese)

[16]楊曉冰, 陳宏, 滑留帥, 等. 固原地區(qū)秦川牛利雜群體生長曲線擬合分析[J].中國牛業(yè)科學(xué), 2011, 37(6): 4-7.

YANG X B, CHEN H, HUA L S, et al. Growth curve analysis of Limousin×Qinchuan cattle in Guyuan Region[J].ChinaCattleScience,2011, 37(6): 4-7. (in Chinese)

[17]楊敏, 朱化彬, 殷宗俊, 等. 皮南肉牛選育初期的生長曲線分析[C]//第三屆中國牛業(yè)發(fā)展大會論文集. 穆棱: 中國畜牧業(yè)協(xié)會, 2008.

YANG M, ZHU H B, YIN Z J, et al. Skin south beef cattle breeding in the early growth curve analysis[C]//The Third China Industry Development Conference. Muling: China Animal Agriculture Association, 2008. (in Chinese)

[18]張麗, 張良志, 張愛玲. 南陽黃牛體重曲線擬合的比較分析[J].廣東海洋大學(xué)學(xué)報, 2008, 28(1): 98-100.

ZANG L, ZHANG L Z, ZHANG A L, et al. The comparison of Nanyang cattle weight curve fitting analysi[J].JournalofGuangdongOceanUniversity, 2008, 28(1): 98-100. (in Chinese)

[19]馬光輝, 張文龍, 黃錫霞, 等. 新疆地區(qū)西門塔爾牛不同模型泌乳曲線擬合效果分析[J].新疆農(nóng)業(yè)科學(xué), 2013, 50(10): 1931-1936, 1930.

MA G H, ZHANG W L, HUANG X X, et al. Analysis of fitting effect on different models of lactation curve for Simmental Cattle in Xinjiang[J].XinjiangAgriculturalSciences, 2013, 50(10): 1931-1936, 1930. (in Chinese)

[20]楊運清, 楊玉林. 動物生長模型參數(shù)的常規(guī)估計法[J].東北農(nóng)業(yè)大學(xué)學(xué)報, 1994, 25(2): 132-137.

YANG Y Q, YANG Y L. Routine parameter estimation methods for animal growth models[J].JournalofNortheastAgriculturalUniversity, 1994, 25(2): 132-137. (in Chinese)

[21]胡良平. 《Windows SAS6.12 & 8.0實用統(tǒng)計分析教程》出版[J].中華醫(yī)學(xué)雜志, 2001, 81(8): 475.

HU L P. 《Windows SAS6.12 & 8.0 Practical statistical analysis tutorial》published[J].NationalMedicalJournalofChina, 2001, 81(8): 475. (in Chinese)

[22]于志, 高艷霞, 曹玉鳳, 等. 中國荷斯坦?fàn)倥：陀膳ＩL發(fā)育規(guī)律的研究[J].中國畜牧獸醫(yī), 2014, 41(6): 121-125.

YU Z, GAO Y X, CAO Y F, et al. Study on growth and development pattern of Chinese Holstein Calf and Heifer[J].ChinaAnimalHusbandry&VeterinaryMedicine, 2014, 41(6): 121-125. (in Chinese)

[23]張文舉. 徽成盆地黃牛生長曲線方程的研究[J].甘肅畜牧獸醫(yī), 1998, 28(1): 3-6.

ZHANG W J. Cattle growth curve equation of the badge into basin research[J].GansuAnimalandVeterinarySciences,1998, 28(1): 3-6. (in Chinese)

[24]LARID A K. Postnatal growth of birds and mammals[J].Growth, 1966, 30(3): 349-363.

[25]熊飛, 苗永旺, 李大林, 等. 檳榔江水牛體重及體尺生長規(guī)律的研究[J].云南農(nóng)業(yè)大學(xué)學(xué)報, 2011, 26(4):472-478.

XIONG F, MIAO Y W, LI D L, et al. Growth regularities of body weight and size of Binglangjiang buf-falo[J].JournalofYunnanAgriculturalUniversity,2011, 26(4): 472-478. (in Chinese)

[26]DARMANI KUHI H, KEBREAB E, LOPEZ S, et al. An evaluation of different growth functions for describing the profile of live weight with time (age) in meat and egg strains of chicken[J].PoultSci, 2003, 82(10): 1536-1543.

[27]陳江凌, 汪聰勇, 耿繁軍. 種公牛早期生長曲線分析[J].安徽農(nóng)業(yè)科學(xué), 2014, 42(14): 4352-4353.

CHEN J L, WANG C Y, GENG F J. Analysis of growth curve of growing bulls[J].JournalofAnhuiAgriculturalSciences, 2014, 42(14): 4352-4353. (in Chinese)

[28]姚亞鈴. 湘西黃牛生長發(fā)育規(guī)律及雜交利用效果研究[D]. 長沙: 湖南農(nóng)業(yè)大學(xué), 2006.

YAO Y L. Study on Xiangxi Cattle growth developmental mechanism and utilization of hybrid effect[D]. Changsha: Hunan Agricultural University, 2006. (in Chinese)

[29]馬茂華, 孔令韶. 新疆呼圖壁種牛場地區(qū)13種優(yōu)勢植物水溶性鹽分的含量特點及數(shù)量分析[J].生態(tài)學(xué)報, 1996, 16(5): 509-516.

MA M H, KONG L S. Quantitative analysis and characteristics of water-soluble salt contents in 13 dominant plant species in the area of Hutubi cattle farm, Xinjiang[J].ActaEcologicaSinica, 1996, 16(5): 509-516. (in Chinese)

[30]張世元, 張偉. 西門塔爾后備牛的培育[J].新疆畜牧業(yè), 2011(4): 32-35.

ZHANG S Y, ZHANG W. Simon tal replacement cattle breeding[J].JournalofAnimalHusbandryinXinjiang, 2011(4): 32-35. (in Chinese)

[31]付雪峰, 季生林, 黃錫霞, 等. 舍飼新疆褐牛犢牛體重與體尺指標(biāo)的相關(guān)及回歸分析[J].新疆農(nóng)業(yè)科學(xué), 2009, 46(4): 844-848.

FU X F, JI S L, HUANG X X, et al. Correlation and regression analysis of body weight and body size measurements for Xinjiang Brown cattle calf under pen feeding[J].XinjiangAgriculturalSciences, 2009, 46(4): 844-848. (in Chinese)