锦江母牛体质量与体尺指标的相关与回归关系

2017-03-01高雨飞彭志鹏欧阳克蕙瞿明仁黎力之熊小文温庆琪许兰娇

江苏农业科学 2017年2期

关键词：管围锦江体尺

许超，高雨飞，彭志鹏，欧阳克蕙，瞿明仁，黎力之，熊小文，温庆琪，许兰娇

(江西农业大学江西省动物营养重点实验室/营养饲料开发工程研究中心，江西南昌 330045)

锦江母牛体质量与体尺指标的相关与回归关系

许超，高雨飞，彭志鹏，欧阳克蕙，瞿明仁，黎力之，熊小文，温庆琪，许兰娇

(江西农业大学江西省动物营养重点实验室/营养饲料开发工程研究中心，江西南昌 330045)

以江西省高安市锦江牛保种区2015年测量获得的132条相关数据为基础材料，对锦江牛体质量与年龄、体高、十字部高、体斜长、胸围、腹围、管围的相关系数进行分析，并进一步分析了估测锦江牛成年母牛体质量的回归模型。结果表明：锦江牛体质量与年龄、体高、十字部高、体斜长、胸围、腹围、管围的相关系数分别为0.496、0.704、0.699、0.661、0.893、0.838、0.605；得到了3个估测体质量的回归模型，估测值与实测值之间的相关程度分别为 0.893、0.907、0.913。

锦江牛；体质量；相关系数；回归模型；体尺指标

锦江牛是优良的地方品种之一，产于我国南方亚热带丘陵地区，主要分布于江西省高安市、上高县。该品种具有体型偏小、耐粗饲、抗病力强、性成熟早、适应性强的遗传特点，是进行黄牛杂交改良的优秀母本，也是培育当地肉牛新品种不可多得的育种素材。本研究旨在探讨锦江牛成牛母牛体质量与体尺指标的相关性，进一步估测体质量的回归分析，以期为该品种牛的实际生产及选育工作提供依据。

1 材料与方法

1.1 材料来源与样本数量

试验采用江西省高安市锦江牛产区2015年8月测量获得的132条相关数据。体质量(kg)为第2天清晨测量的牛空腹质量，体斜长(cm)为肩端到坐骨端的距离，体高(cm)为鬐甲最高点至地面的垂直距离，十字部高(cm)为两腰角中央到地面的垂直距离，胸围(cm)为沿肩胛骨后角处量取的体躯垂直周径，管围(cm)为牛左前肢管骨上1/3最细处测量的水平周径。

1.2 统计方法

采用SPSS 19.0软件的Pearson和Linear过程[1]对试验数据进行处理。

1.2.1 Pearson相关系数分析计算模型

式中：R为Pearson相关系数，X为自变量，Y为因变量。R值介于-1～1之间，当R值为正数时为正相关，表示因变量随自变量的增大而增大；当R值为负数时为负相关，表示因变量随自变量的增大而减小；当R值等于0时，表示因变量与自变量之间没有相关性。

1.2.2 多元线性回归分析计算模型

Y=b0+b1x1+b2x2+…+bnxn。

式中：Y为因变量，b0为常数，b1、b2…bn为回归系数，x1、x2…xn为回归系数对应的自变量。

2 结果与分析

2.1 相关性分析结果

以年龄、体高、十字部高、体斜长、胸围、腹围、管围作为自变量，以体质量作为因变量。各指标之间的相关系数见表1。

由表1可知，年龄、体高、十字部高、体斜长、胸围、腹围、管围与体质量之间的相关系数分别为0.496、0.704、0.699、0.661、0.893、0.838、0.605。经检验，年龄、体尺指标与体质量之间的相关性均极显著(P<0.01)。其中，胸围与体质量之间的相关性最大；在体尺指标中，体高与十字部高的相关性最大，管围与其他体尺指标的相关性较小。

2.2 通径分析结果

各体尺指标与体质量的相关系数只能反映2个性状间的表型相关，如果进一步分析各性状指标间的直接或间接相互作用程度还需进行通径分析[2]。

由通径分析(表2)可知，锦江牛胸围通过体斜长和腹围对体质量的间接通径系数分别为0.177×0.596=0.105、0.206×0.867=0.179，总间接作用为0.105+0.179=0.284。由通径分析的理论可知，其相关系数为0.609+0.177×0.596+0.206×0.867=0.893，该结果与表1一致。直接作用中，胸围(0.609)最大，腹围(0.206)次之，体斜长(0.177)最小；间接作用中，胸围(0.891)最大，腹围(0.300)次之，体高(0.209)最小。

2.3 体质量与体尺指标的回归分析结果

分别采用Linear过程“Enter”法和“Stepwise”法建立多元线性回归方程，统计结果见表3、表4。由表3可知，3个回归模型分别为：体质量=-291.738+3.560×胸围；体质量=-353.458+3.085×胸围+1.186×体斜长；体质量=-359.525+2.428×胸围+1.052×体斜长+0.707×腹围。

