Odhadování křivky přežití pomocí survreg()
Teď můžeme odhadnout přežití pacientek s rakovinou prsu v datech GBSG2 pomocí Weibullova modelu.
Weibullovo rozdělení má dva parametry, které určují tvar křivky přežití.
Balíček survival a data GBSG2 jsou v tomto cvičení již načteny.
Toto cvičení je součástí kurzu
Survival Analysis in R
Pokyny k cvičení
- Odhadni Weibullův model pro pacientky s rakovinou prsu.
- Z modelu vypočítej odhadovanou křivku přežití pomocí funkce
predict()s argumentemtype = "quantile". - Vytvoř
data.frames časovými body a odpovídajícími pravděpodobnostmi přežití. - Prohlédni si prvních několik řádků výsledku pomocí funkce
head().
Interaktivní cvičení na vyzkoušení si v praxi
Vyzkoušejte si toto cvičení dokončením tohoto ukázkového kódu.
# 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(___)