Оценка кривой выживаемости с помощью 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(___)

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