[R] lmer question

From: Markus Bean <markusbean_at_yahoo.com>
Date: Mon 18 Apr 2005 - 16:10:01 EST


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?



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 16:15:43 2005

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