Построение кривых выживаемости объекта coxPH с помощью ggsurvplot - каков самый простой способ?

Question

Я хотел бы отобразить кривые выживаемости из ранее подобранной модели coxPH в R.

R-Версия: 3.6.2

Версии пакета:

• выживаемость: 3,1-8
• выживший: 0.4.6

У меня есть большой скрипт двухлетней давности, в котором есть разные фрагменты кода, которые выглядят так:

require('survival')
require('survminer')
res.cox <- coxph(Surv(time, status) ~ age + sex + wt.loss, data =  lung)
ggsurvplot(survfit(res.cox, data = lung))

К сожалению, когда я снова начал использовать код сейчас (сейчас, два года спустя), это возвращает сообщение об ошибке:

Error: object of type 'symbol' is not subsettable

Я думаю, проблема может быть связана с тем, что формула вызова функции coxPH не копируется в объект Surffit. Я обнаружил, что следующее работает и создает хороший сюжет.

require('survival')
require('survminer')
ggsurvplot(survfit(coxph(Surv(time, status) ~ age + sex + wt.loss, data =  lung), data = lung))

Мне интересно, это ошибка или особенность? Поскольку у меня есть много фрагментов кода, затронутых этим, я хотел спросить, есть ли простое решение, чтобы заставить его снова работать без необходимости копировать и вставлять вызов функции coxph.

Спасибо за вашу помощь,

Моенту