#import lme4 package and lmerTest package
library(lmerTest)
#import TVbo data from lmerTest package
data(TVbo)
m <- lmer(log(Coloursaturation)~TVset*Picture+
(1|Assessor)+(1|Assessor:TVset), data=TVbo)
step(m, test.effs="TVset", reduce.fixed=TRUE, reduce.random=TRUE)
Run the code above in your browser using DataLab