glmnet.recovery.RdSimulate from glmnet penalized regression model
glmnet.recovery(y, X, ...)simulated response
Other recovery:
ols.recovery(),
pls.recovery(),
simple.recovery(),
spls.recovery()
set.seed(1)
X = data.frame(matrix(rnorm(100 * 20), 100, 20))
y = X[,1] + rnorm(100)
# simulate from elastic-net regression:
ysim = glmnet.recovery(y=y, X=X)