if (FALSE) {
data("pbc",package='survival')
pbc$status[pbc$status>=1]=pbc$status[pbc$status>=1]-1
pbc$id=NULL
fctrs<-c('trt','ascites','spiders','edema','hepato','stage')
for(f in fctrs)pbc[[f]]=as.factor(pbc[[f]])
pbc=na.omit(pbc)
orsf=ORSF(data=pbc,ntree=30)
times=seq(365, 365*4,length.out = 10)
predict(orsf,newdata=pbc[c(1:5),],times=times)
}
Run the code above in your browser using DataLab