Porovnání Weibullova a log-normálního modelu I
V tomto cvičení přidáme správné pravděpodobnosti přežití do datového rámce, který pak použijeme k vykreslení křivek přežití. surv_wide je datový rámec ve širokém formátu obsahující informace o hormonální terapii a křivky přežití pro Weibullův a log-normální model.
Balíčky survival, survminer a reshape2 a data GBSG2 jsou v tomto cvičení již načtené.
Toto cvičení je součástí kurzu
Survival Analysis in R
Pokyny k cvičení
- Převeď datový rámec do dlouhého formátu pomocí funkce
melt(). - Pomocí
surv_long$surv_idpřidej do datového rámce správné pravděpodobnosti přežitísurv.
Interaktivní cvičení na vyzkoušení si v praxi
Vyzkoušejte si toto cvičení dokončením tohoto ukázkového kódu.
# Melt the data.frame into long format.
surv_long <- melt(surv_wide, ___ = c("horTh", "dist"), ___ = "surv_id", ___ = "time")
# Add column for the survival probabilities
surv_long$surv <- surv[as.numeric(surv_long$___)]
# Add columns upper, lower, std.err, and strata contianing NA values
surv_long[, c("upper", "lower", "std.err", "strata")] <- NA