# Pon los datos en el formato data.frame
diu <- data.frame(tiempo=c(10, 13, 18, 19, 23, 30, 36, 38, 54,
                  56, 59, 75, 93, 97, 104, 107, 107, 107),
                  estatus=c(1, 0, 0, 1, 0, 1, 1, 0, 0,
                            0, 1, 1, 1, 1, 0, 1, 0, 0))

library(survival)

# Encuentra el estimador Kaplan-Meier
diu.surv <- survfit(Surv(tiempo, estatus), data=diu)

# Muestra el estimador y otras cantidades importantes
summary(diu.surv)

# Obten una grafica del estimador
plot(diu.surv, conf.int=FALSE, mark.time=FALSE)
title(main="Estimador Kaplan-Meier\npara el estudio del DIU",
     xlab="Tiempo de discontinuidad",
     ylab="Funcion de supervivencia estimada")