Question

我想使用ez软件包来引导重复测量的方差分析对象。

library(ez)    

ezANOVA(data = ANT, 
        dv = rt,
        wid = subnum,
        between = group,
        within = flank, detailed=FALSE)$ANOVA

输出：

       Effect DFn DFd          F            p p<.05       ges
2       group   1  18   16.22307 7.894348e-04     * 0.2428108
3       flank   2  36 1740.31250 1.524670e-36     * 0.9841982
4 group:flank   2  36   13.01761 5.571290e-05     * 0.3178194

样本数据：

data(ANT)
summary(ANT)

我知道使用lme4::lmer()和parameters::parameters(model, bootstrap = TRUE)可以得到相同的结果。但是，我希望尽可能多地留在ANOVA“ style” / ezANOVA环境中。

Answer 1

我现在找到的最接近的答案是使用mixed包中的afex函数。 nsim是要引导的样本数，seed的行为类似于set.seed()，以提高可重复性。

library(afex)

mixed(
  formula = rt ~ group*flank + (1 | subnum),
  data=ANT,
  method = "PB",
  args_test = list(nsim = 1000, seed = 123), 
  test_intercept = TRUE,
  type = "II")

结果：

Model: rt ~ group * flank + (1 | subnum)
Data: ANT
       Effect df       Chisq p.value
1       group  1   10.93 ***   .0010
2       flank  2 1852.17 ***   .0010
3 group:flank  2   16.36 ***   .0010
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘+’ 0.1 ‘ ’ 1

但是，我很高兴听到其他选择。

自举重复测量（ezANOVA）

1 个答案: