Skip to content

plotting glmmPQL function

2 messages · rhelpless, Dieter Menne

rhelpless
# 
hello all,

i'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...  i know this should be basic, but i can't seem to
figure it out.  many thanks in advance.

j
Dieter Menne
#
rhelpless wrote:
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.


Dieter