如何在 SPSS 中比较所有可能的组组合与 EMMEANS？
How to compare all possible group combinations with EMMEANS in SPSS?
假设您有一个 2x2 设计，并且您正在使用 SPSS 中的 ANOVA 测试这 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
我如何计算这些比较？
UNIANOVA 中的 EMMEANS 不提供像这样的交互中单元格之间的所有成对比较。还有一些其他程序，例如 GENLIN，确实提供这些，但使用大样本卡方统计而不是 t 或 F 统计。在 UNIANOVA 中，您可以使用 LMATRIX 子命令获得这些，或者您可以使用 EMMEANS 的一些技巧。
对于 EMMEANS 的诡计，创建一个具有四个水平的单一因子来索引 2x2 单元格布局，然后将其作为单向模型处理。其主要效果与 2x2 布局的整体 3 自由度模型相同，当然，带有 COMPARE 的 EMMEANS 在这方面效果很好。
无需创建新变量，您可以将 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 但无需创建新变量的另一个技巧是仅使用交互项指定 DESIGN 并抑制截距。这使得参数估计只有四个单元格均值：
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.
在这种情况下，方差分析 table 中显示的一个效应是 4 df 效应，针对 0 测试所有均值，因此它没有意义，但您想要的比较很容易获得。请注意，此技巧仅适用于未重新参数化为满秩的过程。
假设您有一个 2x2 设计，并且您正在使用 SPSS 中的 ANOVA 测试这 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
我如何计算这些比较？
UNIANOVA 中的 EMMEANS 不提供像这样的交互中单元格之间的所有成对比较。还有一些其他程序，例如 GENLIN，确实提供这些，但使用大样本卡方统计而不是 t 或 F 统计。在 UNIANOVA 中，您可以使用 LMATRIX 子命令获得这些，或者您可以使用 EMMEANS 的一些技巧。
对于 EMMEANS 的诡计，创建一个具有四个水平的单一因子来索引 2x2 单元格布局，然后将其作为单向模型处理。其主要效果与 2x2 布局的整体 3 自由度模型相同，当然，带有 COMPARE 的 EMMEANS 在这方面效果很好。
无需创建新变量，您可以将 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 但无需创建新变量的另一个技巧是仅使用交互项指定 DESIGN 并抑制截距。这使得参数估计只有四个单元格均值：
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.
在这种情况下，方差分析 table 中显示的一个效应是 4 df 效应，针对 0 测试所有均值，因此它没有意义，但您想要的比较很容易获得。请注意，此技巧仅适用于未重新参数化为满秩的过程。