if (FALSE) {
n <- 100; p <- 30
beta <- c(runif(10, -2, 2), rep(0, p - 10))
sim <- simulXy(n = n, p = p, beta = beta, seed = 1, family = gaussian())
fit <- islasso.path(y ~ ., data = sim$data, family = gaussian())
plot(fit, yvar = "coefficients", gof = "AICc", label = TRUE)
plot(fit, yvar = "se", gof = "AICc")
plot(fit, yvar = "gradient", gof = "AICc")
plot(fit, yvar = "gof", gof = "AICc")
}
Run the code above in your browser using DataLab