体育科研中运用多元回归分析时应注意的两个问题

   【说明】目前在体育科研中，有些人在做多元回归分析时，还是经常出现错误。因此，很有必要将我十几年前在全国体育统计学会报告会上发表的论文，供大家参考。

   回归分析是研究变量之间的相关关系的一种数理统计方法。在体育领域中存在着大量多因素（变量）的问题。所以，近年来体育科研中多元回归和逐步回归分析方法运用得越来越多了。

   回归分析主要可用于解决下列问题：

  （一）确定几个变量之间是否存在相关关系。如果存在，则找出它们之间合适的数学表达式──回归方程。如在体质研究、运动员选材研究中为了研究各项形态、机能、素质指标之间的相关关系，常常通过大样本统计数据计算各指标之间的相关系数ｒij，并建立各指标之间的回归方程。如1979年全国体质调研后建立了18～25岁男青年的身高、体重对肺活量的二元回归方程：

   Ｙ＝－2708＋27.89Ｘ1＋35.56Ｘ2
（肺活量）        （身高）   （体重）

  （二）根据回归方程用一个或几个自变量的值预测或控制另一个因变量的取值。对于一些不容易测的指标，可以通过建立回归方程用几个比较容易测的指标来推测它的估计值。如果自变量选得恰当，这样的回归方程是很有实用价值的。如北京体院高强等人研究建立了推测快肌纤维百分比的三元回归方程：

    Ｙ＝5.90＋35.53 Ｘ1＋18.2 Ｘ2一6.41 Ｘ3
 （快肌％）      （相对肌力） （MPE ％）  （IEMG）

    只要进行肌肉力量和肌电图的几项测试，就可以实现肌纤维成分的无损伤测定。

  （三）进行因素分析。当许多个变量（因素）都对某一个变量有影响时，可以通过回归分析找出其中哪几个是重要因素，哪几个是次要因素。如对田径十项全能世界级优秀运动员的十项全能总分和十个单项成绩进行逐步回归计算，从中可筛选出四个比较重要的项目是：400米、铅球、110米栏、跳高。 
  有的研究还根据多元回归方程的系数ｂi。或标准回归系数ｂi’ 来确定各指标的“权重”。
  但是，近来在运用多元回归和逐步回归的研究论文中也出现了一些对回归分析的运用条件注意不够的问题。

一、样本含量ｎ和自变量个数k之间的关系问题

  复相关系数 Ｒ 是检验多元回归方程效果的重要指标。一般讲，Ｒ越接近１即表示回归方程的效果越好。因此当计算的结果Ｒ值接近ｌ时，有些作者就立即认定计算结果十分理想。如“广东省少体校游泳运动员因素分析与运动模型，’（1983年全国体育统计报告会论文）一文中，男15～16岁组自由泳成绩与形态机能指标的十元回归方程Ｒ＝0.999，十三元回归方程Ｒ＝１。又如“对男女优秀篮球运动员五大关节活动幅度与运动成绩逐步回归的探讨”（体育科学86年１期）一文中，辽宁男篮队员运动成绩与四项关节活动幅度建立的四元回归方程Ｒ＝0.995，辽宁女篮队员运动成绩与七项关节活动幅度建立的七元回归方程Ｒ＝１。作者都认为回归方程十分理想，并根据回归方程作了进一步的分析和结论。但是，再看回归方程中各个指标与因变量（Ｙ）的相关系数则发现有许多指标和Ｙ是不相关的（p>0.05）。可见，Ｒ接近１的原因并不是这些指标和Ｙ相关程度极高，而是作者用于计算的样本 n 太小（游泳一文ｎ＝15，篮球一文ｎ＝10）。
  在《体育统计方法》及一些数理统计书中，都指出：复相关系数 R 与多元回归方程中自变量的个数 k 及样本含量 n 有关。当 n 相对于 k 并不很大时，常有较大的 R，当 n=k+1时，即使这 k 个自变量与 Y 并不相关，也会 R＝1。

  如设变量 Y与另外 K 个变量 Ｘ1, Ｘ2……Ｘk的内在联系是线性的，它的第 α 次测试数据是：
    （Ｙa,Ｘa1,Ｘa2……，Ｘak）       a＝1,2,…,n
    则这一组数据可以有如下的结构式：
