생존분석(survival analysis, kaplan-meier plot)

library(survival)

os=survfit(Surv('생존기간','censored여부(0/1)')~'그룹',data = 'data가있는 데이터프레임')

그림그리기

plot(os,main='제목',xlab='Days', ylab='Survival probability',lwd=2,col=c(1,2),

mark.time=T)

p.value얻기

pv=survdiff(Surv('생존기간','censored여부(0/1)')~'그룹',data = 'data가있는 데이터프레임')

pchisq(pv$chisq,length(pv$n)-1,lower.tail=F)