[R-sig-ME] offset in glmmPQL
Ben Bolker
bolker at ufl.edu
Fri Feb 19 16:47:20 CET 2010
It's surprising that your example doesn't work, because this example
(which is a hidden example in the ?glmmPQL help page - it doesn't appear
when you type ?glmmPQL, but it gets run when you say example(glmmPQL))
does ...
library(MASS)
library(nlme)
summary(glmmPQL(y ~ trt + week + offset(week), random = ~ 1 | ID,
family = binomial, data = bacteria))
Can you provide a reproducible example of what doesn't work?
Maarten de Groot wrote:
> Dear all,
>
> I have a small question about including an offset value into a glmmPQL
> with a poisson distribution. I want to compare the number of songs of an
> animal between treatments. Because the experiment time is not always the
> same I want to include this as an offset in my formula
> ("mod<-glmmPQL(nosongs~offset(experiment.time)+treatment,random=~1|block/malenr,family=poisson,data=mcrs)").
>
>
> However R seems not to recognize the command"offset" when using a glmmPQL.
>
> Hopefully some one can help me.
>
> Kind regards,
>
> Maarten
>
> _______________________________________________
> R-sig-mixed-models at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
--
Ben Bolker
Associate professor, Biology Dep't, Univ. of Florida
bolker at ufl.edu / people.biology.ufl.edu/bolker
GPG key: people.biology.ufl.edu/bolker/benbolker-publickey.asc
More information about the R-sig-mixed-models
mailing list