## To Adjust or Not to Adjust? Sensitivity Analysis of M -Bias and Butterfly-Bias

Peng Ding and Luke Miratrix Department of Statistics, Harvard University lmiratrix@stat.harvard.edu

Emails: pengding@fas.harvard.edu and

## Abstract

' M -Bias,' as it is called in the epidemiologic literature, is the bias introduced by conditioning on a pretreatment covariate due to a particular ' M -Structure' between two latent factors, an observed treatment, an outcome, and a 'collider.' This potential source of bias, which can occur even when the treatment and the outcome are not confounded, has been a source of considerable controversy. We here present formulae for identifying under which circumstances biases are inflated or reduced. In particular, we show that the magnitude of M -Bias in linear structural equation models tends to be relatively small compared to confounding bias, suggesting that it is generally not a serious concern in many applied settings. These theoretical results are consistent with recent empirical findings from simulation studies. We also generalize the M -Bias setting (1) to allow for the correlation between the latent factors to be nonzero, and (2) to allow for the collider to be a confounder between the treatment and the outcome. These results demonstrate that mild deviations from the M -Structure tend to increase confounding bias more rapidly than M -Bias, suggesting that choosing to condition on any given covariate is generally the superior choice. As an application, we re-examine a controversial example between Professors Donald Rubin and Judea Pearl.

Key Words : Causality; Collider; Confounding; Controversy; Covariate.

## 1 Introduction

The hallmark of an observational study is selection bias (Heckman, 1979, Copas and Li, 1997, Hern´ an, Hernandez-Diaz, and Robins, 2004). Many statisticians believe that 'there is no reason to avoid adjustment for a variable describing subjects before treatment' in observational studies (Rosenbaum, 2002, pp 76), because 'typically, the more conditional an assumption, the more generally acceptable it is' (Rubin, 2009). This advice, recently dubbed the 'pretreatment criterion' (VanderWeele and Shpitser, 2011), is widely used in empirical studies, as more covariates generally seem to make the ignorability assumption, i.e., the assumption that conditionally on the observed pretreatment covariates, treatment assignment is independent of the potential outcomes (Rosenbaum and Rubin, 1983), more plausible. And, as the validity of causal inference in observational studies relies strongly on this (untestable) assumption (Rosenbaum and Rubin, 1983), it seems reasonable to make all efforts to render it plausible.

However, other researchers (Pearl, 2009b,c, Shrier, 2008, 2009, Sj¨lander, o 2009), mainly from the causal diagram community, do not accept this view because of the possibility of a so-called M -Structure, illustrated in Figure 1(c). In sharp contrast to Rubin and Rosenbaum's advice, Pearl (2009b) and Pearl (2009c) warn practitioners that spurious bias may arise due to adjusting for a collider M in an M -Structure, even if it is a pretreatment covariate . This form of bias, typically called M -bias, a special version of so-called 'collider bias,' has since generated considerable controversy and confusion.

We attempt to resolve some of these debates by an analysis of M -bias under the causal diagram or directed acyclic graph (DAG) framework. For readers unfamiliar with the terminologies from the DAG (or Bayesian Network) literature, more details can be found in Pearl (1995) or Pearl (2009a). We here use only a small part of this larger framework. Arguably the most important structure in the DAG, and certainly the one at root of almost all controversy, is the 'V-Structure' illustrated in Figure 1(b). Here, U and W are marginally independent with a common outcome M , which shapes a 'V' with the vertex M being called a 'collider.' From a data-generation viewpoint, one might imagine Nature generating data in two steps: She first picks independently two values for U and W from two distributions, and then she combines them (possibly along with some additional random variable) to create M . Given this, conditioning on M can cause a spurious correlation between U and W , which is known as the collider bias (Greenland, 2002), or, in epidemiology, Berkson's Paradox (Berkson, 1946). Conceptually, this correlation happens because if one cause of an observed outcome is known to have not occurred, the other cause becomes more likely. Consider an automatic-timer sprinkler system where the sprinkler being on is independent of whether it is raining. Here, the weather gives no information on the sprinkler. However, given wet grass, if one observes a sunny day, one will likely conclude that the sprinklers have recently run. Correlation has been induced.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0324v1-figure-b73e3f93-89b9-454e-b731-3867ed554d08.jpg)

Where things get interesting is when this collider is made into a pre-treatment variable. Consider Figure 1(c), an extension of Figure 1(b). Here U and W are now also causes of the treatment T and the outcome Y , respectively. Nature, as a last, third step generates T as a function of U and some randomness, and Y as a function of W and some randomness. This structure is typically used to represent a circumstance where a researcher observes T , Y , and M in nature and is attempting to derive the causal impact of T on Y . U and W are unobserved, or latent. Clearly, the causal effect of T on Y is zero, which is also equal to the marginal association between T and Y . If a researcher regressed Y on T , he or she would obtain a zero in expectation, which is correct for estimating the causal effect. But perhaps there is a concern that M , a pretreatment covariate, may be a confounder that is masking a treatment effect. Typically, one would then 'adjust' for M to take this possibility into account, e.g., by including M in a regression or by matching units on similar values of M . If we do this in this circumstance, however, then we will not find a zero causal effect, in expectation. This is the so-called ' M -Bias,' and this special structure is called the ' M -Structure' in the DAG literature.

Previous qualitative analysis for binary variables shows that collider bias generally tends to be small (Greenland, 2002), and simulation studies (Liu, Brookhart, Schneeweiss, Mi, and Setoguchi, 2012) again demonstrate that M -Bias is small in many realistic settings. While mathematically describing the magnitudes of M -Bias in general models is intractable, it is possible to derive exact formulae of the biases as functions of the correlation coefficients in linear structural equation models (LSEMs). The LSEM has a long history in statistics (Wright, 1921, 1934) to describe dependence among multiple random variables. Sprites (2002) uses linear models to illustrate M -Bias in observational studies, and Pearl (2013) also utilize the transparency of such linear models to examine various types of causal phenomena, biases, and paradoxes. We here extend these works and provide exact formulae for biases, allowing for a more detailed quantitative analysis of M -bias.

While M -Bias does exist when the true underlying data generating process (DGP) follows the exact M -Structure, it might be rather sensitive to various deviations from the exact M -Structure. Furthermore, some might argue that an exact M -Structure is unlikely to hold in practice. Gelman (2011), for example, doubts the exact independence assumption required for the M -Structure in the social sciences by arguing that there are '(almost) no true zeros' in this discipline. Indeed, since U and W are often latent characteristics of the same individual, the independence assumption U W is a rather strong structural assumption. Furthermore, it might be plausible that the pretreatment covariate M is also a confounder between, i.e., has some causal impact on both, the treatment and outcome. We extend our work by accounting for these departures from a pure M -Structure, and find that even slight departures from the M -Structure can dramatically change the forms of the biases.

