鲁西黄牛成年母牛体重与体尺指标的相关回归分析

2 0 0 7 年增干u 第二稻中国牛业发展大会论文集 鲁西黄牛成年母牛体重与体尺指标 的相关回归分析 周正奎1 ，李姣1 ，姬爱国”，黄萌“，许尚忠”，陈金宝1 ，高雪1 ，任缸艳1 ( 1 冲国农业科学院北京畜牧兽医研究所中国农业科学院肉牛研究中- 心，北京1 0 0 0 9 4 2 西北农林科技大学动物科技学院，陕西杨凌7 1 2 1 0 0 ) 摘要：班山东省鲁西黄牛产1 五2 0 0 6 年测量的1 1 0 条相关数据为基础，分析了体重与体高。体长、胸围、管围的 相关系数，构建了估测鲁西黄牛成年母牛体重的回归模型结果表明：体重与体高，体长，胸围管围的相关系 数分别为：0 7 5 1 、0 8 1 4 、09 0 9 、0 8 5 5 ；得到了2 个估测体重的回归模型，估测值与实剥值之间的相关程度分别 为0 9 5 4 和0 9 5 3 ，并根据2 个模型用V i s u a l B a s i c i i - 言蝙写了鲁西黄牛体重计算程序。 关键词：鲁西黄牛；体重；相关系数；回归模型；计算程序 中图分类号：$ 8 2 3文献标识码：A 准确度量牛的体重是制定繁殖育种和选种选配计 划及饲养标准、测定生产性能的重要依据，是研究生 长发育、衡量畜体变化以及评定养牛场技术管理水平 和经营效果的主要指标。在农村和一些牧场由于条件 所限，直接称量牛的体重困难较大。据此，寻求准确 而简易的估测牛体重的方法，在生产和科研中具有重 要的实践意义。 鲁西黄牛是我国中原良种黄牛四大品种之一，是 我国优质肉牛重要种质资源，主要分布于山东省西南 部的菏泽和济宁两地区境内，鲁西黄牛体躯高大、骨 骼细致、背腰平直、肌肉发达、肉质细嫩。耐粗饲， 适应性强。近几年来随着市场需求的扩大，鲁西黄 牛的饲养正向周边延伸。本文旨在分析鲁西黄牛成年 母牛体尺指标与体重的相关性，并进一步进行了估测 体重的回归分析，以期指导本品种牛的生产实际及选 育工作。 1 材料和方法 1 1 材料来源与样本数量 试验数据来自山东省粱 山县科龙畜牧产业有限公司和鄄城县鲁西黄牛原种场 2 0 0 6 年7 月份测量的共计1 1 0 条相关数据。体高( c m ) 为 投稿日期：2 耵7 - 0 6 - 1 4 通讯作者 基金项目：国家科技支撑计划重大项目( N o 嬲B A D 0 1 A 1 0 ) 和国家 高技术研究与发展计划( 8 6 3 ) 项目( N 0 2 0 0 6 A A I O Z l 9 7 ) 瓷f f 0 。 髫甲撮高点到地面的垂直距离：胸围( e r a ) 是肩胛后角 处体躯的垂直周径；体长( c m ) 即为体斜长，肩端至坐 骨端的距离；管围( c I l l ) 是前肢胫部上1 3 处的周径，在 前管的最细处量取。为便于描述方便，将体重、体高、 体长、胸围、管围5 个指标进行简写，符号分别为T Z 、 T G 、T C 、X W 、G W 。 1 ，2 统计方法所有数据采用S A S8 0 程序R e g r e s s i o n L i n e a r 过程进行处理。 1 2 1 相关系数的计算公式 R ： 苎坚二型兰二里 。4 Z ( X - X ) 2 E ( y - 歹) 2 式中，R 为相关系数，x 为自变量，y 为依变量， 相关系数印为两个变量的乘积和除以各自变量的平方 和。R 的数值介于1 到1 之间，当R 为正数时为正相关， 表示依变量随自变量的增大而增大，当R 值为负数时 为负相关，表示依变量随自变量的增大而减小，R 等 于0 时，表示没有相关性。 1 2 2多元线性回归分析计算模型 Y = b o + b l x l + b 2 x 2 + b 。k ，其中，Y 为依变量，b o 为常 数，b l 、b 2 、b n 为回归系数，x l 、1 2 x n 为回归 系数对应的自变量。 1 3 程序的编写环境使用M i c r o s o f tV i s u a lB a s i c6 0 企业版。应用程序向导生成一个具有计算功能的应用 中固耄牧京是留 第二稻中国牛业发展大会论文集 程序，最后编译为可执行程序。 2 结果 2 1 相关性分析设T G 、T C 、X W 、G W 作为自变量。 T Z 作为依变量。 由表1 可以看出，T G 、T C 、x W 、G W 与T Z 之间 的相关系数分别为0 7 5 1 、O 8 1 4 、0 9 0 9 、0 8 5 5 。 裹1 相关系数统计结果 样本含量体重相关系数T G T CX W G W 1 1 0 R 0 7 5 1n 8 t 4 0 9 0 9 0 ，8 5 5 2 2 回归分析分别采用R e g r e s s i o nL i n e a t 过程“F u l l m o d e l ”法和“F o r w a r ds e l e c t i o n ”法建立多元线性回 归方程，统计结果见表2 、表3 。 襄2 回归模型系数表 由表2 可知，这2 个回归模型分别为： 回归模型1 ： Y ( T Z k 培0 2 0 0 32 7 + o 6 0 42 8T G + 1 5 5 46 2T C + 3 4 9 29 7X W + 1 52 2 52 8G W 回归模型2 ； Y ( 1 Z 净7 7 8 0 9 3 4 64 - 】6 7 1T C4 - 36 5 68 0 X W + 1 6 3 8 6 5 7G W 表3 回归模型拟合度综述衰 从表3 可以看到这两个回归模型拟合度的决定系 数的值分别为0 9 1 01 和0 9 0 85 ，说明线性度很好。 2 3 计算程序图1 为鲁西黄牛成年母牛体重计算程 序界面，界面包括体尺数据输入框，模型类别选择框， 分别输入测得的数据并选择其中的模型后点击计算， 即可在体重输出窗口显示计算结果。 固中固毫教鸯志 图I 计算程序界面 2 0 0 7 年增刊 3 讨论与结论 3 1 本文共得到2 个估测鲁西黄牛成年母牛体重的回 归模型。表明了鲁西黄牛成年母牛体重与体高、体长、 胸围、管围存在明显的线性关系，此结论与国内外相 关研究相吻合。判断回归模型好坏的标准是多方面的， 就数量统计角度而言，根据回归模型来计算拟台值与 实测值之间的决定系数( R 2 ) ，通过比较R 2 值的大小以 择其优劣更为适宣。由表3 可以看出，第一个回归模型 的R 2 值大于第二个回归模型的R 2 值，因此以选择第一 个回归模型更为适宜，但如果同时考虑体尺测量的繁 琐程度可以考虑第二个回归模型。 3 2 在实测体重与估测体重数据比较中我们发现4 1 0 蚝 体重是一个界值。实测体重小于4 1 0 虹的成牛母牛， 两个模型的体重估测值有略高于实测值的趋势，而大 于4 1 0k g 的成年母牛，估测值则有略低于实涣4 值的趋 势。这可能与鲁西黄牛的饲养水平和营养因素有关。 3 3 本文根据两个回归模型用V i s u a lB a s i c6 0 软件编写 了鲁西黄牛成年母年的体重计算程序。此程序为免安装 的可执行程序，具有操作简单界面友好，实用性强， 替代了人工繁琐的计算过程提高计算的精确性。 参考文献： 【l 】王占红何永涛郭维春，等，沿江牛成年母牛体重与体尺指 标的相关与回归分析啊中国畜牧杂志，2 0 0 7 4 3 ( I ) ：1 2 - 1 3 【2 l 杨公社，毛玉眭用二元回归方程式估铡寨川牛成年母牛的体 重【J 1 黄牛杂志，1 9 8 3 ，3 ：3 9 - 4 7 跚武彬，庞之洪陈幼春，等黄牛体尺差异性分析方法的研究 【刀畜牧兽医学报，1 9 9 0 ，2 1 ( 2 ) ：1 2 1 - 1 2 6 【4 】张显华，王戚，王振民黄牛体尺的聚类与主成分分析【J 1 黄 牛杂志，1 9 9 7 ，2 3 ( 2 ) ：9 - 1 0 ， 2 0 0 7 年增刊 【5 】张爱玲，张丽娟，耿社民，等秦川母牛不同年龄阶段体尺和 体重的主成分分析丁西北农林科技大学学报( 自然科学版) ， 2 0 0 3 3 1 ( 2 1 ：2 9 - 3 2 f 6 】拱楠侯军S A Sf o rW i n d o w s ( v 8 ) 统计分折系统教程新编 第二届中鹭牛业发展大会论文集 嗍北京：清华大学出版社 2 4 0 2 9 4 口】王根林，易建明t 粱学武，等 出版社，2 0 0 0 ，7 1 7 2 ， 北京交通大学出版社2 0 0 4 养牛学【h q 北京：中国农业 R e g r e s s i o n A n a l y s i s o f B o d y W e i g h t a n d B o d y M e a s u r e m e n t s o f A d u l t C o w so f L u x i Y e l l o w C a t t l e Z H O U Z h e n g K n i l ，L I J i a 0 1 ，J I A i G u o ”，H U A N G M e n g “，X U S h a n g z b o n g ” C I - m NJ i n b a 0 1 ，G A OX u e l ，R E NH o n g - y a n l ( 1 1 a s t i t u t co f A n i m a lS c i e n c e ，C h i n e s eA c a d e m yo f A g r i c u l t u r a lS c i e n c e s ，B e e fC a t t l eR e s e a r c hC e n t e r , C A A S 。 B e i j i n g1 0 0 0 9 4 ，C h i n a ； 2 C o H e g eo f A n i m a lS c i e n c ea n dT e c h n o l o g y , N o r t h w e s tA g r i c u l t u r ea n dF o r e s t r y U n i v e r s i t y S h a n n x iy a n g l i n g7 1 2 1 0 0 ，C h i n a ) A b s t r a c t ：B a s eO l l1 1 0L u x ic o wm e a s u r e m e n t so fS h a n d o n gp r o v i n c ei n2 0 0 6t h er e l a t e dc o e f f i c i e n tw e r ea n a l y z e d b e t w e e nb o d yw e i g ha n db o d yh e i 曲Lb o d yl e n g t h ，h e a r tg i r t h ，c i r c u mf e r e n c eo fC a n f i O R E s t a b l i s ht h ef o r e c a s tL u x i C O Ww e i g hr e g r e s s i o nm o d e lw a se s t a b l i s h e d T h er e s u l ts h o w e dt h a tt h eb o d yw e i 【g ha n db o d yh e i g 虮b o a yl e n g t h h e a r tg i r t h , c i r c u mf e r e n c eo fc a n l l o l lr e l a t e dc o e f - t i c i e n ta l e0 7 5 l 、0 8 1 4 、0 9 0 9 、0 8 5 5 ：G e tt w of o r e c a s tw e i 血 r e g r e s s i o n m o d e l s w e r er e