Re: [R] lmer question

From: Douglas Bates <bates_at_stat.wisc.edu>
Date: Mon 18 Apr 2005 - 23:31:06 EST

Markus Bean wrote:
> Hi --
>
> I'm using lmer for binomial data. I am trying to
> replicate estimates provided by Agresti (2002,
> Categorical data analysis, Wiley) using abortion data
> in table 10.13 (estimates provided in table 12.3 p.
> 505).
>
> I fit the same model using these three commands:
>
> a1 <- lmer(resp ~ sex + option1 + option2 + (1|id),
> data=abort,family=binomial, method = c("AGQ"))
>
> a2 <- lmer(resp ~ sex + option1 + option2 + (1|id),
> data=abort,family=binomial, method = c("Laplace"))
>
> a3 <- lmer(resp ~ sex + option1 + option2 + (1|id),
> data=abort,family=binomial, method = c("PQL"))
>
> All three methods provide the exact same estimates
> (which should not be the case), and the estimates are
> incorrect. I know the data are correctly entered
> because I obtain correct estimates with gllamm in
> Stata.
>
> I am I doing something wrong here in my commands, or
> is the lmer module not implementing AGQ and Laplace
> properly with this version?

The latter. Actually it is not implementing them at all in this version. I will upload a new release of the lme4 package later this week and insert appropriate messages.



R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Mon Apr 18 23:37:15 2005

This archive was generated by hypermail 2.1.8 : Fri 03 Mar 2006 - 03:31:15 EST