Weibull- en lognormaalmodel vergelijken I
In deze oefening willen we de juiste overlevingskansen toevoegen aan een data frame. Dit data frame wordt gebruikt om de overlevingscurves te plotten. surv_wide is een breed data frame met informatie over hormonale therapie en de overlevingscurves voor de Weibull- en lognormaalmodellen.
De pakketten survival, survminer en reshape2 en de GBSG2-gegevens zijn al voor je geladen in deze oefening.
Deze oefening maakt deel uit van de cursus
Survivalanalyse in R
Oefeninstructies
- Zet het data frame om naar lang formaat met melt (melt de data).
- Gebruik
surv_long$surv_idom de juiste overlevingskansensurvtoe te voegen aan het data frame.
Praktische interactieve oefening
Probeer deze oefening eens door deze voorbeeldcode in te vullen.
# 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