Warning: fopen(/home/virtual/epih/journal/upload/ip_log/ip_log_2024-04.txt): failed to open stream: Permission denied in /home/virtual/lib/view_data.php on line 83 Warning: fwrite() expects parameter 1 to be resource, boolean given in /home/virtual/lib/view_data.php on line 84 Applied mediation analyses: a review and tutorial
Skip Navigation
Skip to contents

Epidemiol Health : Epidemiology and Health

OPEN ACCESS
SEARCH
Search

Articles

Page Path
HOME > Epidemiol Health > Volume 39; 2017 > Article
Review
Applied mediation analyses: a review and tutorial
Theis Lange1,2orcid, Kim Wadt Hansen3orcid, Rikke Sørensen4, Søren Galatius3
Epidemiol Health 2017;39:e2017035.
DOI: https://doi.org/10.4178/epih.e2017035
Published online: August 6, 2017

1Department of Biostatistics, University of Copenhagen, Copenhagen, Denmark

2Center for Statistical Science, Peking University, Copenhagen, Denmark

3Department of Cardiology, University Hospital Bispebjerg, Copenhagen, Denmark

4Department of Cardiology, Rigshospitalet, Copenhagen, Denmark

Correspondence: Theis Lange Department of Biostatistics, University of Copenhagen, 5 Østre Farimagsgade, Copenhagen 1014, Denmark E- mail: thlan@sund.ku.dk
• Received: May 19, 2017   • Accepted: August 6, 2017

©2017, Korean Society of Epidemiology

This is an open-access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

  • 23,456 Views
  • 993 Download
  • 78 Web of Science
  • 80 Crossref
  • 76 Scopus
  • In recent years, mediation analysis has emerged as a powerful tool to disentangle causal pathways from an exposure/treatment to clinically relevant outcomes. Mediation analysis has been applied in scientific fields as diverse as labour market relations and randomized clinical trials of heart disease treatments. In parallel to these applications, the underlying mathematical theory and computer tools have been refined. This combined review and tutorial will introduce the reader to modern mediation analysis including: the mathematical framework; required assumptions; and software implementation in the R package medflex. All results are illustrated using a recent study on the causal pathways stemming from the early invasive treatment of acute coronary syndrome, for which the rich Danish population registers allow us to follow patients’ medication use and more after being discharged from hospital.
As any statistician collaborating with applied researchers (be they doctors, psychologists, or social scientists) knows, once a cause-and-effect relationship has been established, the next point on the agenda is very often “How does this effect come about? What are the underlying mechanisms?” Indeed, such questions have been driving science forward for centuries. Mediation analysis in the general meaning of the term refers to a collection of tools and ways of thinking designed to help applied researchers identify, formalize, and quantify possible mechanisms (i.e., causal pathways) linking a cause to an effect. To name just one example, the search for the mechanisms linking exposure to contaminations with subsequent disease, which was already under way in the 16th century, culminating in Louis Pasteur’s identification of bacteria as the “mediating factor,” was based on reasoning about mediation.
In contrast, statistical mediation analysis, which will be the object of interest in this tutorial, is concerned with quantifying specific causal pathways described by one or more measurements of specific variables that are either assumed or have been shown to be affected by the exposure and themselves affect the outcome. Statistical mediation analysis is broadly said to have been initiated in the seminal 1986 paper by Baron & Kenny [1]. As will be demonstrated in this tutorial, one of the main contributions of statistical mediation analysis is to translate the loose or intuitive concepts of, for example, Pasteur’s “mediating factors” into statements expressed as statistical models using mathematical formalism. Another important contribution, which will also be thoroughly discussed in this tutorial, is the derivation of the assumptions that must be satisfied before causal pathways can meaningfully be estimated from data.
For the remainder of this tutorial, mediation analysis will be taken to mean statistical mediation analysis only. The reader is expected to be familiar with statistical modelling and inference, as well as the distinctions between statistical associations and causal effects (i.e., why observational studies are harder to interpret than randomized studies). A prior knowledge of theoretical causal inference in general, or mediation analysis in particular, is neither assumed nor required.
The rest of this tutorial is structured as follows. Section 1 introduces the case that will be used to illustrate theoretical concepts throughout. In section 2, the mathematical framework for mediation analysis is introduced and the required assumptions are presented and discussed. Methods for estimation in real-life settings are presented in section 3, and the tools applied to the illustrative case are described in section 4. Finally, some current methodological challenges within mediation analysis (in particular, sensitivity analyses and multiple mediators) are discussed in section 5. Due to the nature of the illustrative case, special focus will be given to handling complex mediators. The simplementation will largely build on the recently released R package medflex [2].
Acute coronary syndrome (ACS) presents as a cardiac emergency caused by sudden obstruction of a coronary artery, most frequently due to thrombus formation in an existing atherosclerotic lesion in the vessel wall. In the acute phase, treatment aims to prevent sudden cardiac death and complications by halting the progression of thrombus formation, managing symptoms, and identifying and treating coronary obstructions; the latter goal involves early cardiac catheterization. Once stabilized, patients receive secondary preventive medication and undergo risk factor modification to prevent future cardiovascular events, including death.
Using Danish register data, we have previously established [3]that in a population of patients with a first hospitalization for ACS, the use of an early invasive treatment strategy was associated with a lower short-term risk of cardiac death and readmission for myocardial infarction than a conservative approach. It has been speculated that some or, in selected subgroups, all of the long-term benefit provided by an invasive treatment strategy is mediated through better secondary preventive medical therapy. In this case study, we will explore the relationship between an early invasive treatment strategy, secondary preventive medication, and death from all causes.
Following previous research and current guidelines [4], we define an early invasive strategy as cardiac catheterization performed within 72 hours of index hospitalization, thus assuming an intention to treat with reperfusion therapy, if appropriate based on coronary anatomy. In contrast, we define a conservative approach to be when an angiographic assessment was performed more than 72 hours after the index hospitalization or not at all.
The general recommendations for secondary preventive medical therapy in the Danish national guidelines for treating ACS include: lifelong aspirin, a P2Y12 receptor inhibitor (clopidogrel, prasugrel, or ticagrelor) for 12 months, lifelong statin therapy, and treatment with a β-blocker for at least 2 years. We defined a person as adhering to secondary medication for a given drug if a prescription was filled within 30 days of discharge or if the patient was in possession of a sufficient quantity of the drug to cover the initial 30 days after discharge (see Hansen et al. [5] for further references). As the recommended secondary mediation includes 4 drugs, we have 4 variables that may function as potential mediators. The outcome was death from all causes during follow-up. As this was an observational study, we attempted to control confounding by including information on age, sex, calendar year, net household income, educational level, cohabitation status, myocardial infarction, cardiac arrhythmia, heart failure, pulmonary oedema, cardiogenic shock, valvular heart disease, cerebrovascular disease, cancer, chronic obstructive pulmonary disease, diabetes with complications, acute and chronic renal failure, sepsis, pneumonia, anaemia, respiratory insufficiency, prior revascularization, prior in-hospital bleeding, and the use of antihypertensive medications, aspirin, lipid-lowering drugs, vitamin K antagonists, glucose-lowering drugs, loop diuretics, or chronic obstructive pulmonary disease medication.
As all Danes are given a unique identification number at birth, which is recorded in all subsequent dealings with the health care system, we were able to identify all patients hospitalized with a first ACS. We excluded patient who were discharged on the day of admission to ensure that there had indeed been time to conduct proper electrocardiographic monitoring and sequential measurements of cardiac troponins. Only patients aged 30-90 years were included. In patients younger than 30 years, it was assumed that atherosclerosis may not be have been the predominant underlying cause of ACS; similarly, patients older than 90 years were excluded, as they were deemed too frail for invasive procedures. Finally, we restricted the study cohort to those patients who had not died or emigrated within 30 days of discharge. We were forced to use the condition of 30-day survival to ensure that patients actually had the time to initiate secondary therapy. This shortcoming is discussed further in section 5.
A total of 49,640 patients (mean age, 66.3 years; standard deviation, 12.8 years; 35% females; 83% myocardial infarction) were included. Forty-six percent had received an early invasive treatment strategy. The number of deaths during follow-up (median: 3.6 years) was 10,847 (21.9%). Concomitant use of all 4 drugs (aspirin, a P2Y12 receptor inhibitor, statin, and β-blocker) after discharge was observed in 56% of patients (68 vs. 45% in the early and conservative invasive groups, respectively). Receiving an early invasive treatment strategy was associated with a lower incidence rate of all-cause death (3.1 vs. 8.1 deaths per 100 person-years; adjusted hazard ratio [HR], 0.71; 95% confidence interval [CI], 0.67 to 0.74; p<0.001) compared to a conservative approach.
This case will be used throughout this tutorial. It should be stressed that the analyses presented in this tutorial are meant as a pedagogical tool for explaining mediation analyses. Accordingly, a full discussion of the medical implications as well as some case-specific limitations are not included. The interested reader is referred to Hansen et al. [5] for an in-depth discussion.
The first step of any mediation analysis is to describe pre-existing beliefs about the causal structure in which the mediation analysis is to be conducted. Directed acyclic graphs (DAGs) are the method of choice for doing so. An example is given in Figure 1, where the assumed causal structure of the illustrative case is presented. An arrow in a DAG implies that we believe a possible causal connection exists between the 2 variables in question. A causal connection, such as AB, can loosely be interpreted as: “if we actively change the variable A, the distribution of B might change.” Note that the real assumption in the DAG AB is that an intervention on B will not change A. For a more detailed introduction to DAGs see Pearl [6]. The defining feature of a mediator is that it is positioned between the exposure and the outcome when following the directions of the arrows in the DAG. The DAG must also include all likely common causes of any pair constructed from an exposure, mediator, and outcome. In Figure 1, early intervention is assumed to affect the secondary treatment strategy (here defined as drug initiation within 30 days), which in turn affects mortality; this effect is called an indirect effect. There might also be other mechanisms, not involving the secondary treatment strategy, through which early intervention can affect mortality. These are subsumed within the direct arrow from early invasive strategy to death, which therefore can be thought of as mediation through all mediators except for the secondary treatment strategy; this is called a direct effect.
From an intuitive point of view, mediation analysis boils down to describing what would happen if a) the indirect pathway was the only causal pathway between exposure and outcome and b) the indirect pathway could be deactivated completely. However, this intuition is not sufficient to mathematically define the corresponding parameters to be estimated. We therefore introduce so-called counterfactual variables [6]. Building on the variables defined in Figure 1, these are:
  • Yi (a, m) is the outcome achieved for person i if, possibly contrary to fact, exposure had been set to a and mediator to m.

  • Mi (a) is the mediator achieved for person i if, possibly contrary to fact, exposure had been set to a.

