# NOT RUN {
###Campylobacter infections in Canada (see help("campy"))
campyfit <- tsglm(ts=campy, model=list(past_obs=1, past_mean=c(7,13)))
campyfit_intervmultiple <- interv_multiple(fit=campyfit, taus=80:120,
deltas=c(0,0.8,1), external=FALSE, B=2,
signif_level=0.05) #runs several hours!
plot(campyfit_intervmultiple)
# }
Run the code above in your browser using DataLab