Оценка медианной выживаемости с помощью модели Вейбулла
Теперь мы можем оценить выживаемость пациенток с раком молочной железы из набора данных GBSG2 с помощью модели Вейбулла (функция survreg()).
Напомним, что модель Вейбулла оценивает гладкую функцию выживаемости, а не ступенчатую, как метод Каплана–Мейера.
Функция predict() с параметром type = "quantile" позволяет вычислить квантили функции распределения. Воспользуемся ею для вычисления медианной выживаемости.
Пакет survival и данные GBSG2 уже загружены в этом упражнении.
Это упражнение является частью курса
Анализ выживаемости в R
Инструкции к упражнению
- Постройте модель Вейбулла для пациенток с раком молочной железы.
- Вычислите медианную выживаемость по этой модели с помощью функции
predict()с параметромtype = "quantile".
Интерактивное практическое упражнение
Попробуйте выполнить это упражнение, дополнив этот пример кода.
# Weibull model
wb <- ___(Surv(time, cens) ~ 1, data = GBSG2)
# Compute the median survival from the model
predict(wb, type = "quantile", p = ___, newdata = data.frame(1))