成年酉州乌羊体尺长与体质量的关系

王高富　刘良佳　任航行　蒋婧　余琴　黎年富　夏元友　周鹏

摘要：为研究成年酉州乌羊体尺长与体质量的关系，测定16只酉州乌羊成年公羊和74只成年母羊的体斜长、体高、胸围、管围和体质量等指标，运用SPSS软件分别分析公母羊体尺长与体质量间的相关系数、体尺长对体质量的直接和间接影响以及对体质量的决策程度，分别建立公羊和母羊体尺长与体质量的最优回归模型。结果表明，胸围和体斜长是酉州乌羊母羊体质量的最主要决策变量，胸围是公羊体质量的最主要决策变量。母羊和公羊体尺对体质量的最优回归模型分别是Y=-37569+0262X1+0702X3，Y=-44157+1035X3，其中X1表示体斜长，X3表示胸围。

关键词：酉州乌羊；体斜长；胸围；体质量；回归分析

中图分类号： S8268+92文献标志码： A

文章编号：1002-1302（2017）21-0165-03

收稿日期：2016-06-15

基金项目：重庆市农发资金项目（编号：15404）；重庆市基本科研业务费计划项目（编号：15435）。

作者简介：王高富（1980—），男，重庆人，硕士，副研究员，研究方向为羊遗传育种。E-mail：wanggaofs20031216@163com。

通信作者：周鹏，高级畜牧师，研究方向为动物遗传育种。E-mail：cqzp2006@163com。

酉州乌羊主要分布于武陵山区腹地的酉阳土家族苗族自治县，是重庆市特有的本地山羊品种，2009年正式被评定为国家级遗传资源1]，该羊肉质细嫩，全身皮肤、眼、鼻、嘴、肛门、阴门等可视黏膜均为乌色2]，肉、骨、内脏均可入药3-4]，被当地群众称为“药羊”3]。但酉州乌羊群体数量小，现存栏数量不足5 000只，生长速度慢，成年个体体质量小，对其进行选育提高已经迫在眉睫。

体质量性状是肉羊的重要经济性状，是肉羊育种过程中主要选择性状之一。而体尺性状反映肉羊体格大小、体躯结构和组织器官发育情况，它与肉羊的繁殖、生长、抗病力以及对外界生活条件的适应能力等均有一定的相关性。因为体尺性状数量较多，研究这些性状与体质量的关系，找出较为主要的体尺性状作为育种指标具有一定的实用价值。为了使酉州乌羊的保种选育工作更加高效快捷，本研究拟测定成年酉州乌羊的4个体尺和体质量性状，分析它们之间的关系，以期建立体质量与体尺之间的回归模型。

1材料与方法

11试验对象、测定指标及方法

试验在重庆市酉州乌羊资源保护场进行，随机选择常规饲养条件下的成年公羊16只，成年母羊74只进行体尺和体质量测定，试验羊体况中上等，健康无病。测定指标包括体斜长（X1）、体高（X2）、胸围（X3）、管围（X4）和成年体质量（Y）；测定方法按照《畜禽遗传资源调查技术手册》执行5]，主要测定工具是测杖、卷尺、电子秤等。

12统计方法

试验数据录入Excel数据库中，用SPSS统计软件分析体质量与体尺间的相关系数、通径系数、体尺对体质量的决策程度，并建立体质量与体尺的最优回归模型6]。

2结果与分析

21各性状的表型参数

由表1所测各性状的平均数、标准差及变异系数可知，酉州乌羊成年公母羊体质量的变异系数在各性状中是最大的，分别达到3153%和2403%，这表明酉州乌羊的体质量变化较大，具有较大的选育提高空间。另外，成年公羊的体斜长、体高、胸围和管围等体尺指标的变异系数均大于成年母羊，并且变异系数均較大，说明酉州乌羊体尺性状也有较大的选育提高空间。

22体尺、体质量表型相关分析

由表2、表3可以看出，成年公母羊体尺与体质量都呈极显著正相关（P<001），且公羊各性状间的相关系数均大于母羊。其中公羊胸围与体质量的相关系数最大，达到0972，体斜长次之，为0920，其他依次为管围和体高；母羊胸围与体质量的相关系数最大，为0893，其他依次为体斜长、体高和管围。而公母羊在各体尺之间亦有较强的正相关，除母羊体斜长与管围呈显著正相关（P<005），其他性状间都达到了极显著相关（P<001），这说明成年酉州乌羊的体尺、体质量是相互依存的不同性状， 对其中任何性状进行选择都应该综合考虑。

23不同性别酉州乌羊体质量的正态性检验

