[R-sig-ME] glmer != lmer(,..family=binomioal)

Dieter Menne dieter.menne at menne-biomed.de
Mon Dec 3 14:54:31 CET 2007


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)




More information about the R-sig-mixed-models mailing list