I am using the lmer function from the lme4 package. I wrote the following statement, specifying the method to be adaptive Gaussian quadrature. I am getting an error saying "method = "AGQ" not yet implemented for supernodal representation". Please help. > fit <- lmer(response~beta1+(1|patient),family=binomial,method="AGQ",data=data.2) Boikanyo.