#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
lsmeans(m1, test.effs="Gender:Information")
#import TVbo data from lmerTest package
data(TVbo)
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