[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