This paper theoretically compares the bias from conditioning on an M to not under several scenarios and finds that M -Bias is indeed small relative to other concerns unless there is a strong correlation structure for the variables. We further show that these findings extend to a binary treatment regime as well. This argument proceeds in several stages. First, in Section 2, we examine a pure M -Structure and introduce our LSEM framework. We then discuss the cases when the latent variables U and W may be correlated and M may also be a confounder between the treatment T and the outcome Y . In Section 3, we generalize the results in Section 2 to a binary treatment. In Section 4, we illustrate the theoretical findings using a controversial example between Professors Donald Rubin and Judea Pearl (Rubin, 2007, Pearl, 2009c). Section 5 discusses the relevance of our findings by examining M -Bias in actual practice and by comparing asymptotic to finite sample properties. We conclude with a brief discussion and present all technical details in the Appendix.

## 2 M -Bias and Butterfly-Bias in LSEMs

We begin by examining pure M -Bias in a LSEM. As our primary focus is bias, we assume data are ample and that anything estimable is estimated with nearly perfect precision. In particular, when we say we obtain a result from a regression, we implicitly mean we obtain that result in expectation; in practice an estimator will be near the given quantities. We do not compare relative uncertainties of different estimators given the need to estimate more or fewer parameters. There are likely degrees-of-freedom issues that would implicitly advocate using estimators with fewer parameters, but in the circumstances considered here these concerns are likely to be minor as all the models have few parameters.

A causal DAG can be viewed as a hierarchical DGP. In particular, any variable on the graph can be viewed as a function of its parents and some additional noise, i.e., if R had parents A B , , and C , we would have

$$R = f ( A, B, C, \varepsilon _ { R } ) \text{ with } \varepsilon _ { R } \bigsqcup ( A, B, C ).$$

Generally noise terms such as e R are considered to be independent from each other, but they can also be given an unknown correlation structure corresponding to earlier variables not explicitly included in the diagram. This is typically represented by drawing the dependent noise terms jointly from some multivariate distribution. This framework is quite general; we can represent any distribution that can be factored as a product of conditional distributions corresponding to a DAG (which is one representation of the Markov Condition, a fundamental assumption for DAGs).

LSEMsare special cases of the above with additional linearity and additivity constraints. For simplicity, and without loss of generality, we also rescale all primary variables ( U W M T Y , , , , ) to have zero mean and unit variance. For example, consider this data generating process corresponding to Figure 1(a):

$$\begin{cases} \begin{array} { c c } M, \epsilon _ { T }, \epsilon _ { Y } & \text{ind} \\ T & = & a M + \sqrt { 1 - a ^ { 2 } } \epsilon _ { T }, \\ Y & = & b T + c M + \sqrt { 1 - b ^ { 2 } - c ^ { 2 } } \epsilon _ { Y }, \end{array} \end{cases}$$

where we use A ∼ [ 0 1 , ] to denote a random variable with mean zero and variance one.

In the causal DAG literature, we think about causality as reaching in and fixing a given node to a set value, but letting Nature take her course otherwise. For example, if we were able to set T at t , the above data generation process would be transformed to:

$$\begin{cases} \begin{array} { c c } M, \varepsilon _ { T }, \varepsilon _ { Y } & \text{ind} \\ T & = & t, \\ Y & = & b t + c M + \sqrt { 1 - b ^ { 2 } - c ^ { 2 } } \varepsilon _ { Y }. \end{array} \end{cases}$$

The previous cause, M , of T has been broken, but the impact of T on Y remains intact. This changes the distribution of Y but not M . More importantly, this results in a distribution distinct from that of conditioning on T = t . Consider the case of positive a b , , and c . If we observe a high T , we can infer a high M (as T and M are correlated) and a high Y due to both the bT and cM terms in Y 's equation. However, if we set T to a high value, M is unchanged. Thus, while we will still have the large bT term for Y , the cM term will be 0 in expectation. Thus, the expected value for Y will be less.

This setting as compared to conditioning is represented with the 'do' operator. Given the 'do' operator, we define a local causal effect of T on Y at T = t as:

$$\tau _ { t } = \frac { \partial \mathbb { E } \{ Y \ | \, \text{do} ( T ) = t \} } { \partial t }.$$

For linear models, the local causal effect is a constant, and thus we do not need to specify t . We use 'do' here purely to indicate the different distributions. For a more technical overview, see Pearl (1995) or Pearl (2009a). Our results, with more formality, can easily be expressed in this more technical notation.

Figure 2: M -Structure with Possibly Correlated Hidden Causes

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0324v1-figure-64cabb63-b689-49d1-ae49-8994d1f27f8a.jpg)

If we extend the M -Structure in Figure 1(c) by allowing possible correlation between the two hidden causes U and W , we obtain the DAG in Figure 2. This in turn gives the following DGP:

$$\begin{cases} \begin{array} { c } \varepsilon _ { M }, \varepsilon _ { T }, \varepsilon _ { Y } & \text{ind} \\ ( U, W ) & \end{array} \begin{array} { c } \text{$\varepsilon_{\text{$M}$} } & [ 0, 1 ], \\ ( U, W ) & \end{array} \begin{array} { c } \text{$\varepsilon_{\text{$O}$}, 1, 1, \rho ], \\ M & = & b U + c W + \sqrt{1-b^{2}-c^{2}} \varepsilon _ { M }, \\ T & = & a U + \sqrt{1-a^{2}} \varepsilon _ { T }, \\ Y & = & d W + \sqrt{1-d^{2}} \varepsilon _ { Y }, \end{array} \begin{array} { c } \text{$\varepsilon_{\text{$O}$}, 1, 1, \rho ], \\ \end{array} \begin{array} { c } \text{$\varepsilon_{\text{$O}$}, 1, 1, \rho ], \\ \end{array} \begin{array} { c } \text{$\varepsilon_{\text{$O}$}, 1, 1, \rho ], \\ \end{array} \begin{array} { c } \text{$\varepsilon_{\text{$O}$}, 1, 1, \rho ], \\ \end{array} \begin{array} { c } \text{$\varepsilon_{\text{$O}$}, 1, 1, \rho ], \\ \end{array} \begin{array} { c } \text{$\varepsilon_{\text{$O}$} }. \end{array}.$$

where we use ( A B , ) ∼ [ 0 0;1 1 , , , r ] to denote a bivariate random vector with means zero, variances one and correlation coefficient r .

Here, the true causal effect of T on Y is zero, namely, t t = 0 for all t . The unadjusted estimator for the causal effect obtained by regressing Y onto T is the same as the covariance between T and Y :

$$\text{Bias} _ { u n a d j } = \text{Cov} ( T, Y ) = \text{Cov} ( a U, d W ) = a d \text{Cov} ( U, W ) = a d \rho.$$

The adjusted estimator (see Lemma 2 in Appendix A for a proof) obtained by regressing Y onto ( T M , ) is

$$\text{Bias} _ { a d j } = \frac { a d \rho ( 1 - b ^ { 2 } - c ^ { 2 } - b c \rho ) - a b c d } { 1 - ( a b + a c \rho ) ^ { 2 } }.$$

The results above and some of the results discussed later in this paper can be obtained directly from traditional path analysis (Wright, 1921, 1934). However, we provide elementary proofs, which can easily be extended to binary treatment, in the Appendix. If we allowed for a treatment effect, our results would remain essentially unchanged; the only difference would be due to restrictions on the correlation terms needed to maintain unit variance for all variables.

