Helper function to estimate autoregressive mode
fit_ar(counts, control.dates = NULL, order.max = 5, aic = FALSE, plot = FALSE)Output from `compute_excpected`
Dates to use to estimate covariance
Maximum order of autoregressive process
Logical that determines if AIC is used
logical that determines if an autocorrelation plot is generated for exploration purposes