Estimateur de Kaplan-Meier

Soit S(t) la probabilité qu’un membre d’une population donnée ait une durée de vie supérieure à t. Pour un échantillon de taille N dans cette population, les durées observées jusqu’à chaque décès des membres de l’échantillon N sont :

${\displaystyle t_{1}\leq t_{2}\leq t_{3}\leq \cdots \leq t_{N}.}$

À chaque n_i correspond un t_i, n_i étant le nombre de personnes « à risque » juste avant le temps t_i, et d_i le nombre de décès au temps t_i.

On note que les intervalles entre chaque évènement ne sont pas uniformes. Par exemple, une petite quantité de données peut commencer avec 10 cas. Supposons que le sujet 1 décède au jour 3, les sujets 2 et 3 au jour 11 et le sujet 4 disparaît du suivi (donnée censurée) au jour 9. Les données pour les 2 premiers sujets seraient les suivantes :

${\displaystyle i}$	1	2
${\displaystyle t_{i}}$	3	11
${\displaystyle d_{i}}$	1	2
${\displaystyle n_{i}}$	10	8

L’estimateur de Kaplan-Meier est l’estimation du maximum de vraisemblance non-paramétrique de S(t). C’est un produit de la forme :

${\displaystyle {\hat {S}}(t)=\prod \limits _{t_{i}<t}{\frac {n_{i}-d_{i}}{n_{i}}}.}$

Lorsqu’il n’y aucune censure, n_i est le nombre de survivants juste avant le temps t_i.
Lorsqu’il y a censure, n_i est le nombre de survivants moins le nombre de pertes (cas censurés). Ce sont seulement ces cas survivants qui continuent à être observés (qui n’ont pas encore été censurés) qui sont « à risque » de décès (observé).

Ici une autre définition possible parfois utilisée :

${\displaystyle {\hat {S}}(t)=\prod \limits _{t_{i}\leq t}{\frac {n_{i}-d_{i}}{n_{i}}}.}$

Les deux définitions diffèrent uniquement aux moments des évènements observés. La dernière définition est « continue à droite » tandis que la première est « continue à gauche ». Soit T la variable aléatoire qui mesure le temps d’échec et soit F(t) sa Fonction de répartition cumulative. On note que :

${\displaystyle S(t)=P[T>t]=1-P[T\leq t]=1-F(t).\,}$

En conséquence, la définition continue à gauche de ${\displaystyle \scriptstyle {\hat {S}}(t)}$ peut être préférée pour rendre l’estimation compatible avec une estimation continue à droite de F(t).

L’estimateur de Kaplan-Meier est une statistique, et certains estimateurs sont utilisés pour approcher sa variance. Un de ces estimateurs les plus courants est la formule de Greenwood :

${\displaystyle {\widehat {\mathrm {Var} }}({\widehat {S}}(t))={\widehat {S}}(t)^{2}\sum \limits _{t_{i}<t}{\frac {d_{i}}{{n_{i}}({n_{i}-d_{i}})}}.}$

En 1983, Edward L. Kaplan raconte la genèse de l'estimateur de Kaplan-Meier[3].

Le tout débute en 1952, révèle Kaplan, quand Paul Meier (alors en stage post-doctoral à l'université Johns-Hopkins, au Maryland), après avoir pris connaissance de l'article de Greenwood, publié en 1926, sur la durée du cancer, veut proposer un puissant estimateur de survie appuyé sur les résultats d'essais cliniques. En 1953, le mathématicien Kaplan (travaillant alors aux Laboratoires Bell, au New Jersey) veut proposer un estimateur de la durée qu'auront les tubes à vide utilisés pour amplifier et retransmettre les signaux dans le système de câbles téléphoniques sous-marins. Kaplan soumet son projet d'article au professeur John W. Tukey, qui œuvrait aussi pour les Laboratoires Bell et qui venait d'être maître de thèse de Meier[4] à Princeton, au New Jersey. Chacun des deux jeunes chercheurs avait soumis son manuscrit au Journal of the American Statistical Association, qui leur recommandait d'entrer en contact l'un avec l'autre, pour fusionner les deux articles. Alors, Kaplan et Meier entreprennent, par correspondance (courrier postal), de réconcilier leurs points de vue. Durant les quatre ans que dure cette phase, leur seule crainte est qu'un tiers publie avant eux un article proposant une solution équivalente.

L'article Nonparametric estimation from incomplete observations est finalement publié en 1958 (Journal of the American Statistical Association, vol. 53, p. 457–481)[1].

Plusieurs langages de programmation et logiciels statistiques proposent des implémentations de l'estimateur de Kaplan-Meier. On peut notamment citer :

• SAS, au travers de la procédure proc lifetest[5].
• R, via le package survival[6]^,[7]^,[8].
• Stata, via la commande sts[9]^,[10].
• Python, avec les packages lifelines[11] et scikit-survival[12].

• (en) Cet article est partiellement ou en totalité issu de l’article de Wikipédia en anglais intitulé « Kaplan–Meier estimator » (voir la liste des auteurs).

• Portail des probabilités et de la statistique