Exercício ignorando a censura
Você dá uma festa e, à 1h da manhã, os convidados de repente começam a dançar. Você fica curioso para analisar por quanto tempo seus convidados vão dançar e começa a coletar dados. O problema é que você se cansa e vai dormir depois de um tempo.
Você obtém os seguintes dados de tempos de dança com censura à direita em dancedat:
nameé o nome do seu amigo.timeé o tempo de dança com censura à direita.obs_endindica se você observou o fim da dança do seu amigo (1) ou se você foi dormir antes de ele parar de dançar (0).
Você começa a analisar os dados de manhã, mas está cansado e, a princípio, ignora o fato de que tem observações censuradas. Depois você se lembra deste curso da DataCamp e faz do jeito certo.
O pacote survival já está carregado para você neste exercício.
Este exercício faz parte do curso
Análise de Sobrevivência em R
Instruções do exercício
- Estime a função de sobrevivência fingindo que todas as observações censuradas são observações reais.
- Estime a função de sobrevivência deste conjunto de dados via Kaplan-Meier.
- Plote as curvas de sobrevivência correta e incorreta juntas usando
ggsurvplot_combine()e compare. Observe como elas diferem.
Exercício interativo prático
Experimente este exercício completando este código de exemplo.
# 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 = ___))