The subscript i will be omitted when referring to a randomly picked person. The counterfactual variable Yi (1, m) corresponds to the death time observed in a double-intervention randomized trial where early intervention had been used and secondary medication set to m. Likewise, the counterfactual variable Mi (1) is the secondary medication observed in a single-intervention randomized trial where early intervention had been used. One can combine the two counterfactuals, yielding so-called nested counterfactuals defined as Y (a, M (a*)). When a=a*, the nested counterfactual simply corresponds to the observations one would observe if early intervention had been set to a. In the mediation analysis literature, the effect one would observe in a simple randomized trial is referred to as the total effect of treatment (i.e., early intervention), and it is defined as a comparison of the distribution of Y (1, M (1)) with that of Y (0, M (0)). The comparison could be done as a comparison of average values, but with a survival outcome, it would be more common to compare the 2 arms of the trial using a Cox model, leading to a causal HR quantifying the effect of treatment. The books by Pearl [6] and Hernán & Robins [7] provide a thorough introduction to why 1 arm of a randomized trial can be used to estimate the distribution of the counterfactual variable Y (1, M (1)), which is a quantity defined for the whole population, not only the people in the A=1 arm.
Realizing that the traditional 2-arm randomized controlled trial can be viewed as a double-intervention trial where, for instance, in the a=0 arm, treatment is set to 0 and the mediator to the value it would naturally take for that person when treatment is set to 0, leads to the following definition of the so-called natural direct and indirect effects. For ease of presentation alone, we compared the counterfactuals using their average values, but other scales such as odds ratios (ORs) could equally well have been used.
 Total effect of treatment:
  =E[Y (1, M (1))]−E[Y (0, M (0))]
 =(E[Y (1, M (1))]−E[Y (1, M (0))])+(E[Y (1, M (0))]−E[Y (0, M (0))])
 =natural indirect effect + natural direct effect
Written in words, the natural indirect effect is the effect you see by changing the mediator, as if you had changed the treatment without actually changing the treatment itself. Likewise, the natural direct effect is the effect you see by changing the treatment, but keeping the mediator fixed at whatever level it would be had you not changed the treatment. Thus, by introducing the nested counterfactual E[Y (a, M (a*))] for aa* we can give a precise mathematical definition of mediation. This definition was originally introduced by Pearl [8] and much work has since been published on identification, estimation, and applications, culminating in the recent book by Vanderweele [9], where a comprehensive list of references can be found. As the definition of natural direct and indirect effects at its core builds on comparing distributions of nested counterfactuals, these effects can just as easily be expressed on other scales than the averages. For a survival outcome, it would, for instance, be more common to decompose the HR as follows:
HR of a= 1 vs. a= 0=hazard for Y (1, M (1))hazard for Y (0, M (0))= hazard for Y (1, M (1))hazard for Y (1, M (0))×hazard for Y (1, M (0))hazard for Y (0, M (0))=natural indirect HR× natural direct HR
From the derivations above, it is apparent that the key to employing natural direct and indirect effects is to identify and estimate the distribution, or aspects of the distribution, of the nested counterfactuals Y (a, M (a*)) for potentially different a and a*. As in the DAG in Figure 1, we will allow non- randomized study settings as well. The following assumptions are sufficient to identify natural direct and indirect effects on any scale from independent observations of the triplet (C, A, Y) [9].
No uncontrolled confounding
Assume that the variables collected in C are sufficient for controlling confounding for a) the exposure-outcome relationship, b) the exposure-mediator relationship, and c) the mediator-outcome relationship (in [c], A is included in the set of control variables). Mathematically, the conditions are:
(1a)
Y (a, M (a))  A | C
(1b)
M (a)  A | C
(1c)
Y (a, m)  M | (A, C)
Positivity
Assume that for any values of confounders, all exposure values have non-zero probability and likewise that for any values of confounders and exposure, all mediator values have non-zero probability. Mathematically, the conditions are:
(2a)
P (A=a|C=c)> 0 for all a, c
(2b)
P (M=m|C=c, A=a)> 0 for all a, c, m
and equivalent using densities when A or M are continuous.
Consistency
Assume that the nested counterfactuals will actually take the observed values when the treatment and mediator are actively set to the values they would naturally have had in the absence of an intervention. Mathematically, the condition is:
(3a)
P (Y (A, M)=Y)= 1 and P (M (A)=M)= 1
Identification of natural effects
Assume that the counterfactual out come, Y (a, m) is independent of the counterfactual mediator, M (a*) when ever a and a* are different. Mathematically, the condition is:
(4a)
Y (a, m)  M (a*) | C for any m and aa* 
While these assumptions are structural and therefore not possible to verify using observed data alone, the identification assumption (4a) is by far the most difficult to comprehend. This assumption imposes independence between 2 distinct counterfactual worlds (because a and a* are assumed to be different). From an applied perspective, assumption (4a) can be replaced by assuming that there are no confounders of the mediator-outcome relationship that are themselves affected by exposure. Or, perhaps more practically, one can assume that the indirect and direct effects are created by distinct and causally unrelated mechanisms.
To see why these conditions suffice, we will next derive an explicit formula for E[g (Y (a, M (a*))]. The arbitrary measurable function g : R → R is included to demonstrate that it is the full distribution of the nested counterfactual that we have identified, not only the mean. For ease of exposition only, we will assume that C and M are discrete, with state space C and M, respectively.
E[g(Y (a, M(a )))]= cC E[g(Y (a, M(a ))) | C = c]P(C = c)=cC,mME[g(Y (a, M(a ))) | M(a ) = m, C = c]P(M(a ) = m | C = c)P(C = c)=cC,mME[g(Y (a, m)) | M(a ) = m, C = c]P(M(a ) = m | C = c)P(C = c)=icC,mME[g(Y (a, m)) | C = c]P(M(a ) = m | C = c)P(C = c)=iicC,mME[g(Y (a, m)) | A = a, M = m, C = c]P(M(a ) = m | A = a , C = c)P(C = c)=iiicC,mME[g(Y ) | A = a, M = m, C = c]P(M = m | A = a , C = c)P(C = c)
where equality i follows from (4a), equality ii from (1a-1c), and equality iii from (3a). The final expression only depends on the observed data and can therefore be estimated from the observed data. It appears as if the positivity assumption is not needed; however, it is precisely the positivity assumption that guarantees that, in large samples, all quantities in the final expression can be non-parametrically estimated. If one is only interested in a given function g and contrasts such as E[g (Y (1, M (1)))] − E[g (Y (0, M (1)))], then the identification assumption can be reduced to certain no-interaction assumptions.
When g is reduced to the identity function, the formula is known as Pearl’s mediation formula or just the mediation formula. While the mediation formula in principle allows non-parametric estimation of any mediation analysis, it can rarely be applied directly, as one will suffer the curse of dimensionality when trying to estimate E [g(Y) | A=a, M=m, C=c] in all strata. As in all other branches of statistics, the curse of dimensionality is countered by introducing parametric modelling assumptions. This will be the theme of the next section.
Several suggestions have been made for operationalizing the estimation of natural direct and indirect effects, such as the SPSS/SAS macros developed by Valeri & Vanderweele [10] and the R packages mediation and medflex. The topic of this section will be the class of natural effect models (NEMs) originally introduced by Lange et al. [11] and Vansterlandt et al. [12] and implemented in the R package medflex [2].
The idea underlying NEM is to phrase a mediation analysis as a multiple regression problem, thereby a) parameterizing the quantities of interest, b) allowing the choice of outcome model to follow the convention for that type of outcome (i.e., a Cox model for survival outcomes), and c) harvesting the extensive existing software implementing regression type models. In particular, a NEM is a regression model for the nested counterfactual. Expressed as a generalized linear model (GLM), it becomes:
g(E[Y (a, M (a*))])=α0 + α1a + α2a*
If, for instance, g is the logit function, then α1 would be the marginal natural direct effect log-OR. NEMs can also be formulated conditionally on measured covariates and do not have to be only simple additive and linear effect models:
g(E[Y (a, M (a*)) | C=c])=W (a, a*, c)
where W is a deterministic function. The class of NEM models, along with the corresponding estimation techniques [11,12], also applies to survival models, such as the Cox model. Estimation algorithms for NEMs are all derived under the assumptions listed in the preceding section, and build on the trick that first we duplicate the original data set, then we create an artificial exposure, A*, which takes on different values in the 2 replications of each observation. Finally, we use an auxiliary model to link the artificial observations (i.e., those where A ≠ A*) to the mediators, which is either done through weighting or through imputation. Once this is done, the NEM can be estimated by simply applying standard software applied to the extended data set and using both A and A*, possibly along with C, as the model specification. This entire procedure has been automated and implemented in the R package medflex for any GLM type outcome model. In the following, we describe in detail how to estimate a NEM for a survival outcome with a multidimensional mediator. This is the situation we have in the illustrative case.
  • - Using the original data alone, fit a parametric survival model to the outcome conditioned on confounders, exposure, and mediator. This could, for instance, be a Weibull based parametric time-to-event model.

  • - Construct a new data set by repeating each observation in the original data set twice and including an additional variable a*, which is equal to the original exposure for the first replication and equal to the opposite of the actual exposure for the second replication. In addition, add an identification variable to indicate which data rows originate from the same subject.

  • - Use the predict functionality, possibly along with the Weibull distribution function, to impute possible survival times for the rows where A ≠ A*. In the imputations, the value of the exposure is set to a*, while mediators and confounders are set to their observed values; that is, impute values for the survival times Yi(a*, M (a)).

  • - Fit a Cox model to the extended data set including A, A*, and C, but not the mediator. The coefficient of A will be the natural indirect log-HR and the coefficient of a* will be the natural direct log-HR.

  • - Repeat steps 3 and 4 ten times and combine the obtained log- HRs as with ordinary multiple imputation; that is, take the average of the log-HR estimates.

  • - CIs for the natural effect estimates, as well as derived quantities such as mediated proportions, can be obtained by bootstrapping, which involves repeating steps 1-5 a total of 1,000 times, each time creating a new data set by random sampling with replacement from the original data set.

