translate aov to lmer

65 Views Asked by At

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)
1

There are 1 best solutions below

8
Stéphane Laurent On

It's

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