The count data themselves are obviously non-negative, & due in part to large numbers of 0s, are highly overdispersed (variance >> mean), so was leaning towards using negative binomial (NB) GLM, or perhaps zero-inflated/adjusted (hurdle) approach. Because the depth parameter itself varies as function of space & time (my other factors of interest), I’m having some difficulty understanding how to incorporate it into a model. Can I use depth as a continuous variable or should I aggregate into categories, & in either case, how to deal w/ the fact that amongst other categorical factors of interest, years (some years wetter/drier than others, etc.) & regions (again, plots in south being deeper than north) the distribution of depths vary so much, providing weak overlap?

Then there is the issue of non-independence in the dataset that could be contributing to the overdispersion as well. As I’ve described above, the unit of interest is a given week’s survey at a given plot. This design confounds the spatial replication (n=3 plots/region ea wk) w/ the temporal replication (n=many wks of survey data for ea plot ea year). Is there a way to incorporate a “repeated measures” type of analysis for ea plot’s data in a given year into a Poisson, NB, or appropriate ZIP/ZINB, ZAP/ZANB model? While birds do move a fair amount from week to week, seems it's important to indirectly account for structural characteristics of the plots (reduced variability intra-plot vs inter-plot amongst wks of data) that such an analysis could provide.

I'll be running the models in SAS & sincerely appreciate any recommendations on these matters.