Index Date Imputation for Survival Analysis in Externally Controlled Trials with Delayed Treatment Initiation

Résumé technique : Imputation de la date d'index pour l'analyse de survie dans les essais à groupe témoin externe avec initiation de traitement différée

Énoncé du problème
Les essais à groupe témoin externe (ECTs) comparent les résultats d'un essai à bras unique à des groupes témoins externes dérivés d'essais historiques, de registres ou d'études observationnelles. Un défi méthodologique critique surgit lorsque le suivi de l'essai à bras unique est indexé sur l'initiation du traitement, alors que les données du groupe témoin externe sont indexées sur un jalon clinique antérieur (par exemple, le diagnostic, la rechute ou l'achèvement d'une thérapie initiale).

Ce désalignement crée deux problèmes principaux :

1. Biais de temps immortel (Immortal Time Bias) : Les patients de l'essai à bras unique doivent survivre du jalon antérieur jusqu'au point d'initiation du traitement pour être inclus dans l'analyse. Les contrôles externes, indexés sur le jalon antérieur, ne partagent pas cette exigence de survie, ce qui fausse les ensembles de risque et gonfle les estimations de survie pour le bras de traitement.
2. Distribution tronquée du temps d'initiation : Dans l'essai à bras unique, les temps d'initiation du traitement ne sont observés que pour les patients qui survivent suffisamment longtemps pour initier le traitement. Par conséquent, la distribution empirique des temps d'initiation observés est une distribution sélectionnée par la survie (tronquée), et non la distribution marginale des temps d'initiation dans la population cible. Les méthodes de rééchantillonnage standard (par exemple, la correspondance de la distribution du temps de prescription ou Prescription-Time Distribution Matching) qui utilisent la distribution empirique des temps observés ne parviennent pas à tenir compte de cette troncature, ce qui conduit à une imputation biaisée des dates d'index pour les contrôles externes.

Méthodologie : Imputation de la date d'index (IDI)
Les auteurs proposent l'Imputation de la date d'index (IDI), une approche tenant compte de la troncature, conçue pour reconstruire des dates d'index comparables pour les contrôles externes dans les contextes d'initiation de traitement différée. La procédure cible l'Effet Moyen du Traitement sur les Traités (ATT) au sein de la population post-initiation (patients éligibles pour l'intervention après une phase de pré-initiation).

La procédure IDI se compose de quatre étapes :

