The formal approach to quantitative causal inference in epidemiology: misguided or misrepresented?

Two recent articles, one by Vandenbroucke, Broadbent and Pearce (henceforth VBP) and the other by Krieger and Davey Smith (henceforth KDS), criticize what these two sets of authors characterize as the mainstream of the modern ‘causal inference’ school in epidemiology. The criticisms made by these authors are severe; VBP label the field both ‘wrong in theory’ and ‘wrong in practice’, and KDS—at least in some settings—feel that the field not only ‘bark[s] up the wrong tree’ but ‘miss[es] the forest entirely’. More specifically, the school of thought, and the concepts and methods within it, are painted as being applicable only to a very narrow range of investigations, to the exclusion of most of the important questions and study designs in modern epidemiology, such as the effects of genetic variants, the study of ethnic and gender disparities and the use of study designs that do not closely mirror randomized controlled trials (RCTs). Furthermore, the concepts and methods are painted as being potentially highly misleading even within this narrow range in which they are deemed applicable. We believe that most of VBP’s and KDS’s criticisms stem from a series of misconceptions about the approach they criticize. In this response, therefore, we aim first to paint a more accurate picture of the formal causal inference approach, and then to outline the key misconceptions underlying VBP’s and KDS’s critiques. KDS in particular criticize directed acyclic graphs (DAGs), using three examples to do so. Their discussion highlights further misconceptions concerning the role of DAGs in causal inference, and so we devote the third section of the paper to addressing these. In our Discussion we present further objections we have to the arguments in the two papers, before concluding that the clarity gained from adopting a rigorous framework is an asset, not an obstacle, to answering more reliably a very wide range of causal questions using data from observational studies of many different designs

Causal mediation analysis with multiple mediators.

In diverse fields of empirical research-including many in the biological sciences-attempts are made to decompose the effect of an exposure on an outcome into its effects via a number of different pathways. For example, we may wish to separate the effect of heavy alcohol consumption on systolic blood pressure (SBP) into effects via body mass index (BMI), via gamma-glutamyl transpeptidase (GGT), and via other pathways. Much progress has been made, mainly due to contributions from the field of causal inference, in understanding the precise nature of statistical estimands that capture such intuitive effects, the assumptions under which they can be identified, and statistical methods for doing so. These contributions have focused almost entirely on settings with a single mediator, or a set of mediators considered en bloc; in many applications, however, researchers attempt a much more ambitious decomposition into numerous path-specific effects through many mediators. In this article, we give counterfactual definitions of such path-specific estimands in settings with multiple mediators, when earlier mediators may affect later ones, showing that there are many ways in which decomposition can be done. We discuss the strong assumptions under which the effects are identified, suggesting a sensitivity analysis approach when a particular subset of the assumptions cannot be justified. These ideas are illustrated using data on alcohol consumption, SBP, BMI, and GGT from the Izhevsk Family Study. We aim to bridge the gap from "single mediator theory" to "multiple mediator practice," highlighting the ambitious nature of this endeavor and giving practical suggestions on how to proceed

gesttools: General Purpose G-Estimation in R

In this paper we present gesttools, a series of general purpose, user friendly functions with which to perform g-estimation of structural nested mean models (SNMMs) for time-varying exposures and outcomes in R. The package implements the g-estimation methods found in Vansteelandt and Sjolander (2016) and Dukes and Vansteelandt (2018), and is capable of analysing both end of study and time-varying outcome data that are either binary or continuous, or exposure variables that are either binary, continuous, or categorical. It also allows for the fitting of SNMMs with time-varying causal effects, effect modification by other variables, or both, as well as support for censored data using inverse weighting. We outline the theory underpinning these methods, as well as describing the SNMMs that can be fitted by the software. The package is demonstrated using simulated, and real-world inspired datasets

Intensive care unit (ICU)-acquired bacteraemia and ICU mortality and discharge:Addressing time-varying confounding using appropriate methodology

Background: Studies often ignore time-varying confounding or may use inappropriate methodology to adjust for time-varying confounding. Aim: To estimate the effect of intensive care unit (ICU)-acquired bacteraemia on ICU mortality and discharge using appropriate methodology. Methods: Marginal structural models with inverse probability weighting were used to estimate the ICU mortality and discharge associated with ICU-acquired bacteraemia among patients who stayed more than two days at the general ICU of a London teaching hospital and remained bacteraemia-free during those first two days. For comparison, the same associations were evaluated with (i) a conventional Cox model, adjusting only for baseline confounders and (ii) a Cox model adjusting for baseline and time-varying confounders. Findings: Using the marginal structural model with inverse probability weighting, bacteraemia was associated with an increase in ICU mortality (cause-specific hazard ratio (CSHR): 1.29; 95% confidence interval (CI): 1.02-1.63)and a decrease in discharge (CSHR: 0.52; 95% CI: 0.45-0.60). By 60 days, among patients still in the ICU after two days and without prior bacteraemia, 8.0% of ICU deaths could be prevented by preventing all ICU-acquired bacteraemia cases. The conventional Cox model adjusting for time-varying confounders gave substantially different results [for ICU mortality, CSHR: 1.08 (95% CI: 0.88-1.32); for discharge, CSHR: 0.68 (95% CI: 0.60-0.77)]. Conclusion: In this study, even after adjusting for the timing of acquiring bacteraemia and time-varying confounding using inverse probability weighting for marginal structura