In our illustrative case, the outcome of interest is death from any cause during follow-up. Censoring was almost exclusively administrative, as emigration is rare in this population. The underlying time scale is in years, starting at 30 days after hospital discharge. The mediator is the 4-dimensional variable indicating whether each of the 4 recommended secondary treatments was followed.
As exposure to aspirin, P2Y12 inhibitors, statins, and β-blockers could not reasonably be said to constitute distinct causal pathways, but were highly interdependent, mediation was only assessed through the combined 4-dimensional mediator. Accordingly, the counterfactual mediator was M (a) є {0, 1}4, where a=0 indicates conservative invasive treatment and a=1 indicates early invasive treatment. The nested counterfactual was death time in years, starting 30 days after hospital discharge. To accommodate censoring, the nested counterfactual outcome technically had 2 dimensions, namely, an event time and an event indicator Y (a, M (a*))=((T, δ) (a, M (a*)). For ease of exposition, we will only refer to the underlying event time and event indicator when required by context. Our effect measure of interest was natural direct and indirect HRs decomposing the total effect, which had a HR of 0.71 (95% CI, 0.67 to 0.74). Accordingly, the final natural effects model (i.e., the one fitted in step 4 of our suggested approach) should be a Cox proportional hazard model.
Table 1 presents simple descriptive statistics of the data. Because of the very large sample size, all associations between the treatment strategy and confounders, as well as mediators, were highly significant.
Figure 2 presents Kaplan-Meier curves for the 2 treatment strategy groups. It also presents the curves obtained by fitting a parametric survival model with a Weibull error distribution using only the treatment groups as covariates. It clearly shows the large differences in raw survival between the groups, and, more importantly, the figure also demonstrates that the simple 3-parameter model does a good job of capturing the distributions (Figure 2).
Step 1: For the actual mediation analysis, we fit a parametric survival model with a Weibull error distribution to the survival times using the treatment group, mediators, and a long list of potential confounders as explanatory variables (Appendix 1). The Appendix 2 presents both the employed R code and the full table of estimated parameters. Note that for technical reasons relating to R, it is better to use a copy of the exposure variable when fitting this model.
Step 2: The data set can be duplicated and the auxiliary variable inserted by copying the original data set twice. In both copies a new variable is created (called, say, exposure Star). In the first copy, the new variable is set to the values of the actual received treatment (exposure Star=exposure) while in the second copy, the new variable is set to the opposite value of the actual received treatment. Finally, the two copies are appended producing a single data set with twice as many rows as in the original data set. See the R code in Appendix 2 for coding advice.
Step 3: We now set the temporary exposure variable, used when fitting the imputation model in step 1, to the values in the just created exposure variable (i.e., exposure Star). As the employed survival model is parametric, we can randomly draw survival times conditional on observed mediators and the just-defined temporary exposure variable. This corresponds to randomly drawing the nested counterfactuals variables Ya*,a, where a* and a are potentially different. To avoid extrapolating the imputation model outside what is supported by the data, any imputed survival time longer than the decided maximum follow-up (7 years in this case) will be artificially censored at the time of maximal follow-up. For the R implementation, see lines 21 to 32 of the Appendix 2.
Step 4-5: For each draw of the imputed nested counterfactuals, we fit a Cox model including the exposure that was actually received (the coefficient of this variable will estimate the natural indirect log-HR), the created exposure variable (the coefficient of this variable will estimate the natural direct log-HR), and all confounders. Note that the mediators are not included in this model. This is repeated 10 times, and the 10 resulting model fits are combined using standard formulas for multiple imputation. The resulting estimates are reported in Appendix 1. The fitted model is a natural effects Cox model.
Step 6: Finally, CIs are established by 5,000 bootstrap repetitions of steps 1-5. From the bootstrapped replications, we also directly obtain CIs for derived quantities, as the total effect (the sum of natural direct and indirect log-HRs) and the mediated proportion (natural indirect log-HR divided by total effect log-HR). The results are reported in Table 2.
From the Table 2, it is observed that after controlling for confounders, the use of early invasive treatment was associated with a reduction in 1-year mortality of 30% (OR, 0.70). The effect of early invasive treatment has 2 components: an indirect effect through secondary preventive medication, reducing risk by 10% (OR, 0.90), and the effect through all other pathways, reduces risk by a further 23%. An equivalent statement is that between a quarter and a third of the beneficial effect of early invasive treatment was achieved through the use of the 4 discharge medications. Arguably, this part of the survival gain could be achieved without adopting a full early invasive strategy, but instead by increasing the use of the 4 discharge medications to the levels seen in patients who underwent an early invasive treatment.
R package medflex to avoid own coding
As demonstrated above and in the Appendix 2, the natural effect Cox models requires some independent coding from the researcher. This is in contrast to most other types of outcomes (continuous, binary, counts, etc.) where estimating the natural effect models has been completely automated in the R package medflex [2]. At present, the medflex package does not support survival models; however, such functionality is expected to be introduced in upcoming versions of the package. To illustrate the use of the package, we will reanalyse the illustrative case using 1-year survival, which is essentially fully observed in the data and can therefore be analysed using a natural effects logistic model, which is fully supported by the medflex package. Across the sample, we have 3,610 deaths within 1 year, corresponding to 7.3%. The medflex package will carry out the same steps as described above, but in an automated manner. Accordingly, the first step is to specify an imputation model and feed this model to the function neModel. Here, we also specify the number of mediators. The code is included below, where BinaryMort is the outcome and dhrkag3 the exposure. The mediators are asatreat30, adptreat30, statintreat30, and betatreat30 and the number of mediators is specified using the nMed argument. All other variables are confounders; see Appendix 1 for definitions of the variables.
  • fitAux ← glm (BinaryMort ~ dhrkag3 + asatreat30 + adptreat30

  • + statintreat30 + betatreat30 + i_alder + factor (sex)

  • + factor (indkgrp) + factor (uddankat) + boralene + factor (fi_aar)

  • + mi + card + cochf + puled + shock + cervas + mal + diabet

  • + arf + crf + anemi + pneumoni + sepsis + klap + bleed

  • + Antihyp_12mb + Lipidlow_12mb + ASA_12mb + VitKant_12mb

  • + Diureti_loop_12mb + COPD_12mb + tidl_reva, data=workData, family=“binomial”)

  • extendedData ← neImpute(fitAux, nMed=4)

The last step is then to specify the natural effects model within the neModel function and extract the estimates for natural direct and indirect effects. The code is presented below.
  • fit NEM_binaryOutcome ← neModel (BinaryMort ~ dhrkag30 + dhrkag31

  • + i_alder + factor (sex) + factor (indkgrp) + factor (uddankat)

  • + boralene + factor (fi_aar) + mi + card + cochf + puled + shock

  • + cervas + mal + diabet + arf + crf + anemi + pneumoni + sepsis

  • + klap + bleed + Antihyp_12mb + Lipidlow_12mb + ASA_12mb

  • + VitKant_12mb + Diureti_loop_12mb + COPD_12mb + tidl_reva, expData=extendedData, family=“binomial”, se=“robust”

Summary (neEffdecomp (fitNEM_binaryOutcome))
The neImpute function creates the 2 new auxiliary exposure variables dhrkag30 and dhrkag31, which correspond to the natural direct and indirect effects, respectively. As the natural effects model is in this case a logistic regression, the estimates are presented as ORs in the Table 3.
From the Table 3, it is observed that the mediated proportion is similar to what was found in the Cox-based analysis, but with wider CIs. As we are using fewer events in the analysis, the wider CIs are to be expected. The effect estimates are similar to the Cox-based analysis, but numerically smaller; however, as one is a HR and the other is an OR, they cannot be directly compared.
As outlined in the preceding sections, mediation analysis with a single well-defined mediator (possibly multi-dimensional) and associated simple causal structure has by now been very well researched. This includes theoretical considerations and software implementations. For the applied researcher, a review of existing software solutions written by Starkopf et al. [13] is under review and available as a working paper upon request. On the purely applied side, we still need to see more applications, mainly to address subject matter problems, but also to establish common best practices for conducting mediation analyses.
This by no means implies that there are no unsolved methodological questions within mediation analysis. We see some of the most pressing problems as:
First, how to handle measurement error for the mediator. Currently, the best suggestion is to conduct sensitivity analyses assessing the potential impact of such measurement errors. This is of course good, but it would be more fruitful to have methods that could handle measurement errors directly. Mplus has capabilities in this direction [14], but they come at the cost of numerous parametric assumptions, and worse, a reduced causal interpretation, because effects are expressed on a latent and somewhat arbitrarily defined scale. Second, further methods to handle causally ordered mediators and/or mediators measured repeatedly over time. Important work in this regard was recently published [15]. Moreover, in the context of a survival outcome, the problem is further complicated by the fact that death has a truncation effect on the mediator process [16].
Third, existing software for mediation analysis should be extended to make it easier to conduct sensitivity analyses.
It is our hope that this tutorial has shown the potential of mediation analyses in discovering the causal mechanisms underlying a given cause-and-effect relation, and has demonstrated the relative ease with which mediation analyses can be conducted using standard software.

The authors have no conflicts of interest to declare for this study.

Figure 1.
Generic directed acyclic graph for mediation analysis (A) and for the illustrative example (B). C, confounder; A, exposure; M, mediator; Y, outcome.
epih-39-e2017035f1.gif
Figure 2.
Kaplan-Meier curves (full line) along with survival curves (finely dashed lines) from a fitted parametric model with a Weibull error distribution. The lower curves are for the conservative strategy group, while the upper is for the early invasive strategy.
epih-39-e2017035f2.gif
Table 1.
Descriptive statistics
Conservative Early invasive p-value
n 26,858 22,782
Mean age (yr) 69.0 63.0 <0.001
Male (%) 59.4 70.7 <0.001
Table 2.
Summary of mediation analysis
HR 95% CI
Lower limit Upper limit
Effect1
 Natural indirect 0.90 0.88 0.92
 Natural direct 0.77 0.76 0.79
 Total 0.70 0.69 0.70
Mediated proportion 0.30 0.25 0.34

HR, hazard ratio; CI, confidence interval.

1 The effects are HRs for all-cause mortality except for the mediated proportion. These results are based on a natural effects Cox model conditional on all recorded baseline confounders.

Table 3.
Natural direct and indirect ORs when only looking at 1-year survival
OR 95% CI
Lower limit Upper limit
Effect
 Natural indirect 0.84 0.78 0.89
 Natural direct 0.66 0.62 0.70
 Total 0.55 0.51 0.60
Mediated proportion 0.30 0.22 0.38

OR, odds ratio; CI, confidence interval.

  • 1. Baron RM, Kenny DA. The moderator-mediator variable distinction in social psychological research: conceptual, strategic, and statistical considerations. J Pers Soc Psychol 1986;51:1173-1182.ArticlePubMedPDF
  • 2. Steen J, Loeys T, Moerkerke B, Vansteelandt S. Medflex: an R package for flexible mediation analysis using natural effect models. J Stat Softw 2017;doi: https://doi.org/10.18637/jss.v076.i11.Article
  • 3. Hansen KW, Sorensen R, Madsen M, Madsen JK, Jensen JS, von Kappelgaard LM, et al. Effectiveness of an early versus a conservative invasive treatment strategy in acute coronary syndromes: a nationwide cohort study. Ann Intern Med 2015;163:737-746.ArticlePubMed
  • 4. Roffi M, Patrono C, Collet JP, Mueller C, Valgimigli M, Andreotti F, et al. 2015 ESC guidelines for the management of acute coronary syndromes in patients presenting without persistent ST-segment elevation: Task Force for the Management of Acute Coronary Syndromes in Patients Presenting without Persistent ST-Segment Elevation of the European Society of Cardiology (ESC). Eur Heart J 2016;37:267-315.ArticlePubMedPDF
  • 5. Hansen KW, Sorensen R, Madsen M, Madsen JK, Jensen JS, von Kappelgaard LM, et al. Effectiveness of an early versus a conservative invasive treatment strategy in acute coronary syndromes: a nationwide cohort study. Ann Intern Med 2015;163:737-746.ArticlePubMed
  • 6. Pearl J. Causality: models, reasoning, and inference. 2nd ed. Cambride: Cambridge University Press; 2009.
  • 7. Hernán MA, Robins JM. Causal inference. London: Chapman & Hall; 2012.
  • 8. Pearl J. Direct and indirect effects. In: Jack Breese J, Koller D, eds. American Association for Artificial Intelligence. Uncertainty in artificial intelligence. Proceedings of the Seventeenth Conference. 2001 Aug 2-5. Seattle, USA. San Francisco: Morgan Kaufmann Pub. 2001. 411-420.
  • 9. Vanderweele TJ. Explanation in-causal inference: methods for mediation and interaction. New York: Oxford University Press; 2015.
  • 10. Valeri L, Vanderweele TJ. Mediation analysis allowing for exposure-mediator interactions and causal interpretation: theoretical assumptions and implementation with SAS and SPSS macros. Psychol Methods 2013;18:137-150.ArticlePubMedPMCPDF
  • 11. Lange T, Vansteelandt S, Bekaert M. A simple unified approach for estimating natural direct and indirect effects. Am J Epidemiol 2012;176:190-195.ArticlePubMedPDF
  • 12. Vansteelandt S, Bekaert M, Lange T. Imputation strategies for the estimation of natural direct and indirect effects. Epidemiol Methods 2012;1:131-158.ArticlePDF
  • 13. Starkopf L, Andersen MP, Gerds TA, Torp-Pedersen C, Lange T. Comparison of five software solutions to mediation analysis; 2017 [cited 2017 Sep 14]. Available from: https://ifsv.sund.ku.dk/biostat/annualreport/images/0/0a/Research_Report_17-01.pdf.
  • 14. Muthen B, Muthen LK, Asparouhov T. Regression and mediation analysis using Mplus. Los Angeles: Muthen & Muthen; 2016.
  • 15. Steen J, Loeys T, Moerkerke B, Vansteelandt S. Flexible mediation analysis with multiple mediators. Am J Epidemiol 2017;186:184-193.ArticlePubMedPDF
  • 16. Strohmaier S, Haase N, Wetterslev J, Lange T. A simple to implement algorithm for natural direct and indirect effects in survival studies with a repeatedly measured mediator; 2017 [cited 2017 Aug 31]. Available from: http://www.stat-center.pku.edu.cn/Stat/Index/research_show/id/227.
Appendix 1.
Full model fits
When estimating the natural effect Cox model in the illustrative case, a Weibull parametric survival model was used as an intermediate step in order be able to impute the nested counterfactual. The model contains exposure, mediators, and all considered confounders.
After imputation, a natural effects Cox model can be estimated by fitting a Cox model to the extended data set. The table below presents estimates for a single imputation.
Table S1a.
Model t for Weibul based parametric survival model
Variable (in R notation) Estimate SE p-value
(Intercept) 13.302 0.101 0.000
dhrkag3TEMP 0.250 0.024 0.000
asa_treat30 0.080 0.015 0.000
adp_treat30 0.095 0.011 0.000
statin_treat30 0.204 0.013 0.000
beta_treat30 0.098 0.012 0.000
i_alder -0.064 0.001 0.000
factor(sex)2 0.232 0.021 0.000
factor(indkgrp)2 0.159 0.028 0.000
factor(indkgrp)3 0.316 0.038 0.000
factor(uddankat)2 0.030 0.022 0.165
factor(uddankat)3 0.146 0.034 0.000
boralene -0.155 0.021 0.000
factor(fi_aar)2006 0.003 0.029 0.918
factor(fi_aar)2007 -0.027 0.032 0.399
factor(fi_aar)2008 -0.146 0.038 0.000
factor(fi_aar)2009 -0.209 0.042 0.000
factor(fi_aar)2010 -0.173 0.048 0.000
factor(fi_aar)2011 -0.280 0.055 0.000
mi -0.598 0.032 0.000
card -0.026 0.027 0.339
cochf -0.283 0.026 0.000
puled -0.234 0.077 0.003
shock -0.107 0.170 0.530
cervas -0.341 0.034 0.000
mal -0.985 0.039 0.000
diabet -0.388 0.036 0.000
arf -0.405 0.066 0.000
crf -0.395 0.051 0.000
anemi -0.285 0.040 0.000
pneumoni -0.380 0.028 0.000
sepsis -0.155 0.074 0.036
klap -0.260 0.037 0.000
bleed -0.102 0.048 0.034
Antihyp_12mb -0.021 0.022 0.352
Lipidlow_12mb 0.093 0.024 0.000
ASA_12mb -0.101 0.022 0.000
VitKant_12mb -0.048 0.037 0.188
Diureti_loop_12mb -0.442 0.023 0.000
COPD_12mb -0.362 0.023 0.000
tidl_reva 0.278 0.064 0.000
Log(scale) -0.022 0.009 0.012
Table S1b.
Full model t for natural eect model based on a single imputation
Description R variable names log-HR SE p-value
Indirect dhrkag3 -0.10 0.01 0.00
Direct dhrkag3STAR -0.26 0.01 0.00
Age i_alder 0.07 0.00 0.00
Gender (female) sex2 -0.21 0.01 0.00
Income (middle) factor(indkgrp)2 -0.16 0.02 0.00
Income (high) factor(indkgrp)3 -0.32 0.03 0.00
Education (middle) factor(uddankat)2 -0.05 0.01 0.00
Education (high) factor(uddankat)3 -0.12 0.02 0.00
Live alone boralene 0.17 0.01 0.00
Year (2006) factor(fi_aar)2006 -0.03 0.02 0.19
Year (2007) factor(fi_aar)2007 0.00 0.02 0.92
Year (2008) factor(fi_aar)2008 0.14 0.03 0.00
Year (2009) factor(fi_aar)2009 0.20 0.03 0.00
Year (2010) factor(fi_aar)2010 0.16 0.03 0.00
Year (2011) factor(fi_aar)2011 0.25 0.03 0.00
Myocardial infarction mi 0.49 0.02 0.00
Cardiac arrhythmia card 0.05 0.02 0.01
chronic obstructive pulmonary disease cochf 0.30 0.02 0.00
Pulmonary oedema puled 0.30 0.06 0.00
Cardiogenic shock shock 0.22 0.11 0.05
Cerebrovascular disease cervas 0.37 0.02 0.00
Diabetes with complications diabet 0.36 0.03 0.00
acute renal failure arf 0.44 0.05 0.00
Chronic renal failure crf 0.44 0.04 0.00
Anaemia anemi 0.34 0.03 0.00
Pneumonia pneumoni 0.46 0.02 0.00
Sepsis sepsis 0.20 0.05 0.00
Valvular heart disease klap 0.26 0.03 0.00
Prior in-hospital bleeding bleed 0.17 0.03 0.00
Use of antihyp. medication last 12M Antihyp_12mb 0.06 0.01 0.00
Use of lipid-lowering drugs last 12M Lipidlow_12mb -0.04 0.02 0.02
Asprin ASA_12mb 0.16 0.01 0.00
Use of vitamin K antagonists last 12M VitKant_12mb 0.09 0.02 0.00
Use of glucose-lowering drugs lst 12M?? Diureti_loop_12mb 0.50 0.02 0.00
Use of loop diuretics or COPD last 12M COPD_12mb 0.39 0.02 0.00
Prior revascularization tidl_reva -0.24 0.04 0.00
Appendix 2.
Employed R code
The full R code used for estimating the natural effect Cox model is presented below.
epih-39-e2017035-app2.pdf

Figure & Data

References

    Citations

    Citations to this article as recorded by  
    • Temporary relationship between sleep duration and depression and its impact on future risk of cardiovascular disease
      Xue Tian, Xue Xia, Yijun Zhang, Qin Xu, Yanxia Luo, Anxin Wang
      Journal of Affective Disorders.2024; 350: 559.     CrossRef
    • Potential Mediators for Treatment Effects of Novel Diabetes Medications on Cardiovascular and Renal Outcomes: A Meta‐Regression Analysis
      José M. Rodriguez‐Valadez, Malak Tahsin, Umesh Masharani, Meyeon Park, M. G. Myriam Hunink, Joseph Yeboah, Lihua Li, Ellerie Weber, Asem Berkalieva, Luuk Avezaat, Wendy Max, Kirsten E. Fleischmann, Bart S. Ferket
      Journal of the American Heart Association.2024;[Epub]     CrossRef
    • Metabolic Syndrome Is Associated With an Increased Risk of Rheumatoid Arthritis: A Prospective Cohort Study Including 369,065 Participants
      Peiyang Luo, WanLi Xu, Ding Ye, Weiwei Chen, Jiacheng Ying, Bin Liu, Jiayu Li, Xiaohui Sun, Zhixing He, Chengping Wen, Yingying Mao
      The Journal of Rheumatology.2024; 51(4): 360.     CrossRef
    • Conditional Separable Effects
      Mats J. Stensrud, James M. Robins, Aaron Sarvet, Eric J. Tchetgen Tchetgen, Jessica G. Young
      Journal of the American Statistical Association.2023; 118(544): 2671.     CrossRef
    • Mediation Analysis in Manual Therapy Research
      Thiago da Silva Rocha Paz, Pedro Teixeira Vidinha Rodrigues, Bruno Moreira Silva, Arthur de Sá Ferreira, Leandro Alberto Calazans Nogueira
      Journal of Chiropractic Medicine.2023; 22(1): 35.     CrossRef
    • An exact regression‐based approach for the estimation of natural direct and indirect effects with a binary outcome and a continuous mediator
      Mariia Samoilenko, Geneviève Lefebvre
      Statistics in Medicine.2023; 42(3): 353.     CrossRef
    • Cardiovascular disease, muscle function, and long-term falls risk: The Perth Longitudinal Study of Ageing Women
      Abadi K. Gebre, Marc Sim, Jack Dalla Via, Alexander J. Rodríguez, Kun Zhu, John T. Schousboe, Jonathan M. Hodgson, Catherine P. Bondonno, Richard L. Prince, Joshua R. Lewis
      Archives of Gerontology and Geriatrics.2023; 107: 104911.     CrossRef
    • Is the socioeconomic inequality in stroke prognosis changing over time and does quality of care play a role?
      Vibe Bolvig Hyldgård, Rikke Søgaard, Jan Brink Valentin, Theis Lange, Dorte Damgaard, Søren Paaske Johnsen
      European Stroke Journal.2023; 8(1): 351.     CrossRef
    • Neighbourhood‐related socioeconomic perinatal health inequalities: An illustration of the mediational g‐formula and considerations for the big data context
      Lizbeth Burgos Ochoa, Lindsey van der Meer, Adja J. M. Waelput, Jasper V. Been, Loes C. M. Bertens
      Paediatric and Perinatal Epidemiology.2023; 37(4): 341.     CrossRef
    • Mediation analysis of circadian preferences and other behavioral and sociodemographic factors associated with subjective well-being in undergraduate students
      Julia Ribeiro da Silva Vallim, Gabriela Sant’Ana Lima, Heloisa Mayumi Suyama Tsuji, Beatriz Duarte Palma Xylaras, Sergio Tufik, Marcelo Demarzo, Vânia D’Almeida
      Chronobiology International.2023; 40(4): 376.     CrossRef
    • Temporal relationship between hyperuricemia and hypertension and its impact on future risk of cardiovascular disease
      Xue Tian, Shuohua Chen, Penglian Wang, Qin Xu, Yijun Zhang, Xiaoli Zhang, Shouling Wu, Yanxia Luo, Anxin Wang
      European Journal of Internal Medicine.2023; 111: 82.     CrossRef
    • Childhood adversity and risk of later labor market marginalization in young employees in Sweden
      Emma Björkenstam, Magnus Helgesson, Ellenor Mittendorfer-Rutz
      European Journal of Public Health.2023;[Epub]     CrossRef
    • Adolescent BMI trajectory and associations with adult metabolic syndrome and offspring obesity
      Darren M. Frank, Patrick T. Bradshaw, Mahasin Mujahid, Elissa Epel, Barbara A. Lararia
      Obesity.2023; 31(7): 1924.     CrossRef
    • TCA cycle metabolites associated with adverse outcomes after acute coronary syndrome: mediating effect of renal function
      Raul Sanchez-Gimenez, Óscar M. Peiró, Gil Bonet, Anna Carrasquer, George A. Fragkiadakis, Mònica Bulló, Christopher Papandreou, Alfredo Bardaji
      Frontiers in Cardiovascular Medicine.2023;[Epub]     CrossRef
    • Mediation analysis of the testosterone treatment effect to prevent type 2 diabetes in the Testosterone for Prevention of Type 2 Diabetes Mellitus trial
      Kristy P Robledo, Ian C Marschner, David J Handelsman, Karen Bracken, Bronwyn G A Stuckey, Bu B Yeap, Warrick Inder, Mathis Grossmann, David Jesudason, Carolyn A Allan, Gary Wittert
      European Journal of Endocrinology.2023; 189(1): 50.     CrossRef
    • The role of negative affects as mediators in the relationship between stress and mental health in Ecuadorian adolescents
      Rodrigo Moreta-Herrera, Doris Zumba-Tello, Jaisalmer de Frutos-Lucas, Stefanía Llerena-Freire, Alexandra Salinas-Palma, Almudena Trucharte-Martínez
      Health Psychology Report.2023;[Epub]     CrossRef
    • An Introduction to Causal Mediation Analysis With a Comparison of 2 R Packages
      Sangmin Byeon, Woojoo Lee
      Journal of Preventive Medicine and Public Health.2023; 56(4): 303.     CrossRef
    • The mediating role of negative symptoms in “secondary factors” determining social functioning in chronic schizophrenia
      Na Hu, Wei Li, Hu Deng, Jiaqi Song, Hanxue Yang, Jiabao Chai, Wenqian Huang, Hong Wang, Xuanzi Zhou, Pan Zhang, Sushuang He, Yonghua Cui, Tengteng Fan, Ying Li
      Frontiers in Psychiatry.2023;[Epub]     CrossRef
    • Social deprivation and kidney failure due to an undiagnosed nephropathy
      Hamza Sakhi, Mathilde Beaumier, Cécile Couchoud, Mathilde Prezelin-Reydit, Jennifer Radenac, Thierry Lobbedez, Denis Morin, Vincent Audard, Valérie Chatelet
      Nephrology Dialysis Transplantation.2023; 39(1): 133.     CrossRef
    • Arsenic Exposure and Epigenetic Aging: The Association with Cardiovascular Disease and All-Cause Mortality in the Strong Heart Study
      Enoch X. Jiang, Arce Domingo-Relloso, Ahlam Abuawad, Karin Haack, Maria Tellez-Plaza, M. Danielle Fallin, Jason G. Umans, Lyle G. Best, Ying Zhang, Allison Kupsco, Daniel W. Belsky, Shelley A. Cole, Ana Navas-Acien
      Environmental Health Perspectives.2023;[Epub]     CrossRef
    • Aortic Stiffness, Inflammation, and Incidence of Cardiovascular Events in Elderly Participants From the General Population
      Benjamin Nilsson Wadström, Margaretha Persson, Gunnar Engström, Peter M. Nilsson
      Angiology.2022; 73(1): 51.     CrossRef
    • Effect of social deprivation on peritoneal dialysis uptake: A mediation analysis with the data of the REIN registry
      Mathilde Beaumier, Eve Calvar, Ludivine Launay, Clémence Béchade, Antoine Lanot, Nicole Schauder, Fatouma Touré, Mathilde Lassalle, Cécile Couchoud, Valérie Châtelet, Thierry Lobbedez
      Peritoneal Dialysis International: Journal of the International Society for Peritoneal Dialysis.2022; 42(4): 361.     CrossRef
    • Are sex differences in the outcome of peritoneal dialysis explained by nurse assistance? A cohort study with data from the Registre de Dialyse Péritonéale de Langue Française
      Isabelle Danneville, Mathilde Beaumier, Valérie Chatelet, Annabel Boyer, Antoine Lanot, Clémence Bechade, Thierry Lobbedez
      Nephrology Dialysis Transplantation.2022; 37(8): 1520.     CrossRef
    • A Systematic Review of the Reporting Quality of Observational Studies That Use Mediation Analyses
      Rodrigo R. N. Rizzo, Aidan G. Cashin, Matthew K. Bagg, Sylvia M. Gustin, Hopin Lee, James H. McAuley
      Prevention Science.2022; 23(6): 1041.     CrossRef
    • Illustrations of serial mediation using PROCESS, Mplus and R
      Laura Lemardelet, Pier-Olivier Caron
      The Quantitative Methods for Psychology.2022; 18(1): 66.     CrossRef
    • Social Context and Rejection Expectations Modulate Neural and Behavioral Responses to Social Feedback
      Pauline Petereit, Sarah Jessen, Tatiana Goregliad Fjaellingsdal, Ulrike M. Krämer
      Journal of Cognitive Neuroscience.2022; 34(5): 823.     CrossRef
    • Translating questions to estimands in randomized clinical trials with intercurrent events
      Mats J. Stensrud, Oliver Dukes
      Statistics in Medicine.2022; 41(16): 3211.     CrossRef
    • Associations between Brain Microstructure and Phonological Processing Ability in Preschool Children
      Ying Zhou, Guangfei Li, Zeyu Song, Zhao Zhang, Huishi Huang, Hanjun Li, Xiaoying Tang
      Children.2022; 9(6): 782.     CrossRef
    • Impact of diet on ten-year absolute cardiovascular risk in a prospective cohort of 94 321 individuals: A tool for implementation of healthy diets
      Emilie W. Kjeldsen, Jesper Q. Thomassen, Katrine L. Rasmussen, Børge G. Nordestgaard, Anne Tybjærg-Hansen, Ruth Frikke-Schmidt
      The Lancet Regional Health - Europe.2022; 19: 100419.     CrossRef
    • MarZIC: A Marginal Mediation Model for Zero-Inflated Compositional Mediators with Applications to Microbiome Data
      Quran Wu, James O’Malley, Susmita Datta, Raad Z. Gharaibeh, Christian Jobin, Margaret R. Karagas, Modupe O. Coker, Anne G. Hoen, Brock C. Christensen, Juliette C. Madan, Zhigang Li
      Genes.2022; 13(6): 1049.     CrossRef
    • Racial disparities in epithelial ovarian cancer survival: An examination of contributing factors in the Ovarian Cancer in Women of African Ancestry consortium
      Holly R. Harris, Kristin A. Guertin, Tareq F. Camacho, Courtney E. Johnson, Anna H. Wu, Patricia G. Moorman, Evan Myers, Traci N. Bethea, Elisa V. Bandera, Charlotte E. Joslin, Heather M. Ochs‐Balcom, Lauren C. Peres, Will T. Rosenow, Veronica W. Setiawan
      International Journal of Cancer.2022; 151(8): 1228.     CrossRef
    • Cardiorespiratory fitness as a mediator in the relationship between lung function and blood pressure in adults
      F.R. Almeida, T.L.V.D.P. Ostolin, V.R. Almeida, B.B. Gonze, E.F. Sperandio, M.S.M.P. Simões, I. Godoy, S.E. Tanni, M. Romiti, R.L. Arantes, V.Z. Dourado
      Brazilian Journal of Medical and Biological Research.2022;[Epub]     CrossRef
    • The association between prenatal greenspace exposure and Autism spectrum disorder, and the potentially mediating role of air pollution reduction: A population-based birth cohort study
      Lief Pagalan, Tim F. Oberlander, Gillian E. Hanley, Laura C. Rosella, Celeste Bickford, Whitney Weikum, Nancy Lanphear, Bruce Lanphear, Michael Brauer, Matilda van den Bosch
      Environment International.2022; 167: 107445.     CrossRef
    • Mediator Effect of Cardiorespiratory Fitness on the Association between Physical Activity and Lung Function in Adults: Cross-Sectional Results from the Epimov Study
      Thatiane Lopes Valentim Di Paschoale Ostolin, Bárbara de Barros Gonze, Evandro Fornias Sperandio, Rodolfo Leite Arantes, Marcello Romiti, Victor Zuniga Dourado
      International Journal of Environmental Research and Public Health.2022; 19(15): 9377.     CrossRef
    • Is there sex disparity in vascular access at dialysis initiation in France? A mediation analysis using data from the Renal Epidemiology and Information Network registry
      Mathilde Beaumier, Maxence Ficheux, Cécile Couchoud, Mathilde Lassalle, Ludivine Launay, Cécile Courivaud, Aurélien Tiple, Thierry Lobbedez, Valérie Chatelet
      Clinical Kidney Journal.2022; 15(11): 2144.     CrossRef
    • Plasma trimethylamine-N-oxide, its precursors and risk of cardiovascular events in patients with acute coronary syndrome: Mediating effects of renal function
      Raul Sanchez-Gimenez, Óscar M. Peiró, Gil Bonet, Anna Carrasquer, Georgios A. Fragkiadakis, Mònica Bulló, Christopher Papandreou, Alfredo Bardaji
      Frontiers in Cardiovascular Medicine.2022;[Epub]     CrossRef
    • Mediators linking obesity to childhood asthma
      Yang‐Ching Chen, Yen‐Tsung Huang, Wen‐Harn Pan, Yungling L. Lee
      Pediatric Allergy and Immunology.2022;[Epub]     CrossRef
    • Pleiotropic modifiers of age-related diabetes and neonatal intestinal obstruction in cystic fibrosis
      Melis A. Aksit, Hua Ling, Rhonda G. Pace, Karen S. Raraigh, Frankline Onchiri, Anna V. Faino, Kymberleigh Pagel, Elizabeth Pugh, Adrienne M. Stilp, Quan Sun, Elizabeth E. Blue, Fred A. Wright, Yi-Hui Zhou, Michael J. Bamshad, Ronald L. Gibson, Michael R.
      The American Journal of Human Genetics.2022; 109(10): 1894.     CrossRef
    • Time of coronary revascularization: methodology of a mediation analysis study
      Boris Sobolev, Lisa Kuramoto
      CMAJ Open.2022; 10(4): E1052.     CrossRef
    • Insulin resistance mediates obesity-related risk of cardiovascular disease: a prospective cohort study
      Xue Tian, Shuohua Chen, Penglian Wang, Qin Xu, Yijun Zhang, Yanxia Luo, Shouling Wu, Anxin Wang
      Cardiovascular Diabetology.2022;[Epub]     CrossRef
    • Obstructive sleep apnea, chronic obstructive pulmonary disease and NAFLD: an individual participant data meta-analysis
      Ingrid Jullian-Desayes, Wojciech Trzepizur, Jérôme Boursier, Marie Joyeux-Faure, Sébastien Bailly, Meriem Benmerad, Marc Le Vaillant, Sandrine Jaffre, Thierry Pigeanne, Acya Bizieux-Thaminy, Marie-Pierre Humeau, Claire Alizon, François Goupil, Charlotte C
      Sleep Medicine.2021; 77: 357.     CrossRef
    • The -514C>T polymorphism in the LIPC gene modifies type 2 diabetes risk through modulation of HDL-cholesterol levels in Mexicans
      M. T. Guerra-García, H. Moreno-Macías, A. Ochoa-Guzmán, M. L. Ordoñez-Sánchez, R. Rodríguez-Guillen, P. Vázquez-Cárdenas, V. M. Ortíz-Ortega, M. Peimbert-Torres, C. A. Aguilar-Salinas, M. T. Tusié-Luna
      Journal of Endocrinological Investigation.2021; 44(3): 557.     CrossRef
    • Items for consideration in a reporting guideline for mediation analyses: a Delphi study
      Aidan G Cashin, James H McAuley, Sallie Lamb, Sally Hopewell, Steven J Kamper, Christopher M Williams, Nicholas Henschke, Hopin Lee
      BMJ Evidence-Based Medicine.2021; 26(3): 106.     CrossRef
    • Joint effects of alcohol use, smoking and body mass index as an explanation for the alcohol harm paradox: causal mediation analysis of eight cohort studies
      Sebastián Peña, Pia Mäkelä, Tiina Laatikainen, Tommi Härkänen, Satu Männistö, Markku Heliövaara, Seppo Koskinen
      Addiction.2021; 116(8): 2220.     CrossRef
    • A Regression-Based Approach for Mediation Analysis with Censored Data
      N A Salehah, Suhartono
      Journal of Physics: Conference Series.2021; 1752(1): 012009.     CrossRef
    • Benefit of a more extended pelvic lymph node dissection among patients undergoing radical prostatectomy for localized prostate cancer: A causal mediation analysis
      Marian S. Wettstein, Luke A. David, Clinsy Pazhepurackel, Aatif A. Qureshi, Alex Zisman, Michael Nesbitt, Karim Saba, Jaime O. Herrera‐Caceres, Christian D. Fankhauser, Ardalan Ahmad, Robert J. Hamilton, Daniel Eberli, Alexandre R. Zlotta, Tullio Sulser,
      The Prostate.2021; 81(5): 286.     CrossRef
    • The impact of adjusting for pure predictors of exposure, mediator, and outcome on the variance of natural direct and indirect effect estimators
      Awa Diop, Geneviève Lefebvre, Caroline S. Duchaine, Danielle Laurin, Denis Talbot
      Statistics in Medicine.2021; 40(10): 2339.     CrossRef
    • The Triglyceride-Glucose Index and Obesity-Related Risk of End-Stage Kidney Disease in Austrian Adults
      Josef Fritz, Wolfgang Brozek, Hans Concin, Gabriele Nagel, Julia Kerschbaum, Karl Lhotta, Hanno Ulmer, Emanuel Zitt
      JAMA Network Open.2021; 4(3): e212612.     CrossRef
    • Plasma Concentrations of Magnesium and Risk of Dementia: A General Population Study of 102 648 Individuals
      Jesper Qvist Thomassen, Janne S Tolstrup, Børge G Nordestgaard, Anne Tybjærg-Hansen, Ruth Frikke-Schmidt
      Clinical Chemistry.2021; 67(6): 899.     CrossRef
    • Use of acute care services by adults with a migrant background: a secondary analysis of a EurOOHnet survey
      Ellen Keizer, Oliver Senn, Morten Bondo Christensen, Linda Huibers
      BMC Family Practice.2021;[Epub]     CrossRef
    • Adipose Tissue Redox Microenvironment as a Potential Link between Persistent Organic Pollutants and the 16-Year Incidence of Non-hormone-Dependent Cancer
      Vicente Mustieles, Francisco M. Pérez-Carrascosa, Josefa León, Theis Lange, Jens-Peter Bonde, Celia Gómez-Peña, Francisco Artacho-Cordón, Rocío Barrios-Rodríguez, Rocío Olmedo-Requena, José Expósito, José J. Jiménez-Moleón, Juan P. Arrebola
      Environmental Science & Technology.2021; 55(14): 9926.     CrossRef
    • Discussion on “Causal mediation of semicompeting risks” by Yen‐Tsung Huang
      Mats J. Stensrud, Jessica G. Young, Torben Martinussen
      Biometrics.2021; 77(4): 1160.     CrossRef
    • Physical Activity From Adolescence Through Midlife and Associations With Body Mass Index and Endometrial Cancer Risk
      Pedro F Saint-Maurice, Joshua N Sampson, Kara A Michels, Steven C Moore, Erikka Loftfield, Kathleen McClain, Michael B Cook, Britton Trabert, Charles E Matthews
      JNCI Cancer Spectrum.2021;[Epub]     CrossRef
    • Effects of Patient and Surgery Characteristics on Persistent Postoperative Pain
      Xinlei Mi, Baiming Zou, Parisa Rashidi, Raheleh Baharloo, Roger B. Fillingim, Margaret R. Wallace, Paul L. Crispen, Hari K. Parvataneni, Hernan A. Prieto, Chancellor F. Gray, Tiago N. Machuca, Steven J. Hughes, Gregory J.A. Murad, Elizabeth Thomas, Atif I
      The Clinical Journal of Pain.2021; 37(11): 803.     CrossRef
    • Parametric-Regression–Based Causal Mediation Analysis of Binary Outcomes and Binary Mediators: Moving Beyond the Rareness or Commonness of the Outcome
      Mariia Samoilenko, Geneviève Lefebvre
      American Journal of Epidemiology.2021; 190(9): 1846.     CrossRef
    • Non-Persistence with Medication as a Mediator for the Social Inequality in Risk of Major Adverse Cardiovascular Events in Patients with Incident Acute Coronary Syndrome: A Nationwide Cohort Study
      Christina Boesgaard Graversen, Jan Brink Valentin, Mogens Lytken Larsen, Sam Riahi, Teresa Holmberg, Søren Paaske Johnsen
      Clinical Epidemiology.2021; Volume 13: 1071.     CrossRef
    • Association of race and ethnicity with clinical phenotype, genetics, and survival in pediatric acute myeloid leukemia
      Shannon E. Conneely, Casey L. McAtee, Rohit Gupta, Joseph Lubega, Michael E. Scheurer, Rachel E. Rau
      Blood Advances.2021; 5(23): 4992.     CrossRef
    • The conduct and reporting of mediation analysis in recently published randomized controlled trials: results from a methodological systematic review
      Tat-Thang Vo, Cecilia Superchi, Isabelle Boutron, Stijn Vansteelandt
      Journal of Clinical Epidemiology.2020; 117: 78.     CrossRef
    • Performance of methods to conduct mediation analysis with time‐to‐event outcomes
      Lizbeth Burgos Ochoa, Judith J.M. Rijnhart, Brenda W. Penninx, Klaas J. Wardenaar, Jos W.R. Twisk, Martijn W. Heymans
      Statistica Neerlandica.2020; 74(1): 72.     CrossRef
    • Mediation analysis for estimating cardioprotection of longitudinal RAS inhibition beyond lowering blood pressure and albuminuria in type 1 diabetes
      Jingchuan Guo, Ashley I. Naimi, Maria M. Brooks, Matthew F. Muldoon, Trevor J. Orchard, Tina Costacou
      Annals of Epidemiology.2020; 41: 7.     CrossRef
    • Contribution of income and job strain to the association between education and cardiovascular disease in 1.6 million Danish employees
      Elisabeth Framke, Jeppe Karl Sørensen, Per Kragh Andersen, Annemette Coop Svane-Petersen, Kristina Alexanderson, Jens Peter Bonde, Kristin Farrants, Esben Meulengracht Flachs, Linda L Magnusson Hanson, Solja T Nyberg, Ebbe Villadsen, Mika Kivimäki, Reiner
      European Heart Journal.2020; 41(11): 1164.     CrossRef
    • Reducing socio-economic inequalities in all-cause mortality: a counterfactual mediation approach
      Jessica E Laine, Valéria T Baltar, Silvia Stringhini, Martina Gandini, Marc Chadeau-Hyam, Mika Kivimaki, Gianluca Severi, Vittorio Perduca, Allison M Hodge, Pierre-Antoine Dugué, Graham G Giles, Roger L Milne, Henrique Barros, Carlotta Sacerdote, Vittorio
      International Journal of Epidemiology.2020; 49(2): 497.     CrossRef
    • IV estimation without distributional assumptions
      Theis Lange, Aksel K. G. Jensen
      Biometrical Journal.2020; 62(3): 688.     CrossRef
    • Use of antimalarial drugs is associated with a lower risk of preeclampsia in lupus pregnancy: A prospective cohort study
      Miguel Ángel Saavedra, Dafhne Miranda‐Hernández, Alejandra Lara‐Mejía, Antonio Sánchez, Sara Morales, Claudia Cruz‐Reyes, Pilar Cruz‐Domínguez, Gabriela Medina, Luis Javier Jara
      International Journal of Rheumatic Diseases.2020; 23(5): 633.     CrossRef
    • Education and incident type 2 diabetes: quantifying the impact of differential exposure and susceptibility to being overweight or obese
      Jimmi Mathisen, Aksel K. G. Jensen, Ingelise Andersen, Gregers S. Andersen, Ulla A. Hvidtfeldt, Naja H. Rod
      Diabetologia.2020; 63(9): 1764.     CrossRef
    • How polluted is your fat? What the study of adipose tissue can contribute to environmental epidemiology
      Vicente Mustieles, Juan P Arrebola
      Journal of Epidemiology and Community Health.2020; 74(5): 401.     CrossRef
    • Analytical approaches to testing pathways linking greenspace to health: A scoping review of the empirical literature
      Angel M. Dzhambov, Matthew H.E.M. Browning, Iana Markevych, Terry Hartig, Peter Lercher
      Environmental Research.2020; 186: 109613.     CrossRef
    • Association Between Chronic Hepatitis C Virus Infection and Myocardial Infarction Among People Living With HIV in the United States
      Jessica Williams-Nguyen, Stephen E Hawes, Robin M Nance, Sara Lindström, Susan R Heckbert, H Nina Kim, W Chris Mathews, Edward R Cachay, Matt Budoff, Christopher B Hurt, Peter W Hunt, Elvin Geng, Richard D Moore, Michael J Mugavero, Inga Peter, Mari M Kit
      American Journal of Epidemiology.2020; 189(6): 554.     CrossRef
    • Associations Between Complementary and Alternative Medicine and Health Care Expenditures Among Adults With Chronic Back Pain
      R. Henry Olaisen
      Medical Care.2020; 58(8): 689.     CrossRef
    • Inflammation, non-endothelial dependent coronary microvascular function and diastolic function—Are they linked?
      Hannah E. Suhrs, Jakob Schroder, Kira B. Bové, Naja D. Mygind, Daria Frestad, Marie M. Michelsen, Theis Lange, Ida Gustafsson, Jens Kastrup, Eva Prescott, Salvatore De Rosa
      PLOS ONE.2020; 15(7): e0236035.     CrossRef
    • Direct and indirect pathways between low income status and becoming a high-cost health care user in Ontario, Canada: a mediation analysis of health risk behaviors
      Luke Mondor, Tristan Watson, Kathy Kornas, Catherine Bornbaum, Walter P. Wodchis, Laura C. Rosella
      Annals of Epidemiology.2020; 51: 28.     CrossRef
    • Social bonds do not mediate the relationship between early adversity and adult glucocorticoids in wild baboons
      Stacy Rosenbaum, Shuxi Zeng, Fernando A. Campos, Laurence R. Gesquiere, Jeanne Altmann, Susan C. Alberts, Fan Li, Elizabeth A. Archie
      Proceedings of the National Academy of Sciences.2020; 117(33): 20052.     CrossRef
    • No effect of plasma trimethylamine N-Oxide (TMAO) and plasma trimethyllysine (TML) on the association between choline intake and acute myocardial infarction risk in patients with stable angina pectoris
      Anthea Van Parys, Vegard Lysne, Jannike Øyen, Jutta Dierkes, Ottar Nygård
      Human Nutrition & Metabolism.2020; 21: 200112.     CrossRef
    • Direct effect and indirect effect on an outcome under nonlinear modeling
      Kai Wang
      The International Journal of Biostatistics.2020;[Epub]     CrossRef
    • The effect of early-life and adult socioeconomic position on development of lifestyle-related diseases
      Linda Ejlskov, Henrik Bøggild, Claus D Hansen, Jesper Wulff, Steen M Hansen, Liis Starkopf, Theis Lange, Thomas Gerds, Christian Torp-Pedersen
      European Journal of Public Health.2019; 29(3): 562.     CrossRef
    • The Role of Genetic Variant rs13266634 in SLC30A8/ZnT8 in Postoperative Hyperglycemia After Major Abdominal Surgery
      Deborah Brunke-Reese, Paddy Ssentongo, Anna E Ssentongo, Brett E Phillips, Eric M Pauli, Arthur Berg, Shannon L Kelleher, David I Soybel
      The Journal of Clinical Endocrinology & Metabolism.2019; 104(9): 3877.     CrossRef
    • Hormone replacement therapy, mammographic density, and breast cancer risk: a cohort study
      Shadi Azam, Theis Lange, Stephanie Huynh, Arja R. Aro, My von Euler-Chelpin, Ilse Vejborg, Anne Tjønneland, Elsebeth Lynge, Zorana J. Andersen
      Cancer Causes & Control.2018; 29(6): 495.     CrossRef
    • Disparidade racial na sobrevivência em 10 anos para o câncer de mama: uma análise de mediação usando abordagem de respostas potenciais
      Mário Círio Nogueira, Maximiliano Ribeiro Guerra, Jane Rocha Duarte Cintra, Camila Soares Lima Corrêa, Vívian Assis Fayer, Maria Teresa Bustamante-Teixeira
      Cadernos de Saúde Pública.2018;[Epub]     CrossRef
    • Lorcaserin and metabolic disease: weight‐loss dependent and independent effects
      H. Bays, C. Perdomo, E. Nikonova, R. Knoth, M. Malhotra
      Obesity Science & Practice.2018; 4(6): 499.     CrossRef
    • Angiopoietin-2: A Potential Mediator of the Glycocalyx Injury in Adult Nephrotic Patients
      Maria Chaves, Matheus Mendes, Maximilian Schwermann, Raquel Queiroz, Regina Coelho, Francisco Salmito, Gdayllon Meneses, Alice Martins, Ana Moreira, Alexandre Libório
      Journal of Clinical Medicine.2018; 7(11): 401.     CrossRef

    Figure

    Epidemiol Health : Epidemiology and Health