cox回归图的风险数

Question

如果我有一个以上的自变量，我可以为考克斯图制作"number at risk“表吗？如果可能，我在哪里可以找到我在数据上使用的相关代码(我搜索了，但找不到)：

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‘未找到’，尽管只有一个变量和函数survfit有效。

Answer 1

考克斯图的

“风险数”表

这不是考克斯图，而是Kaplan-Meier图。您正在尝试绘制Cox模型，而您想要的是使用survfit拟合KM曲线，然后绘制结果拟合：

library("survival")
library("survminer")
fit <- survfit(Surv(time,status) ~ ph.ecog + sex , data = lung)
ggsurvplot(fit, data = lung, risk.table = TRUE)

​

​

既然你现在提到你有连续的预测器，也许你可以考虑你期望的风险表或KM图显示什么。下面是一个将连续度量(Age)打包的示例：

library("survival")
library("survminer")
#> Loading required package: ggplot2
#> Loading required package: ggpubr
#> Loading required package: magrittr
lung$age_bin <- cut(lung$age, quantile(lung$age))
fit <- survfit(Surv(time,status) ~ age_bin + sex , data = lung)
ggsurvplot(fit, data = lung, risk.table = TRUE)

​

​