## veteran data
## randomized trial of two treatment regimens for lung cancer
data(veteran, package = "randomForestSRC")
veteran_rf <- rfsrc(Surv(time, status) ~ ., data = veteran, ntree = 100)
veteran_prtl <- plot.variable(veteran_rf, surv.type = "surv",
partial = TRUE, time=30, xvar.names = "age",
show.plots=FALSE)
gg_dta <- gg_partial(veteran_prtl)
plot(gg_dta)
Run the code above in your browser using DataLab