1. Estimation de la distribution du temps d'initiation et des scores de propension :

   • La méthode estime la distribution marginale des temps d'initiation du traitement ($F_R$) dans la population cible du bras unique. Crucialement, elle corrige le biais de troncature inhérent aux données observées du bras unique (où l'initiation n'est observée que si $T > R$).
   • Sous une hypothèse de travail selon laquelle le temps d'initiation $R$ est indépendant des covariables de base $X$, $F_R$ est estimée à l'aide d'un estimateur pondéré par l'inverse de la survie non paramétrique.
   • Simultanément, un score de propension pour l'appartenance à la cohorte ($p(X) = P(G=1|X)$) est estimé à l'aide d'une régression logistique. Pour garantir que le score de propension reflète la population cible au moment de l'entrée dans le parcours (plutôt que la cohorte observée sélectionnée par la survie), les patients du bras unique sont pondérés par l'inverse de leur probabilité de satisfaire à la condition d'éligibilité (survie jusqu'à $R$) au moment de l'initiation observée.

2. Imputation des temps de pseudo-initiation :

   • Pour chaque patient du groupe témoin externe, un temps de pseudo-initiation ($\hat{R}_j$) est tiré de la distribution marginale estimée $\hat{F}_R$.
   • Le temps de suivi est ré-indexé à partir du jalon initial jusqu'à $\hat{R}_j$.
   • Les patients du groupe témoin externe ne sont conservés dans l'analyse que s'ils satisfont à l'ensemble de risque post-initiation et à la condition d'éligibilité à $\hat{R}_j$ (par exemple, ils doivent avoir survécu et être restés sans progression jusqu'au temps imputé).

3. Estimation des contrastes de survie marginaux :

   • La cohorte ajustée du groupe témoin externe est analysée à l'aide des poids ATT dérivés de l'étape 1 pour aligner la distribution des covariables avec la population cible du bras unique.
   • Les fonctions de survie sont estimées à l'aide d'estimateurs de Kaplan-Meier pondérés ou de modèles de Cox pondérés sur les temps de suivi ré-indexés.

4. Inférence :

   • Un bootstrap non paramétrique est utilisé pour tenir compte de la variabilité introduite par l'imputation stochastique des dates d'index et l'estimation des paramètres de nuisance (scores de propension, poids de troncature).

Principales contributions

• Estimation tenant compte de la troncature : Contrairement aux méthodes existantes d'émulation de date d'index (par exemple, PTDM) qui échantillonnent à partir de la distribution empirique des temps d'initiation observés, l'IDI estime explicitement la distribution marginale du temps d'initiation en corrigeant le fait que les temps observés sont tronqués par la survie et l'éligibilité.
• Cadre intégré : L'IDI couple l'alignement temporel (imputation des dates d'index) avec l'ajustement des covariables (pondération par score de propension ou appariement) pour traiter à la fois le désalignement de l'origine temporelle et le biais de confusion au niveau de la population.
• Définition formelle de l'estimande post-initiation : L'article définit formellement l'estimande cible comme la fonction de survie indexée au temps de décision $R$ pour les patients qui restent éligibles à $R$, se distinguant ainsi des stratégies assignées à la ligne de base.

Résultats
Les auteurs ont évalué l'IDI via des études de Monte Carlo et une analyse d'émulation à partir des données de l'essai oncologique randomisé ECOG-ACRIN 5508.

• Études de simulation :

  • Dans des scénarios de désalignement de la date d'index et de troncature basée sur la survie, les analyses naïves et la pondération par score de propension sans alignement de la date d'index ont montré un biais substantiel et une faible couverture des intervalles de confiance.
  • L'IDI combinée à la pondération par score de propension a considérablement réduit le biais et maintenu des taux de couverture nominaux.
  • Sensibilité à la troncature : Une comparaison ciblée avec une imputation de type PTDM a montré que, bien que les deux méthodes se soient comportées de manière similaire sous une légère troncature, la PTDM est devenue de plus en plus biaisée à mesure que la sévérité de la troncature augmentait. La PTDM attribuait systématiquement des dates de pseudo-index plus précoces aux contrôles, prolongeant artificiellement leur survie. L'IDI est restée approximativement non biaisée à travers tous les niveaux de troncature en estimant correctement la distribution marginale de l'initiation.
  • La méthode a démontré sa robustesse même lorsque l'hypothèse de travail d'indépendance entre le temps d'initiation et les covariables était violée, bien que la variabilité ait augmenté.

• Essai émulé (ECOG-ACRIN 5508) :

  • Les auteurs ont émulé un essai à groupe témoin externe comparant la maintenance par pemetrexed (bras unique) versus la maintenance par bevacizumab (groupe témoin externe), en utilisant l'essai randomisé comme référence.
  • L'analyse externe naïve (indexée à l'enregistrement de l'essai) favorisait fortement le bras de traitement (Hazard Ratio [HR] = 0,447), s'écartant significativement de la référence randomisée (HR = 0,860).
  • Après application de l'IDI avec pondération par score de propension, l'HR estimé s'est rapproché de la référence (HR = 0,828), et les différences de probabilité de survie à un jalon donné se sont également mieux alignées avec les résultats randomisés. Cela a démontré la capacité de l'IDI à réduire la divergence causée par le désalignement de la date d'index et les différences d'éligibilité.

Signification et affirmations
L'article affirme que l'IDI fournit une stratégie pratique pour aligner les dates d'index dans les analyses de survie impliquant une initiation de traitement différée, un scénario courant en oncologie et dans les maladies rares où les essais randomisés sont irréalisables.

• Nécessaire mais non suffisant : Les auteurs soulignent que l'IDI traite les composantes d'alignement temporel et de troncature de l'analyse de groupe témoin externe, qui sont nécessaires pour une interprétation causale, mais ne sont pas suffisantes. Une utilisation valide nécessite toujours un chevauchement adéquat des covariables, une ascertaination comparable des résultats et la capacité d'évaluer l'éligibilité au moment de la date d'index imputée.
• Distinction par rapport aux méthodes existantes : Ce travail se distingue des méthodes ciblant des stratégies assignées à la ligne de base (par exemple, le clonage-censure-pondération) et de l'émulation standard de la date d'index qui ignore la troncature.
• Limites : Les auteurs reconnaissent des limites, notamment le recours à une hypothèse de travail d'indépendance pour la distribution du temps d'initiation (bien qu'une robustesse ait été observée), le potentiel de réduction de la taille de l'échantillon effectif dû au filtrage de l'éligibilité au moment de l'imputation, et l'exigence que des données d'éligibilité longitudinales (par exemple, le statut de progression) soient disponibles dans la cohorte du groupe témoin externe pour émuler pleinement l'ensemble de risque post-initiation.

En résumé, l'article présente l'IDI comme une correction méthodologique pour un biais spécifique survenant dans les essais à groupe témoin externe, démontrant par simulation et émulation qu'elle peut récupérer des estimations non biaisées des effets de traitement là où les approches standard échouent en raison de l'origine temporelle mal alignée et des temps d'initiation sélectionnés par la survie.