生存曲線の信頼区間

PART2 CHAPTER5のシミュレーション。
 

library(survival)
data<- matrix(c(4.07,6.54,1.39,6.17,5.89,4.76,3.67,1,0,1,0,1,1,0),nc=2)
res<- survfit(Surv(data[,1],data[,2])~rep(1,nrow(data)),type="kaplan-meier")
plot(res,frame=FALSE,xlab="time",ylab="Over all survival")