Оцінювання кривої виживання за допомогою survreg()

Тепер ми можемо оцінити виживання пацієнток із раком молочної залози в даних GBSG2, використовуючи модель Вейбулла.

Розподіл Вейбулла має два параметри, які визначають форму кривої виживання.

Пакет survival і дані GBSG2 уже завантажено для вас у цій вправі.

Ця вправа є частиною курсу

Аналіз виживаності в R

Переглянути курс

Інструкції до вправи

Оцініть модель Вейбулла для пацієнток із раком молочної залози.
Обчисліть оцінену криву виживання з моделі за допомогою функції predict() з параметром type = "quantile".
Створіть data.frame з точками часу та відповідними ймовірностями виживання.
Перегляньте перші кілька рядків результату за допомогою функції head().

Інтерактивна практична вправа

Спробуйте виконати цю вправу, доповнивши цей зразок коду.

# Weibull model
wb <- ___(___(___, ___) ~ ___, data = ___)

# Retrieve survival curve from model probabilities 
surv <- seq(.99, .01, by = -.01)

# Get time for each probability
t <- predict(___, type = ___, p = ___, newdata = data.frame(1))

# Create data frame with the information
surv_wb <- data.frame(time = ___, surv = ___)

# Look at first few lines of the result
head(___)

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