Nothing
N_TEST <- 20
test_that("Custom threshold is consistent with asymptotic critical value", {
for (i in 1:N_TEST) {
ts <- list(testData(model="blocks"),
testData(model="fms"),
testData(model="mix"),
testData(model="stairs10"),
testData(model="teeth10"))
for (x in ts) {
n <- length(x)
alpha <- runif(1, 0, 1)
G_left <- max(floor(runif(1, 20, 40)), ceiling(runif(1, 0.05*n, 0.2*n)))
G_right <- max(floor(runif(1, 20, 40)), ceiling(runif(1, 0.05*n, 0.2*n)))
# Just re-define the critical value, but in 'custom' branch of code
threshold.custom <- mosum.criticalValue(n, G_left, G_right, alpha)
threshold.function1 <- function(G, n, alpha) { mosum.criticalValue(n,G,G,alpha) }
threshold.function2 <- function(G_l, G_r, n, alpha) { mosum.criticalValue(n,G_l,G_r,alpha) }
mcpts1 <- mosum(x, G_left, G.right=G_right, alpha=alpha)
mcpts2 <- multiscale.bottomUp(x, G=c(G_left,G_right), alpha=alpha)
mcpts3 <- multiscale.localPrune(x, G=c(G_left,G_right), alpha=alpha)
mcpts4 <- mosum(x, G_left, G.right=G_right, alpha=alpha,
threshold="custom",
threshold.custom=threshold.custom)
mcpts5 <- multiscale.bottomUp(x, G=c(G_left,G_right), alpha=alpha,
threshold="custom",
threshold.function=threshold.function1)
mcpts6 <- multiscale.localPrune(x, G=c(G_left,G_right), alpha=alpha,
threshold="custom",
threshold.function=threshold.function2)
expect_equal(mcpts1$cpts, mcpts4$cpts)
expect_equal(mcpts2$cpts, mcpts5$cpts)
expect_equal(mcpts3$cpts, mcpts6$cpts)
}
}
})
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.