R

Question

我在受试者之间进行了R中的重复测量ANOVA，给出了与SPSS相同的结果。

设计如下：积极性（2个级别）x新近度（2个级别）x组（3个级别）

因此，我有一个三方互动。现在，我想分别针对三组进行所有可能的事后t检验。使用以下代码并使用lsmeans和afex包这是可行的：

anova_all <- aov_ez("code", "err", stim.lng, +
within=.("recency", "positivity"), between=.("group"))
ref <- lsmeans(anova_all,~recency*positivity|group)
comps <- contrast(ref,method="pairwise", adjust="none")
summary(comps)

代码运行，但与SPSS输出相比，我获得了不同的结果（就p值和t值而言）。这让我感到困惑，我相信我应该得到与SPSS类似的结果，因为只进行了配对样本t检验。

我猜关键点可能是两种软件的不同自由度。例如，当我在9个受试者的特定组中进行配对t检验时，在SPSS中，自由度为9，而在R中为89-102（参见附图）。为什么这么高的R？是否还有其他用SP复制SPSS的结果？

这里是我的设计和SPSS和R在特定患者组中的三种不同输出，用于所有可能的成对事后检验：

R的输出：

SPSS的输出：