Ｙ1＝β0＋β1Ｘ11＋β2Ｘ12＋……＋βkＸ1k＋ξ1
Ｙ2＝β0＋β1Ｘ21＋β2Ｘ22＋……＋βkＸ2k＋ξ2
         ……
Ｙn＝β0＋β1Ｘn1＋β2Ｘn2＋……＋βkＸnk＋ξn

  其中：β0,β1……是 K+1 个待估计参数，Ｘ1，Ｘ2，……，Ｘk是 K个可以精确测得的变量，ξ1，ξ2,……ξn是 n个相互独立且服从正态分布 N（θ，δ）的随机变量，这就是多元回归的数学模型。把它写成矩阵形式：
                  Ｙ＝Ｘβ＋ξ

  用最小二乘法原理，建立正规方程，可解出 bo，bl……，bk。它们是参数 β0,β1……βk的最小二乘估计，则多元回归方程为：
    Ｙ＝bo＋b1Ｘ1＋b2Ｘ2＋……＋bkＸk
    计算bo, b1, b2,……，bk时要计算下列四个矩阵：X,A,C,B
    其中：X是 n 组测试数据Ya的结构矩阵，
          A＝X'X 是正规方程组的系数 L ij矩阵，
          C＝A^-1 是系数矩阵A的逆矩阵，  
          B 是正规方程组的常数项 Liy矩阵。
   多元回归的系数 b=A^-1 B．可见，多元回归方程的系数 bi 的解是和所测得的原始数据的组数（样本含量）n 紧密相系的。
   可以推导出：当 n=k + 1时，结构矩阵 X 是方阵，这时解出的bi，即使 Xi和 Y 全都不相关，仍可以得到：
  Q=E(Y-Ý )＝0 亦即：U=Σbi Lij＝Lyy 所以 R＝sqrt(U/Lyy)＝1。

  为了用计算实例来说明以上的理论推导，我们随机抽取了 K=34，n＝10的一个样本。其中大部分 X和Y 是不相关的（即r iy的值，经检验为 P>0.05）。选取其中9个和Y不相关的X作九元回归计算，riy最小的为-0.011,最大的为0.471,均为P>0.05,计算结果见表1，虽然每个X的偏回归平方和Pi作F检验时，都是P>0.05，但复相关系数R＝0.9999975。
  把多元回归的X增加到15，20和34个时，R值分别为：0.9999891，1.000009,1.0000830可见当样本含量n太小与自变量个数K的比例不恰当时，解出的bi是不能正确地反映各自变量和Y之间的内在关系的。如果用这样的多元回归方程为依据去分析和作结论，就一定会产生谬误。

   表1:     K＝9   n＝10的多元回归计算结果

    R=0 .9999975      Lyy＝10.05696      U＝10.05692

   有的数理统计专著上提出：n 应是 K 的5～10倍，这是我们在计算时应特别注意的问题。
   体育科研中，特别是研究优秀运动员的时候，样本含量一般都比较小。而希望分析与专项成绩有关的指标时，又常常取许多个Xi，所以，K和n的比例常是一个大问题。
   我们认为：在运用多元回归分析时，首先应尽可能地扩大样本含量n，特别要注意n与K的比例。其次，在算出多元回归方程后，不仅要看R值的大小（或对方程作方差分析后的F值大小），还必须仔细看一下每个自变量的Pi的F检验结果，以决定这个方程是否具有实用价值。

