Порівняння моделі Вейбулла та оцінки Каплана—Майєра
Давайте побудуємо криву виживання, яку отримуємо з моделі Вейбулла для даних GBSG2!
Пакети survival і survminer та дані GBSG2 вже завантажено для вас у цій вправі.
Ця вправа є частиною курсу
Аналіз виживаності в R
Інструкції до вправи
- Обчисліть модель Вейбулла для даних
GBSG2. - Обчисліть криву виживання з моделі.
- Побудуйте криві виживання, які ви отримали з двох оцінок.
Інтерактивна практична вправа
Спробуйте виконати цю вправу, доповнивши цей зразок коду.
# Weibull model
wb <- ___(___(___, ___) ~ ___, ___)
# Retrieve survival curve from model
surv <- seq(.99, .01, by = -.01)
# Get time for each probability
t <- predict(___, type = ___, p = ___, newdata = data.frame(___))
# Create data frame with the information needed for ggsurvplot_df
surv_wb <- data.frame(time = ___, surv = ___,
upper = NA, lower = NA, std.err = NA)
# Plot
ggsurvplot_df(fit = ___, surv.geom = ___)