library("tscount")
###Campylobacter infections in Canada (see help("campy"))
interventions <- interv_covariate(n=length(campy), tau=c(84, 100),
delta=c(1, 0)) #detected by Fokianos and Fried (2010, 2012)
#Linear link function with Negative Binomial distribution:
campyfit <- tsglm(campy, model=list(past_obs=1, past_mean=13),
xreg=interventions, distr="nbinom")
broom::glance(campyfit)
Error: glance doesn't know how to deal with data of class tsglm
Which means working without broom / creating the broom tool for tsglm.
Bad news:
Which means working without
broom
/ creating thebroom
tool fortsglm
.