Message-ID: <loom.20071203T135043-886@post.gmane.org>
Date: 2007-12-03T13:54:31Z
From: Dieter Menne
Subject: 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)