> When trying to fit a generalized linear mixed model using glmmPQL: > >> fit0 <- glmmPQL(ifelse(response=="A",1,0)~gender,data=set1, > random=~1|subject,family=binomial) > iteration 1 > Error in solve.default(pdMatrix(a, fact = TRUE)) : > Lapack routine dgesv: system is exactly singular > > Could this be occuring because the paired responses for each subject are > always the same? If this were the case, what would be the best way to > handle this situation? > > I replaced the response variable with fake data that are not always the > same for each pair and then I don't get this (or any) error message. > >
Sorry, I forgot that a variable like gender will be the same for both responses. Rick B. ______________________________________________ [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
