Tutorials › Marginal Structural Models › Why Outcome-Blind Diagnostics Have Limits

Diagnostics and Sensitivity to Misspecification · Tutorial 285 of 501

Why Outcome-Blind Diagnostics Have Limits

Learn what treatment-model and balance diagnostics can establish, what causal assumptions they cannot verify, and how to interpret apparently reassuring results.

Advanced 13 min read

What You'll Learn

  • Distinguish observed treatment-model fit from the causal assumptions needed for identification.
  • Derive the population covariate-balance identity under a correct propensity score.
  • Construct observationally indistinguishable causal models with different treatment effects.
  • Identify ways summary balance can conceal important distributional differences.
  • Explain why treatment versions and unmeasured confounding are not ruled out by outcome-blind checks.
  • Compare balance diagnostics with outcome-informed checks and sensitivity analysis.

What Can an Outcome-Blind Diagnostic Establish?

Treatment-model diagnostics and covariate-balance checks are useful because they can reveal problems in the observed treatment and covariate data. They can show, for example, that estimated treatment probabilities are extreme or that selected measured covariates remain imbalanced after weighting. But these checks use treatment and covariate information, not the counterfactual outcomes that define a causal effect. Their scope is therefore limited: reassuring diagnostics do not establish that the causal assumptions required by a marginal structural model hold.

To make that distinction precise, consider a single binary treatment \(A\), measured pretreatment covariates \(L\), and potential outcomes \(Y(1)\) and \(Y(0)\). The average treatment effect is \(E[Y(1)-Y(0)]\). The standard identification assumptions include consistency, conditional exchangeability, and positivity. Consistency says that the observed outcome equals the potential outcome under the treatment actually received. Conditional exchangeability says that \((Y(1),Y(0))\) is independent of \(A\) given \(L\). Positivity says that \(0<P(A=1\mid L=l)<1\) for covariate values in the target population.

An outcome-blind diagnostic is any check calculated without using \(Y\) or potential outcomes. It may assess a treatment model or compare weighted distributions of observed covariates. Such a check can provide evidence about those observed quantities. It cannot, by itself, confirm that \(L\) contains all common causes of treatment and outcome, that treatment is well-defined, or that the causal assumptions hold.

Why Correct Propensities Balance Covariates

Write \(e(L)=P(A=1\mid L)\) for the true propensity score. The identity below explains why inverse-probability weights can balance measured covariates when the propensity score is correct. It also clarifies what a balance check examines: expectations of functions of the measured \(L\), not the unobserved causal assumptions.

Result: Suppose \(0<e(L)<1\) almost surely and \(g(L)\) is integrable. Then
$$ E\left[\frac{A g(L)}{e(L)}\right]=E[g(L)] \quad\text{and}\quad E\left[\frac{(1-A)g(L)}{1-e(L)}\right]=E[g(L)]. $$

Proof. By iterated expectation and the definition of \(e(L)\),

$$ E\left[\frac{A g(L)}{e(L)}\right] =E\left[E\left\{\left.\frac{A g(L)}{e(L)}\right|L\right\}\right] =E\left[\frac{g(L)}{e(L)}E(A\mid L)\right] =E[g(L)]. $$

The second identity follows in the same way, since \(E(1-A\mid L)=1-e(L)\):

$$ E\left[\frac{(1-A)g(L)}{1-e(L)}\right] =E\left[\frac{g(L)}{1-e(L)}E(1-A\mid L)\right] =E[g(L)]. $$

Integrability ensures that these expectations exist, and positivity ensures that the denominators are nonzero. Taking \(g(L)\) to be individual covariates, or selected functions such as squares and interactions, gives corresponding population balance statements. This result is about a measured-data distribution. It does not say that balanced groups have exchangeable potential outcomes, nor that all functions of \(L\) have been checked. \(\square\)

In practice, estimated weights and finite samples replace the population quantities in the identity. A diagnostic therefore reports evidence about a fitted procedure and a set of observed summaries. It is not a direct test of the identity’s causal interpretation or a proof that the true propensity score was used.

Worked Example: Perfect Measured Balance With Hidden Confounding

Let \(U\) be an unmeasured binary variable with \(P(U=1)=0.5\). Treatment is assigned with \(P(A=1\mid U=1)=0.9\) and \(P(A=1\mid U=0)=0.1\). Let measured \(L\) be an independent Bernoulli variable with probability \(0.5\), and let \(Y(1)=Y(0)=U\).

