R语言之survival-生存分析

library(survival)library(survminer)两分组最佳cutoffcutoff<-surv_cutpoint(survival_OS_data, #数据集                      time="OS.time",#“ ”里写数据集时间变量                      event="OS",##“ ”里数据集结局变量名称                      variables=c("TP53"));summary(cutoff) #输出结果

好了，确定了最佳cutoff。其实我分了三组（利用x-tile软件分组，确定两组最佳cutoff后就可以根据cutoff值进一步利用x-tile软件划分分三个组），这里只显示高低两组。
下面作图

fit <- survfit(Surv(OS.time, OS)~group,                 data=survival_OS_data_hl)
TCGS_OS <- ggsurvplot(fit,