# # # # # # # # # #
# Short and noisy time series
set.seed(12)
N <- 300
noise <- rnorm(N)
signal <- c(rep(0,200),rep(10,100))
x <- signal+noise
# SCARM Filter
signal.extr <- scarm.filter(x)
par(mfrow=c(2,1))
ts.plot(x, col="darkgrey", main="time series data and signal extraction")
lines(signal.extr$signal)
ts.plot(signal.extr$adapted.width, main="adapted window width")
Run the code above in your browser using DataLab