The above can also be expressed in the potential outcomes framework (Neyman, 1923/1990, Rubin, 1974). In particular, for a given unit let Nature draw e M , e T , e Y , U , and V as before. Let T be the 'natural treatment' for that unit, i.e., what treatment it would receive sans intervention. Then calculate Y t ( ) for any t of interest using the 'do' operator. These are what we would see if we set T = t . How Y t ( ) changes for a particular unit defines that unit's collection of potential outcomes. Then E { Y t ( ) } for some t is the expected potential outcome over the population for a particular t . We can examine the derivative of this function as above to get a local treatment effect. This connection is exact: the findings in this paper are the same as what one would find using this DGP and the potential outcomes framework. We here examine regression as the estimator. Note that matching would produce identical results as the amount of data grew (assuming the data generating process ensures common support, etc.).

Exact M -Bias. The M -Bias originally considered in the literature is the special case where the correlation coefficient between U and W is r = 0. In this case, the unadjusted estimator is unbiased and the absolute bias of the adjusted estimator is | abcd | / { 1 -( ab ) 2 } . With moderate correlation coefficients a b c d , , , the denominator 1 -( ab ) 2 is close to one, and the bias is close to -abcd . Since abcd is a product of four correlation coefficients, it can be viewed as a 'higher order bias.' For example, if a = b = c = d = 0 2, then 1 . -( ab ) 2 = 0 9984 . ≈ 1, and the bias of the adjusted estimator is -abcd / { 1 -( ab ) 2 } = -0 0016 . ≈ 0; if a = b = c = d = 0 3, then 1 . -( ab ) 2 = 0 9919 . ≈ 1, and the bias of the adjusted estimator is -abcd / { 1 -( ab ) 2 } = -0 0082 . ≈ 0 . Even moderate correlation results in little bias.

In Figure 3, we plot the bias of the adjusted estimator as a function of the correlation coefficients, and let these coefficients change to see how the bias changes. In the first subfigure, we assume all the correlation coefficients have the same magnitude ( a = b = c = d ), and we plot the absolute bias of the adjusted estimator versus a . The constraints on variance and correlation only allow for some combinations of values for a b c , , and d which limits the domain of the figures. In this case, for example, | a | ≤ √ 2 2 due to the requirement that / b 2 + c 2 = 2 a 2 ≤ 1. Other figures have limited domains due to similar constraints. In the second subfigure of Figure 3, we assume that M is more predictive to the treatment T than to the outcome Y , with a = b = 2 c = 2 d . In the third subfigure of Figure 3, we assume that M is more predictive to the outcome Y , with 2 a = 2 b = c = d . The biases are generally very small within wide ranges of the feasible regions of the correlation coefficients. However, the biases do blow up when the correlation coefficients are extremely large. Near the boundary of the feasible regions in Figure 3, the M -Structure is approximately deterministic, which is rare in social sciences. Pearl (2009b) does not exclude the worst cases, and thus he considers M -Bias as a severe problem.

Figure 3: M -Bias with Independent U and W . The three subfigures correspond to the cases when ( U W , ) are equally/more/less predictive to the treatment than to the outcome. In each subfigure, we show the proportions of the areas where the adjusted estimator has a bias smaller than 0 01. .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0324v1-figure-1ef59469-1e57-4d68-b7a1-eaf9f4aa4d32.jpg)

In Figure 4(a), we assume a = b and c = d and examine a broader range of relationships. Here, the grey area satisfies | Bias ad j | &lt; min ( | a | , | c | ) / 20. For example, when the absolute values of the correlation coefficients are smaller than 0 5 . (the square with dashed boundary in Figure 4(a)), the corresponding area is almost grey, implying small bias.

Due to the four dimensional sensitivity parameters ( a b c d , , , ) , a full exploration and graphical illustration over all possible values of the sensitivity parameters is formidably hard. In the absence of prior knowledge about the DAG, our sensitivity analysis here is based on some simplifications (e.g., a = b and c = d ), which may reflect some real situations. Using the bias formulae in this paper, we can easily conduct sensitivity analysis for other parameter combinations, depending on our practical problem and background knowledge about the DAG.

As a side note, Pearl (2013) noticed a surprising fact: the stronger the correlation between T and M , the larger the absolute bias of the adjusted estimator, since the absolute bias is monotone increasing in | ab | . From the second and the third subfigure of Figure 3, we see that when M is more predictive of the treatment, the biases of the adjusted estimator indeed tends to be larger.

(a) Pure M -Bias. Within the grey region, the absolute bias of the adjusted estimator is less than 1 / 20 of the minimum of | a | (= | b | ) and | c | (= | d | ) .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0324v1-figure-e80d396f-6729-47be-a9f0-cc99a4d529c5.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0324v1-figure-b5ec7a96-16bc-4b6d-8904-c44729585ab6.jpg)

(b) M -Bias with Correlated U and W . Within the grey region, the adjusted estimator is superior.

Figure 4: M -Bias under Different Scenarios

glyph[negationslash]

Correlated Latent Variables. When the latent variables U and W are correlated with r = 0, both the unadjusted and adjusted estimators may be biased. The ques- tion then becomes: which is worse? The ratio of the absolute biases is

$$\left | \frac { \text{Bias} _ { a d j } } { \text{Bias} _ { u n a d j } } \right | = \left | \frac { \rho ( 1 - b ^ { 2 } - c ^ { 2 } - b c \rho ) - b c } { \rho \left \{ 1 - ( a b + a c \rho ) ^ { 2 } \right \} } \right |,$$

which does not depend on d (the relationship between W and Y ). For example, if the correlation coefficients a b c , , , r all equal 0 2, the ratio above is 0 714; in this . . case the adjusted estimator is superior to the unadjusted one by a factor of 1 4. . Figure 4(b) compares this ratio to 1 for all combinations of r and a (= b = ) c . Generally, the adjusted estimator has smaller bias except when a b , , and c are quite large.

In Figure 5, we again assume a = b = c = d and investigate the absolute biases as functions of a for fixed r at 0 1 . , 0 2 . , and 0 4 . . When the correlation coefficients a (= b = c = ) d are not dramatically larger than r , the adjusted estimator has smaller bias than the unadjusted one.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0324v1-figure-2ee88e68-9f8c-4b3c-98ba-c140d89caa90.jpg)

glyph[negationslash]

Figure 5: M -Bias with correlated U and W when r = 0 with a = b = c = d . In each subfigure, we show the proportion of the areas where the adjusted estimator has a smaller bias than the unadjusted estimator.

glyph[negationslash]

The Disjunctive Cause Criterion. In order to remove biases in observational studies, VanderWeele and Shpitser (2011) propose a new 'disjunctive cause criterion' for selecting confounders, which requires controlling for all the covariates that are either causes of the treatment, causes of the outcome, or causes of both. According to the 'disjunctive cause criterion,' when r = 0, we should control for

( U W , ) if possible. Unfortunately, neither of ( U W , ) is observable. However, controlling the 'proxy variable' M for ( U W , ) may reduce bias when r is relatively large. In the special case with b = 0, the ratio of the absolute biases is