二、自变量之间相关的问题

  在体育科研中常常运用多元回归或逐步回归来进沂因素分析。即从许多个自变量Xi中分析哪几个是影响Y的主要因素，哪几个是影响Y的次要因素。
  按照回归分析的原理，在多元回归中，回归平方和 U 是所有自变量对 Y 变差（波动）的总贡献。在原有的自变量中减去一个自变量后，U 就会减小，减小的数值越大，就说明这个自变量对 Y 的作用越大。所以，把取消一个自变量后，回归平方和 U 减少的值称为：Y 对这个自变量的偏回归平方和Pi。但是，通常情况下由于自变量之间互相存在着相关关系，所以 ΣPi≠U，只有当各自变量之间的相关系数 rij=0 即正规方程的系数矩阵为对角阵（除主对角线元素L11,L22，……，Lkk外其余元素全为0）时，才有U＝ΣPi。
  通过计算偏回归平方和Pi后可以作比较，Pi值越大的自变量对Y的作用就越大。如果Pi经F检验为P＞0.05，就可认为它对Y的影响不显著。但是，这样用Pi来进行比较的前提是：各自变量之间的相关程度应该很小，即r ij的绝对值都不大，其分析结果才正确。
  可是，体育科研中对许多指标作回归分析时，往往自变量之间的相关程度都很高。如以100米跑成绩为Y,分析各项形态、机能、素质指标和100米跑成绩之间的关系时，形态指标中的身高、指距、下肢长……这类长度指标之间的 rij都在0.8一0.9左右，体重和胸围、上臂紧张围、大腿围……等围度类指标之间的rij约在 0.8-0.9之间，素质指标中30米跑、立定跳远、纵跳……之间的rij也在0.7左右。如果计算者不注意自变量之间的相关性的问题，按计算的结果下结论，也会出现错误的。
  因为，假设在一个多元回归中，X1和X2对Y的影响都很大，但是它们之间高度相关。由于X1对Y的影响中很大部分可以从X2对Y的影响中反映出来，这时如从回归方程中减去X1，U并不会减小很多，即P1的值一定不大，同样的理由，P2也不会很大。因此，这时看P1、P2的值就不能正确地反映出X1, X2对Y的重要性。
  我们在一次四元回归计算中遇到这样一个实例（n=31），四个变量和Y之间的相关系数如表2。由表2可见，X1 与 X2的相关系数r12=0.949呈高度相关，而其他Xi之间不仅相关系数rij的绝对值较小．而且P>0.05,即相互之间相关不显著。再看表上r3y＝0.245（P>0.05），即X3与Y相关不显著。但做四元回归计算出的Pi（见表3），P2＜P3而且对P2作F检验，F＝0.528( P>0.05)，似乎可得出结论X2对 Y的影响比X3小，而且由于经F检验不显著应把X2剔除。

     表2

   注：** 为P＜0.01    * 为P＜0.05

     表3

  当我们考虑到 X1和X2之间高度相关，而把 X1剔除后再作三元回归计算（见表3），P2就大大地增大为3.364034不仅F检验结果P<0.01，而且因为P2>P3说明X2比X3的作用更大些。同时从表3上还可以看出，在四元回归中P1、P2都小于1,远比三元回归中的P2小。这个实例可以让我们更加理解回归分析原理中关于：自变量之间的相关性会给因素分析带来麻烦的阐述。
  我们认为：在做多元回归分析时，如遇到自变量之间相关程度较高，最好先做R型聚类分析，然后从各类中选出一个典型指标，再做多元回归计算。但是有时两个相关的同类指标按专业知识分析都不应剔除，如在选材的形态指标中下肢长B和下肢长C虽然高度相关，但这两个指标在选材中各有不同的意义，不能替代，都必须保留。这时如用多元回归计算Pi进行分析时，就应特别注意了。

   多元回归分析是体育科研中一个常用的计算方法，以上是实际使用中较易出差错的两个问题，提出来与大家共同商讨，不当之处请批评。