require(survival)
data("bladder4state")
obj<- multidf(gap1=bladder4state$y1, event1=bladder4state$d1,
gap2=bladder4state$y2, status=bladder4state$d2,
size=bladder4state$size)
obj2<-obj[[1]]
KM(time = obj2$time, status = obj2$status, t = 20)
fit <- survfit(Surv(obj2$time, obj2$status) ~ 1, data = obj2)
summary(fit, time = 20)$surv
Run the code above in your browser using DataLab