# NOT RUN {
## import lme4 package and lmerTest package
library(lmerTest)
## specify lmer model
m1 <- lmer(Informed.liking ~ Gender*Information +(1|Consumer), data=ham)
## calculate least squares means for interaction Gender:Information
lsmeansLT(m1, test.effs="Gender:Information")
m <- lmer(Coloursaturation ~ TVset*Picture + (1|Assessor), data=TVbo)
plot(lsmeansLT(m))
lsmeansLT(m, test.effs="TVset")
# }
Run the code above in your browser using DataLab