表1 体质量与体尺指标的相关性分析结果

注：“*”“**”分别表示纵横指标间在0.05、0.01水平下相关性显著。

表2 锦江牛体重与体尺指标通径分析结果

表3 回归模型系数

表4 回归模型拟合度

注：因变量为体质量。类别1的回归模型组成：常量、胸围；类别2的回归模型组成：常量、胸围、体斜长；类别3的回归模型组成：常量、胸围、体斜长、腹围。

由表4可知，3个回归模型的复相关系数(R)分别为 0.893、0.907、0.913，拟合度的确定系数(R2)分别为0.797、0.823、0.833。可见，回归模型的线性度较好，尤其是第2、3个回归模型拟合度的确定系数(R2)更高，线性度更好。

2.4 回归系数的显著性检验

由回归系数显著性检验结果(表5)可知，3个回归模型中回归系数的F值分别为510.223、299.316、212.626，其相应的显著性检验值均小于0.001，差异极显著(P<0.01)，表明3个回归模型中回归系数是合理的。

3 结论与讨论

3.1 锦江牛母牛体质量与体尺指标的相关性

胸围是对中国黄牛及其杂交后代影响最大的体尺因素[3]，本研究对体高、十字部高、体斜长、胸围、腹围、管围6个体尺指标进行相关分析并得到了相似结果。管围是评估公牛体质量回归模型中是常见的自变量，与公牛体质量呈明显正相关关系，但对于母牛其相关性却很弱[4-6]。本研究中，锦江牛母牛管围与体质量的相关性很弱，这与王占红等的研究结论[7]一致。

表5 回归系数显著性检验结果

注：显著性检验概率值为双尾概率，99%置信区间。

3.2 通径分析

通径分析表明，胸围、腹围、体斜长对体质量的直接作用力和间接作用力均很大，且三者与体质量的相关性极显著，因此对锦江牛母牛选育时应首先考虑胸围、腹围、体斜长3个因素。本研究所测体尺指标对体质量回归方程的最大决定系数R2=0.833，表明方程中包含了估计体质量的主要体尺因素；1-R2=0.167，表明还有一些对体质量有影响的因素尚未考虑，仍需进一步研究。

3.3 锦江牛母牛体质量与体尺指标的回归分析

本研究共得到3个估测锦江成年母牛体质量的回归模型，表明锦江牛成年母牛体质量与胸围、体斜长、腹围存在明显的线性关系，此结论与国内外相关研究相吻合。判断回归模型优劣的标准是多方面的，就数理统计角度而言，根据回归模型计算拟合值与实测值之间的决定系数(R2)，通过比较R2值的大小来择其优劣更为适宜。由表3可知，第2个回归模型的R2值0.823与第3个回归模型的R2值0.833大致相等，两者均大于第1个回归模型的R2值0.793。但如果同时考虑体尺测量的繁琐程度可选取第1个回归模型，即：体质量=-291.738+3.560×胸围。

3.4 结论

锦江牛母牛体高、十字部高、体斜长、胸围、腹围、管围6个体尺指标与体质量的相关系数分别为0.704、0.699、0.661、0.893、0.838、0.605，其中胸围对体质量的影响最大。锦江牛母牛体尺指标与体质量的多元线性回归模型为：(1)体质量=-291.738+3.560×胸围(R=0.893，P<0.01)；(2)体质量=-353.458+3.085×胸围+1.186×体斜长(R=0.907，P<0.01)；(3)体质量=-359.525+2.428×胸围+1.052×体斜长+0.707×腹围(R=0.913，P<0.01)。

[1]苏金明，傅荣华，周建斌，等. 统计软件SPSS for Windows实用指南[M]. 北京：电子工业出版社，2000：418-439.

[2]田亚磊，高腾云，张花菊，等. 鲁山牛腿山羊体重与体尺性状相关性分析[J]. 江苏农业科学，2009(6)：288，291.

[3]焦平林. 秦川淮北及秦淮F1黄牛体尺体重回归分析[J]. 黄牛杂志，1998，24(4)：6-7.

[4]李永禄，邱怀，冯仰廉，等. 养牛学[M]. 北京：中国农业出版社，1987：235-236.

[5]Faroppa V. Standard of the breed and regulations[J]. La Razza Bovina Piermontese，1993(4)：9-19.

[6]曹红鹤. 意大利皮埃蒙特肉牛线性体型评分方法[J]. 黄牛杂志，1999，25(4)：17-19.

[7]王占红，何永涛，郭维春，等. 沿江牛成年母牛体重与体尺指标的相关与回归分析[J]. 中国畜牧杂志，2007，43(1)：12-13.

10.15889/j.issn.1002-1302.2017.02.045

2015-11-23

公益性行业(农业)科研专项(编号：201303143)；国家现代农业产业技术体系(编号：nycytx-38)。

许超(1988—)，男，江苏宿迁人，硕士研究生，主要从事反刍动物营养研究。E-mail：229462924@qq.com。

欧阳克蕙。E-mail：ouyangkehui@sina.com。