# \donttest{
if (requireNamespace("shapr") & requireNamespace("mlbench")) {
library(shapr)
# Boston housing dataset
library(mlbench)
data(BostonHousing2)
dat <- BostonHousing2
y <- dat$cmedv
x <- subset(dat, select = -c(cmedv, medv, town, chas))
# Fit a glmnet model using nested CV
# Only 3 outer CV folds and 1 alpha value for speed
fit <- nestcv.glmnet(y, x, family = "gaussian", n_outer_folds = 3, alphaSet = 1)
# Generate SHAP values using shapr
sh <- nestcv.explain(fit, pred_nestcv_glmnet)
# Plot overall variable importance
plot_shap_bar(sh, x)
# Plot beeswarm plot
plot_shap_beeswarm(sh, x, size = 1)
}
# }
Run the code above in your browser using DataLab