repaired icens_prodlim.c in version 1.1.0 todo:
library(prodlim) dat <- SimSurv(100,cens.baseline=1/100,surv.baseline=1/100) fit <- prodlim(Hist(time,status)~1,data=dat) quantile(fit) fit <- prodlim(Hist(time,status)~X2,data=dat) quantile(fit)
library(prodlim) library(pseudo) dat=SimSurv(4000) x=prodlim(Hist(time,status)~1,data=dat) system.time({u=jackknife(x)}) system.time({v=pseudosurv(time=dat$time,event=dat$status,tmax=unique(dat$time))})