Directed acyclic graphs (DAGs; e.g., Greenland, et al, 1999) are a part of a formalism of causal inference from the counterfactual interpretation of causality camp. In these graphs the presence of an arrow from variable $A$ to variable $B$ asserts that variable $A$ directly causes (some change in risk of) variable $B$, and the absence of such an arrow asserts that variable $A$ does not directly cause (some change in risk of) variable $B$.

As an example, the statement "tobacco smoke exposure directly causes a change in risk of mesothelioma" is represented by the black arrow from "tobacco smoke exposure" to "mesothelioma" in the not a DAG causal diagram below.

Likewise, the statement "asbestos exposure directly causes a change in risk of mesothelioma" is represented by the black arrow from "asbestos exposure" to "mesothelioma" in the not a DAG causal graph below.

I use the term not a DAG to describe the below causal graph because of the red arrow, which I intend to assert something like "asbestos exposure causes a change in the direct causal effect of tobacco smoke exposure on risk of mesothelioma" (asbestos does physical damage to the cells of the lung that, in addition to directly causing a change in risk of mesothelioma, also renders the cells more susceptible to the carcinogenic harms of tobacco smoke exposure with the result that exposure to both asbestos and tobacco result in an increase in risk that is more than the sum of the two separate risks), and this does not quite fit with the formal meaning of causal arrows in DAGs I described at the start of my question (i.e. because the red arrow does not terminate in a variable).

Not a Directed Acyclic Graph: tobacco smoke exposure causes increased mesothelioma risk; asbestos exposure causes increased mesothelioma risk; asbestos exposure causes an increase in the causal effect of tobacco on mesothelioma risk.

How does one correctly represent interaction effects within the visual formalism of a DAG?


