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
_______________________________________________
R-sig-ecology mailing list
R-sig-ecology@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-ecology
