劉更壽
(青海省大通種牛場(chǎng),青海 西寧 810102)
?
2.5歲大通牦牛母牛體重與體尺指標(biāo)的相關(guān)回歸分析
劉更壽
(青海省大通種牛場(chǎng),青海 西寧 810102)
為了分析2.5歲大通牦牛母牛體尺與體重的相關(guān)性,以便在實(shí)際工作中應(yīng)用,在11月份對(duì)31頭2.5歲大通牦牛母牛進(jìn)行稱重和體尺測(cè)量,根據(jù)得出的數(shù)據(jù)分析其體重與體高、體斜長、胸圍、管圍的相關(guān)系數(shù),同時(shí)估測(cè)大通牦牛2.5歲母牛的體重回歸模型。結(jié)果顯示,大通牦牛2.5歲母牛體重與體高、體斜長、胸圍、管圍的相關(guān)系數(shù)分別為0.309、0.713、0.846和0.595,經(jīng)檢驗(yàn)體高、管圍與體重的相關(guān)系數(shù)達(dá)到了顯著水平(p<0.05),體斜長、胸圍與體重的相關(guān)系數(shù)達(dá)到了極顯著水平(p<0.01),得到了3個(gè)估測(cè)體重的回歸模型,估測(cè)值與實(shí)測(cè)值之間相關(guān)程度分別為0.923、0.864和0.893。大通牦牛2.5歲母牛體重與體高、體斜長、胸圍、管圍之間顯著相關(guān),說明體重與體尺之間存在著明顯的線性關(guān)系。
大通牦牛;母牛;體重;體尺;相關(guān)性
大通牦牛是中國農(nóng)科院蘭州畜牧與獸藥研究所和青海省大通種牛場(chǎng)連續(xù)20年執(zhí)行農(nóng)業(yè)部“六五”、“七五”、“八五”、“九五”重點(diǎn)項(xiàng)目培育成功的牦牛新品種[1]。其肉用性能好,遺傳性較穩(wěn)定,是牦牛的理想遺傳資源[2]。在其它牛及動(dòng)物上對(duì)體重和體尺研究方面的報(bào)道較多[3、4],但是對(duì)大通牦牛體尺與體重相關(guān)性分析還未見報(bào)道,然而研究牦牛體重與體尺的相關(guān)性有較強(qiáng)實(shí)踐意義。在以往的畜牧獸醫(yī)工作中,無論是品種資源調(diào)查、選種選配,或者是計(jì)算日糧、考慮用藥量等等,都需要了解牛的體重。實(shí)稱牛的體重一般不易做到,尤其是在牧區(qū)更難做到,根據(jù)牛的體尺來估算牛的體重確實(shí)是可靠易行的簡(jiǎn)便方法。在實(shí)際工作中發(fā)現(xiàn),如把每個(gè)年齡段的樣本合在一起分析,體重與體尺的相關(guān)系數(shù)差異性不顯著,說明不同的年齡和性別體重與體尺的相關(guān)系數(shù)也有差異,為此按年齡和性別分別進(jìn)行分析,旨在為以后的實(shí)際工作中估測(cè)體重把誤差控制在最小的范圍內(nèi)。
1.1材料來源
試驗(yàn)牦牛是青海省大通種牛場(chǎng)2011年11月份測(cè)定31頭大通牦牛2.5歲母牛。體重(kg)是空腹進(jìn)行稱重,為了便于描述,將體重、體斜長、體高、胸圍、管圍5個(gè)指標(biāo)用下列大寫字母代替,分別為TZ、TXC、TG、XW、GW。
1.2統(tǒng)計(jì)方法
所得數(shù)據(jù)采用Spss17.0程序Pearson和Linear過程[5]進(jìn)行處理。
1.2.1Pearson相關(guān)系數(shù)分析計(jì)算模型
式中,R為Pearson相關(guān)系數(shù),數(shù)值介于-1~1之間。當(dāng)R值為正數(shù)是為正相關(guān),表示依變量隨自變量的增大而增大;當(dāng)R值為負(fù)數(shù)時(shí)為負(fù)相關(guān),表示依變量隨自變量的增大而減?。划?dāng)R值等于0時(shí)表示依變量與自變量之間沒有相關(guān)性。X為自變量,Y為依變量。
1.2.2多元線性回歸分析計(jì)算模型
Y=b0+b1X1+b2X2+…+bn+Xn,其中Y為依變量,b0為常數(shù)、b1,b2+…bn,為回歸系數(shù),X1,X2,…Xn,為回歸系數(shù)對(duì)應(yīng)的自變量。
2.1相關(guān)分析
設(shè)TXC、TG、XW、GW的值作為自變量,TZ的值作為依變量。
2.2回歸分析
分別采用Linear過程“Enter”法和“Stepwise”建立多元線性回歸方程,統(tǒng)計(jì)結(jié)果見表2和表3。
表1相關(guān)系數(shù)統(tǒng)計(jì)結(jié)果
樣本含量體重相關(guān)系數(shù)TGTXCXWGW30RTZ0.3090.7130.8640.595
注:經(jīng)檢驗(yàn),體高、管圍與體重的相關(guān)系數(shù)達(dá)到顯著水平(p<0.05),體斜長、胸圍與體重的相關(guān)系數(shù)達(dá)到極顯著水平(p<0.01)。
表2回歸模型系數(shù)
模型模型分組非標(biāo)準(zhǔn)化系數(shù)1b0(常數(shù)項(xiàng))-181.595b1(TG回歸系數(shù))-1.046b2(TXC的回歸系數(shù))0.691b3(XW的回歸系數(shù))1.910b4(GW的回歸系數(shù))6.9012b0(常數(shù)項(xiàng))-179.180b1(XW的回歸系數(shù))2.4623b0(常數(shù)項(xiàng))-225.149b1(XW回歸系數(shù))2.132b2(TXC的回歸系數(shù))6.094
表3回歸模型擬合度綜述①
模型類別模型的相關(guān)系數(shù)(R)決定系數(shù)(R2)1②0.9230.8522③0.8640.7463④0.8930.797
從表2可知,這3個(gè)回歸模型分別為:
第一個(gè)回歸模型:TZ=-181.595+ -1.046×TG+0.691×TXC+1.910×XW+6.901×GW
第二個(gè)回歸模型:TZ=-179.180+2.462×XW
第三個(gè)回歸模型:TZ=-225.194+2.132×XW+6.094×TXC
這3個(gè)回歸模型的相關(guān)系數(shù)R分別為0.923、0.846和0.746,擬合度的決定系數(shù)R2值分別為0.852、0.746和0.797,說明線性度較好。
3.1通過分析2.5歲大通牦牛母牛體高、體斜長、胸圍、管圍4個(gè)指標(biāo)與體重之間的相關(guān)系數(shù),分別為0.309、0.713、0.864和0.595,說明大通牦牛2.5歲母牦牛體高、體斜長、胸圍、管圍之間存在顯著正相關(guān)。這與窮達(dá)、薛長安[6,7]的研究結(jié)果一致。
3.2該研究表明得到了3個(gè)估測(cè)大通牦牛2.5歲母牦牛體重的回歸模型。表明大通牦牛2.5歲母牦牛體重與體高、體斜長、胸圍、管圍存在明顯的線性關(guān)系,在實(shí)際應(yīng)用中可根據(jù)各回歸模型的R2值大小以及體尺測(cè)量的繁瑣程度來對(duì)3個(gè)模型進(jìn)行選擇。其中最簡(jiǎn)單的模型是TZ=-179.180+2.462×XW。
3.3大通牦牛2.5歲母牛體重與體尺指標(biāo)的相關(guān)回歸分析模型在實(shí)際生產(chǎn)中具有重要的意義。牦牛是放牧型動(dòng)物,活動(dòng)范圍較為寬廣,對(duì)其稱重較為不便,然而在畜牧獸醫(yī)工作中,無論是品種資源調(diào)查、選種選配,或者是計(jì)算日糧、考慮用藥量等都需要了解牛的體重,測(cè)量牦牛體尺,通過體重與體尺的回歸模型測(cè)算牦牛體重是比較方便可行的方法。
3.4此次的研究結(jié)果表明,2.5歲大通牦牛母牛體重與體尺之間存在明顯的線性關(guān)系,但是根據(jù)得出的兩個(gè)回歸模型對(duì)樣本進(jìn)行估算時(shí),各樣本間誤差的離散程度較大,這說明還有某些因素對(duì)體重與體尺造成一定的影響。筆者認(rèn)為,測(cè)量體尺時(shí)由于牦牛膽小,對(duì)人有過激的反應(yīng),致使牦牛站立的姿勢(shì)不一致,這就對(duì)測(cè)量體高和胸圍會(huì)造成較大的誤差,所以對(duì)牦牛體尺的測(cè)量有待進(jìn)一步的研究,包括測(cè)量的方法和部位,以達(dá)到把系統(tǒng)誤差控制在最小的范圍內(nèi)。
[1]陸仲璘,何曉林,閻萍.世界上第一個(gè)牦牛培育新品種-“大通牦牛”簡(jiǎn)介[P]中國草食動(dòng)物“大通牦?!毙缕贩N培育及其培育技術(shù)論文集2005:12-14.
[2]楊博輝,姚軍,王敏強(qiáng),等.大通牦牛肌肉纖維組織學(xué)特性研究[J]中國草食動(dòng)物,2001(5):34-35.
[3]馬存壽.雜種牛荷斯坦成年母牛體重與體尺指標(biāo)的相關(guān)回歸分析[J]青海畜牧獸醫(yī)雜志,2012(1):28-29.
[4]田亞磊,朱東亮,張聰,等.太行黑山羊體重、體重的相關(guān)性分析[J]山西農(nóng)業(yè)科學(xué),2009(5):31-32.
[5]凌莉.生物統(tǒng)計(jì)學(xué)基礎(chǔ)[M]北京:2010:130-135.人民衛(wèi)生出版社.
[6]窮達(dá),彭措巴姆.西藏嘉黎成年母牦牛體重與體尺指標(biāo)的相關(guān)回歸分析[J].畜牧與飼料科學(xué),2011,32(6):11-12.
2016-06-01
S823.2
A
1005-2739(2016)05-0057-03