Ejercicio ignorando la censura
Organizas una fiesta y, a la 1 a. m., de repente tus invitados empiezan a bailar. Te pica la curiosidad por analizar cuánto tiempo bailarán y empiezas a recoger datos. El problema es que te cansas y te vas a dormir al cabo de un rato.
Obtienes los siguientes tiempos de baile con censura por la derecha en dancedat:
namees el nombre de tu amigo o amiga.timees el tiempo de baile con censura por la derecha.obs_endindica si observaste el final del baile de tus amigos (1) o si te fuiste a dormir antes de que dejaran de bailar (0).
Empiezas a analizar los datos por la mañana, pero estás cansado y, al principio, ignoras que tienes observaciones censuradas. Luego recuerdas este curso en DataCamp y lo haces correctamente.
El paquete survival ya está cargado para este ejercicio.
Este ejercicio forma parte del curso
Análisis de supervivencia en R
Instrucciones del ejercicio
- Estima la función de supervivencia fingiendo que todas las observaciones censuradas son observaciones reales.
- Estima la función de supervivencia de este conjunto de datos mediante Kaplan-Meier.
- Dibuja las curvas de supervivencia correcta y errónea juntas usando
ggsurvplot_combine()y compáralas. Fíjate en cómo difieren.
Ejercicio interactivo práctico
Prueba este ejercicio y completa el código de muestra.
# Create dancedat data
dancedat <- data.frame(
name = c("Chris", "Martin", "Conny", "Desi", "Reni", "Phil",
"Flo", "Andrea", "Isaac", "Dayra", "Caspar"),
time = c(20, 2, 14, 22, 3, 7, 4, 15, 25, 17, 12),
obs_end = c(1, 1, 0, 1, 1, 1, 1, 1, 0, 0, 0))
# Estimate the survivor function pretending that all censored observations are actual observations.
km_wrong <- survfit(___(time) ~ 1, data = dancedat)
# Estimate the survivor function from this dataset via kaplan-meier.
km <- survfit(___(___, ___) ~ ___, data = dancedat)
# Plot the two and compare
ggsurvplot_combine(list(correct = ___, wrong = ___))