我被困在这里了如何将以下命令行转换为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)