Porównanie modelu Weibulla i estymatora Kaplana-Meiera
Narysujmy krzywą przeżycia uzyskaną z modelu Weibulla dla danych GBSG2!
Pakiety survival i survminer oraz dane GBSG2 są już wczytane w tym ćwiczeniu.
To ćwiczenie jest częścią kursu
Analiza przeżycia w R
Instrukcje do ćwiczenia
- Oblicz model Weibulla dla danych
GBSG2. - Wyznacz krzywą przeżycia na podstawie modelu.
- Narysuj krzywe przeżycia uzyskane z obu estymatorów.
Interaktywne ćwiczenie praktyczne
Spróbuj tego ćwiczenia, uzupełniając ten przykładowy kod.
# 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 = ___)