杜泊羊（Dorper sheep）原产于南非，简称杜泊羊，属于世界上著名的肉用绵羊品种之一。该品种生长发育快，具有典型的肉用体型、肉用品质好、体质结实、繁殖性能好，无明显繁殖季节限制，可常年发情和配种，且对炎热、干旱和寒冷等气候有良好适应性[1]。家畜体尺数据直接反映家畜体格大小、体躯结构及发育等情况，与家畜生理机能、生产性能、抗病力以及对外界环境条件的适应能力等密切相关。国内外学者针对有关杜泊羊生产性能、生长发育、饲养管理及改良地方品种等方面研究较多[2-12]，但对体质量与体尺指标相关性研究尚缺乏报道。本试验对以湖羊为受体采用胚胎移植技术扩繁的杜泊羊体质量与体尺指标进行测定，初步掌握杜泊羊生长发育情况，并进行相关性分析、通径分析、逐步回归分析以及对表型相关系数分解，探究体尺指标与体质量关系，建立体质量与体尺指标之间最优回归方程，为肉羊新品种选育工作提供参考。1材料与方法1.1试验动物试验动物杜泊羊均为2017年9月底在甘肃中盛华美羊业西峰羊场（以湖羊为受体，从天津奥群公司购买的纯种杜泊A级冷冻胚胎）采用胚胎移植技术所产，6月龄左右，健康状况良好。1.2测定指标随机抽取60只（公、母羊各半）健康羊只进行体质量（y）、体尺（xi）测定。体质量称量用电子秤，精度为0.01 kg。体尺指标测量参考家畜育种学[13]，测量工具为测杖和卷尺，测量指标包括体高（x1）、体斜长（x2）、胸围（x3）、胸深（x4）、胸宽（x5）和管围（x6）。1.3数据统计与分析数据采用Excel 2007软件整理，采用SPSS 19.0软件进行相关性分析、通径分析和逐步回归分析，并对相关系数进行分解，建立体质量与体尺指标之间的最优回归方程，结果以“平均值±标准差”表示。2结果与分析2.1体质量与体尺指标的表型统计量（见表1）由表1可知，体质量变异系数最大，为10.8%；其次为胸宽、胸深、管围、胸围、体高，其变异系数分别为9.4%、5.9%、5.8%、5.2%、4.5%；体斜长的变异系数最小，为3.7%。10.13557/j.cnki.issn1002-2813.2021.13.001.T001表1杜泊羊体质量与体尺指标的表型统计量项目体质量y/kg体高/cm体斜长/cm胸围/cm胸深/cm胸宽/cm管围/cm变异系数/%10.84.53.75.25.99.45.8平均数50.4863.6369.3895.1526.7521.8011.03标准差5.442.882.554.941.572.050.642.2杜泊羊体质量与体尺指标的表型相关性分析（见表２）由表2可知，体质量与体高、体斜长、胸围和管围呈极显著相关（P0.01），与胸深呈显著相关（P0.05），体质量与胸围的表型相关系数最大，为0.560；与体斜长、管围、体高、胸深、胸宽的表型相关系数为0.505、0.372、0.331、0.263、0.066。其中，体高与体斜长和胸围、体斜长与胸围，胸围与胸深均呈极显著相关（P0.01）。体斜长与胸深，胸围与胸宽和管围均呈显著相关（P0.05）。10.13557/j.cnki.issn1002-2813.2021.13.001.T002表2杜泊羊体质量与体尺指标的表型相关系数指标体高体斜长胸围胸深胸宽管围体质量体高1体斜长0.390**1胸围0.455**0.489**1胸深0.2310.271*0.360**1胸宽0.0420.0590.266*0.1041管围0.1650.1910.323*0.0180.0171体质量0.331**0.505**0.560**0.263*0.0660.372**1注：*表示相关性显著（P0.05），**表示相关性极显著（P0.01）。2.3体质量与体尺指标之间的通径分析为进一步揭示体质量与体尺指标间的相关关系，对体质量和体尺指标进行通径分析，用体尺指标与体质量的表型相关系数建立以通径系数为变量的多元一次方程组，方程组如下：P1+0.390P2+0.455P3+0.231P4+0.042P5+0.165P6=0.331 (1)0.390P1+P2+0.489P3+0.271P4+0.059P5+0.191P6=0.505 (2)0.455P1+0.489P2+P3+0.360P4+0.266P5+0.323P6=0.560 (3)0.231P1+0.271P2+0.360P3+P4+0.104P5+0.018P6=0.263 (4)0.042P1+0.059P2+0.266P3+0.104P4+P5+0.017P6=0.066 (5)0.165P1+0.191P2+0.323P3+0.018P4+0.017P5+P6=0.372 (6)用矩阵法对多元一次方程组进行求解，得到杜泊羊体尺指标对体质量的通径系数分别为：P1=0.021；P2=0.278；P3=0.339；P4=0.063；P5=-0.051；P6=0.206。2.4体尺指标对体质量的相关系数分解（见表３）根据相关系数的组成效应，可将杜泊羊各体尺指标（xi）与体质量（y）的表型相关系数（皮尔森相关系数）分解为体尺对体质量的直接作用和该体尺指标通过其他体尺性状对体质量的间接作用两部分。由表3可知，杜泊羊的胸围对体质量的直接作用最大；胸宽对体质量的直接作用最小。胸围对体重的综合间接作用最大，其次为体高、体斜长、胸深、管围。各体尺指标对体质量的直接作用大小为：胸围体斜长管围胸深体高胸宽；各体尺指标对体质量的间接作用大小为：胸围体高体斜长胸深管围胸宽。10.13557/j.cnki.issn1002-2813.2021.13.001.T003表3杜泊羊体尺对体质量相关系数的分解自变量对体质量的相关系数直接作用间接作用体高体斜长胸围胸深胸宽管围合计体高0.3310.021—0.1970.2550.0610.0030.0610.577体斜长0.5050.2780.129—0.2740.0710.0040.0710.549胸围0.5600.3390.1510.247—0.0950.0180.1200.631胸深0.2630.0630.0760.1370.202—0.0070.0070.429胸宽0.066-0.0510.0140.0300.1490.027—0.0060.226管围0.3720.2060.0550.0960.1810.0050.001—0.3382.5体尺指标对体质量的逐步回归分析2.5.1体质量和各体尺指标分析模型（见表4~表6）由表4~6可知，体高、胸深、胸宽、管围对体质量影响不显著，应当从自变量中剔除。体斜长和胸围对体质量的影响显著。剔除自变量体高、胸深、胸宽和管围后，对因变量体质量与体斜长、胸围进行多元逐步回归分析。10.13557/j.cnki.issn1002-2813.2021.13.001.T004表4体质量和各体尺指标整体回归分析模型模型RR2调整R2估计标准误更改统计量R2Fdf1df2显著性10.654a0.4270.3634.3460.4276.5956530.000注：a预测变量：常量、胸宽、管围、胸深、体高、体斜长、胸围；下表同。10.13557/j.cnki.issn1002-2813.2021.13.001.T005表5体质量和各体尺指标整体回归分析的多元回归方差分析模型平方和自由度均方F值P值1747.3166124.5536.5950.000a10.13557/j.cnki.issn1002-2813.2021.13.001.T006表6体质量和各体尺指标总体回归分析的多元回归系数参数估计值模型指标非标准化系数标准系数t显著性B的95%置信区间B标准误差下限上限1常量-51.07219.490-2.6200.011-90.164-11.980体高0.0410.2270.0210.1780.859-0.4160.497体斜长0.5920.2620.2782.2570.0280.0661.119胸围0.3740.1550.3392.4170.0190.0640.685胸深0.2170.3920.0630.5530.583-0.5701.004胸宽-0.1360.289-0.051-0.4710.639-0.7170.444管围1.7610.9470.2061.8590.069-0.1393.6612.5.2体质量和各体尺指标逐步回归分析模型（见表7~表9）逐步回归分析即按自变量对因变量作用的显著程度，从大到小逐个引入回归方程，即将自变量的t值从大到小逐个引入回归方程。由表6可知，各个自变量的显著性由大到小依次为胸围体斜长管围（P0.05）胸深（P0.05）胸宽（P0.05）体高（P0.05），因此建立体质量和胸围、体斜长之间的回归分析模型2（见表7）。10.13557/j.cnki.issn1002-2813.2021.13.001.T007表7杜泊羊体质量和胸围、体斜长指标逐步回归分析模型模型RR2调整R2估计标准误更改统计量R2Fdf1df2显著性10.560a0.3140.3024.5480.31426.5361580.00020.620b0.3840.3634.3460.0706.5161570.013注：a.预测变量：常量、胸围，b.预测变量：常量、胸围、体斜长；下表同。由模型2所得到的方程为体质量和体尺指标之间的最优回归方程为：y=-37.724＋0.649x2＋0.454x3 (R＝0.620，Ｐ0.05）。标准化后的回归方程为：y＝0.304x2＋0.412x3（见表9）。方差分析的F值为17.788（见表8），达到极显著水平（Ｐ0.01）。10.13557/j.cnki.issn1002-2813.2021.13.001.T008表8体质量和胸围、体斜长指标逐步回归分析的多元回归方差分析模型平方和自由度均方F值P值1548.7891548.78926.5360.000a2671.8452335.92317.7880.000b10.13557/j.cnki.issn1002-2813.2021.13.001.T009表9体质量和胸围、体斜长指标逐步回归分析的多元回归系数参数估计值模型指标非标准化系数标准系数t显著性B的95%置信区间B标准误差下限上限1常量-8.30811.427-0.7270.470-31.18314.566胸围0.6180.1200.5605.1510.0000.3780.8582常量-37.72415.876-2.3760.021-69.515-5.934胸围0.4540.1310.4123.4570.0010.1910.717体斜长0.6490.2540.3042.5530.0130.1401.1573讨论3.1体质量和体尺指标的表型统计量变异系数又称离散系数，是反映标志变动相对程度的指标。体质量的变异系数（10.8%）大于10%，说明6月龄杜泊羊的体质量变异幅度较大，体质量具有较大的选择与提高潜力。其余6个体尺指标的变异系数均小于10%，说明杜泊羊其对应的6个指标变异幅度较小。3.2体质量与体尺指标的相关性分析本试验通过对6月龄杜泊羊体质量和其体尺指标进行相关性分析，发现杜泊羊体尺与体质量都呈正相关，其中体质量与体高、体斜长、胸围和管围相关性极显著（P0.01），与胸深相关性显著（P0.05）。研究表明，体高、体斜长、胸围、管围和胸深对体质量的综合作用较大，表现为该体尺指标对体质量产生直接作用或通过辅助其他性状对体质量产生间接作用。在体高、体斜长、胸围、胸深、胸宽和管围6个自变量之间亦有较强的正相关，其中体高与体斜长和胸围、体斜长与胸围，胸围与胸深均相关性极显著（P0.01）；体斜长与胸深，胸围与胸宽和管围均相关性显著（P0.05）。与许鑫等[14]对杜泊羊、王欣荣等[15]对藏羊、罗建川等[16]对贵州杂交绵羊、顾亚玲等[17]对萨福克羔羊、王志明等[18]对甘肃高山细毛羊、李少斌等[19]对甘肃高山细毛羊杂交一代以及喇永富等[20]对湖羊体质量与其体尺指标分析结果相类似，说明杜泊羊的各个体尺指标是相互作用、相互联系的，选种时需要对杜泊羊的各个体尺指标都要加以考虑，从而最终能达到最佳的选种效果。3.3体质量与其体尺指标间的通径分析及其相关系数的分解通径系数表示自变量对因变量的直接作用大小，可以真实地反映两个变量间的直接相关程度，它是标准化的偏回归系数[21-24]。通过对体质量和各体尺指标之间的通径分析以及其相关系数的分解可以看出：6月龄杜泊羊体质量与体高、体斜长、胸围和管围均表现出极显著相关性。本研究中，胸围、体斜长对体质量既有较强的直接作用，又有较强的间接作用；管围主要是通过直接作用对体质量产生影响，对体质量的间接作用较弱；体高主要是通过间接作用对体质量产生影响，对体质量的直接作用较弱；胸深、胸宽对体质量的直接作用、间接作用均较弱。体高、体斜长、胸深、胸宽和管围主要是通过胸围对体质量产生间接作用，胸围主要是通过体斜长对体质量产生间接作用。在各体尺指标中对体质量的直接作用大小为：胸围体斜长管围胸深体高胸宽；对体质量的间接作用大小为：胸围体高体斜长胸深管围胸宽。上述分析表明，杜泊羊的生长发育是1个复杂的有机整体，各指标间有着种种密切而复杂的联系。影响体质量的因素不仅多，而且影响途径也不相同。所以在选种过程中除加强对体质量影响直接作用强的因素如杜泊羊的胸围、体斜长选择的同时，也应该关注其他因素对体质量的作用。3.4体尺指标对体质量的逐步回归分析通过对6月龄杜泊羊体质量与其体尺指标之间的逐步回归分析，得到体质量和体尺指标之间标准化后的回归方程为y＝0.304x2＋0.412x3（P0.01），说明胸围和体斜长是影响杜泊羊体质量的主要因素，且用这2个体尺指标来估计杜泊羊体质量是可信的，也说明所选性状对体质量有比较大的决定作用，拟合的方程具有一定的参考价值。其他体尺指标的P值均大于0.05，说明在多元回归中，这些体尺指标对体质量的影响不显著。由方程可知，当体斜长保持不变时，胸围每增加1 cm，杜泊羊的体质量增加0.412 kg；当胸围保持不变时，体斜长每增加1 cm，杜泊羊的体质量增加0.304 kg。本试验回归结果入选指标与许鑫等[14]对特克赛尔羊、白萨福克羊及杜泊羊，王伟峰[25]对白萨福克羊，王高富等[26]对酉州乌羊，索效军等[27]对甘肃高山细毛羊的研究结果相同。已有相关研究[15,27-30]对藏羊、萨福克、乌骨山羊、隆林山羊、辽宁绒山羊等的回归方程中，除了胸围、体斜长入选外，有些也入选了体高、胸宽和胸深等因素，可能是由于试验羊的品种、生活环境、饲养管理水平等诸多因素不同而造成的。4结论本试验条件下，除胸宽外，体高、体斜长、胸围、胸深和管围与6月龄杜泊羊体质量均表现出显著相关性；各体尺指标中对其体质量的直接作用大小为：胸围体斜长管围胸深体高胸宽；对其体质量的间接作用大小为：胸围体高体斜长胸深管围胸宽；胸围、体斜长是影响其体质量的主要因素，其体质量和体尺指标之间的最优回归方程为y=-37.724＋0.649x2＋0.454x3（R＝0.620，P0.05），标准化后的回归方程为y＝0.304x2＋0.412x3（F＝17.788，P0.01）。