Число в группе риска для регрессионного графика Кокса

Question

могу ли я составить таблицу «число в группе риска» для графика Кокса, если у меня более одной независимой переменной? если это возможно, где я могу найти соответствующий код (я искал, но не смог найти) код, который я использовал для своих данных:

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

таблица «число под угрозой» для графика Кокса

Это не график Кокса, это график Каплана-Мейера. Вы пытаетесь построить модель Кокса, когда вам нужно подогнать кривые КМ, используя survfit, а затем построить результирующее соответствие:

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

Так как теперь вы упоминаете, что у вас есть постоянные предикторы, возможно, вы могли бы подумать о том, что вы ожидаете показать на таблице риска или на графике КМ. Вот пример создания непрерывной меры (возраста):

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)