Lee, I suggest you look into generalized estimating equations, implemented by package geepack in R. I believe you can incorporate an AR1 correlation structure within a Poisson framework via the argument 'corstr="ar1"' in function geeglm.
best, Nathaniel
