data(dataKi67)
# Kaplan Meier estimations of graft survivals in Wintzer et al. (1991)
plot(dataKi67$year[dataKi67$study.num==1],
dataKi67$surv[dataKi67$study.num==1],
xlab="Post transplantation time (years)",
ylab="Graft survival", ylim=c(0.6,1), xlim=c(0, 4), type="n")
# Goup 1
lines(c(0, dataKi67$year[dataKi67$study.num==1 & dataKi67$classe==1]),
c(1, dataKi67$surv[dataKi67$study.num==1 & dataKi67$classe==1]),
type="b", col=1, lty=1, lwd=2)
# Goup 2
lines(c(0, dataKi67$year[dataKi67$study.num==1 & dataKi67$classe==2]),
c(1, dataKi67$surv[dataKi67$study.num==1 & dataKi67$classe==2]),
type="b", col=2, lty=2, lwd=2)
# legend
legend("bottomleft", c("group #1 (log Ki67 < 2.49)",
"group #2 (log Ki67 > 2.49)"), col=c(1, 2), lty=c(1, 2), lwd=c(2, 2))
Run the code above in your browser using DataLab