if (requireNamespace("survival", quietly = TRUE)) {
shapboost <- SHAPBoostSurvival$new(
max_number_of_features = 1,
evaluator = "coxph",
metric = "c-index",
verbose = 0,
xgb_params = list(
objective = "survival:cox",
eval_metric = "cox-nloglik"
)
)
X <- as.data.frame(survival::gbsg[, -c(1, 10, 11)])
y <- as.data.frame(survival::gbsg[, c(10, 11)])
subset <- shapboost$fit(X, y)
}
Run the code above in your browser using DataLab