SPSS对一组数据进行正态性检验有2种方法，即JP3]Kolmogorov-Smirnov Test和Shapiro-Wilk Test。Kolmogorov-Smimov Test检验结果较精确，适用于大样本的检测。而Shapiro-Wilk Test适用于小样本的检验。本研究母羊样本数量73只，公羊样本数量16只，均属于小样本，因此对因变量Y进行正态性检验后利用Shapiro-Wilk Test的输出结果。母羊Shapiro-Wilk统计量0968，P值=0061>005（表4），公羊统计量0953，P值=0535>005（表5），所以公母羊体质量均服从正态分布，可以进行回归分析。

24通径分析

各体尺性状对体质量的相关仅反映两性状间的表型相关，而不能反映各体尺性状与体质量的直接关系大小。因此，必须深入分析各因变量（体尺）对自变量（体质量）的直接影响，即通径分析，进一步揭示各性状间的相关原因。根据通径分析原理，利用各性状间的表型相关系数建立体尺对体质量的通径系数正规方程组。

241母羊正规方程组方程组如下：

P1+0563P2+0728P3+0240P4=0755；

0563P1+P2+0626P3+0329P4=0624；

0728P1+0626P2+P3+0503P4=0893；

0240P1+0329P2+0503P3+P4=0402。

解方程得：P1=017；P2=0085；P3=0709；P4=-0025。endprint

通过对通径系数进行显著性检验可知：P3达到了极显著的水平（P<001），P1达到显著水平（P<005），P2、P4均未达到显著水平，表明胸围和体斜长对母羊体质量的直接影响较大，其中以胸围对母羊体质量的直接影响作用最大（P3=0709），其次是胸体斜长（P1=017），体高和管围对体质量的直接影响较小，这与表型相关的分析结果一致。

242公羊正规方程组方程组如下：

P1+0871P2+0928P3+0866P4=0920；

0871P1+P2+0873P3+0824P4=0874；

0928P1+0873P2+P3+0890P4=0972；

0866P1+0824P2+0890P3+P4=0912。

解方程得：P1=0046；P2=0055；P3=0697；P4=-0206。

通过对通径系数进行显著性检验可知：P3达到了极显著的水平（P<001），P1、P2、P4均未达到显著水平，表明胸围对公羊体质量的直接影响较大，体斜长、体高和管围对体质量的直接影响较小。

25相关系数分解

通径分析在多元回归的基础上将相关系数riy分解为直接通径系数（某一自变量对因变量的直接作用）和间接通径系数（该自变量通过其他自变量对因变量的间接作用）。通径分析的理论已证明，任一自变量（x）与因变量（y）之间的简单相关系数（riy）=xi与y之间的直接通径系数（Piy）和所有xi与Y的间接通径系数，任一自变量xi对y的间接通径系数=相关系数（riy）×通径系数（Pjy）。

251母羊各体尺对体质量的相关系数分解

由表6可知，在4个自变量对体质量的直接影响方面，胸围（X3）的直接作用最大，体斜长次之，管围直接作用最小。通过分析各个间接通径系数发现，体斜长通过胸围对体质量的间接作用最大，其间接通径系数达到0516 2，其次是体高和管围通过胸围对体质量的间接作用，其间接通径系数分别达到了0443 8和0356 6，而体高和管围对体质量的直接影响却较小。由此可见，体高（X2）和管围（X4）对体质量的影响主要取决于间接作用；胸围（X3）对体质量不但有较强的直接影响，而且辅助其他性状对体质量产生较大的间接作用。因此，胸围和体斜长对酉州乌羊母羊体质量的增加具有重要作用，而体高和管围对体质量改变影响不大，在选育时可以不必过多考虑。

252公羊各体尺对体质量的相关系数分解

由表7可知，4个自变量对体质量的直接影响中，胸围X3的直接作用最大。通过分析各个间接通径系数发现，体斜长通过胸围对体质量的间接作用最大，其间接通径系数达到0646 8，其次是管围和体高通过胸围对体质量的间接作用，其间接通径系数分别达到了0620 3和0608 5，而体斜长、体高和管围对体质量的直接影响却较小。由此可见，体斜长（X1）、体高（X2）和管围（X4）对体质量的影响主要取决于间接作用；胸围（X3）对体质量不但有较强的直接影响，而且还辅助其他性状对体质量产生较大的间接作用。因此，胸围对酉州乌羊公羊体质量的增加具有重要作用。

26决策程度分析

袁志发等提出决策系数的概念，指出在复杂的路径信息中，选择什么样的路径对因变量Y最好，可由决策系数来决定7]。决策系数计算公式为R2（i）=2Piriy-Pi2。

由决策系数公式计算出酉州乌羊母羊各变量决策系数：

R2（1）=2×017×0755-0172=0227 8；

R2（2）=2×0085×0624-00852=0098 9；

R2（3）=2×0709×0893-07092=0763 6；

R2（4）=2×（-0025）×0402-（-0025）2=-0207 3。

决策系数排序为R2（3）﹥R2（1）﹥R2（2）﹥R2（4），且R2（4）﹤0。由决策系数可知，酉州乌羊母羊胸围（X3）和体斜长（X1）是体质量的最主要决策变量，其次是体高（X2）和管围（X4）。

