杀杀
cox 森林图绘制
cox的结果多为数据,绘制成表格较为不便,因此可以用森林图进行美化展示
R语言中森林图的实现代码较为简单:
library(forestplot)
library(survival)
library(survminer)
library(ggplot2)
model2 <- coxph(Surv(os_m, os_event) ~ age+gender+ki67+kps, data = testdata)
ggforest(model2, #coxph得到的Cox回归结果
data = testdata, #数据集
main = 'Hazard ratio of ALLdata', #标题
cpositions = c(0.05, 0.15, 0.35), #前三列数值的距离
fontsize = 1.4, #字体大小
refLabel = 'Reference', #相对变量的标签
noDigits = 3 #HR、95%CI的小数位
)
最终绘制图形

网友评论