# NOT RUN {
n <- 200
signal <- sin(2 * pi * 1:n / n)
signal[51:100] <- signal[51:100] + 5
signal[151:200] <- signal[151:200] + 5
y <- rnorm(n) + signal
est <- BinSegBstrap(y = y)
plot(y)
lines(signal)
lines(est$est, col = "red")
n <- 100
signal <- sin(2 * pi * 1:n / n)
signal[51:100] <- signal[51:100] + 5
y <- rnorm(n) + signal
test <- BstrapTest(y = y)
est <- estimateSingleCp(y = y)
plot(y)
lines(signal)
lines(est$est, col = "red")
# }
Run the code above in your browser using DataLab