[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