Квантили кривой выживаемости на основе модели Вейбулла
Теперь можно оценить выживаемость пациенток с раком молочной железы из набора данных GBSG2 с помощью модели Вейбулла.
Функция predict() с параметром type = "quantile" позволяет вычислить квантили функции распределения. Как мы уже знаем из этого курса, функция выживаемости равна 1 минус функция распределения
(\(S = 1 - D\)), поэтому квантили функции выживаемости легко вычислить с помощью функции predict().
Пакет survival и данные GBSG2 уже загружены в этом упражнении.
Это упражнение является частью курса
Анализ выживаемости в R
Инструкции к упражнению
- Постройте модель Вейбулла для пациенток с раком молочной железы.
- Найдите момент времени, в котором вероятность пережить этот момент составляет 70 процентов (используйте функцию
predict()с параметромtype = "quantile").
Интерактивное практическое упражнение
Попробуйте выполнить это упражнение, дополнив этот пример кода.
# Weibull model
wb <- ___(Surv(___, ___) ~ 1, data = GBSG2)
# 70 Percent of patients survive beyond time point...
predict(wb, type = ___, p = ___, newdata = data.frame(1))