自举重复测量 (ezANOVA)答案

【问题标题】：Bootstrapping repeated measure (ezANOVA)自举重复测量 (ezANOVA)
【发布时间】：2019-12-02 13:18:06
【问题描述】：

我想使用 ez 包引导重复测量 anova 对象。

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“风格”/ezANOVA 环境中。

【问题讨论】：

标签： r lme4 anova

【解决方案1】：

我目前找到的最接近的答案是使用 afex 包中的 mixed 函数。 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

但是，我很乐意听到其他选项。

【讨论】：