if (FALSE) {
set.seed(1)
n <- 100
p <- 100
p1 <- 20 #number of nonzero coefficients
coef.veri <- sort(round(c(seq(.5, 3, l=p1/2), seq(-1, -2, l=p1/2)), 2))
sigma <- 1
coef <- c(coef.veri, rep(0, p-p1))
X <- matrix(rnorm(n*p), n, p)
mu <- drop(X%*%coef)
y <- mu + rnorm(n, 0,sigma)
lambda <- 2
o <- islasso(y ~ ., data = data.frame(y = y, X),
family = gaussian(), lambda = lambda)
plot(o)
}
Run the code above in your browser using DataLab