$$\left | \frac { \text{Bias} _ { a d j } } { \text{Bias} _ { u n a d j } } \right | = \frac { 1 - c ^ { 2 } } { 1 - ( a c \rho ) ^ { 2 } } \leq 1 ;$$

in another special case with c = 0, the ratio of the absolute biases is

$$\left | \frac { \text{Bias} _ { a d j } } { \text{Bias} _ { u n a d j } } \right | = \frac { 1 - b ^ { 2 } } { 1 - ( a b ) ^ { 2 } } \leq 1.$$

Therefore, if either U or W is not causative to M , the adjusted estimator is always better than the unadjusted one.

Butterfly-Bias: M -Bias with Confounding Bias. Models, especially in the social sciences, are approximations. They rarely hold exactly. In particular, for any covariate M of interest, there is likely to be some concern that M is indeed a confounder, even if it is also a possible source of M -Bias. If we let M both be a confounder as well as the middle of an M -Structure we obtain a 'Butterfly-Structure' (Pearl, 2013) as shown in Figure 6. In this circumstance, conditioning will help with confounding bias, but hurt with M -Bias. Ignoring M will not resolve any confounding, but will avoid M -Bias. The question then becomes that of determining which is the lesser of the two evils.

Figure 6: Butterfly-Structure

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0324v1-figure-ce28232f-9677-4788-8b35-1f701ad7fc06.jpg)

We can examine this trade-off for a LSEM corresponding to Figure 6. The DGP is given by the following equations:

$$\begin{cases} \begin{array} { c c c } U, W, \varepsilon _ { M }, \varepsilon _ { T }, \varepsilon _ { Y } & \underset { \sim } { \text{$ind} } \\ M & = & b U + c W + \sqrt { 1 - b ^ { 2 } - c ^ { 2 } } \varepsilon _ { M }, \\ T & = & a U + e M + \sqrt { 1 - a ^ { 2 } - e ^ { 2 } } \varepsilon _ { T }, \\ Y & = & d W + f M + \sqrt { 1 - d ^ { 2 } - f ^ { 2 } } \varepsilon _ { Y }. \end{array} \end{cases}$$

Again, the true causal effect of T on Y is zero. The unadjusted estimator obtained by regressing Y onto T is the covariance between T and Y :

$$\text{Bias} _ { u n a d j } = \text{Cov} ( T, Y ) = a b f + c d e + e f.$$

It is not, in general, zero, implying bias. The adjusted estimator (see Lemma 3 in Appendix for a proof) obtained by regressing Y onto ( T M , ) has bias

$$\text{Bias} _ { a d j } = - \frac { a b c d } { 1 - ( a b + e ) ^ { 2 } }.$$

If the values of e and f are relatively high (i.e., M has a strong effect on both T and Y ), the confounding bias is large and the unadjusted estimator will be severely biased. For example, if a b c d e , , , , , and f all equal 0 2, the bias of the unadjusted . estimator is 0 056, but the bias of the adjusted estimator is only . -0 0017, an order . of magnitude smaller. Generally, the largest term for the unadjusted bias is the second-order term of e f , while the adjusted bias only has, ignoring the denominator, a fourth-order term of abcd . This suggests adjustment is generally preferable and that M -bias is in some respect a 'higher order bias.'

Detailed comparison of the ratio of the biases is difficult, since we can vary six parameters ( a b c d e , , , , , f ) . In Figure 7(a), we assume all the correlation coefficients have the same magnitude, and plot bias for both estimators as a function of the correlation coefficient within the feasible region, defined by the restrictions -√ 2 2 / &lt; a &lt; ( -1 + √ 5 ) / 2, due to the restrictions

$$b ^ { 2 } + c ^ { 2 } < 1, \, a ^ { 2 } + e ^ { 2 } < 1, \, d ^ { 2 } + f ^ { 2 } < 1, \, \text{and} \, | a ^ { 2 } + e | < 1.$$

Within 74 9% of the feasible region, the adjusted estimator has smaller bias than the . unadjusted one. The unadjusted estimator only has smaller bias than the adjusted estimator when the correlation coefficients are extremely large. In Figure 7(b), we assume a = b = c = d and e = f , and compare | Bias ad j | and | Bias unad j | within the feasible region of ( a e , ) defined by (2). We can see that the adjusted estimator is superior to the unadjusted one for 71% (colored in grey in Figure 7(b)) of the feasible region. In the area satisfying | e | &gt; | a | in Figure 7(b), where the connection between M to T and Y is stronger than the other connections, the area is almost entirely grey suggesting that the adjusted estimator is preferable. This is sensible because here the confounding bias has larger magnitude than the M -Bias. In the area satisfying | a | &lt; | e | , where M -bias is stronger than confounding bias, the unadjusted estimator is superior for some values, but still tends to be inferior when the correlations are roughly the same size.

(a) Absolute biases of both estimators with a = b = c = d = e = f .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0324v1-figure-1d175990-39f9-4cc6-95c6-83dc88387b8b.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0324v1-figure-c8aea970-1691-485c-998f-6ad385bf064a.jpg)

(b) Comparison of the absolute biases with a = b = c = d and e = f . Within 71 2% (in grey) of the feasible region, the . adjusted estimator has smaller bias than the unadjusted one.

Figure 7: Butterfly-Bias

## 3 Extensions to a Binary Treatment

One might worry that the conclusions in the previous section are not applicable for a binary treatment. It turns out, however, that they are. In this section, we extend the results in Section 2 to binary treatments by representing the treatment through a latent Gaussian variable as shown in Figure 8.

Correlated Latent Variables. We extend Figure 2 to Figure 8(a). Here, T ∗ is the old T . The generating equations for T and T ∗ become

$$T = I ( T ^ { * } \geq \alpha ), \text{ and } T ^ { * } = a U + \sqrt { 1 - a ^ { 2 } } \epsilon _ { T }.$$

Other variables and noise terms remain the same. Although it might be relaxed, we make reference to the Normally assumption of the error terms for mathematical simplicity. The intercept a determines the proportion of the individuals receiving the treatment: F ( -a ) = P T ( = ) 1 , where F ( ) · is the cumulative distribution function of a standard Normal distribution. When a = 0, the number of individuals

U

a

ρ

M

T

W

d

Y

T*

(a) Correlated Hidden Causes

Figure 8: DAGs for Binary Treatment

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0324v1-figure-5da20aa4-efc5-4864-9be3-fe4a844b8ea1.jpg)

exposed to the treatment and control are balanced; when a &lt; 0, more individuals are exposed to the treatment; when a &gt; 0, the reverse.

The true causal effect of T on Y is again zero. Let f ( ) = · F ′ ( ) · and h a ( ) ≡ f a ( ) / { F ( a ) F ( -a ) } . Then Lemma 6 in Appendix shows that the unadjusted estimator has bias

$$\text{Bias} _ { u n a d j } = a d \rho \eta ( \alpha ),$$

and the adjusted estimator has bias

