Sample 1

fit<-survfit(Surv(PFS_MONTHS,PFS_STATUS)~gene1[,a],gene1)
p <- ggsurvplot(fit, risk.table = T, pval = T, pval.method = TRUE,
                palette = surv_color,
                xlab = "Time (Months)",
                surv.median.line = "hv",
                ylab = "PFS(%)",
                legend.labs = c("WILD", "MUT"))
  tiff(paste("results_AD/gene_", i,"_surv", ".tiff"), compression="lzw",units="in",res=600,pointsize=8,width=6, height=7)