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.
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