$$\text{Bias} _ { a d j } = \frac { a d \eta ( \alpha ) \{ \rho ( 1 - b ^ { 2 } - c ^ { 2 } - b c \rho ) - b c \} } { \rho \{ 1 - ( a b + a c \rho ) ^ { 2 } \phi ( \alpha ) \eta ( \alpha ) \} }.$$

When r = 0, the unadjusted estimator is unbiased, but the adjusted estimator has bias

$$- \frac { a b c d \eta ( \alpha ) } { 1 - ( a b ) ^ { 2 } \phi ( \alpha ) \eta ( \alpha ) }.$$

glyph[negationslash]

When r = 0, the ratio of the absolute biases is

$$\left | \frac { \text{Bias} _ { a d j } } { \text{Bias} _ { u n a d j } } \right | = \left | \frac { \rho ( 1 - b ^ { 2 } - c ^ { 2 } - b c \rho ) - b c } { \rho \{ 1 - ( a b + a c \rho ) ^ { 2 } \phi ( \alpha ) \eta ( \alpha ) \} } \right |.$$

The patterns for a binary treatment do not differ much from a continuous treatment. As before, if the correlation coefficient is moderately small, the M -Bias also tends to be small. As shown in Figure 9 (analogous to Figure 5), when | r | is comparable to | a | (= | b | = | c | = | d | ) , the adjusted estimator is less biased than the unadjusted estimator. Only when | a | is much larger than | r | is the unadjusted estimator superior.

b

c

Figure 9: M -Bias with correlated ( U W , ) and binary treatment. Compare to Figure 5

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0324v1-figure-31c63d2e-82d0-4fbd-8796-97d8e0d26f6c.jpg)

Butterfly-Bias with a Binary Treatment. We can extend the LSEM ButterflyBias setup to binary treatment just as we extended the M -Bias setup. Compare Figure 8(b) to Figure 6. T becomes T ∗ and T is built from T ∗ as above. The structural equations for T and T ∗ for butterfly bias in the binary case are then

$$T = I ( T ^ { * } \geq \alpha ), \text{ and } T ^ { * } = a U + e M + \sqrt { 1 - a ^ { 2 } - e ^ { 2 } } \varepsilon _ { T }.$$

The other equations and variables are the same as before.

Although the true causal effect of T on Y is zero, Lemma 7 in Appendix shows that the unadjusted estimator has bias

$$\text{Bias} _ { u n a d j } = ( c d e + a b f + e f ) \eta ( \alpha ),$$

and the adjusted estimator has bias

$$\text{Bias} _ { a d j } = - \frac { a b c d \eta ( \alpha ) } { 1 - ( a b + e ) ^ { 2 } \phi ( \alpha ) \eta ( \alpha ) }.$$

Therefore, the ratio of the absolute biases is

$$\left | \frac { \text{Bias} _ { a d j } } { \text{Bias} _ { u n a d j } } \right | = \left | \frac { a b c d \eta ( \alpha ) } { ( c d e + a b f + e f ) \{ 1 - ( a b + e ) ^ { 2 } \phi ( \alpha ) \eta ( \alpha ) \} } \right |.$$

Complete investigation of the ratio of the biases is intractable with seven varying parameters ( a b c d e , , , , , f , a ) . However, in the very common case with a = 0, which gives equal-sized treatment and control groups, we again find trends similar to the continuous treatment case. See Figure 10. As before, only in the cases with very small e (= f ) but large a (= b = c = ) d , does the unadjusted estimator tend to be superior. Within a reasonable region of a , these patterns are quite similar.

(a) Absolute biases with a = b = c = d = e = f .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0324v1-figure-1c0a5a1c-d414-4f0f-8309-816695a9afd7.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0324v1-figure-0c216afe-b2cd-4742-a9e1-2a937e53b2c7.jpg)

(b) Comparison of the absolute bias with a = b = c = d and e = f . Within 74 4% . (in grey) of the feasible region, the adjusted estimator is better than the unadjusted estimator.

Figure 10: Butterfly-Bias with a Binary Treatment

## 4 Illustration: The Rubin-Pearl Controversy

Pearl (2009c) cites Rubin (2007)'s example about the causal effect of smoking habits ( T ) on lung cancer ( Y ), and argues that conditioning on the pretreatment covariate 'seat-belt usage' ( M ) would introduce spurious associations, since M could be reasonably thought of as an indicator of a person's attitudes toward societal norms ( U ) as well as safety and health related measures ( W ). Assuming all the analysis is already conditioned on other observed covariates, we focus our discussion on the five variables ( U W M T Y , , , , ) , of which the dependence structure is illustrated by Figure 11. Since the patterns with a continuous treatment and a binary treatment are similar, we focus our discussion on LSEMs.

As Pearl (2009c) points out,

If we have good reasons to believe that these two types of attitudes are marginally independent, we have a pure M-structure on our hand.

In the case with r = 0, conditioning on M will lead to spurious correlation between T and Y under the null, and will bias the estimation of the causal effect of T on Y . However, Pearl (2009c) also recognizes that the independence assumption seems very strong in this example, since U and W are both background variables about the habit and personality of a person. Pearl (2009c) further argues:

But even if marginal independence does not hold precisely, conditioning on 'seat-belt usage' is likely to introduce spurious associations, hence bias, and should be approached with caution.

Although we believe most things should be approached with caution, our work, above, suggests that even mild perturbations of an M -Structure can switch which of the two approaches, conditioning or not conditioning, is likely to remove more bias. In particular, Pearl (2009c) is correct in that the adjusted estimator indeed tends to introduce more bias than the unadjusted one when an exact M -Structure holds and thus the general advice 'to condition on all observed covariates' may not be sensible in this context. However, in the example of Rubin (2007), the exact independence between a person's attitude toward societal norms U and safety and health related measures W is questionable, since we have good reasons to believe that other hidden variables such as income and family background will affect both U and W simultaneously, and thus Pearl's fears may be unfounded.

To examine this further, we consider two possible deviations from the exact M -Structure, and investigate the biases of the unadjusted and adjusted estimators for each.

- (a) (Correlated U and W ) Assume the DGP follows the DAG in Figure 11, with an additional correlation between the attitudes U and W as shown in Figure 2. If we then assume that all the correlation coefficients have the same positive magnitude, earlier results demonstrate that the adjusted estimator is preferable as it strictly dominates the unadjusted estimator except for extremely large values of the correlation coefficients.
- Furthermore, in Rubin (2007)'s example, attitudes toward societal norms U are more likely to affect the 'seat-belt usage' variable M than safety and health

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0324v1-figure-9552d8ea-1cb7-4855-be6d-3f588674ee79.jpg)

An Arrow from W to T

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0324v1-figure-5b022bd6-b194-4069-9b88-2ddf6017d0ee.jpg)

- (a) M -Bias, with possible deviations: correlated U and V , and an additional arrow from W to T
- (b) Biases of resulting unadjusted and adjusted estimators with a = b = c = d = g .

Figure 11: Sensitivity Analysis of Pearl (2009c)'s Critique on Rubin (2007)

related measures W , which further strengthens the case for adjustment. If we were willing to assume that c is zero but r is not, equation (1) in Section 2 again shows that the adjusted estimator is superior.

