TCGA_生存分析

#R 生存曲线
setwd("D:/TCGA_ProstateCancer_Survival/Survival_analysis")
library(survival)
library(survminer)
rt <- read.table("survivalInput.txt",header=T,sep="\t")
rt$futime <- rt$futime/365 
a <- (rt$expression < median(rt$expression))
diff <- survdiff(Surv(futime,fustat)~a,data=rt) 
pValue <- 1 - pchisq(diff$chisq,df=1)
pValue <- round(pValue,5)
fit <- survfit(Surv(futime,fustat) ~ a,data=rt) 
summary(fit) 
ggsurvplot(fit,pval=TRUE,linetype="solid",palette=c("red","blue"),legend.title="",legend=c(0.95,0.95),legend.labs=c("AR High-expression","AR low-expression"),conf.int=T,conf.int.style="ribbon",conf.int.alpha=0.1,font.legend=c(6))
ggsave(file="a.pdf",width = 10, height = 10)