\donttest{
set.seed(1)
n <- 480
p_true <- 6
p <- 50
x <- mvtnorm::rmvnorm(n, rep(0, p), diag(p))
theta_0 <- rbind(
  runif(p_true, -5, -2),
  runif(p_true, -3, 3),
  runif(p_true, 2, 5),
  runif(p_true, -5, 5)
)
theta_0 <- cbind(theta_0, matrix(0, ncol = p - p_true, nrow = 4))
y <- c(
  x[1:80, ] %*% theta_0[1, ] + rnorm(80, 0, 1),
  x[81:200, ] %*% theta_0[2, ] + rnorm(120, 0, 1),
  x[201:320, ] %*% theta_0[3, ] + rnorm(120, 0, 1),
  x[321:n, ] %*% theta_0[4, ] + rnorm(160, 0, 1)
)
small_lasso_data <- cbind.data.frame(y, x)
result_no_vp <- fastcpd.lasso(
  small_lasso_data,
  beta = "BIC",
  cost_adjustment = NULL,
  pruning_coef = 0
)
summary(result_no_vp)
result_20_vp <- fastcpd.lasso(
  small_lasso_data,
  beta = "BIC",
  cost_adjustment = NULL,
  vanilla_percentage = 0.2,
  pruning_coef = 0
)
summary(result_20_vp)
}
