Littérature scientifique sur le sujet « Immortal time bia »

Negli studi clinici spesso è di interesse confrontare la sopravvivenza di pazienti appartenenti a due o più gruppi di trattamento. In alcune situazioni, la classificazione non è effettuata all’inizio del follow-up ma cambia nel tempo. Ad esempio, tutti i pazienti sono sottoposti ad un trattamento iniziale ed alcuni lo continuano mentre altri cambiano dopo un certo periodo di tempo. In questo caso il trattamento è rappresentato da una variabile binaria tempo-dipendente. Un contesto tipico è il confronto tra chemioterapia e trapianto di cellule staminali nella Leucemia Linfoblastica Acuta. In questa situazione, il metodo Kaplan-Meier non è utilizzabile in quanto affetto da immortal time bias. Due approcci non-parametrici alternativi sono stati proposti in letteratura. Andersen et al. (1983) suggeriscono di classificare i pazienti ad un tempo “landmark” che corrisponde al punto iniziale della stima della curva di sopravvivenza, includendo solo i pazienti ancora a rischio al landmark. Il secondo metodo, proposto da Simon e Makuch (1984), consiste nell’aggiornamento dinamico dei “risk sets” dei due gruppi di trattamento tempo-dipendenti. Entrambi i metodi sono stati presentati in maniera euristica e senza specificare le quantità teoriche che corrispondono agli stimatori proposti. Perciò, l’interpretazione delle curve stimate dai due metodi non è mai stata chiarita. Quando l’interesse non è rivolto alla sopravvivenza globale ma alla predizione profilo-specifica, ovvero tenendo conto delle caratteristiche individuali dei soggetti, occorre utilizzare metodi di regressione parametrici o semi-parametrici. Il modello di Cox è quello più popolare ma in presenza di effetti tempo-dipendenti e/o di covariate tempo-dipendenti non può essere utilizzato per ottenere delle curve. Tra le possibili alternative sono stati considerati il modello parametrico di Hanley e Miettinen (2009) e il modello di regressione semi-parametrico basato sul landmark di Van Houwelingen (2007). Il primo è basato sulla stima della funzione azzardo nel tempo applicando una regressione logistica ad un dataset esteso creato dalla suddivisione del tempo di sopravvivenza osservato di ciascun soggetto in un certo numero di unità di tempo e trattando il numero di eventi in ogni singolo intervallo di tempo come una variabile casuale Binomiale. Il secondo metodo scaturisce dall’idea di utilizzare il modello di Cox su molteplici partizioni del dataset ciascuna creata partendo da un tempo landmark progressivo e includendo solo i soggetti a rischio al landmark; la classificazione del trattamento per questi pazienti è fissata a quel tempo consentendo di aggiornare dinamicamente il valore delle covariate tempo-dipendenti in ciascun modello e permettendo ai coefficienti stimati di variare nel tempo. Gli scopi del presente lavoro sono la revisione e lo sviluppo di metodi per: 1) descrivere la sopravvivenza in funzione di un covariata binaria tempo-dipendente sia da una prospettiva fissa sia dinamicamente nel tempo; 2) la valutazione dell’impatto su queste quantità dei fattori prognostici, in particolare il tempo di attesa al trapianto, utilizzando dei parametri interpretabili; 3) lo sviluppo di predizioni profilo-specifiche. Nella prima parte del lavoro si intende chiarire il significato delle le quantità teoriche stimate dai metodi landmark e Simon e Makuch. In aggiunta, si presenta un approccio innovativo basato su domande controfattuali e predizione dinamica, verificando la validità dei risultati attraverso delle simulazioni. Nella seconda parte, si presentano i modelli di regressione di Hanley-Miettinen e del landmark e si mostra come utilizzarli per ottenere la stima dell’effetto del tempo i attesa al trapianto e per produrre delle predizioni profilo-specifiche su dati reali inerenti a pazienti affetti da Leucemia Linfoblastica Acuta, confrontando la performance dei modelli attraverso delle simulazioni.
In clinical studies it is often of interest to compare the survival experience of patients in two or more treatment groups. In some situations the categorization is not fixed at baseline but changes during the follow-up, where patients, for example, start from an initial treatment and either continue it or switch to an alternative one after some time (waiting time). Thus, treatment is a binary non reversible time-dependent variable. A typical problem is comparing outcomes of chemotherapy vs stem-cell transplantation in Acute Lymphoblastic Leukemia (ALL) where patients are treated initially with chemotherapy and during the follow-up they can receive bone marrow transplant. In this context, the standard Kaplan-Meier method is unreliable since it is affected by the immortal time bias. Two alternative non-parametric approaches were proposed in the literature. Andersen et al. (1983) suggests to classify patients at a landmark time which corresponds to the starting point for the estimation of the Kaplan-Meier survival curve, involving only patients still at risk at the landmark. The second, proposed by Simon and Makuch (1984), consists in dynamically updating in time the risk set of the two time-dependent treatment groups. Both methods were presented mostly relying on heuristic bases and without specifying the theoretical quantities corresponding to the proposed estimators. Thus, the interpretations of the curves estimated by the two methods was never clarified. When the focus is not on the overall survival experience but rather on profile-specific prediction, i.e. accounting for the individual characteristics of the subjects, one must resort to semi-parametric or parametric regression models. The Cox model is the most popular one but in the presence of time-varying effects and/or time-dependent covariates it cannot be used to obtain survival curves. Among the possible alternatives we considered the full parametric model by Hanley and Miettinen (2009) and the semi-parametric landmark regression model by Van Houwelingen (2007). The first is based on estimating the hazard function over time by applying a logistic regression to an expanded dataset created by splitting the observed survival time of each subject into a number of time-units and to treat the number of events in every single interval as a Binomial random variable. The second originates from the idea of fitting the Cox model to multiple subsets of data, each one created starting from a sliding landmark time point and including only the subjects at risk at the landmark; the treatment classification for these patients is frozen at that time allowing to dynamically update the time-dependent covariates in each model and to let the parameter estimates to vary in time. The aims of the dissertation are reviewing and developing methods for: 1) the description of the survival experience according to a binary time-dependent treatment indicator both from a fixed perspective and dynamically update in time; 2) the assessment of the impact on these quantities of prognostic factors, in particular the waiting time to transplant, through interpretable parameters; 3) the development of profile-specific predictions. In the first part of this work we wish to clarify the theoretical quantities estimated by the landmark and Simon-Makuch methods. In addition, we present a novel approach based on counterfactual questions and dynamic prediction, checking the validity of our findings using simulations. In the second part, we review the Hanley-Miettinen and landmark regression models and we show how to use them to properly estimate the effect of waiting time to transplant and to make profile-specific dynamic predictions on a real dataset on ALL, comparing the performance of the two models using simulations.