- (b) (An arrow from W to T ) Pearl's example seems a bit confusing on further inspection, even if we accept his independence assumption U W . In particular, one's 'attitudes towards safety and health related measures' likely impact one's decisions about smoking. Therefore, we might reasonably expect an arrow from W to T . In Figure 11(a), we remove the correlation between U and W , but we allow an arrow from W to T , i.e., the generating equation for T becomes T = aU + gW + √ 1 -a 2 -g 2 e T . Lemma 8 in Appendix gives the associated formulae for biases of the adjusted and unadjusted estimators. Figure 11(b) shows that, assuming a = b = c = d = g (i.e., equal correlations), the adjusted estimator is uniformly better.

## 5 Two Further Issues

One controversy about M -Bias is whether M -Structure is rare or not in practice, and we go through several examples to discuss this issue. In the second part of this section, we make a distinction between asymptotic and finite sample properties of M -Bias.

Is M -Structure Rare? Although Pearl (2009c) argues that M -Bias is a structural property, Rubin (2009) claims that M -bias is a rare phenomenon such as 'trying to balance a multidimensional cone on its point with no external supports in some visible directions.' As mentioned in the introduction, Gelman (2011) argues that, in social sciences, 'true zeros' are rare and consequently the independence structure in the exact M -Structure is also rare. Section 4 revisited the controversial example between Professors Pearl and Rubin, and Figure 11(a) illustrated two possible deviations from the exact M -Structure. Both deviations suggested conditioning is a superior choice. In the following, we review three other examples of M -Bias in the current literature, and investigate the plausibility of the exact M -Structure.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0324v1-figure-7e8eaa04-9f01-4976-8d25-a0b308fbba09.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0324v1-figure-6dc6c5f8-516b-4573-8c77-42f08f2bdb9d.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0324v1-figure-927d238f-88ba-48d2-85ca-6905c7ebede4.jpg)

(a) Glymour (2006)

(b) Kelcey and Carlisle (2011)

(c) Liu et al. (2012)

Figure 12: Three M -Structures with Possible Deviations in Dashed Arrows

As shown in Figure 12(a), Glymour (2006) postulates a possible M -Structure with exposure 'low income,' outcome 'diabetes,' and M variable 'mother had diabetes,' where 'family income during childhood' affects both exposure and M , and 'mother's genetic diabetes risk' affect both outcome and M . However, this M -Structure is subject to several plausible deviations: 'mother's genetic diabetes risk' may affect 'family income during childhood; 'mother had diabetes' may affect 'low education;' and 'family income during childhood' may affect 'diabetes.'

Kelcey and Carlisle (2011) have ( T Y M U W , , , , ) as 'teacher's instructional reading practice,' 'student's reading comprehension achievement,' 'teacher's reading knowledge,' 'professional development in reading,' and 'teacher's general knowledge.' See Figure 12(b). However, this M -Structure is dubious because of the possible correlation between the latent ( U W , ) and the confounding effect of M on the relationship between T and Y .

Figure 12(c) is a possible M -Structure investigated by Liu et al. (2012), where ( T Y M U W , , , , ) are 'use of selective serotonin reuptake inhibitors (SSRI),'

'lung cancer,' 'coronary artery disease,' 'depression,' and 'ever smoker.' Although it is plausible that 'coronary artery disease' is not the confounder between 'use of SSRI' and 'lung cancer,' it is very likely that 'depression' affects both 'ever smoker' and 'lung cancer.'

In summary, while all the examples above were quite useful to illustrate M -Bias in theoretical research, it is unwise to believe that these M -Structures are exact based on our background knowledge. Therefore, we suggest researchers conduct sensitivity analysis, such as illustrated earlier, according to their scientific knowledge about the structure of the DAG and the associated parameters.

Asymptotic versus Finite Sample Properties. The discussion in the previous sections are mainly based on asymptotic theory assuming large samples. As argued by Pearl (2009b), this approach allows for investigating the existence of bias in a certain DAG, and asymptotic analysis helps reveal the structural property of a DAG. A referee pointed out that the asymptotic theory is quite different from the more practical finite sample theory. In finite sample data analysis, practitioners, often interested in interval estimation and hypothesis testing, are typically more interested in whether associated confidence intervals cover the true causal parameters at nominal rates, and whether tests for null hypotheses about the causal effect have valid size. These questions are related to the asymptotic property of the DAGs, but also depend on the sample size, the procedure for constructing confidence interval, and choice of test statistic. Theoretical discussion of the finite sample theory is unfortunately more difficult. Simulation study, however, is an alternative tool for these questions. Some studies exist. In particular, Liu et al. (2012) simulate large cohort studies under an M -Structure corresponding to their science question of interest, and find that the impact of M -Bias was small for most of their 178 scenarios unless the association between M and the unmeasured confounders is very large.

## 6 Discussion

For objective causal inference, Rubin and Rosenbaum suggest balancing all the pretreatment covariate in observational studies to parallel with the design of randomized experiments (Rubin, 2007, 2008, 2009, Rosenbaum, 2002), which is called the 'pretreatment criterion' (VanderWeele and Shpitser, 2011). However, Pearl and other researchers (Pearl, 2009b,c, Shrier, 2008, 2009, Sj¨lander, 2009) criticize the o 'pretreatment criterion' by pointing out that this criterion may lead to biased inference in presence of a possible M -Structure even if the treatment assignment is unconfounded. We investigate this controversy in detail for LSEMs, ideally providing a template for future research about more general DAGs (e.g., nonparametric and nonlinear models). While we agree that Pearl's warning is very insightful, our asymptotic theory shows that, at least for LSEMs, this conclusion is quite sensitive to various deviations from the exact M -Structure, e.g., to circumstances where latent causes may be correlated or the M variable may also be a confounder between the treatment and the outcome. We also go through several candidate M -Structures in the existing literature, and find that exact M -Structure is likely to be rare with various deviations typically being more plausible. Overall, this coupled with our asymptotic theory suggests that for linear systems, except in some extreme cases, adjusting for all the pretreatment covariates is in fact a reasonable choice.

## Acknowledgment

The authors thank all the participants in the 'Causal Graphs in Low and High Dimensions' seminar at Harvard Statistics Department in Fall, 2012, and thank Professor Peter Spirtes for sending us his slides (Sprites, 2002). Comments from the associate editor and two reviewers greatly improved the quality of our paper.

## Appendix: Lemmas and Proofs

Lemma 1 In the linear regression model Y = b 0 + b T T + b MM + e with e ( T M , ) and E ( e ) = 0 , we have

$$\beta _ { T } \ = \ \frac { C o v ( Y, T ) V a r ( M ) - C o v ( Y, M ) C o v ( M, T ) } { V a r ( T ) V a r ( M ) - C o v ^ { 2 } ( M, T ) }.$$

Proof. Solve for ( b T , b M ) using the following moment conditions

$$\begin{cases} \text{ Cov} ( Y, T ) & = & \beta _ { T } \text{Var} ( T ) + \beta _ { M } \text{Cov} ( M, T ), \\ \text{ Cov} ( Y, M ) & = & \beta _ { T } \text{Cov} ( M, T ) + \beta _ { M } \text{Var} ( M ). \end{cases}$$

