Ошибка "объект 'статус' не найден" при использовании coxph в R
Всякий раз, когда я добираюсь до функции "coxph" в моем коде, она выдает ошибку "Ошибка в Surv(время, статус): объект" статус "не найден". Но я не получаю такой же ошибки в другом месте. Что я могу сделать?
Это мой текущий код:
install.packages("survminer")
library(survminer)
library(survival)
data(lung)
require("survival")
km_fit <- survfit(Surv(time, status) ~ sex, data=lung)
ggsurvplot(km_fit, linetype = "strata", conf.int = TRUE, pval = TRUE, palette = "Dark2")
cox <- coxph(Surv(time, status), age+sex+ph.ecog+ph.karno+pat.karno+meal.cal+wt.loss, data=lung)
summary(cox)
1 ответ
Решение
Ошибка в вашем коде заключена в запятую вместо использования ~. Если вы замените его, как в следующем примере, он работает:
coxph(Surv(time, status) ~ age+sex+ph.ecog+ph.karno+pat.karno+meal.cal+wt.loss, data=lung)