Kwantielen van de overlevingscurve uit een Weibull-model
We kunnen nu de overleving van borstkankerpatiënten in de GBSG2-gegevens schatten met een Weibull-model.
De functie predict() met type = "quantile" laat je de kwantielen van de verdelingsfunctie berekenen. Zoals je tot nu toe in deze cursus hebt geleerd, is de overlevingsfunctie 1 min de verdelingsfunctie
(\(S = 1 - D\)), dus we kunnen de kwantielen van de overlevingsfunctie eenvoudig berekenen met de functie predict().
Het pakket survival en de GBSG2-gegevens zijn voor je geladen in deze oefening.
Deze oefening maakt deel uit van de cursus
Survivalanalyse in R
Oefeninstructies
- Schat een Weibull-model voor de borstkankerpatiënten.
- Bepaal het tijdstip waarop de kans om langer dan dat tijdstip te overleven 70 procent is (gebruik de functie
predict()mettype = "quantile").
Praktische interactieve oefening
Probeer deze oefening eens door deze voorbeeldcode in te vullen.
# Weibull model
wb <- ___(Surv(___, ___) ~ 1, data = GBSG2)
# 70 Percent of patients survive beyond time point...
predict(wb, type = ___, p = ___, newdata = data.frame(1))