Kvantily křivky přežití z Weibullova modelu
Teď můžeme odhadnout přežití pacientek s rakovinou prsu v datech GBSG2 pomocí Weibullova modelu.
Funkce predict() s argumentem type = "quantile" umožňuje vypočítat kvantily distribuční funkce. Jak jsme se v tomto kurzu dozvěděli, funkce přežití je 1 − distribuční funkce (\(S = 1 - D\)), takže kvantily funkce přežití snadno získáme právě pomocí funkce predict().
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.
- Najdi časový bod, ve kterém je pravděpodobnost přežití delšího než tento bod rovna 70 procentům (použij funkci
predict()s argumentemtype = "quantile").
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 <- ___(Surv(___, ___) ~ 1, data = GBSG2)
# 70 Percent of patients survive beyond time point...
predict(wb, type = ___, p = ___, newdata = data.frame(1))