## Examples from randomForestSRC::rfsrc
##
data(pbc, package = "randomForestSRC")
pbc_rf <- rfsrc(Surv(days, status) ~ ., pbc, nsplit = 10)
pbc_vs <- var.select(pbc_rf)
xvar <- pbc_vs$topvars[1:6]
pbc_prtl_time <- mclapply(c(1,3,5), function(tm){
plot.variable(pbc_rf, surv.type = "surv",
time = tm,
xvar.names = xvar, partial = TRUE,
show.plots = FALSE)
})
save(pbc_prtl_time, file="data/pbc_prtl_time.rda", compress="xz")
# Format for 1 year survival partial plot for ggRandomForests
gg_dta <- gg_partial(pbc_prtl_time[[1]])
# Without specifying which plot, plot.gg_partial returns a list of all four
# figures.
plot(gg_dta)
plot.gg_partial(gg_dta[[1]])
plot.gg_partial(gg_dta[[4]])
# Or a panel of the batch
plot(gg_dta, panel=TRUE)
Run the code above in your browser using DataLab