De overlevingscurve schatten met survreg()

We kunnen nu de overleving van de borstkankerpatiënten in de GBSG2-gegevens schatten met een Weibull-model.

De Weibull-verdeling heeft twee parameters die de vorm van de overlevingscurve bepalen.

Het pakket survival en de gegevens GBSG2 zijn voor je ingeladen in deze oefening.

Deze oefening maakt deel uit van de cursus

Survivalanalyse in R

Cursus bekijken

Oefeninstructies

Schat een Weibull-model voor de borstkankerpatiënten.
Bereken de geschatte overlevingscurve uit het model met de functie predict() en type = "quantile".
Maak een data.frame met de tijdspunten en bijbehorende overlevingskansen.
Bekijk de eerste rijen van het resultaat met de functie head().

Praktische interactieve oefening

Probeer deze oefening eens door deze voorbeeldcode in te vullen.

# 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(___)

Code bewerken en uitvoeren