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")