require(survival) # print(KM <- survfit(Surv(futime, fustat) ~ rx, data=ovarian)) layout(t(1:2)) plot(KM, main="治療種類別カプラン=マイヤプロット", lty=1:2) legend("bottomright", lty=1:2, legend=names(table(ovarian$rx))) plot(KM, main="治療種類別二重対数プロット", fun=function(y) {log(-log(y))}, lty=1:2) legend("bottomright", lty=1:2, legend=names(table(ovarian$rx))) # res <- coxph(Surv(futime,fustat)~rx+age+factor(resid.ds),data=ovarian) summary(res)