Re: [R] Error: cannot use PQL when using lmer

From: Ben Bolker <bolker_at_ufl.edu>
Date: Sun, 06 Jul 2008 02:09:11 +0000 (UTC)

 <hpdutra <at> yahoo.com> writes:

>
>
> > library(MASS)
> > attach(bacteria)
> > table(y)
> y
> n y
> 43 177
> > y<-1*(y=="y")
> > table(y,trt)
> trt
> y placebo drug drug+
> 0 12 18 13
> 1 84 44 49
> > library(lme4)
> > model1<-lmer(y~trt+(week|ID),family=binomial,method="PQL")
> Error in match.arg(method, c("Laplace", "AGQ")) :
> 'arg' should be one of “Laplace”, “AGQ”
>

  What is your question?
  Doug Bates warned a few weeks ago that the newer version of lmer would no longer use PQL for GLMMs (he found that it was unreliable, even as a starting method for Laplace fits). I think you can still get the older version if you want it, or you can use glmmPQL from the MASS package (glmmPQL has some advantages anyway).

   It might be better to forward further discussion to r-sig-mixed.

   Ben Bolker



R-help_at_r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Received on Sun 06 Jul 2008 - 02:14:11 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Sun 06 Jul 2008 - 08:32:05 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive