数据库代码公开：单基因单肿瘤uniCox生存分析

#以MKI67为例
geneName = "MKI67"
#以LUAD为例
CancerName = "LUAD"
#读取基因表达量
data = read.table(paste0("expTime.",geneName,".txt"), 
                  header=T, 
                  sep="\t", 
                  check.names=F,
                  row.names = 1)
#生存单位改成年
data$OS.time = data$OS.time/365     
#提取目标肿瘤的数据
data = data[data$CancerType == CancerName,]
#单因素COX分析
outTab=data.frame()
cox = coxph(Surv(OS.time, OS) ~ data[,geneName], 
          data = data)
coxSummary = summary(cox)
outTab = rbind(outTab,
             cbind(id=paste0("TCGA-",CancerName,"-OS"),
                   HR=coxSummary$conf.int[,"exp(coef)"],
                   HR.95L=coxSummary$conf.int[,"lower .95"],
                   HR.95H=coxSummary$conf.int[,"upper .95"],
                   pvalue=coxSummary$coefficients[,"Pr(>|z|)"])