Predict survival based on rpart tree object
Combines the rpart result with a stratified Kaplan-Meier (prodlim) to predict survival
pecRpart(formula, data, ...)
formula |
passed to rpart |
data |
passed to rpart |
... |
passed to rpart |
list with three elements: ctree and call
library(prodlim) if (!requireNamespace("rpart",quietly=TRUE)){ library(rpart) library(survival) set.seed(50) d <- SimSurv(50) nd <- data.frame(X1=c(0,1,0),X2=c(-1,0,1)) f <- pecRpart(Surv(time,status)~X1+X2,data=d) predictSurvProb(f,newdata=nd,times=c(3,8)) }
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.