Skip to content
Prev 450 / 20628 Next

glmer != lmer(,..family=binomioal)

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)