> Up until now my procedure is to evaluate two models, viz. > > m0 <- lmer(response~(1|subject),data) > m1 <- lmer(response~gender+(1|subject),data) > > The likelihood-ratio test approach directly compares these two. Since these models differ in their fixed effects, you need REML=FALSE for the LRT to be meaningful.