I'm getting a strange error from glmmPQL. Consider the following sample code: set.seed(8) N. <- 1000 z <- rnorm(N.) pr.good <- exp(-1e-4*exp(2+2*z)) quantile(pr.good) DF. <- data.frame(yield=rbinom(N., N., pr.good)/N., Offset=rep(-10, N.), nest=1:N.) fit <- glmmPQL(fixed=1-yield~offset(Offset), random=~1|nest, family=binomial(link="cloglog"), data=DF., weights=rep(N., N.)) Using R 1.9.0pat, Windows 2000, with MASS updated just a few hours ago, I get the following error: iteration 1 iteration 2 Error in logLik.reStruct(object, conLin) : NA/NaN/Inf in foreign function call (arg 3) Any suggestions? Best Wishes, Spencer Graves