I'm stuck here. How to translate the following command lines to lmer syntax.
model <- aov(
response ~ treatment*time + Error(id/time),
data = selfesteem2)
summary(model)
I've tried this among others without success.
anova(lmerTest::lmer(
response ~ treatment*time + (1|id:time),
data = selfesteem2))
The data:
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)
It's