[R] plotting glmmPQL function

Dieter Menne dieter.menne at menne-biomed.de
Tue Sep 2 12:13:15 CEST 2008

rhelpless wrote:
> 'm an R newbie struggling a bit with the glmmPQL function (from the nlme
> pack).  i think i've managed to run the model successfully, but can't seem
> to plot the resulting function. plot(glmmPQL(...)) plots the residuals
> rather than the model

Use predict on the fit; you should probably convert the plus-minus infinity
range by using logistic transformation for plotting.

glmm = glmmPQL(y ~ trt + I(week > 2), random = ~ 1 | ID,
                family = binomial, data = bacteria)
bacteria$pred = predict(glmm)

Creating a special rhelpless mail account is not considered good etiquette
here, so better use your real name if you want a reply.


View this message in context: http://www.nabble.com/plotting-glmmPQL-function-tp19266698p19267495.html
Sent from the R help mailing list archive at Nabble.com.

More information about the R-help mailing list