Comment obtenir l'estimateur de Kaplan-Meier pour une analyse de survie avec R ?

L’estimateur de Kaplan-Meier donne la fonction de survie non paramétrique.
Pour l’obtenir sous R on peut utiliser le package survival.

On se place ici dans un cas très simple où il n’y a ni censure ni troncature.
Pour bien comprendre le code, je vous conseille vivement de regarder la documentation du package en question!!



#survival analysis

ls()
rm(list=ls())

library(survival)

#on crée un jeu de données correspondant à des durées (étudiées dans l'analyse de survie)

z<-c(14,14,14,17,17,17,17,17,17,17,17,17,17,17,17,17,17,17,17,17,17,17,20,20,20,20,20,20,20,20,20,20,20,20,20,20,20,20,20,20,20,20,20,20,20,20,20,20,23)
d<-data.frame(delay=z)

#on crée une colonne status, ici tous les individus sont "morts" pendant l'expérience
# mort au sens de l'analyse de survie

d$status<-1

s<-survfit(Surv(d$delay,d$status)~1)
plot(s,main="survival function")