set.seed (1111)
n <- 200
p <- 5
X <- matrix (rnorm (n * p), ncol = p)
X[,2] <- X[,1] + rnorm (n, sd = 0.1)
X[,3] <- X[,1] + rnorm (n, sd = 0.1)
true.beta <- c (1, 2, 0, 0, -1)
y <- drop (X %*% true.beta) + rnorm (n)
control.obj <- lqa.control (max.steps = 200, conv.eps = 1e-3,
conv.stop = FALSE)
obj <- lqa (y ~ X, family = gaussian (), penalty = lasso (1.5),
control = control.obj)
obj$coef
Run the code above in your browser using DataLab