data(DataSet8.1)
DataSet8.1$block <- factor(x = DataSet8.1$block)
DataSet8.1$set <- factor(x = DataSet8.1$set)
DataSet8.1$trt <- factor(x = DataSet8.1$trt)
library(lmerTest)
Exam8.1Lmer <- lmer(y ~ set + trt %in% set + (1|set/block), DataSet8.1)
summary(Exam8.1Lmer)
anova(Exam8.1Lmer)
library(emmeans)
emmeans(object = Exam8.1Lmer, specs = ~trt|set)
contrast(emmeans(object = Exam8.1Lmer, specs = ~trt|set), method = "pairwise", by = "set")
Run the code above in your browser using DataLab