First calculate the marginal treatment probability: \(P(A=1)=0.5(0.9)+0.5(0.1)=0.5\). Since \(L\) is independent of \(A\), treatment groups have the same distribution of \(L\); weighting is not needed to achieve measured balance. The marginal treatment probability is also bounded away from zero and one.

Now calculate the observed outcome means. Among treated individuals, \(P(U=1\mid A=1)=0.5(0.9)/0.5=0.9\), so \(E(Y\mid A=1)=0.9\). Among untreated individuals, \(P(U=1\mid A=0)=0.5(0.1)/0.5=0.1\), so \(E(Y\mid A=0)=0.1\). The observed difference is \(0.8\).

But \(Y(1)=Y(0)=U\), so the individual treatment effect and the average treatment effect are both zero. The apparent association is entirely due to \(U\), which affects treatment and outcome but is absent from the balance check. The example does not show that balance is useless; it shows that balance in measured \(L\) does not establish exchangeability.

Different Causal Worlds Can Look the Same to Diagnostics

The limitation is deeper than an inability to measure a particular confounder. Even the complete observed distribution of treatment and covariates can be compatible with different causal effects. An outcome-blind diagnostic, which uses less information than the complete observed distribution, cannot distinguish such cases.

Result: There can be two causal data-generating mechanisms with the same distribution of \((A,L)\) but different average treatment effects. Consequently, no diagnostic that depends only on \((A,L)\) can determine which effect is correct.

Construction. Let \(L\) be independent of \(A\) in both mechanisms, with \(P(A=1)=0.5\). In the first mechanism, treatment is randomized, and the potential outcomes have \(P(Y(1)=1)=0.9\) and \(P(Y(0)=1)=0.1\). The average treatment effect is \(0.9-0.1=0.8\). In the second mechanism, let \(U\) be Bernoulli with probability \(0.5\), let \(Y(1)=Y(0)=U\), and set \(A=U\) with probability \(0.9\) and \(A=1-U\) with probability \(0.1\). Its average treatment effect is zero. Both mechanisms have the same distribution of \((A,L)\), namely \(A\) is Bernoulli with probability \(0.5\) and independent of \(L\). Thus every diagnostic computed only from \((A,L)\) has the same distribution under both mechanisms, although their causal effects differ. \(\square\)

The construction is not a claim that all treatment-model diagnostics are uninformative. It establishes a boundary on what they can verify: causal identification requires assumptions linking treatment assignment to potential outcomes, and those assumptions are not determined by the distribution of treatment and measured covariates alone.

Worked Example: The Same Observed Treatment–Outcome Table, Different Effects

Here is a stronger illustration: even the observed treatment–outcome data can agree across causal mechanisms. In the first mechanism, let \(A\) be randomized with probability \(0.5\), and let the binary potential outcomes satisfy \(P(Y(1)=1)=0.9\) and \(P(Y(0)=1)=0.1\). The average treatment effect is \(0.8\). The observed joint probabilities are \(P(A=1,Y=1)=0.45\), \(P(A=1,Y=0)=0.05\), \(P(A=0,Y=1)=0.05\), and \(P(A=0,Y=0)=0.45\).

In the second mechanism, let \(U\) be Bernoulli with probability \(0.5\), set \(Y(1)=Y(0)=U\), and let \(A=U\) with probability \(0.9\), flipping that assignment with probability \(0.1\). Then \(P(A=1,Y=1)=P(U=1,A=1)=0.5(0.9)=0.45\); \(P(A=1,Y=0)=0.5(0.1)=0.05\); \(P(A=0,Y=1)=0.5(0.1)=0.05\); and \(P(A=0,Y=0)=0.5(0.9)=0.45\). These are exactly the same observed probabilities as in the first mechanism.

The second mechanism has zero causal effect because the potential outcomes are equal for each individual. Treatment depends on \(U\), so exchangeability fails if \(U\) is omitted. Both mechanisms also have \(P(A=1)=0.5\). Adding an independent measured \(L\) preserves the same treatment frequencies and measured balance in both. Thus even agreement of the full observed table cannot, without assumptions, distinguish the causal effects; treatment-and-covariate diagnostics alone have still less information.

Balance Summaries Can Conceal Distributional Differences

Balance is always balance with respect to something: a chosen set of functions, summaries, or distances between distributions. Equality of means does not imply equality of distributions. A check that reports only mean differences can therefore miss differences in spread, tails, or shape. More extensive balance checks can examine more features, but any finite collection of summaries still describes only what was checked.

Worked Example: Equal Means, Unequal Distributions

Suppose two weighted treatment groups have the following distributions of a covariate \(L\):

\(L\)Treated-group probabilityControl-group probability
-10.250.40
00.500.20
10.250.40

