Kwantyle krzywej przeżycia z modelu Weibulla
Teraz możemy oszacować przeżycie pacjentek z rakiem piersi na podstawie danych GBSG2, korzystając z modelu Weibulla.
Funkcja predict() z argumentem type = "quantile" pozwala obliczyć kwantyle funkcji rozkładu. Jak już wiesz z tego kursu, funkcja przeżycia to 1 minus funkcja rozkładu
(\(S = 1 - D\)), dlatego za pomocą funkcji predict() możemy łatwo wyznaczyć kwantyle funkcji przeżycia.
Pakiet survival oraz dane 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.
- Wyznacz punkt czasowy, w którym prawdopodobieństwo przeżycia dłużej niż ten punkt wynosi 70% (korzystając z funkcji
predict()z argumentemtype = "quantile").
Interaktywne ćwiczenie praktyczne
Spróbuj tego ćwiczenia, uzupełniając ten przykładowy kod.
# Weibull model
wb <- ___(Surv(___, ___) ~ 1, data = GBSG2)
# 70 Percent of patients survive beyond time point...
predict(wb, type = ___, p = ___, newdata = data.frame(1))