[R] glmmPQL:

[Spencer Graves]

      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