[R-sig-eco] correlated residuals

[Alexandre Villers]

Good morning,

I'm trying to fit a very simple model: the proportion of occupied nest 
boxes according to food availability on a time serie (37 years). The 
problem is that the glm (binomial error family with logit  link) 
presents a marked overdispersion ( ~ 24). Checking the residuals showed 
me they are correlated with a lag of 3 years (preys have cyclic 
population dynamic) and accounting for this temporal correlation will 
certainly improve the model.
I know how to cope with this with glmmPQL (with a corARMA correlation 
structure) but I "cheat" a bit given I don't really need random effect.

Is there a way to specify such correlation structure with gee or geepack 
? Any link to a pdf or webpage explaining the "how to" would be 
appreciated...

Best regards

Alex