# NOT RUN {
# }
# NOT RUN {
# Fit example model:
tv <- lmer(Sharpnessofmovement ~ TVset * Picture +
(1 | Assessor:TVset) + (1 | Assessor:Picture) +
(1 | Assessor:Picture:TVset) + (1 | Repeat) + (1 | Repeat:Picture) +
(1 | Repeat:TVset) + (1 | Repeat:TVset:Picture) + (1 | Assessor),
data = TVbo)
# Backward reduce the model:
(st <- step(tv)) # takes ~10 sec to run
# Pairwise comparisons of LS-means for Picture and TVset:
plot(st, which=c("Picture", "TVset"), pairwise = TRUE)
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab