scholarly journals Marginal Structural Models with Counterfactual Effect Modifiers

2018 ◽  
Vol 14 (1) ◽  
Author(s):  
Wenjing Zheng ◽  
Zhehui Luo ◽  
Mark J van der Laan

Abstract In health and social sciences, research questions often involve systematic assessment of the modification of treatment causal effect by patient characteristics. In longitudinal settings, time-varying or post-intervention effect modifiers are also of interest. In this work, we investigate the robust and efficient estimation of the Counterfactual-History-Adjusted Marginal Structural Model (van der Laan MJ, Petersen M. Statistical learning of origin-specific statically optimal individualized treatment rules. Int J Biostat. 2007;3), which models the conditional intervention-specific mean outcome given a counterfactual modifier history in an ideal experiment. We establish the semiparametric efficiency theory for these models, and present a substitution-based, semiparametric efficient and doubly robust estimator using the targeted maximum likelihood estimation methodology (TMLE, e.g. van der Laan MJ, Rubin DB. Targeted maximum likelihood learning. Int J Biostat. 2006;2, van der Laan MJ, Rose S. Targeted learning: causal inference for observational and experimental data, 1st ed. Springer Series in Statistics. Springer, 2011). To facilitate implementation in applications where the effect modifier is high dimensional, our third contribution is a projected influence function (and the corresponding projected TMLE estimator), which retains most of the robustness of its efficient peer and can be easily implemented in applications where the use of the efficient influence function becomes taxing. We compare the projected TMLE estimator with an Inverse Probability of Treatment Weighted estimator (e.g. Robins JM. Marginal structural models. In: Proceedings of the American Statistical Association. Section on Bayesian Statistical Science, 1-10. 1997a, Hernan MA, Brumback B, Robins JM. Marginal structural models to estimate the causal effect of zidovudine on the survival of HIV-positive men. Epidemiology. 2000;11:561–570), and a non-targeted G-computation estimator (Robins JM. A new approach to causal inference in mortality studies with sustained exposure periods - application to control of the healthy worker survivor effect. Math Modell. 1986;7:1393–1512.). The comparative performance of these estimators is assessed in a simulation study. The use of the projected TMLE estimator is illustrated in a secondary data analysis for the Sequenced Treatment Alternatives to Relieve Depression (STAR*D) trial where effect modifiers are subject to missing at random.

2016 ◽  
Vol 12 (1) ◽  
pp. 233-252 ◽  
Author(s):  
Wenjing Zheng ◽  
Maya Petersen ◽  
Mark J. van der Laan

Abstract In social and health sciences, many research questions involve understanding the causal effect of a longitudinal treatment on mortality (or time-to-event outcomes in general). Often, treatment status may change in response to past covariates that are risk factors for mortality, and in turn, treatment status may also affect such subsequent covariates. In these situations, Marginal Structural Models (MSMs), introduced by Robins (1997. Marginal structural models Proceedings of the American Statistical Association. Section on Bayesian Statistical Science, 1–10), are well-established and widely used tools to account for time-varying confounding. In particular, a MSM can be used to specify the intervention-specific counterfactual hazard function, i. e. the hazard for the outcome of a subject in an ideal experiment where he/she was assigned to follow a given intervention on their treatment variables. The parameters of this hazard MSM are traditionally estimated using the Inverse Probability Weighted estimation Robins (1999. Marginal structural models versus structural nested models as tools for causal inference. In: Statistical models in epidemiology: the environment and clinical trials. Springer-Verlag, 1999:95–134), Robins et al. (2000), (IPTW, van der Laan and Petersen (2007. Causal effect models for realistic individualized treatment and intention to treat rules. Int J Biostat 2007;3:Article 3), Robins et al. (2008. Estimaton and extrapolation of optimal treatment and testing strategies. Statistics in Medicine 2008;27(23):4678–721)). This estimator is easy to implement and admits Wald-type confidence intervals. However, its consistency hinges on the correct specification of the treatment allocation probabilities, and the estimates are generally sensitive to large treatment weights (especially in the presence of strong confounding), which are difficult to stabilize for dynamic treatment regimes. In this paper, we present a pooled targeted maximum likelihood estimator (TMLE, van der Laan and Rubin (2006. Targeted maximum likelihood learning. The International Journal of Biostatistics 2006;2:1–40)) for MSM for the hazard function under longitudinal dynamic treatment regimes. The proposed estimator is semiparametric efficient and doubly robust, offering bias reduction over the incumbent IPTW estimator when treatment probabilities may be misspecified. Moreover, the substitution principle rooted in the TMLE potentially mitigates the sensitivity to large treatment weights in IPTW. We compare the performance of the proposed estimator with the IPTW and a on-targeted substitution estimator in a simulation study.


2014 ◽  
Vol 2 (2) ◽  
pp. 147-185 ◽  
Author(s):  
Maya Petersen ◽  
Joshua Schwab ◽  
Susan Gruber ◽  
Nello Blaser ◽  
Michael Schomaker ◽  
...  

AbstractThis paper describes a targeted maximum likelihood estimator (TMLE) for the parameters of longitudinal static and dynamic marginal structural models. We consider a longitudinal data structure consisting of baseline covariates, time-dependent intervention nodes, intermediate time-dependent covariates, and a possibly time-dependent outcome. The intervention nodes at each time point can include a binary treatment as well as a right-censoring indicator. Given a class of dynamic or static interventions, a marginal structural model is used to model the mean of the intervention-specific counterfactual outcome as a function of the intervention, time point, and possibly a subset of baseline covariates. Because the true shape of this function is rarely known, the marginal structural model is used as a working model. The causal quantity of interest is defined as the projection of the true function onto this working model. Iterated conditional expectation double robust estimators for marginal structural model parameters were previously proposed by Robins (2000, 2002) and Bang and Robins (2005). Here we build on this work and present a pooled TMLE for the parameters of marginal structural working models. We compare this pooled estimator to a stratified TMLE (Schnitzer et al. 2014) that is based on estimating the intervention-specific mean separately for each intervention of interest. The performance of the pooled TMLE is compared to the performance of the stratified TMLE and the performance of inverse probability weighted (IPW) estimators using simulations. Concepts are illustrated using an example in which the aim is to estimate the causal effect of delayed switch following immunological failure of first line antiretroviral therapy among HIV-infected patients. Data from the International Epidemiological Databases to Evaluate AIDS, Southern Africa are analyzed to investigate this question using both TML and IPW estimators. Our results demonstrate practical advantages of the pooled TMLE over an IPW estimator for working marginal structural models for survival, as well as cases in which the pooled TMLE is superior to its stratified counterpart.


2021 ◽  
Vol 21 (1) ◽  
Author(s):  
Amir Almasi-Hashiani ◽  
Saharnaz Nedjat ◽  
Reza Ghiasvand ◽  
Saeid Safiri ◽  
Maryam Nazemipour ◽  
...  

Abstract Objectives The relationship between reproductive factors and breast cancer (BC) risk has been investigated in previous studies. Considering the discrepancies in the results, the aim of this study was to estimate the causal effect of reproductive factors on BC risk in a case-control study using the double robust approach of targeted maximum likelihood estimation. Methods This is a causal reanalysis of a case-control study done between 2005 and 2008 in Shiraz, Iran, in which 787 confirmed BC cases and 928 controls were enrolled. Targeted maximum likelihood estimation along with super Learner were used to analyze the data, and risk ratio (RR), risk difference (RD), andpopulation attributable fraction (PAF) were reported. Results Our findings did not support parity and age at the first pregnancy as risk factors for BC. The risk of BC was higher among postmenopausal women (RR = 3.3, 95% confidence interval (CI) = (2.3, 4.6)), women with the age at first marriage ≥20 years (RR = 1.6, 95% CI = (1.3, 2.1)), and the history of oral contraceptive (OC) use (RR = 1.6, 95% CI = (1.3, 2.1)) or breastfeeding duration ≤60 months (RR = 1.8, 95% CI = (1.3, 2.5)). The PAF for menopause status, breastfeeding duration, and OC use were 40.3% (95% CI = 39.5, 40.6), 27.3% (95% CI = 23.1, 30.8) and 24.4% (95% CI = 10.5, 35.5), respectively. Conclusions Postmenopausal women, and women with a higher age at first marriage, shorter duration of breastfeeding, and history of OC use are at the higher risk of BC.


Biostatistics ◽  
2018 ◽  
Vol 21 (1) ◽  
pp. 172-185 ◽  
Author(s):  
Pål Christie Ryalen ◽  
Mats Julius Stensrud ◽  
Sophie Fosså ◽  
Kjetil Røysland

Abstract In marginal structural models (MSMs), time is traditionally treated as a discrete parameter. In survival analysis on the other hand, we study processes that develop in continuous time. Therefore, Røysland (2011. A martingale approach to continuous-time marginal structural models. Bernoulli 17, 895–915) developed the continuous-time MSMs, along with continuous-time weights. The continuous-time weights are conceptually similar to the inverse probability weights that are used in discrete time MSMs. Here, we demonstrate that continuous-time MSMs may be used in practice. First, we briefly describe the causal model assumptions using counting process notation, and we suggest how causal effect estimates can be derived by calculating continuous-time weights. Then, we describe how additive hazard models can be used to find such effect estimates. Finally, we apply this strategy to compare medium to long-term differences between the two prostate cancer treatments radical prostatectomy and radiation therapy, using data from the Norwegian Cancer Registry. In contrast to the results of a naive analysis, we find that the marginal cumulative incidence of treatment failure is similar between the strategies, accounting for the competing risk of other death.


Epidemiology ◽  
2000 ◽  
Vol 11 (5) ◽  
pp. 550-560 ◽  
Author(s):  
James M. Robins ◽  
Miguel Ángel Hernán ◽  
Babette Brumback

2019 ◽  
Vol 112 (3) ◽  
pp. e178
Author(s):  
Soudeh Ansari ◽  
Michael P. LaValley ◽  
Sara Lodi ◽  
Brooke Hayward ◽  
Gilbert L. Mottla ◽  
...  

2014 ◽  
Vol 2 (1) ◽  
pp. 95-108 ◽  
Author(s):  
Anna L. Decker ◽  
Alan Hubbard ◽  
Catherine M. Crespi ◽  
Edmund Y.W. Seto ◽  
May C. Wang

AbstractWhile child and adolescent obesity is a serious public health concern, few studies have utilized parameters based on the causal inference literature to examine the potential impacts of early intervention. The purpose of this analysis was to estimate the causal effects of early interventions to improve physical activity and diet during adolescence on body mass index (BMI), a measure of adiposity, using improved techniques. The most widespread statistical method in studies of child and adolescent obesity is multivariable regression, with the parameter of interest being the coefficient on the variable of interest. This approach does not appropriately adjust for time-dependent confounding, and the modeling assumptions may not always be met. An alternative parameter to estimate is one motivated by the causal inference literature, which can be interpreted as the mean change in the outcome under interventions to set the exposure of interest. The underlying data-generating distribution, upon which the estimator is based, can be estimated via a parametric or semi-parametric approach. Using data from the National Heart, Lung, and Blood Institute Growth and Health Study, a 10-year prospective cohort study of adolescent girls, we estimated the longitudinal impact of physical activity and diet interventions on 10-year BMI z-scores via a parameter motivated by the causal inference literature, using both parametric and semi-parametric estimation approaches. The parameters of interest were estimated with a recently released R package, ltmle, for estimating means based upon general longitudinal treatment regimes. We found that early, sustained intervention on total calories had a greater impact than a physical activity intervention or non-sustained interventions. Multivariable linear regression yielded inflated effect estimates compared to estimates based on targeted maximum-likelihood estimation and data-adaptive super learning. Our analysis demonstrates that sophisticated, optimal semiparametric estimation of longitudinal treatment-specific means via ltmle provides an incredibly powerful, yet easy-to-use tool, removing impediments for putting theory into practice.


Blood ◽  
2009 ◽  
Vol 114 (22) ◽  
pp. 3365-3365
Author(s):  
Matthieu Resche-Rigon ◽  
Marie Robin ◽  
Regis Peffault de Latour ◽  
Sylvie Chevret ◽  
Gerard P Socie

Abstract Abstract 3365 Poster Board III-253 Introduction: Although allogeneic SCT with RIC has now gained wide acceptance, its eventual benefit again non-transplant approach is largely unknown (outside the setting of large randomized trials). When evaluating the impact on survival of reduced intensity conditioning in malignant hematological diseases, standard estimations based on Cox regression from observational databases could be biased because they ignore covariates that confound treatment decision. In this setting, we applied and compared two different statistical methods that were developed to control for confounding in estimating exposure (or treatment) effect from epidemiological studies. Patients and Methods: The statistical challenge was that allograft tended to be given when a patient was in advanced phase of his/her hematological malignancy, so that treatment was confounded by performance indicators, which in turn lie on the causal pathway between treatment and outcome. Thus, comparison of outcome first used propensity score (PS) analyses that attempt to create a comparison group of non-treated patients that closely resembles the group of treated patients by matching for the likelihood that a given patient has received the treatment. Then, we used marginal structural models (MSMs) that consist in creating, by using inverse probability of treatment weights, a pseudo-population in which the probability of treatment does no longer depend on covariates, and the effect of treatment on outcome is the same as in the original population. Result: Reduced intensity conditioning allograft was performed in 82 patients with chemotherapy-sensitive patients relapsing after autologous transplantation. Patients with myeloma (MM, 23 pts), follicular lymphoma (FL, 28 pts) or Hodgkin disease (HD, 31 pts), were compared to 276 patients who relapsed after autologous transplantation but did not underwent allogeneic stem cell transplantation (142 MM, 115 FL and 19 HD). From original datasets, 21 (91%) matched pairs could be constituted from MM patients, as compared to 19 (68%) of the FL patients, down to 15 (48%) of the HD patients. Based on these PS-matched samples, a significant benefit of reduced intensity conditioning as compared with non allografted patients was observed in MM, with estimated hazard ratio (HR) of death at 0.34 (95% confidence interval, CI: 0.14-0.88), as well as in FL (HR= 0.78, 95%CI: 0.27;2.30) and in HD (HR= 0.24; 95%CI: 0.09-0.62). MSM-based analyses that applied to the reweighted populations confirmed these trends towards survival benefits in FL, though partially erased in MM and HD. Conclusions: We reported the application of marginal structural models, a new class of causal models to estimate the effect of nonrandomized treatments as an alternative to PS based approaches in small samples. We expect that an increasing number of physicians involved in clinical cohorts become familiar with these novel and appealing quantitative methods when assessing innovative treatment effects. Disclosures: No relevant conflicts of interest to declare.


Sign in / Sign up

Export Citation Format

Share Document