Lemma 2 Under the model generated by Figure 2, the regression coefficient of T by regressing Y onto ( T M , ) is

$$\beta _ { T } = \frac { a d \rho ( 1 - b ^ { 2 } - c ^ { 2 } - b c \rho ) - a b c d } { 1 - ( a b + a c \rho ) ^ { 2 } }.$$

Proof. We apply Lemma 1, where all variance terms such as Var ( M ) are 1, and the covariance terms are easily calculated. For example, we have Cov ( Y T , ) = Cov ( dW aU , ) = ad r , and Cov ( Y M , ) = Cov ( dW bU , + cW ) = bd r + cd .

Lemma 3 Under the model generate by Figure 7, the regression coefficient of T from regressing Y onto ( T M , ) is

$$\beta _ { T } = - \frac { a b c d } { 1 - ( a b + e ) ^ { 2 } }.$$

Proof. Similar to the proof of Lemma 2, we apply Lemma 1.

Lemma 4 Assume that ( X 1 , X 2 ) follows a bivariate Normal distribution with means zero, variances one, and correlation coefficient r . Then E ( X 1 | X 2 ≥ z ) -E ( X 1 | X 2 &lt; z ) = r h ( ) z , where h ( ) = z f ( ) z / { F ( ) z F ( - } z ) .

Proof. Since X 1 = rX 2 + √ 1 -r 2 Z with Z ∼ N ( 0 1 , ) and Z X 2, we have

$$\mathbb { E } ( X _ { 1 } \ | \ X _ { 2 } \geq z ) = r \mathbb { E } ( X _ { 2 } \ | \ X _ { 2 } \geq z ) = \frac { r } { \Phi ( - z ) } \int _ { z } ^ { \infty } x \phi ( x ) d x = - \frac { r } { \Phi ( - z ) } \int _ { z } ^ { \infty } d \phi ( x ) = r \frac { \phi ( z ) } { \Phi ( - z ) }.$$

Similarly, we have E ( X 1 | X 2 &lt; z ) = E ( X 1 | -X 2 &gt; -z ) = -r f ( -z ) / F ( ) = z -r f ( ) z / F ( ) z . Therefore, E ( X 1 | X 2 ≥ z ) -E ( X 1 | X 2 &lt; z ) = r f ( ) z { 1 / F ( -z ) + 1 / F ( ) z } = r h ( ) z .

Lemma 5 The covariance between X and B ∼ Bernoulli ( p ) is

$$C o v ( X, B ) = p ( 1 - p ) \{ \mathbb { E } ( X \ | \ B = 1 ) - \mathbb { E } ( X \ | \ B = 0 ) \}.$$

Proof. It follows from the definition of the covariance.

Lemma 6 Under the model generated by Figure 8(a), the regression coefficient of T from regressing Y onto ( T M , ) is

$$\beta _ { T } = \frac { a d \eta ( \alpha ) \{ \rho ( 1 - b ^ { 2 } - c ^ { 2 } - b c \rho ) - b c \} } { \rho \{ 1 - ( a b + a c \rho ) ^ { 2 } \phi ( \alpha ) \eta ( \alpha ) \} }.$$

Proof. We have the following joint Normality of ( Y M T , , ∗ ) :

$$\begin{pmatrix} Y \\ M \\ T ^ { * } \end{pmatrix} \ \sim \ N _ { 3 } \left \{ \begin{pmatrix} 0 \\ 0 \\ 0 \end{pmatrix}, \begin{pmatrix} 1 & b d \rho + c d & a d \rho \\ b d \rho + c d & 1 & a b + a c \rho \\ a d \rho & a b + a c \rho & 1 \end{pmatrix} \right \}.$$

From Lemma 4, we have

$$\mathbb { E } ( M \, | \, T = 1 ) - \mathbb { E } ( M \, | \, T = 0 ) \ & = \ \mathbb { E } ( M \, | \, T ^ { * } \geq \alpha ) - \mathbb { E } ( M \, | \, T ^ { * } < \alpha ) = ( a b + a c \rho ) \eta ( \alpha ), \\ \mathbb { E } ( Y \, | \, T = 1 ) - \mathbb { E } ( Y \, | \, T = 0 ) \ & = \ \mathbb { E } ( Y \, | \, T ^ { * } \geq \alpha ) - \mathbb { E } ( Y \, | \, T ^ { * } < \alpha ) = a d \rho \eta ( \alpha ).$$

Therefore, from Lemma 5, the covariances are Cov ( M T , ) = F ( a ) F ( -a )( ab + ac r h a ) ( ) , and Cov ( Y T , ) = F ( a ) F ( -a ) ad rh a ( ) . According to Lemma 1, the regression coefficient b T is

$$\begin{array} { r c l } \supset & \cdot \cdot \\ \beta _ { T } & = & \frac { \Phi ( \alpha ) \Phi ( - \alpha ) a d \rho \eta ( \alpha ) - ( b d \rho + c d ) \Phi ( \alpha ) \Phi ( - \alpha ) ( a b + a c \rho ) \eta ( \alpha ) } { \Phi ( \alpha ) \Phi ( - \alpha ) - \Phi ^ { 2 } ( \alpha ) \Phi ^ { 2 } ( - \alpha ) ( a b + a c \rho ) ^ { 2 } \eta ^ { 2 } ( \alpha ) } \\ & = & \frac { a d \eta ( \alpha ) \{ \rho ( 1 - b ^ { 2 } - c ^ { 2 } - b c \rho ) - b c \} } { \rho \{ 1 - ( a b + a c \rho ) ^ { 2 } \phi ( \alpha ) \eta ( \alpha ) \} }. \end{array}$$

Lemma 7 Under the model generated by Figure 8(b), the regression coefficient of T from regressing Y onto ( T M , ) is

$$\beta _ { T } = - \frac { a b c d \eta ( \alpha ) } { 1 - ( a b + e ) \phi ( \alpha ) }.$$

Proof. We have the following joint Normality of ( Y M T , , ∗ ) :

$$\begin{pmatrix} Y \\ M \\ T ^ { * } \end{pmatrix} \ \sim \ N _ { 3 } \left \{ \begin{pmatrix} 0 \\ 0 \\ 0 \end{pmatrix}, \begin{pmatrix} 1 & c d + f & c d e + a b f + e f \\ c d + f & 1 & a b + e \end{pmatrix} \right \}.$$

From Lemma 4, we have

E ( M T | = ) 1 -E ( M T | = ) 0 = E ( M T | ≥ a ) -E ( M T | &lt; a ) = ( ab + ) e h a ( )

$$\mathbb { E } ( Y \, | \, T = 1 ) - \mathbb { E } ( Y \, | \, T = 0 ) \ \ = \ \mathbb { E } ( Y \, | \, T ^ { * } \geq \alpha ) - \mathbb { E } ( Y \, | \, T ^ { * } < \alpha ) = ( c d e + a b f + e f ) \eta ( \alpha ).$$