Greenland, S., Pearl, J., and Robins, J. M. (1999). Causal diagrams for epidemiologic research. Epidemiology, 10(1):37–48.

  • $\begingroup$ Is it possible to change the DAG by adding a node for lung cell damage with arrows from tobacco and asbestos pointed in and an arrow to mesothelioma out? That keeps the DAG and provides flexibility to model the individual and compound effects. $\endgroup$
    – John Woolf
    Commented Jun 23, 2015 at 14:35
  • $\begingroup$ Yes, I had thought of something like that, but I wonder if that also implies that tobacco's main effect is not also 'lung cell damage,' when it is. I suppose one could add another variable 'tobacco and asbestos together' caused by tobacco and by asbestos, but this seems clunky... but perhaps it is what one has to go with. $\endgroup$
    – Alexis
    Commented Jun 23, 2015 at 14:44
  • $\begingroup$ What software did you use to create the graph? $\endgroup$ Commented Jun 7, 2018 at 4:32
  • $\begingroup$ @zipzapboing drew it by hand in MS Word, or PageStream or Scribus or something, then edited for color & transparency using GIMP. (It's actually an excerpt from some of my lecture slides.) $\endgroup$
    – Alexis
    Commented Jun 7, 2018 at 6:36

4 Answers 4


Pearl's theory of causality is completely non-parametric. Interactions are not made explicit because of that, neither in the graph nor in the structural equations it represents. However, causal effects can vary (wildly) by assumption.

If an effect is identified and you estimate it from data non-parametrically, you obtain a complete distribution of causal effects (instead of, say, a single parameter). Accordingly, you can evaluate the causal effect of tobacco exposure conditional on asbestos exposure non-parametrically to see whether it changes, without committing to any functional form.

Let's have a look at the structural equations in your case, which correspond to your "DAG" stripped of the red arrow:

Mesothelioma = $f_{1}$(Tobacco, Asbestos, $\epsilon_{m}$)

Tobacco = $f_{2}$($\epsilon_{t}$)

Asbestos = $f_{3}$($\epsilon_{a}$)

where the $\epsilon$ are assumed to be independent because of missing dashed arrows between them.

We have left the respective functions f() and the distributions of the errors unspecified, except for saying that the latter are independent. Nonetheless, we can apply Pearl's theory and immediately state that the causal effects of both tobacco and asbestos exposure on mesothelioma are identified. This means that if we had infinitely many observations from this process, we could exactly measure the effect of setting the exposures to different levels by simply seeing the incidences of mesothelioma in individuals with different levels of exposure. So we could infer causality without doing an actual experiment. This is because there exist no back-door paths from the exposure variables to the outcome variable.

So you would get

P(mesothelioma | do(Tobacco = t)) = P(mesothelioma | Tobacco = t)

The same logic holds for the causal effect of asbestos, which allows you to simply evaluate:

P(mesothelioma | Tobacco = t, Asbestos = a) - P(mesothelioma | Tobacco = t', Asbestos = a)

in comparison to

P(mesothelioma | Tobacco = t, Asbestos = a') - P(mesothelioma | Tobacco = t', Asbestos = a')

for all relevant values of t and a in order to estimate the interaction effects.

In your concrete example, let's assume that the outcome variable is a Bernoulli variable - you can either have mesothelioma or not - and that a person has been exposed to a very high asbestos level a. Then, it is very likely that he will suffer from mesothelioma; accordingly, the effect of increasing tobacco exposure will be very low. On the other hand, if asbestos levels a' are very low, increasing tobacco exposure will have a greater effect. This would constitute an interaction between the effects of tobacco and asbestos.

Of course, non-parametric estimation can be extremely demanding and noisy with finite data and lots of different t and a values, so you might think about assuming some structure in f(). But basically you can do it without that.

  • $\begingroup$ Julian, thank you for the in-depth on Pearl's contrafactual causality logic. Your answer would be improved by speaking to the visual representation of interaction effects that was central to my question (Perhaps by saying "DAGs do not make interaction effects explicit," or what have you.) $\endgroup$
    – Alexis
    Commented Jul 21, 2015 at 16:51
  • 1
    $\begingroup$ Yes, you are right; I've tried to make that clearer in the very beginning. $\endgroup$ Commented Jul 21, 2015 at 19:26
  • $\begingroup$ Julian, that's great. You say that the building blocks of a DAG are not individual links, C=f(D) but more generally sets of links, C=f(D,E,...). My problem is that Pearl is sometimes vague on this, e.g. in "Causality" his Def 2.2.1.: "A causal structure of a set of variables V is a DAG in which each node corresponds to a distinct element of V, and each link represents direct functional relationship among the corresponding variables." Surely the individual links from the parents of C to C mean nothing on their own, they only show participation in a collective influence on C. You'd agree? $\endgroup$ Commented Apr 8, 2019 at 20:07
  • 1
    $\begingroup$ Yes. C = f(D, E) would be shown by two links in the graph from D/E to C. The interpretation of the links is that D and E determine C. So links, of course, are relational. $\endgroup$ Commented Apr 9, 2019 at 10:20

The simple answer is that you already do. Conventional DAGs do not only represent main effects but rather the combination of main effects and interactions. Once you have drawn your DAG, you already assume that any variables pointing to the same outcome can modify the effect of the others pointing to the same outcome. It is a modeling assumption, separate from the DAG, which presumes the lack of an interaction.

In addition, interaction can occur without including an explicit interaction term in your model. If you include main effects only in a model for the risk ratio of Y with respect to treatment T and covariate Q, the estimate of the risk difference will differ depending on the level of Q. In order to accommodate all these possibilities nonparametrically, DAGs make only the weakest assumptions on the functional form of the relationships among the variables, and assuming no interaction is a stronger assumption that allowing for an interaction. This again is to say that DAGs already allow for interaction without any adjustment. See Vanderweele (2009) for a discussion of interaction that uses conventional DAGs but allows for interaction.

Bollen & Paxton (1998) and Muthén & Asparouhov (2015) both demonstrate interactions in path models with latent variables, but these interactions explicitly refer to product terms in a parametric model rather than to interactions broadly. I have also seen diagrams similar to yours where the causal arrow points to a path, but strictly speaking a path is not a unique quantity that a variable can have a causal effect on (even though that may be how we want to interpret our models); it simply represents the presence of a causal effect, not its magnitude.

Bollen, K. A., & Paxton, P. (1998). Interactions of latent variables in structural equation models. Structural Equation Modeling: A Multidisciplinary Journal, 5(3), 267-293.

Asparouhov, T. & Muthén, B. (2020): Bayesian estimation of single and multilevel models with latent variable interactions, Structural Equation Modeling: A Multidisciplinary Journal

VanderWeele, T. J. (2009). On the distinction between interaction and effect modification. Epidemiology, 20(6), 863-871.

  • 4
    $\begingroup$ Also: "Once you have drawn your DAG, you already assume that any variables pointing to the same outcome can modify the effect of the others pointing to the same outcome. It is a modeling assumption, separate from the DAG, which presumes the lack of an interaction." is brilliant... it really helped these ideas land with me. $\endgroup$
    – Alexis
    Commented Jun 16, 2018 at 17:06
  • $\begingroup$ Noah was the cited 2015 paper by Múthven and Asparouhov referring to an then unpublished version of Tihomir Asparouhov & Bengt Muthén (2020): Bayesian estimation of single and multilevel models with latent variable interactions, Structural Equation Modeling: A Multidisciplinary Journal? $\endgroup$
    – Alexis
    Commented Nov 20, 2020 at 3:23
  • $\begingroup$ @Alexis I do believe so. $\endgroup$
    – Noah
    Commented Nov 20, 2020 at 19:23

A new method of representing interactions by creating dedicated nodes was proposed and termed "IDAG" since this question was asked. In my understanding the example sentence from question "asbestos exposure causes a change in the direct causal effect of tobacco smoke exposure on risk of mesothelioma" would be represented as:

Asbestos → Δ Mesothelioma Tobacco

The limitation of this approach is that you would need to represent all other effects on Tobacco and Mesothelioma separately.

For details see:

Anton Nilsson and others, A directed acyclic graph for interactions, International Journal of Epidemiology, Volume 50, Issue 2, April 2021, Pages 613–619, https://doi.org/10.1093/ije/dyaa211


If you want to estimate the non-separable non-linear structural equations directly, there is a growing econometrics literature on this. You do, of course, need to make some assumptions in order to ensure statistical identification (even if you have built a defensible case for causal identification using graphical criteria and causal calculus), but these are not as restrictive as in the linear or parametric nonlinear case. Note that nonparametric quantile regression is equivalent to a non-separable model under certain conditions, so that gives you a fairly feasible implementation option.

Breunig, C. (2020). Specification Testing in Nonparametric Instrumental Quantile Regression. Econometric Theory. https://doi.org/10.1017/S0266466619000288

Dunker, F. (April 16, 2020). Nonparametric instrumental variable regression and quantile regression with full independence. Repository arXiv. https://arxiv.org/pdf/1511.03977.pdf

Chernozhukov, V., Fernández-Val, I., Newey, W., Stouli, S., & Vella, F. (2020). Semiparametric Estimation of Structural Functions in Nonseparable Triangular Models. Quantitative Economics, 11, 503-633. https://qeconomics.org/ojs/index.php/qe/article/viewFile/1328/1320

Babii, A., & Florens, S.P. (Jan. 30, 2020). Are unobservables separable? Repository arXiv. http://arxiv.org/pdf/1705.01654

Su, L., Tu, Y., & Ullah, A. (2015). Testing Additive Separability of Error Term in Nonparametric Structural Models. Econometric Reviews, 34(6-10), 1057-1088. https://doi.org/10.1080/07474938.2014.956621

Lu, X., & White, H. (2014). Testing for separability in structural equations. Journal of Econometrics, 182(1), 14–26. https://doi.org/10.1016/j.jeconom.2014.04.005

  • 1
    $\begingroup$ Thank you, Cameron. Can you edit your answer to specifically address the graphical DAG formalism which is the central substance of my question? Also: welcome to CV! :) $\endgroup$
    – Alexis
    Commented Nov 21, 2020 at 20:15

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.