Software accompanying
H. Cho and P. Fryzlewicz (2021) "Multiple change point detection under serial dependence: wild contrast maximisation and gappy Schwarz algorithm".
- The main routine of WCM.gSa (
wcm.gsa) is contained in main.R. - Codes used for analysing nitrogen oxides concentrations in London (contained in
AirQualityData.csv) are provided in air_data.R. - Codes used for analysing Hadley Centre central England temperature data (contained in
hadcet.csv) are provided in temp_data.R.
To use WCM.gSa, do the following:
- Source main.R into
R. - Read the descriptions for
wcm.gsawithin main.R.
For example,
source("main.R")
set.seed(111)
f <- rep(c(0, 5, 2, 8, 1, -2), c(100, 200, 200, 50, 200, 250))
x <- f + arima.sim(list(ar = c(.75, -.5), ma = c(.8, .7, .6, .5, .4, .3)), n = length(f), sd = 1)
wcm.gsa(x)
If you have any questions, contact haeran.cho@bristol.ac.uk