Überlebenskurve mit survreg() schätzen
Wir können nun das Überleben der Brustkrebspatientinnen in den GBSG2-Daten mit einem Weibull-Modell schätzen.
Die Weibull-Verteilung hat zwei Parameter, die die Form der Überlebenskurve bestimmen.
Das Paket survival und die Daten GBSG2 sind für diese Übung bereits geladen.
Diese Übung ist Teil des Kurses
Survival-Analyse in R
Anleitung zur Übung
- Schätze ein Weibull-Modell für die Brustkrebspatientinnen.
- Berechne die geschätzte Überlebenskurve aus dem Modell mit der Funktion
predict()undtype = "quantile". - Erzeuge ein
data.framemit den Zeitpunkten und den zugehörigen Überlebenswahrscheinlichkeiten. - Schau dir die ersten Zeilen des Ergebnisses mit der Funktion
head()an.
Interaktive Übung
Vervollständige den Beispielcode, um diese Übung erfolgreich abzuschließen.
# Weibull model
wb <- ___(___(___, ___) ~ ___, data = ___)
# Retrieve survival curve from model probabilities
surv <- seq(.99, .01, by = -.01)
# Get time for each probability
t <- predict(___, type = ___, p = ___, newdata = data.frame(1))
# Create data frame with the information
surv_wb <- data.frame(time = ___, surv = ___)
# Look at first few lines of the result
head(___)