Обчисліть модель Вейбулла

Для моделі Вейбулла з коваріатами ми можемо обчислити криву виживання так само, як і для моделі Вейбулла без коваріат. Єдине, що потрібно, — вказати значення коваріат для певної кривої виживання у функції 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)

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