data(DataSet8.4)
DataSet8.4$block <- factor(x = DataSet8.4$block)
DataSet8.4$a <- factor(x = DataSet8.4$a)
DataSet8.4$b <- factor(x = DataSet8.4$b)
library(lmerTest)
Exam8.4lmer <-
lmer(
y ~ a + b %in% a +
(1|block) + (1|block:a) + (1|block:b)
, data = DataSet8.4
)
anova(Exam8.4lmer, ddf = "Kenward-Roger")
library(emmeans)
emmeans(object = Exam8.4lmer, specs = ~a|b)
Run the code above in your browser using DataLab