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