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.
Dieter
--
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.