Построение модели Вейбулла

Для модели Вейбулла с ковариатами кривую выживаемости можно вычислить так же, как и для модели без ковариат. Единственное, что нужно сделать, — указать значения ковариат для заданной кривой выживаемости в функции predict(). Для этого используется аргумент newdata.

Это упражнение является частью курса

Анализ выживаемости в R

Посмотреть курс

Инструкции к упражнению

Постройте модель Вейбулла для набора данных GBSG2 с ковариатой horTh, чтобы проанализировать влияние гормональной терапии на выживаемость пациентов.
Вычислите кривую выживаемости для пациентов, получающих гормональную терапию.
Изучите кривую выживаемости с помощью функции str().

Интерактивное практическое упражнение

Попробуйте выполнить это упражнение, дополнив этот пример кода.

# Weibull model
wbmod <- survreg(Surv(___, ___) ~ ___, data = GBSG2)
coef(wbmod)

# Retrieve survival curve from model
surv <- seq(.99, .01, by = -.01)
t_yes <- predict(___, type = "quantile", p = ___,
  newdata = data.frame(horTh = ___))

# Take a look at survival curve
str(t_yes)

Редактировать и запускать код