data(Prostate)
summary(Prostate)
cor(Prostate$lpsa, Prostate$lcavol)
if (FALSE) {
fit <- islasso(lpsa ~ ., data = Prostate, family = gaussian())
summary(fit, pval = 0.05)
lambda.aic <- aic.islasso(fit, method = "AIC")
fit.aic <- update(fit, lambda = lambda.aic)
summary(fit.aic, pval = 0.05)
}
Run the code above in your browser using DataLab