data(ipd.data)
trt <- survfit(Surv(time,event)~strata(group),
ipd.data,subset=trt==1)
ctrl <- survfit(Surv(time,event)~strata(group),
ipd.data,subset=trt==0)
par(mfrow=c(1,2),pty="s")
plot(ctrl,mark.=FALSE,ylab="S(t)",xli=c(0,12),xlab="time",lt=1:8,main="Control")
legend(bty="n",legend=paste("Trial",1:8),lty=1:8,"topright")
plot(trt,mark.=FALSE,ylab="S(t)",xlab="time",lty=1:8,main="Treated")
legend(bty="n",legend=paste("Trial",1:8),lty=1:8,"topright")
Run the code above in your browser using DataLab