A Comparison of Three Methods of Mendelian Randomization when the Genetic Instrument, the Risk Factor and the Outcome Are All Binary

The method of instrumental variable (referred to as Mendelian randomization when the instrument is a genetic variant) has been initially developed to infer on a causal effect of a risk factor on some outcome of interest in a linear model. Adapting this method to nonlinear models, however, is known to be problematic. In this paper, we consider the simple case when the genetic instrument, the risk factor, and the outcome are all binary. We compare via simulations the usual two-stages estimate of a causal odds-ratio and its adjusted version with a recently proposed estimate in the context of a clinical trial with noncompliance. In contrast to the former two, we confirm that the latter is (under some conditions) a valid estimate of a causal odds-ratio defined in the subpopulation of compliers, and we propose its use in the context of Mendelian randomization. By analogy with a clinical trial with noncompliance, compliers are those individuals for whom the presence/absence of the risk factor X is determined by the presence/absence of the genetic variant Z (i.e., for whom we would observe X = Z whatever the alleles randomly received at conception). We also recall and illustrate the huge variability of instrumental variable estimates when the instrument is weak (i.e., with a low percentage of compliers, as is typically the case with genetic instruments for which this proportion is frequently smaller than 10%) where the inter-quartile range of our simulated estimates was up to 18 times higher compared to a conventional (e.g., intention-to-treat) approach. We thus conclude that the need to find stronger instruments is probably as important as the need to develop a methodology allowing to consistently estimate a causal odds-ratio

Does appropriate empiric antibiotic therapy modify intensive care unit-acquired Enterobacteriaceae bacteraemia mortality and discharge?

BACKGROUND: Conflicting results have been found regarding outcomes of intensive care unit (ICU)-acquired Enterobacteriaceae bacteraemia and the potentially modifying effect of appropriate empiric antibiotic therapy. AIM: To evaluate these associations while adjusting for potential time-varying confounding using methods from the causal inference literature. METHODS: Patients who stayed more than two days in two general ICUs in England between 2002 and 2006 were included in this cohort study. Marginal structural models with inverse probability weighting were used to estimate the mortality and discharge associated with Enterobacteriaceae bacteraemia and the impact of appropriate empiric antibiotic therapy on these outcomes. FINDINGS: Among 3411 ICU admissions, 195 (5.7%) ICU-acquired Enterobacteriaceae bacteraemia cases occurred. Enterobacteriaceae bacteraemia was associated with an increased daily risk of ICU death [cause-specific hazard ratio (HR): 1.48; 95% confidence interval (CI): 1.10-1.99] and a reduced daily risk of ICU discharge (HR: 0.66; 95% CI: 0.54-0.80). Appropriate empiric antibiotic therapy did not significantly modify ICU mortality (HR: 1.08; 95% CI: 0.59-1.97) or discharge (HR: 0.91; 95% CI: 0.63-1.32). CONCLUSION: ICU-acquired Enterobacteriaceae bacteraemia was associated with an increased daily risk of ICU mortality. Furthermore, the daily discharge rate was also lower after acquiring infection, even when adjusting for time-varying confounding using appropriate methodology. No evidence was found for a beneficial modifying effect of appropriate empiric antibiotic therapy on ICU mortality and discharge

Integrated multiple mediation analysis: A robustness–specificity trade-off in causal structure

Recent methodological developments in causal mediation analysis have addressed several issues regarding multiple mediators. However, these developed methods differ in their definitions of causal parameters, assumptions for identification, and interpretations of causal effects, making it unclear which method ought to be selected when investigating a given causal effect. Thus, in this study, we construct an integrated framework, which unifies all existing methodologies, as a standard for mediation analysis with multiple mediators. To clarify the relationship between existing methods, we propose four strategies for effect decomposition: two-way, partially forward, partially backward, and complete decompositions. This study reveals how the direct and indirect effects of each strategy are explicitly and correctly interpreted as path-specific effects under different causal mediation structures. In the integrated framework, we further verify the utility of the interventional analogues of direct and indirect effects, especially when natural direct and indirect effects cannot be identified or when cross-world exchangeability is invalid. Consequently, this study yields a robustness–specificity trade-off in the choice of strategies. Inverse probability weighting is considered for estimation. The four strategies are further applied to a simulation study for performance evaluation and for analyzing the Risk Evaluation of Viral Load Elevation and Associated Liver Disease/Cancer data set from Taiwan to investigate the causal effect of hepatitis C virus infection on mortality

