如果我有多个独立变量,是否可以为Cox图制作“处于危险中的数字”表?如果可以在哪里找到相关代码(我搜索了但找不到)我在数据上使用的代码:
fit <- coxph(Surv(time,event) ~chr1q21_status+CCND1+CRTM1+IRF4,data = myeloma)
ggsurvplot(fit, data = myeloma,
risk.table=TRUE, break.time.by=365, xlim = c(0,4000),
risk.table.y.text=FALSE, legend.labs = c("2","3","4+"))
虽然仅对一个变量及其功能起作用,但得到了此消息对象'ggsurv'找不到'
“ cox图的风险编号”表
这不是考克斯图,这是卡普兰-迈耶图。您想要绘制一个Cox模型时,想要使用survfit
拟合KM曲线,然后绘制结果拟合:
library("survival")
library("survminer")
fit <- survfit(Surv(time,status) ~ ph.ecog + sex , data = lung)
ggsurvplot(fit, data = lung, risk.table = TRUE)