将 aov 翻译为 lmer

Question

我被困在这里了。如何将以下命令行翻译为 lmer 语法。

model <- aov(
  response ~ treatment*time + Error(id/time),
  data = selfesteem2)
summary(model)

我尝试过其他方法但没有成功。

anova(lmerTest::lmer(
  response ~ treatment*time + (1|id:time),
  data = selfesteem2))

---

数据：

library(tidyverse)
data("selfesteem2", package = "datarium")
selfesteem2 <- selfesteem2 %>% # long format
  pivot_longer(c("t1","t2","t3"), values_to = "score", names_to = "time") %>% 
  mutate(time = factor(time)) %>% 
  rename(response = score)

Answer 1

是

anova(lmerTest::lmer(
  response ~ treatment*time + (1|id/time),
  data = selfesteem2))