Two rapid assays for screening of patulin biodegradation

Artículo sobre distintos ensayos para comprobar la biodegradación de la patulinaThe mycotoxin patulin is produced by the blue mould pathogen Penicillium expansum in rotting apples during postharvest storage. Patulin is toxic to a wide range of organisms, including humans, animals, fungi and bacteria. Wash water from apple packing and processing houses often harbours patulin and fungal spores, which can contaminate the environment. Ubiquitous epiphytic yeasts, such as Rhodosporidium kratochvilovae strain LS11 which is a biocontrol agent of P. expansum in apples, have the capacity to resist the toxicity of patulin and to biodegrade it. Two non-toxic products are formed. One is desoxypatulinic acid. The aim of the work was to develop rapid, high-throughput bioassays for monitoring patulin degradation in multiple samples. Escherichia coli was highly sensitive to patulin, but insensitive to desoxypatulinic acid. This was utilized to develop a detection test for patulin, replacing time-consuming thin layer chromatography or high-performance liquid chromatography. Two assays for patulin degradation were developed, one in liquid medium and the other in semi-solid medium. Both assays allow the contemporary screening of a large number of samples. The liquid medium assay utilizes 96-well microtiter plates and was optimized for using a minimum of patulin. The semisolid medium assay has the added advantage of slowing down the biodegradation, which allows the study and isolation of transient degradation products. The two assays are complementary and have several areas of utilization, from screening a bank of microorganisms for biodegradation ability to the study of biodegradation pathways

Outcome in patients perceived as receiving excessive care across different ethical climates: a prospective study in 68 intensive care units in Europe and the USA.

PURPOSE: Whether the quality of the ethical climate in the intensive care unit (ICU) improves the identification of patients receiving excessive care and affects patient outcomes is unknown. METHODS: In this prospective observational study, perceptions of excessive care (PECs) by clinicians working in 68 ICUs in Europe and the USA were collected daily during a 28-day period. The quality of the ethical climate in the ICUs was assessed via a validated questionnaire. We compared the combined endpoint (death, not at home or poor quality of life at 1 year) of patients with PECs and the time from PECs until written treatment-limitation decisions (TLDs) and death across the four climates defined via cluster analysis. RESULTS: Of the 4747 eligible clinicians, 2992 (63%) evaluated the ethical climate in their ICU. Of the 321 and 623 patients not admitted for monitoring only in ICUs with a good (n = 12, 18%) and poor (n = 24, 35%) climate, 36 (11%) and 74 (12%), respectively were identified with PECs by at least two clinicians. Of the 35 and 71 identified patients with an available combined endpoint, 100% (95% CI 90.0-1.00) and 85.9% (75.4-92.0) (P = 0.02) attained that endpoint. The risk of death (HR 1.88, 95% CI 1.20-2.92) or receiving a written TLD (HR 2.32, CI 1.11-4.85) in patients with PECs by at least two clinicians was higher in ICUs with a good climate than in those with a poor one. The differences between ICUs with an average climate, with (n = 12, 18%) or without (n = 20, 29%) nursing involvement at the end of life, and ICUs with a poor climate were less obvious but still in favour of the former. CONCLUSION: Enhancing the quality of the ethical climate in the ICU may improve both the identification of patients receiving excessive care and the decision-making process at the end of life

Support vector machine versus logistic regression modeling for prediction of hospital mortality in critically ill patients with haematological malignancies

Background: Several models for mortality prediction have been constructed for critically ill patients with haematological malignancies in recent years. These models have proven to be equally or more accurate in predicting hospital mortality in patients with haematological malignancies than ICU severity of illness scores such as the APACHE II or SAPS II [1]. The objective of this study is to compare the accuracy of predicting hospital mortality in patients with haematological malignancies admitted to the ICU between models based on multiple logistic regression (MLR) and support vector machine (SVM) based models. Methods: 352 patients with haematological malignancies admitted to the ICU between 1997 and 2006 for a life-threatening complication were included. 252 patient records were used for training of the models and 100 were used for validation. In a first model 12 input variables were included for comparison between MLR and SVM. In a second more complex model 17 input variables were used. MLR and SVM analysis were performed independently from each other. Discrimination was evaluated using the area under the receiver operating characteristic (ROC) curves (+/- SE). Results: The area under ROC curve for the MLR and SVM in the validation data set were 0.768 (+/- 0.04) vs. 0.802 (+/- 0.04) in the first model (p = 0.19) and 0.781 (+/- 0.05) vs. 0.808 (+/- 0.04) in the second more complex model (p = 0.44). SVM needed only 4 variables to make its prediction in both models, whereas MLR needed 7 and 8 variables in the first and second model respectively. Conclusion: The discriminative power of both the MLR and SVM models was good. No statistically significant differences were found in discriminative power between MLR and SVM for prediction of hospital mortality in critically ill patients with haematological malignancies