I used the following syntax for the mixed model and then step but it did not work.
Does it normally work like this or I actually can not use backward elimination with lmer? Thanks!
fullmodel<-lmer(Eeff~NDF+ADF+CP+NEL+DMI+FCM + (1|Study),data=na.omit(phuong))
step(fullmodel, direction = "backward", trace=FALSE )
You could do it, just not with the step function. Since your model is just additive it shouldn’t take that long to do by hand.