Thanks for providing a partially reproducible example. I believe the
error message you cite came from lme. I say this, because I modified
your call to glmmPQL2 to call lme and got the following:
library(nlme)
fit.lme - lme(y ~ trt + I(week 2), random = ~ 1 | ID,
+ data
I just identified an error in my recent post on this subject: There
is a very good reason that Venables Ripley's glmmPQL did NOT include
an argument like the weights.lme in the glmmPQL. included in my
recent post: Their function calls glm first and then provides weights
computed
Dear listers,
On the line of a last (unanswered) question about glmmPQL() of the
library MASS, I am still wondering if it is possible to pass a variance
structure object to the call to lme() within the functions (e.g.
weights=varPower(1), etc...). The current weights argument of glmmPQL is
Have you received a reply to this post? I haven't seen one. I don't
have an answer for you, but if you'd still like help from this list, I
suggest you prepare the simplest possible toy example that you can
conceive and send it to this list, restating your question in terms of
that
Dear listers,
glmmPQL (package MASS) is given to work by repeated call to lme. In the
classical outputs glmmPQL the Variance Structure is given as fixed
weights, Formula: ~invwt. The script shows that the function
varFixed() is used, though the place where 'invwt' is defined remains