data(dataDIVAT4)
divat.surv <- survfit(Surv(time/365.24, status) ~ hospital, data = dataDIVAT4)
plot(divat.surv, lty = 1:6, col=1:6, lwd=2, mark.time=FALSE,
xlab="Post transplantation time (days)", ylab="Patient and graft survival")
Run the code above in your browser using DataLab