由决策系数公式计算出酉州乌羊公羊各变量决策系数：

R2（1）=2×0046×0920-0046 2=0082 5；

R2（2）=2×0055×0874-0055 2=0093 1；

R2（3）=2×0697×0972-0697 2=0869 2；

R2（4）=2×（-0206）×0912-（-0206）2=-0418 2。

决策系数排序为R2（3）﹥R2（2）﹥R2（1）﹥R2（4），且R2（4）﹤0。由决策系数可知，酉州乌羊公羊胸围（X3）是体质量的最主要决策变量，其次是体斜长（X1）、体高（X2）和管围（X4）。

27最优线性回归方程的建立

根据逐步线性回归原理，从所有可供选择的自变量中逐步选择加入或剔除某个自变量，直到建立最优的回归方程。运用SPSS分析程序分别获得酉州乌羊母羊和公羊体质量和体尺之间的最优线性回归方程，回归系数等见表8。

由表8可知母羊各自变量的偏回归系数、方程截距、标准回归系数（即通径系数）、标准误差以及相对应的显著性检验结果，从而可得线性回归方程：

Y=-37569+0262 X1+0702 X3。

对成年母羊回归模型的截距（a）和2个偏回归系数（b1，b3）进行显著性检验，均达到极显著的水平（P<001）；同时，入选的2个体尺性状对体质量的决定系数R2=0821。说明回归方程極也显著存在且所选性状对体质量有较大的决定作用，拟合的方程也具有较强的参考价值。

由表9可知，公羊各自变量的偏回归系数、方程截距、标准回归系数（即通径系数）、标准误差以及相对应的显著性检验结果，从而可得线性回归方程如下：endprint

Y=-44157+1035X3。

对成年公羊回归模型的截距（a）和一个偏回归系数（b3）进行显著性检验，表明达到了极显著的水平（P<001）；同时，入选的此体尺性状对体质量的决定系数R2=0945。说明回归方程极显著存在且所选性状对体质量有较大的决定作用，拟合的方程具有较强的参考价值。

3讨论

31主要选育指标的确定

本研究选择了酉州乌羊的4个成年体尺指标作为研究对象，从相关分析、通径分析及决定程度分析表明，无论是成年公羊还是母羊，其胸围均是影响体质量的一个主要指标，与湖北乌羊和大足黑山羊等研究结果8-9]有所不一致，湖北乌羊公羊体质量最主要影响因素是体长，大足黑山羊体质量最主要影响因素是胸宽和管围，这可能是由于品种和生长环境不同所造成的。因此，在酉州乌羊本品种选育过程中应该以胸围为主并兼顾体斜长和管围作为主要选择指标。

32体质量与体尺性状的作用关系

由相关系数剖分可知，不管母羊或者公羊的体质量和胸围的相关主要由两者直接作用而形成的，间接作用影响较小。而与体质量的相关性次之的为体斜长、体高和管围，对体质量的影响主要通过间接作用产生，直接效应较小。这表明畜体是一个交互作用的有机整体，选择时应综合考虑才能取得较好的效果。

33需要进一步研究的问题

本研究只选择了几个重要的体尺性状，其他体尺性状如胸深、胸宽和腰角宽等是否存在与体质量更强的相关性，还有待进一步探讨。此外，由于本研究公羊样本量较小，变异系数较大，在以后的研究中要扩大公羊样本量，进一步研究体尺与体质量的相关关系。

HS223]参考文献：

1]ZK（#]赵金红，王高富，陈静，等 酉州乌羊胚胎移植效果初报J] 上海畜牧兽医通讯，2012（1）：014

2]周鹏，王高富，任航行，等 酉州乌羊遗传资源的保护对策J] 上海畜牧兽医通讯，2013（5）：75-75

3]夏元友，彭春江，郑义 酉州乌羊品种特性及生产性能初探J] 畜禽业，2005（7）：44-45

4]周鹏，黄勇富，王高富，等 酉州乌羊生理生化指标的研究J] 上海畜牧兽医通讯，2012（4）：32-33

5]陈伟生，徐桂芳 中国家畜地方品种资源图谱（下）M] 北京：中国农业出版社，2004

6]杜家菊，陳志伟 使用SPSS线性回归实现通径分析的方法J] 生物学通报，2010，45（2）：4-6

7]袁志发，周静芋，郭满才，等 决策系数——通径分析中的决策指标J] 西北农林科技大学学报（自然科学版），2001，29（5）：131-133

8]韩燕国，王党伟，何春波，等 湖北乌羊体尺与体重的关系研究J] 中国草食动物科学，2012（专辑）：366-369

9]陈永军，赵中权，张家骅，等 成年大足黑山羊体重与体尺通径分析及最优回归模型的建立J] 草食家畜，2008（3）：71-74endprint