Description Usage Arguments Value Author(s) References See Also Examples
A plot of segments estimated by the backward procedure.
1 2 |
x |
bwd object |
y |
observed data |
... |
graphical parameters |
plot of estimated segments
Seung Jun Shin, Yicaho Wu, Ning Hao
Shin, Wu, and Hao (2018+) A backward procedure for change-point detection with applications to copy number variation detection, arXiv:1812.10107.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 | # simulated data
set.seed(1)
n <- 1000
L <- 10
mu0 <- -0.5
mu <- rep(mu0, n)
mu[(n/2 + 1):(n/2 + L)] <- mu0 + 1.6
mu[(n/4 + 1):(n/4 + L)] <- mu0 - 1.6
y <- mu + rnorm(n)
alpha <- c(0.01, 0.05)
# BWD
obj1 <- bwd(y, alpha = alpha)
# Modified for epidemic changes with a known basline mean, mu0.
obj2 <- bwd(y, alpha = alpha, mu0 = 0)
par(mfrow = c(2,1))
plot(obj1, y)
plot(obj2, y)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.