手把手教你协方差分析的SPSS操作！

一、问题与数据

某研究将73例脑卒中患者随机分为现代理疗组（38例）和传统康复疗法组（35例）进行康复治疗，采用Fugl-Meyer运动功能评分法（FMA）分别记录治疗前、后的运动功能情况，部分数据如下。试问现代理疗和传统康复治疗对脑卒中患者运动功能的改善是否有差异？

二、对数据结构的分析

整个数据资料涉及2组患者（共73例），每名患者有康复治疗前、后2个数据，测量指标为FMA评分。由于治疗前的FMA分数会对治疗后的FMA分数产生影响，因此在比较现代理疗和传统康复疗法对患者运动功能的改善情况时，应把治疗前的FMA评分作为协变量进行调整，若满足协方差分析的应用条件，可采用完全随机设计的协方差分析。

协方差分析可以控制混杂因素对处理效应的影响，提高假设检验的效能和分析结果的精度。其应用条件包括：受试对象的观测指标满足独立性，各处理组的观测指标均来自正态分布总体，且方差相等。需要控制的协变量（自变量）与观测指标（因变量）之间存在线性关系，且每个组用协变量（自变量）与观测指标（因变量）进行直线回归时，回归直线的斜率相同（即各组回归直线平行）。

协方差分析相关的假设检验

1. 各组回归直线是否平行的假设检验；

2. 各组观测指标方差是否相同的假设检验；

3. 协变量（自变量）与观测指标（因变量）之间是否存在线性关系的假设检验；

4. 控制协变量的影响后，各组调整的均数是否相等的假设检验。

三、SPSS分析方法

1、数据录入SPSS

（组别1=现代理疗组，组别2=传统康复疗法组，FMA1=治疗前FMA评分，FMA2=治疗后FMA评分）

2、选择Analyze→General Linear Model→Univariate

3、选项设置

A. 主对话框设置：选择观测指标（FMA2）到Dependent Variable窗口，组别变量到Fixed Factor(s)窗口，协变量（FMA1）到Covariate(s)窗口。

B. Model设置：点击Model按钮→选择Custom选项→将组别和FMA1分别放入Model窗口→将组别和FMA1同时选中（按住Ctrl后分别点击“组别”和“FMA1”），选入Model窗口构成交互项→点击Continue返回主对话框→OK。

放入分组变量与协变量的交互项是为了检验各组回归直线是否平行，若交互项结果满足P>0.05，则尚不能认为各组协变量与观测指标之间的回归直线斜率不等。在各组回归直线平行的条件成立时，才可以考虑进一步使用协方差分析。当处理因素与协变量有交互作用，即各组回归直线平行的条件不成立时（P<>

得到结果后，这一步仅需要查看Tests of Between-Subjects Effects的结果，组别*FMA1一行为各组回归直线是否平行进行假设检验的结果。F=0.703，P（Sig.）=0.405 >0.05，尚不能认为两组治疗前FMA评分与治疗后FMA评分之间回归直线的斜率不等，即满足回归直线平行的条件（这一步是协方差分析的假设检验1）。因此，可以做协方差分析，接下来需要重复上述步骤2、步骤3，并在Model设置中，将分组变量与协变量的交互项从模型中去掉，如下图所示。

C. Options设置：点击Options按钮→在Display部分勾选Descriptive statistics（给出各组及总的例数、均数和标准差）、Homogeneity tests（给出方差齐性检验结果）、Parameter estimates（给出协方差模型的各个参数）→将组别变量放入Display means for窗口（给出各组调整均数的估计值）→点击Continue返回主对话框→OK。

四、结果解读

Descriptive Statistics表格给出了治疗后FMA评分的部分统计信息，包括两组及总的例数（N）、均数（Mean）和标准差（Std. Deviation）。

Levene’s Test of Equality of Error Variances表格给出了方差齐性检验的结果，F=0.199，P（Sig.）=0.657，尚不能认为两组治疗后FMA评分的方差不等，即满足方差齐的条件（这一步是协方差分析的假设检验2）。

Tests of Between-Subjects Effects表格给出了协方差分析结果（不含交互项），其中FMA1一行为协变量与观测指标之间是否存在线性关系的假设检验结果。F=134.213，P（Sig.）<>

组别一行为各组观测指标调整的均数是否相等的假设检验结果。F=7.866，P（Sig.）=0.007 <>

Parameter Estimates表格给出了协方差模型参数估计的结果。本例中的协方差模型为：

Estimated Marginal Means表格给出了协方差分析时观测指标的调整均数，各组调整的均数是利用参数估计的结果计算的。用各组FMA1的总平均数代入上面的协方差模型，即可得到观测指标的调整均数。

五、撰写结论

根据基线运动功能调整后，现代理疗方法对脑卒中患者进行康复治疗的运动功能得分为55（95% CI：52-57），传统康复疗法的运动功能得分为50（95% CI：48-53）。两种方法对脑卒中患者运动功能改善的差异具有统计学意义，现代理疗方法优于传统康复疗法（F=7.866，P=0.007）。