Výpočet Weibullova modelu
U Weibullova modelu s kovariátami můžeme křivku přežití vypočítat stejně jako u modelu bez kovariát. Jediné, co potřebujeme udělat, je zadat hodnoty kovariát pro danou křivku přežití ve funkci predict(). K tomu slouží argument newdata.
Toto cvičení je součástí kurzu
Survival Analysis in R
Pokyny k cvičení
- Vypočítej Weibullův model pro dataset
GBSG2s kovariátouhorTh, abys zjistil/a vliv hormonální terapie na přežití pacientů. - Vypočítej křivku přežití pro pacienty, kteří hormonální terapii dostávají.
- Prohlédni si křivku přežití pomocí funkce
str().
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
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)