Szacowanie mediany przeżycia na podstawie modelu Weibulla
Teraz możemy oszacować przeżycie pacjentek z rakiem piersi ze zbioru danych GBSG2 przy użyciu modelu Weibulla (funkcja survreg()).
Pamiętaj, że model Weibulla szacuje gładką funkcję przeżycia – w odróżnieniu od funkcji schodkowej, którą wyznacza metoda Kaplana-Meiera.
Funkcja predict() z argumentem type = "quantile" pozwala obliczać kwantyle funkcji rozkładu.
Wykorzystamy ją do wyznaczenia mediany przeżycia.
Biblioteka survival oraz zbiór danych GBSG2 są już wczytane w tym ćwiczeniu.
To ćwiczenie jest częścią kursu
Analiza przeżycia w R
Instrukcje do ćwiczenia
- Dopasuj model Weibulla dla pacjentek z rakiem piersi.
- Oblicz medianę przeżycia z tego modelu, używając funkcji
predict()z argumentemtype = "quantile".
Interaktywne ćwiczenie praktyczne
Spróbuj tego ćwiczenia, uzupełniając ten przykładowy kod.
# 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))