Objectif: Définir l’effet des lipides et du traitement de la dyslipidémie sur les cancers de la prostate et de la vessie en utilisant différents devis d’étude et en tenant compte de la présence de plusieurs biais, particulièrement le biais du temps immortel. Devis: Le premier volet utilise un devis rétrospectif de type cas témoins. Un questionnaire semi-quantitatif de fréquence de consommation alimentaire validé a été utilisé. Le génotype COX2 de neuf polymorphisme nucléotidique unique (SNP) a été mesuré avec une plateforme Taqman. Des modèles de régression logistique non conditionnelle ont été utilisés pour comparer le risque de diagnostic d’un cancer de la prostate et l’interaction. Le deuxième volet utilise un devis rétrospectif de type cohorte basée sur les données administratives de la Régie de l’assurance-maladie du Québec (RAMQ). Des modèles de régression de Cox ont été employés pour mesurer l’association entre les statines et l’évolution du cancer de la vessie. Le troisième volet, porte un regard méthodologique sur le biais du temps immortel en examinant sa présence dans la littérature oncologique. Son importance est illustrée avec les données de la cohorte du deuxième volet, et les méthodes de correction possibles son appliquées. Résultats: L’étude du premier volet démontre qu’une diète riche en acides gras oméga-3 d’origine marine était fortement associée à un risque diminué de cancer de la prostate agressif (p<0.0001 pour la tendance). Le ratio de cote pour le cancer de la prostate du quartile supérieur d’oméga-3 était de 0.37 (IC 95% = 0.25 à 0.54). L’effet diététique était modifié par le génotype COX-2 SNP rs4648310 (p=0.002 pour l’interaction). En particulier, les hommes avec faible apport en oméga-3 et la variante rs4648310 avait un risque accru de cancer de la prostate (ratio de cote = 5.49, IC 95%=1.80 à 16.7), effet renversé par un apport en oméga-3 plus grand. L’étude du deuxième volet a observé que l’utilisation de statines est associée à une diminution du risque de progression du cancer de la vessie (risque relatif = 0.44, IC 95% = 0.20 à 0.96, p=0.039). Cette association était encore plus forte pour le décès de toute cause (HR = 0.57, 95% CI = 0.43 to 0.76, p=0.0001). L’effet des statines semble être dose-dépendant. L’étude du troisième volet démontre que le biais du temps immortel est fréquent et important dans les études épidémiologiques oncologiques. Il comporte plusieurs aspects dont certains sont mieux prévenus au stade du choix du devis d’étude et différentes méthodes statistiques permettent un contrôle de ce biais. Conclusion: 1) Une diète riche en oméga-3 aurait un effet protecteur pour le cancer de la prostate. 2) L’utilisation de statines aurait un effet protecteur sur la progression du cancer non invasif de la vessie. Les lipides semblent avoir un effet sur les cancers urologiques.
Purpose: To define the effects of dietary lipids and of treatment of dyslipidemia with statins on prostate and bladder cancers, using different epidemiologic study designs and accounting for biases, particularly immortal time bias. Study Design: The first part used a retrospective a case-control study design. Diet was assessed with a semi-quantitative food frequency questionnaire, and nine COX-2 tag single nucleotide polymorphisms (SNPs) were genotyped. We used logistic regression models to estimate odds ratios (ORs), 95% confidence intervals (CIs), and p-values for association and interaction. The second part used a retrospective cohort study design based on administrative databases of Québec, Canada. Cox regression models were used to measure association between statin use and bladder cancer evolution. The third part focuses on the immortal time bias by describing its presence in the oncologic literature. The importance of this bias is illustrated with data from the cohort used in the second part and statistical correction methods are applied. Results: The first part showed that an increasing intake of omega-3 fatty acids of marine origin was strongly associated with a decreased risk of aggressive prostate cancer (trend p<=0.0001). The OR (95% CI) for prostate cancer comparing the highest to the lowest quartile of omega-3 intake was of 0.37 (0.25 – 0.54). The dietary effect was modified by the rs4648310 COX-2 SNP (interaction p=0.02). This reflected the observation that men with low marine omega-3 intake and the variant rs4648310 SNP had an increased risk of disease (OR = 5.49; 95% CI: 1.80-16.7), which was reversed by increasing intake of marine omega-3. The second part showed that statin use was associated with a decreased risk of bladder cancer progression (HR = 0.44, 95% CI = 0.20 to 0.96, p=0.0388). The inverse association was even stronger for risk of mortality from all causes (HR = 0.57, 95% CI = 0.43 to 0.76, p=0.0001). The statin use effect appears dose-dependent. The third part showed that the immortal time bias is frequent and important in many epidemiological studies in oncology. It has many aspects and some of these are better prevented at time of study design selection. Various statistical methods also allowed control of this bias. Conclusion. 1) Dietary omega-3 appears to decrease prostate cancer risk. 2) Statin use appears to decrease risk of bladder cancer progression. Lipids seem to have an effect on urological cancers.