library(survival)
library(survex)
cph <- coxph(Surv(time, status) ~ ., data = veteran, model = TRUE, x = TRUE, y = TRUE)
rsf_ranger <- ranger::ranger(Surv(time, status) ~ .,
data = veteran,
respect.unordered.factors = TRUE,
num.trees = 100,
mtry = 3,
max.depth = 5
)
cph_exp <- explain(cph)
rsf_ranger_exp <- explain(rsf_ranger,
data = veteran[, -c(3, 4)],
y = Surv(veteran$time, veteran$status)
)
predict(cph_exp, veteran[1, ], output_type = "survival")[, 1:10]
predict(cph_exp, veteran[1, ], output_type = "risk")
predict(rsf_ranger_exp, veteran[1, ], output_type = "chf")[, 1:10]
Run the code above in your browser using DataLab