Question

我进行了双向2x4方差分析，发现交互项是显着的。现在我想对8个细胞进行Tukey的HSD事后检测。我在R中使用了TukeyHSD（）函数，它在Q表中使用k为8进行所有28个可能的单元到单元组合。

比较相同列或行中的单元格与我相关，而对角线比较则不然。所以会有16个比较而不是28个。我看到有几个消息来源称这些无混淆与混淆的比较，他们建议只在Q表（https://ocw.mit.edu/courses/brain-and-cognitive-sciences/9-07-statistical-methods-in-brain-and-cognitive-science-spring-2004/lecture-notes/21_anova4.pdf和{{进行查找的调整k值进行无混淆的比较。 3}}）

我想知道这是否是一个有效的事情，如果有的话，如果有一些方法在R中方便地做。

Answer 1

您可能希望使用 lsmeans 包。你可以这样做：

library(lsmeans)
( a.lsm = lsmeans(model, ~ a | b) )
pairs(a.lsm)

......同样的事情扭转了这些因素。如果您希望将所有比较视为一个系列，请使用

rbind(pairs(a.lsm), pairs(b.lsm))

对同步测试进行多元$ t $校正 - 这是正确的一步校正。