The treated mean is \(0.25(-1)+0.50(0)+0.25(1)=0\). The control mean is \(0.40(-1)+0.20(0)+0.40(1)=0\). A mean-balance diagnostic therefore reports exact equality. The treated second moment is \(0.25(1)+0.50(0)+0.25(1)=0.5\), whereas the control second moment is \(0.40(1)+0.20(0)+0.40(1)=0.8\). Since both means are zero, these are also the variances. The distributions are not equal, despite exact mean balance.

This example concerns the information in a particular summary, not a claim that a specific weighting procedure must produce these probabilities. It demonstrates why diagnostic conclusions should name the features that were checked rather than treating one balance statistic as proof of full distributional equality.

Other Assumptions Outside the Check

Treatment probabilities and covariate balance also cannot verify that the treatment is sufficiently well-defined for consistency to apply. Suppose \(A=1\) denotes “received the intervention,” but that label includes two versions, such as two substantially different doses. If outcomes differ by version, the notation \(Y(1)\) does not identify a unique outcome unless the intervention or version distribution is specified. The frequencies of \(A\) and the balance of \(L\) do not resolve that ambiguity.

Positivity diagnostics have a related boundary. They can reveal observed regions with very few treated or untreated individuals, and extreme estimated weights may warn that comparisons rely on limited data. But empirical overlap in measured \(L\) does not establish positivity conditional on unmeasured variables. Nor does a finite sample with observations in both groups guarantee adequate support for every target contrast.

In longitudinal settings, the same logic applies at each treatment time. Treatment-model checks and weighted balance summaries address recorded treatment and covariate histories. They do not establish that the observed histories contain all relevant confounders, that the sequential exchangeability assumptions hold, or that the treatment regimes are well-defined. The history can be thoroughly inspected and still omit a cause of both later treatment and outcome.

Using Diagnostics Without Overclaiming

The useful interpretation is conditional and specific: “These checks found no serious problem in the treatment probabilities or the measured features we examined.” That is a meaningful statement about the analysis. It is not equivalent to “the causal assumptions have been verified.” The distinction also explains why diagnostics serve complementary roles rather than replacing one another.

  • Treatment-model and balance checks examine estimated treatment probabilities and selected measured covariate distributions.
  • Outcome-informed checks, including the standard approach of negative-control outcomes as model checks, can reveal some discrepancies that outcome-blind checks cannot address. A reassuring result still does not prove exchangeability or consistency.
  • Sensitivity analysis asks how conclusions could change under specified departures from assumptions, such as unmeasured confounding. It quantifies dependence on assumptions rather than verifying them from the observed data.

A common mistake is to reason, “The propensity scores overlap and all measured covariates are balanced, so there is no confounding.” The first clause concerns observed treatment support; the second concerns measured covariates and chosen summaries. Neither entails the third, which is a claim about all common causes of treatment and potential outcomes. A better conclusion is that the diagnostics support the adequacy of the fitted treatment model for the observed features assessed, while causal interpretation remains dependent on substantive assumptions and other evidence.

Key takeaway: Outcome-blind diagnostics are valuable checks on observed treatment assignment and measured covariates. They cannot, on their own, verify exchangeability, consistency, or positivity in unobserved dimensions. Report what was checked, what was found, and which causal assumptions remain unverified.

Check Your Understanding

Answer each question, then compare your reasoning with the solutions.

  1. If inverse-probability weighting exactly balances the mean of every recorded covariate, which causal assumption does that fact still fail to establish?
  2. In the hidden-confounding example, calculate \(P(U=1\mid A=0)\) and the average treatment effect.
  3. For the two distributions in the equal-means example, give the variance in each group and explain why a mean-only check misses their difference.
  4. Why can two causal mechanisms with the same distribution of \((A,L)\) not be distinguished by a diagnostic computed only from \((A,L)\)?

Solutions. (1) It does not establish conditional exchangeability: an unmeasured common cause may remain. It also does not prove consistency or positivity beyond the observed dimensions checked. (2) \(P(U=1\mid A=0)=P(U=1,A=0)/P(A=0)=0.5(0.1)/0.5=0.1\). Since \(Y(1)=Y(0)=U\), \(E[Y(1)-Y(0)]=0\). (3) The means are both zero, so each variance equals its second moment: \(0.5\) in the treated group and \(0.8\) in the control group. A mean-only check compares the same first moment and does not compare the second moments. (4) The diagnostic’s input distribution is identical under the two mechanisms, so any function of \((A,L)\) has the same distribution under both. It contains no information that selects between their different causal effects.

Clinical Biostats