Hi, I'm having trouble with glmmPQL from the MASS package. I'm trying to fit a model with a binary response variable, two fixed and two random variables (nested), with a sample of about 200,000 data points.
Unfortunately, I'm getting an error message that is difficult to understand without knowing the internals of the glmmPQL function. > model <- glmmPQL(primed ~ log(dist) * role , random = ~ dist | > target.utt / prime.utt , family=binomial(link = "logit"), > data=data.utts, niter=5, verbose = TRUE) > Loading required package: nlme > iteration 1 > iteration 2 > iteration 3 > Error in solve.default(pdMatrix(a, fact = TRUE)) : > system is computationally singular: reciprocal condition > number = 8.65949e-32 > In addition: Warning messages: > 1: Singular precision matrix in level -1, block 4 > 2: Singular precision matrix in level -1, block 4 > 3: Singular precision matrix in level -1, block 4 > 4: Singular precision matrix in level -1, block 4 > 5: Singular precision matrix in level -1, block 4 Any suggestions? Will a larger dataset (possible) solve the problem? Thanks David -- David Reitter - ICCS/HCRC, Informatics, University of Edinburgh Blog: http://www.davids-world.com Homepage: http://www.david- reitter.com ______________________________________________ [email protected] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
