Weibull-model berekenen

Voor een Weibull-model met covariaten kunnen we de overlevingscurve net zo berekenen als voor het Weibull-model zonder covariaten. Het enige wat we moeten doen, is de waarden van de covariaten voor een specifieke overlevingscurve opgeven in de functie predict(). Dit kan met het argument newdata.

Deze oefening maakt deel uit van de cursus

Survivalanalyse in R

Cursus bekijken

Oefeninstructies

Bereken een Weibull-model voor de GBSG2-gegevensset met covariaat horTh om het effect van hormonale therapie op de overleving van patiënten te analyseren.
Bereken de overlevingscurve voor patiënten die hormonale therapie krijgen.
Bekijk de overlevingscurve met str().

Praktische interactieve oefening

Probeer deze oefening eens door deze voorbeeldcode in te vullen.

# Weibull model
wbmod <- survreg(Surv(___, ___) ~ ___, data = GBSG2)
coef(wbmod)

# Retrieve survival curve from model
surv <- seq(.99, .01, by = -.01)
t_yes <- predict(___, type = "quantile", p = ___,
  newdata = data.frame(horTh = ___))

# Take a look at survival curve
str(t_yes)

Code bewerken en uitvoeren