Dear lmers,
from reading the docs and the first lines of lmer() code, I got the
understanding that the two formulations below should be equivalent. The first
works as expected, the second gives an envir() error. What's wrong?
Dieter
x=data.frame(resp=rbinom(100,1,0.8),
treat=sample(letters[7:8],100,TRUE),
pat=as.factor(sample(letters[1:4],100,TRUE))
)
glmer(resp~treat+(1|pat),family=binomial,data=x)
lmer(resp~treat+(1|pat),family=binomial,data=x)