我想用rstatix::anova_test()
运行重复测量方差分析。
我的数据是:
tb <- tibble(id = factor(rep(1:100, each = 8)),
area = factor(rep(c(4,2,3,1,5,3,1,5), 100)),
score = round(runif(800, 1, 7), 0)
)
我想用这个公式:anova_test(data = tb, dv = score, wid = id, within = area)
。
对于id
和area
的一些组合,存在多个实例。在这些情况下,我想使用score
的mean
作为dv
的值。
我试着解决这个问题如下:
tbs <- tb %>%
group_by(id, area) %>%
summarise(
mean_score = mean(score, na.rm = TRUE)
)
anova_test(data = tbs, dv = mean_score, wid = id, within = area)
但是,我收到一个错误,指出列id
不存在。显然,列确实存在,并且每行现在都包含id
和area
的唯一组合,那么我在这里做错了什么呢?
在运行anova_test()
之前,我缺少一个%>% ungroup()
。