rmanova<-read.csv("경로/반복측정분산분석.csv")
head(rmanova)

install.packages("survival")
library(survival)

#치료방법별 데이터 생성
r1 = survfit(Surv(time, status == 1) ~ treatment, survival)
r1

summary(r1)

#기간을 한정하여 결과 출력
summary(r1,time=c(0,19))

#치료방법별 생존율
survdiff(Surv(time,status)~treatment, data=survival)

#생존곡선 시각화
install.packages("survminer")
library(survminer)

ggsurvplot(r1,risk.table = T,pval=T)