Skatta överlevnadskurvan med survreg()
Nu kan vi skatta överlevnaden för bröstcancerpatienterna i datamängden GBSG2 med hjälp av en Weibull-modell.
Weibull-fördelningen har två parametrar som bestämmer överlevnadskurvans form.
Paketet survival och datamängden GBSG2 är redan inlästa i den här övningen.
Den här övningen är en del av kursen
Överlevnadsanalys i R
Övningsinstruktioner
- Skatta en Weibull-modell för bröstcancerpatienterna.
- Beräkna den skattade överlevnadskurvan från modellen med funktionen
predict()och argumentettype = "quantile". - Skapa en
data.framemed tidpunkterna och motsvarande överlevnadssannolikheter. - Titta på de första raderna i resultatet med funktionen
head().
Interaktiv övning med praktiskt arbete
Testa den här övningen genom att slutföra den här exempelkoden.
# 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(___)