如何在SPSS中用EMMEANS比较所有可能的组组合？

Question

假设你有一个2x2的设计，你正在使用SPSS中的方差分析来测试这4组之间的差异。这是您的数据的图表：

​

​

在执行方差分析之后，我们可以在组之间进行6种可能的成对比较。它们是：

A-C

B-D

A-D

B-C

A-B

C-D

如果我想执行成对比较，我通常会在UNIANOVA命令后使用以下脚本：

/EMMEANS=TABLES(Var1*Var2) COMPARE (Var1) ADJ(LSD)
/EMMEANS=TABLES(Var1*Var2) COMPARE (Var2) ADJ(LSD)

但是，运行此脚本后，输出仅包含6个可能比较中的4个-缺少两个成对比较，它们是：

A-B

C-D

我如何计算这些比较？

Answer 1

UNIANOVA中的EMMEANS不能提供像这样的交互中细胞之间的所有成对比较。还有一些其他的过程，如GENLIN，确实提供了这些，但使用大样本卡方统计而不是t或F统计。在UNIANOVA中，您可以使用LMATRIX子命令获得它们，也可以使用EMMEANS的一些技巧。

对于EMMEANS的技巧，创建一个具有四个级别的单因素，用于索引单元格的2x2布局，然后将其作为单向模型处理。这方面的主要效果与2x2布局的整体3自由度模型相同，当然EMMEANS和COMPARE在这方面工作得很好。

在不创建新变量的情况下，可以使用LMATRIX：

 /LMATRIX "(1,1) - (2,2)" var1 1 -1 var2 1 -1 var1*var2 1 0 0 -1
 /LMATRIX "(1,2) - (2,1)" var1 1 -1 var1 -1 1 var1*var2 0 1 -1 0

被引用的部分是标签，表示2x2设计中正在比较的单元格。

您可以使用的另一个技巧使指定LMATRIX变得更简单，但不创建新的变量，那就是只使用交互项指定设计，并取消截取。这使得参数估计只有四个单元格表示：

UNIANOVA Y BY var1 var2
/INTERCEPT=EXCLUDE
/DESIGN var1*var1
/LMATRIX "(1,1) - (2,2)" var1*var2 1 0 0 -1
/LMATRIX "(1,2) - (2,1)" var1*var1 0 1 -1 0.

在这种情况下，ANOVA表中显示的一个效应是4df效应，对0进行所有均值测试，因此它不是您感兴趣的，但您想要的比较很容易获得。请注意，此技巧仅适用于不将参数重新参数化为满秩的过程。