Message-ID: <40e66e0b0712030600r1856786fy2a38e858d66b242d@mail.gmail.com>
Date: 2007-12-03T14:00:58Z
From: Douglas Bates
Subject: glmer != lmer(,..family=binomioal)
In-Reply-To: <loom.20071203T135043-886@post.gmane.org>
Can you tell us the version of the lme4 package that you are using, please?
On Dec 3, 2007 7:54 AM, Dieter Menne <dieter.menne at menne-biomed.de> wrote:
> 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)
>
> _______________________________________________
> R-sig-mixed-models at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
>