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