Forest plot for multiple cox regression (Hazard Ratio)

最新推荐文章于 2022-11-27 10:35:52 发布

importAlan

最新推荐文章于 2022-11-27 10:35:52 发布

阅读量719

点赞数

分类专栏： Forest plot 文章标签： r语言

本文链接：https://blog.csdn.net/importAlan/article/details/119146922

版权

本文介绍如何利用R语言进行多重Cox回归分析，并通过Forest plot展示各因素的风险比（Hazard Ratio），帮助理解不同变量对生存时间的影响。

摘要由CSDN通过智能技术生成

library(forestplot)
library(survival)
library(survminer)
setwd("/Users/elijah/Desktop/ICGC/21.multiIndep")
rt=read.table("/Users/elijah/Desktop/ICGC/21.multiIndep/indepInput.txt",header=T,sep="\t",check.names=F,row.names=1)

multiCox=coxph(Surv(futime, fustat) ~ ., data = rt)
summary(multiCox)
multiCoxSum=summary(multiCox)

outTab=data.frame()
outTab=cbind(
  HR=multiCoxSum$conf.int[,"exp(coef)"],
  HR.95L=multiCoxSum$conf.int[,"lower .95"],
  HR.95H=multiCoxSum$conf.int[,"upper .95"],
  pvalue=multiCoxSum$coefficients[,"Pr(>|z|)"],
  N= multiCoxSum$n)
outTab=cbind(id=row.names(outTab),outTab)
outTab=as.data.frame(outTab)
a = 2:ncol(outTab)
for (i in a ){
  outTab[,i]= as.numeric(outTab[,i])
  outTab