【SPSS操作】Non-parametric ANCOVA (Quade‘s) in SPSS

在SPSS中执行Quade的非参数ANCOVA需要手动计算步骤，因为SPSS没有内置的菜单选项。以下是详细的操作指南。

步骤1：准备数据

确保数据包含以下变量：

• 分组变量（Group）：自变量，包含不同的处理组。
• 协变量（X）：需要调整的连续变量。
• 因变量（Y）：结果变量。

步骤2：计算协变量和因变量的秩

1. 对协变量X进行排名：

   • 点击菜单：Transform → Rank Cases。
   • 将变量X放入“Variable(s)”框。
   • 在“Assign Rank 1 to”中选择Smallest Value（默认）。
   • 点击“Ties”选项，选择Mean处理结（ties）。
   • 生成新变量（如Rx）。

2. 对因变量Y进行排名：

   • 重复上述步骤，生成Y的秩变量（如Ry）。

步骤3：调整因变量的秩

1. 用线性回归调整Ry中的协变量效应：
   • 点击菜单：Analyze → Regression → Linear。
   • 因变量选Ry，自变量选Rx。
   • 点击“Save”按钮，勾选Unstandardized Residuals，生成残差变量（如Ry_adj）。

步骤4：非参数检验

1. 使用Kruskal-Wallis检验比较调整后的残差：
   • 点击菜单：Analyze → Nonparametric Tests → Independent Samples。
   • 在“Fields”选项卡中，将Ry_adj放入Test Fields，分组变量Group放入Groups。
   • 在“Settings”选项卡，选择Customize tests，勾选Kruskal-Wallis 1-way ANOVA (k samples)。
   • 点击Run。

结果解读

• 如果Kruskal-Wallis检验的显著性（p值）小于0.05，表明在调整协变量后，组间存在显著差异。

---

注意事项

• 方法局限性：此方法通过残差调整近似Quade检验，可能与原始方法存在差异。
• 严格方法：如需严格实现Quade检验，建议使用R语言（如coin包）或手动编写SPSS语法计算统计量。

语法示例（调整残差法）

RANK VARIABLES=X BY GROUP /RANK INTO Rx.
RANK VARIABLES=Y BY GROUP /RANK INTO Ry.

REGRESSION
  /DEPENDENT Ry
  /METHOD=ENTER Rx
  /SAVE RESID(Ry_adj).

NPAR TESTS
  /K-W=Ry_adj BY GROUP(1, k)  /* k为组数 */
  /STATISTICS DESCRIPTIVES
  /MISSING ANALYSIS.