$$& \mathbb { E } ( M \, | \, T ^ { * } \geq \alpha ) - \mathbb { E } ( M \, | \, T ^ { * } < \alpha ) = ( a b + e ) \eta ( \alpha ), \\ & \mathbb { E } ( Y \, | \, T ^ { * } \geq \alpha ) - \mathbb { E } ( Y \, | \, T ^ { * } < \alpha ) = ( c d e + a b f + e f ) \eta ( \alpha ).$$

From Lemma 5, we obtain their covariances Cov ( M T , ) = F ( a ) F ( -a )( ab + ) e h a ( ) , and Cov ( Y T , ) = F ( a ) F ( -a )( cde + abf + e f ) h a ( ) . According to Lemma 1, the regression coefficient b T is

$$\begin{array} { r c l } \Lambda & = & \cdot \cdot \\ \beta _ { T } & = & \frac { \Phi ( \alpha ) \Phi ( - \alpha ) ( c d + f ) \Phi ( \alpha ) \Phi ( - \alpha ) ( a b + e ) \eta ( \alpha ) } { \Phi ( \alpha ) \Phi ( - \alpha ) - \Phi ^ { 2 } ( \alpha ) \Phi ^ { 2 } ( - \alpha ) ( a b + e ) ^ { 2 } \eta ^ { 2 } ( \alpha ) } \\ & = & - \frac { a b c d \eta ( \alpha ) } { 1 - ( a b + e ) ^ { 2 } \phi ( \alpha ) \eta ( \alpha ) }. \end{array}$$

Lemma 8 Under the model generated by Figure 11(a) with an arrow from W to T, the unadjusted estimator has bias ad r + dg, and the adjusted estimator has bias

$$\frac { d g - ( c d ) ( a b + c g ) } { 1 - ( a b + c g ) ^ { 2 } }.$$

Proof. The unadjusted estimator is Cov ( T Y , ) = ad r + dg . Expanding Lemma 1 gives the above as the regression coefficient of T from regressing Y onto ( T M , ) .

## References

Berkson, J. (1946): 'Limitations of the application of fourfold table analysis to hospital data,' Biometrics Bulletin , 2, 47-53.

Copas, J. B. and H. G. Li (1997): 'Inference for non-random samples (with discussion),' Journal of the Royal Statistical Society: Series B , 59, 55-95.

- Gelman, A. (2011): 'Causality and statistical learning,' American Journal of Sociology , 117, 955-966.
- Glymour, M. M. (2006): 'Using causal diagrams to understand common problems in social epidemiology,' In Methods in Social Epidemiology, Oakes M, and Kaufman J, eds. Jossey-Bass: San Francisco, CA , 393-428.

Greenland, S. (2002): 'Quantifying biases in causal models: classical confounding vs collider-stratification bias,' Epidemiology , 14, 300-306.

Heckman, J. J. (1979): 'Sample selection bias as a specification error,' Econometrica , 153-161.

Hern´ an, M. A., S. Hernandez-Diaz, and J. M. Robins (2004): 'A structural approach to selection bias,' Epidemiology , 15, 615-625.

Kelcey, B. and J. Carlisle (2011): 'The threshold of embedded M collider bias and confounding bias.' Society for Research on Educational Effectiveness Conference , available at http://files.eric.ed.gov/fulltext/ED519118.pdf .

Liu, W., M. A. Brookhart, S. Schneeweiss, X. Mi, and S. Setoguchi (2012): 'Implications of M bias in epidemiologic studies: a simulation study,' American Journal of Epidemiology , 176, 938-948.

- Neyman, J. (1923/1990): 'On the application of probability theory to agricultural experiments. essay on principles. section 9,' Statistical Science , 5, 465-472.
- Pearl, J. (1995): 'Causal diagrams for empirical research,' Biometrika , 82, 669688.

Pearl, J. (2009a): Causality: Models, Reasoning and Inference, 2nd Edition , Cambridge University Press.

Pearl, J. (2009b): 'Letter to the editor,' Statistics in Medicine , 28, 1415-1416.

- Pearl, J. (2009c): 'Myth, confusion, and science in causal analysis,' Technical Report , available at http://ftp.cs.ucla.edu/pub/stat\_ser/r348.pdf .
- Pearl, J. (2013): 'Linear models: A useful microscope for causal analysis,' Journal
- of Causal Inference , 1, 155-170.

Rosenbaum, P. R. (2002): Observational Studies , New York: Springer.

Rosenbaum, P. R. and D. B. Rubin (1983): 'The central role of the propensity score in observational studies for causal effects,' Biometrika , 70, 41-55.

Rubin, D. B. (1974): 'Estimating causal effects of treatments in randomized and nonrandomized studies,' Journal of Educational Psychology , 66, 688-701.

Rubin, D. B. (2007): 'The design versus the analysis of observational studies for causal effects: parallels with the design of randomized trials,' Statistics in Medicine , 26, 20-36.

Rubin, D. B. (2008): 'For objective causal inference, design trumps analysis,' The Annals of Applied Statistics , 2, 808-840.

Rubin, D. B. (2009): 'Should observational studies be designed to allow lack of balance in covariate distributions across treatment groups?' Statistics in Medicine , 28, 1420-1423.

Shrier, I. (2008): 'Letter to the editior,' Statistics in Medicine , 27, 2740-2741.

- Shrier, I. (2009): 'Propensity scores,' Statistics in Medicine , 28, 1315-1318.
- Sj¨ olander, A. (2009): 'Propensity scores and M-structures,' Statistics in Medicine , 28, 1416-1420.

Sprites, P. (2002): 'Presented at: WNAR/IMS Meeting. Los Angeles, CA, June 2002,' .

VanderWeele, T. J. and I. Shpitser (2011): 'A new criterion for confounder selection,' Biometrics , 67, 1406-1413.

- Wright, S. (1921): 'Correlation and causation,' Journal of Agricultural Research , 20, 557-585.
- Wright, S. (1934): 'The method of path coefficients,' The Annals of Mathematical Statistics , 5, 161-215.
Peng Ding Luke Miratrix
To Adjust or Not to Adjust? Sensitivity Analysis of M-Bias and Butterfly-Bias
"M-Bias," as it is called in the epidemiologic literature, is the bias
introduced by conditioning on a pretreatment covariate due to a particular
"M-Structure" between two latent factors, an observed treatment, an outcome,
and a "collider." This potential source of bias, which can occur even when the
treatment and the outcome are not confounded, has been a source of considerable
controversy. We here present formulae for identifying under which circumstances
biases are inflated or reduced. In particular, we show that the magnitude of
M-Bias in linear structural equation models tends to be relatively small
compared to confounding bias, suggesting that it is generally not a serious
concern in many applied settings. These theoretical results are consistent with
recent empirical findings from simulation studies. We also generalize the
M-Bias setting (1) to allow for the correlation between the latent factors to
be nonzero, and (2) to allow for the collider to be a confounder between the
treatment and the outcome. These results demonstrate that mild deviations from
the M-Structure tend to increase confounding bias more rapidly than M-Bias,
suggesting that choosing to condition on any given covariate is generally the
superior choice. As an application, we re-examine a controversial example
between Professors Donald Rubin and Judea Pearl.