Iterative detection procedure for multiple interventions of unknown types occuring at unknown times as proposed by Fokianos and Fried (2010, 2012).

1 2 3 4 5 |

`fit` |
an object of class |

`taus` |
integer vector of times which are considered for the possible intervention to occur. Default is to consider all times. |

`deltas` |
numeric vector that determines the types of intervention to be considered (see Details). |

`external` |
logical value specifying wether the interventions effect is external or not (see Details). |

`B` |
positive integer value giving the number of bootstrap samples for estimation of the p-value. |

`signif_level` |
numeric value with |

`start.control_bootstrap` |
named list that determines how to make initial estimation in the bootstrap, see argument |

`final.control_bootstrap` |
named list that determines how to make final maximum likelihood estimation in the bootstrap, see argument |

`inter.control_bootstrap` |
named list determining how to maximise the log-likelihood function in an intermediate step, see argument |

`parallel` |
logical value. If |

`...` |
additional arguments passed to the function for detection of single intervention effects |

This function performs an iterative procedure for detection of multiple intervention effects. In each step the function `interv_detect`

is applied for each of the possible intervention types provided in the argument `deltas`

. If there is (after a Bonferroni correction) no significant intervention effect the procedure stops. Otherwise the type of intervention with the minimum p-value is chosen. In case of equal p-values preference is given to a level shift (i.e. *δ=1*) and then to the type of intervention with the largest test statistic. The effect of the chosen intervention is removed from the time series. The time series cleaned from the intervention effect is tested for further interventions in a next step.

For each time in `taus`

the test statistic of a score test on an intervention effect occuring at that time is computed, see `interv_test`

. The time with the maximum test statistic is considered as a candidate for a possible intervention effect at that time. The type of the intervention effect is specified by `delta`

as described in `interv_covariate`

. The intervention is included as an additional covariate according to the definition in `tsglm`

. It can have an internal (the default) or external (`external=TRUE`

) effect (see Liboschik et al., 2014).

All p-values given in the output are multiplied by the number of intervention types considered to account for the multiple testing in each step by a Bonferroni correction. Note that this correction can lead to p-values greater than one.

Note that this bootstrap procedure is very time-consuming.

An object of class `"interv_multiple"`

, which is a list with the following components:

`interventions` |
data frame giving the detected interventions, which has the variables |

`fit_H0` |
object of class |

`fit_cleaned` |
object of class |

`model_interv` |
model specification of the model with all detected interventions at their respective times. |

`fit_interv` |
object of class |

`track` |
named list of matrices with the detailed results of the iterative detection procedure. Element |

Tobias Liboschik, Philipp Probst, Konstantinos Fokianos and Roland Fried

Fokianos, K. and Fried, R. (2010) Interventions in INGARCH processes. *Journal of Time Series Analysis* **31(3)**, 210–225, http://dx.doi.org/10.1111/j.1467-9892.2010.00657.x.

Fokianos, K., and Fried, R. (2012) Interventions in log-linear Poisson autoregression. *Statistical Modelling* **12(4)**, 299–322. http://dx.doi.org/10.1177/1471082X1201200401.

Liboschik, T., Kerschke, P., Fokianos, K. and Fried, R. (2013) Modelling interventions in INGARCH processes. *SFB 823 Discussion Paper* **03/13**, http://hdl.handle.net/2003/29878.

`tsglm`

for fitting a GLM for time series of counts.
`interv_test`

for testing for intervention effects and `interv_detect`

for detection of single interventions of given type. `interv_covariate`

for generation of deterministic covariates describing intervention effects.

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 | ```
## Not run:
###Campylobacter infections in Canada (see help("campy"))
#Searching for potential intervention effects (runs several hours!):
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), B=500, signif_level=0.05)
campyfit_intervmultiple
plot(campyfir_intervmultiple)
#Parallel computation for shorter run time on a cluster:
library(parallel)
ntasks <- 3
clust <- makeCluster(ntasks)
setDefaultCluster(cl=clust)
interv_multiple(fit=campyfit, taus=80:120, deltas=c(0,0.8,1), B=500,
signif_level=0.05, parallel=TRUE)
## End(Not run)
``` |

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

All documentation is copyright its authors; we didn't write any of that.