Tutorials › Propensity Score Matching and Instrumental Variables › Worst-Case Bounds For A Binary Outcome

Partial Identification and Bounds · Tutorial 402 of 699

Worst-Case Bounds For A Binary Outcome

You will derive no-assumption bounds for a binary-outcome average treatment effect and learn how to calculate and interpret them.

Advanced 10 min read

What You'll Learn

  • Define the binary-outcome average treatment effect using potential outcomes.
  • Derive bounds for each average potential outcome from the observed treatment groups.
  • Combine those bounds to obtain lower and upper bounds for the average treatment effect.
  • Interpret the bounds in risk-difference units and understand how treatment prevalence affects them.
  • Recognize why observed outcome means do not identify causal effects without additional assumptions.
  • Distinguish population bounds from estimates and confidence intervals based on sample data.

Why Bound the Effect?

When outcomes are binary, the average treatment effect is a difference in probabilities: the average outcome under treatment minus the average outcome without treatment. But for each person, only one of those outcomes is observed. If treatment assignment may depend on the potential outcomes, the observed treatment-group means need not equal the population averages under treatment and no treatment. A useful first step is therefore to ask what the data alone can establish, without assuming that treatment assignment is unconfounded.

The answer is an interval rather than, in general, a single value. The interval uses the treatment prevalence and the observed outcome means in the two treatment groups. It allows the unobserved potential outcomes to take any value consistent with being binary. No random assignment, conditional exchangeability, or instrument is required for these bounds.

Definition: Let \(D\in\{0,1\}\) indicate treatment, and let \(Y(1)\in\{0,1\}\) and \(Y(0)\in\{0,1\}\) denote the potential outcomes under treatment and no treatment. The average treatment effect (ATE) in the target population is \(E[Y(1)-Y(0)]\). Define the treatment prevalence \(p=P(D=1)\), the observed treated-group mean \(m_1=E[Y\mid D=1]\), and the observed untreated-group mean \(m_0=E[Y\mid D=0]\).

The outcome value 1 can represent an event, such as recovery, readmission, or disease. The ATE is then a difference in event probabilities, measured in risk-difference units. For example, an ATE of 0.05 represents an average increase of 5 percentage points in the probability of the outcome under treatment.

For the formulas below, assume treatment is well defined, the observed outcome equals the potential outcome under the treatment actually received, and there is no interference between individuals. These are consistency and no-interference conditions. Also assume \(0<p<1\), so both observed group means are defined. The bounds make no assumption about how treatment assignment relates to either potential outcome, and no assumption about the association between \(Y(1)\) and \(Y(0)\).

Deriving the Bounds

Consistency gives \(Y=Y(1)\) for treated individuals and \(Y=Y(0)\) for untreated individuals. Thus \(m_1\) reveals the mean of \(Y(1)\) among the treated, but it does not reveal the mean of \(Y(1)\) among the untreated. Since the missing potential outcome is binary, its mean must be between 0 and 1. Apply the law of total expectation to \(Y(1)\):

$$ E[Y(1)] = pE[Y(1)\mid D=1]+(1-p)E[Y(1)\mid D=0] = pm_1+(1-p)E[Y(1)\mid D=0]. $$

The unknown conditional mean \(E[Y(1)\mid D=0]\) lies in \([0,1]\). Multiplying by the nonnegative weight \(1-p\) and adding \(pm_1\) gives the bounds for the population mean under treatment. The same reasoning applies to \(Y(0)\): its mean among untreated individuals is observed as \(m_0\), while its mean among treated individuals is unknown but lies in \([0,1]\).

$$ pm_1 \leq E[Y(1)] \leq pm_1+(1-p), \qquad (1-p)m_0 \leq E[Y(0)] \leq (1-p)m_0+p. $$

To bound the ATE, subtract the mean under no treatment from the mean under treatment. A lower bound results by taking the smallest allowed value of \(E[Y(1)]\) and subtracting the largest allowed value of \(E[Y(0)]\). An upper bound uses the largest allowed value of \(E[Y(1)]\) and the smallest allowed value of \(E[Y(0)]\). Therefore every ATE compatible with these restrictions must satisfy:

Result: With a binary outcome, consistency, no interference, and \(0<p<1\), the no-assumption bounds based on \(p\), \(m_1\), and \(m_0\) are
$$ pm_1-\bigl((1-p)m_0+p\bigr) \leq E[Y(1)-Y(0)] \leq \bigl(pm_1+(1-p)\bigr)-(1-p)m_0. $$
Equivalently, the lower endpoint is \(pm_1-(1-p)m_0-p\), and the upper endpoint is \(pm_1+(1-p)-(1-p)m_0\).

This derivation gives bounds that follow from the observed means and the binary range of the potential outcomes. It does not assume that the observed mean difference \(m_1-m_0\) is the ATE. That equality requires additional assumptions connecting treatment assignment to potential outcomes.

Worked Example: A Treatment Prevalence of 40 Percent

Suppose a hypothetical population has treatment prevalence \(p=0.40\). The observed outcome mean is \(m_1=0.30\) among treated individuals and \(m_0=0.20\) among untreated individuals. Find the no-assumption ATE bounds.

The lower endpoint combines the smallest possible treated potential-outcome mean with the largest possible untreated potential-outcome mean:

$$ 0.40(0.30)-\left[0.60(0.20)+0.40\right] =0.12-(0.12+0.40) =-0.40. $$

The upper endpoint combines the largest possible treated potential-outcome mean with the smallest possible untreated potential-outcome mean:

$$ \left[0.40(0.30)+0.60\right]-0.60(0.20) =(0.12+0.60)-0.12 =0.60. $$

The ATE is bounded by \([-0.40,0.60]\), or from a 40-percentage-point decrease to a 60-percentage-point increase in the probability of the outcome. The observed difference is \(m_1-m_0=0.10\), but that 10-percentage-point difference is not, by itself, the ATE. Without further assumptions, the interval includes negative, zero, and positive effects.

Understanding the Interval

The width of the interval has a simple form. Subtract the lower endpoint from the upper endpoint:

$$ \begin{aligned} &\left[pm_1+(1-p)-(1-p)m_0\right] -\left[pm_1-(1-p)m_0-p\right]\\ &=(1-p)+p=1. \end{aligned} $$

Thus, under these no-assumption bounds, the interval has width one risk-difference unit, regardless of the observed means or treatment prevalence. The prevalence changes where the interval sits, but not its width. The endpoints remain within the possible ATE range \([-1,1]\): the lower endpoint is at least \(-1\), and the upper endpoint is at most 1, because \(0\leq m_1,m_0\leq1\).

This fixed width helps explain why the bounds may be informative about the possible scale of the effect without resolving its direction. If the interval includes both negative and positive values, the observed data plus the binary outcome restriction do not determine whether treatment increases or decreases the outcome on average. More restrictions on the relationship between treatment and potential outcomes may narrow the interval, but they are substantive assumptions, not consequences of the observed group means alone.

Worked Example: A High Treatment Prevalence

In a second hypothetical population, \(p=0.75\), \(m_1=0.60\), and \(m_0=0.25\). Calculate the endpoints and interpret them.

For the lower endpoint, the treated contribution is \(pm_1=0.75(0.60)=0.45\). The largest possible contribution from \(Y(0)\) among treated individuals is \(p=0.75\), in addition to the observed untreated contribution \((1-p)m_0=0.25(0.25)=0.0625\). Hence:

$$ 0.45-(0.0625+0.75)=-0.3625. $$

For the upper endpoint, the largest possible contribution to \(E[Y(1)]\) is \(0.45+0.25=0.70\), and the smallest possible \(E[Y(0)]\) is \(0.0625\). Therefore:

$$ 0.70-0.0625=0.6375. $$

The ATE interval is \([-0.3625,0.6375]\), with endpoints approximately a 36.3-percentage-point decrease and a 63.8-percentage-point increase. Its width is \(0.6375-(-0.3625)=1.000\), as predicted. Even though three-quarters of the population is treated, a high treatment prevalence alone does not establish that the treated and untreated groups are comparable.

Worked Example: Equal Observed Means Do Not Imply a Zero ATE

Suppose a hypothetical population has \(p=0.20\), and both observed group means equal 0.50: \(m_1=m_0=0.50\). The observed difference is zero. Calculate the bounds.

The lower endpoint is:

$$ 0.20(0.50)-\left[0.80(0.50)+0.20\right] =0.10-(0.40+0.20) =-0.50. $$

The upper endpoint is:

$$ \left[0.20(0.50)+0.80\right]-0.80(0.50) =(0.10+0.80)-0.40 =0.50. $$

The bounds are \([-0.50,0.50]\), even though the observed means agree exactly. The equality of the two observed means is a statement about outcomes in the groups that received different treatments. It does not tell us what those same individuals' unobserved potential outcomes would have been under the other treatment.

Assumptions, Diagnostics, and Alternatives

The bounds depend on a small set of conditions. First, the outcome must genuinely be binary on the stated scale; if it is not restricted to 0 and 1, the missing-potential-outcome means need different bounds. Second, the treatment and outcome must be defined consistently so that the observed outcome corresponds to the potential outcome under the received treatment. Third, both treatment groups must be represented to use \(m_1\) and \(m_0\) as defined above. These are conditions on the setup and observed data, not tests that prove the bounds' causal interpretation.

Several checks can catch errors in applying the formulas. Confirm that \(p\), \(m_1\), and \(m_0\) are all proportions on the same population and that both observed means lie in \([0,1]\). Check that the lower endpoint does not exceed the upper endpoint and that both endpoints lie within \([-1,1]\). A result outside those ranges usually indicates a substitution or coding error. Also report whether an endpoint is a risk difference, not a relative risk or an odds ratio.

If the treatment prevalence is 0 or 1, one of the observed group means is undefined, so the displayed two-group formula cannot be used as written. The observed arm still identifies its own potential-outcome mean. For \(p=0\), if \(m_0=E[Y\mid D=0]\), then \(E[Y(0)]=m_0\) and \(E[Y(1)]\in[0,1]\), giving an ATE range \([-m_0,1-m_0]\). For \(p=1\), \(E[Y(1)]=m_1\) and \(E[Y(0)]\in[0,1]\), giving \([m_1-1,m_1]\).

The principal alternative is to make additional assumptions to obtain a narrower interval or a point estimate. For example, unconfoundedness asserts that, possibly conditional on measured covariates, treatment assignment is independent of the potential outcomes. Under suitable overlap and consistency conditions, adjustment methods can then identify an ATE. Such an estimate may be more precise, but its causal interpretation rests on the added assumptions. Instrumental-variable analyses use a different set of assumptions and may target a different effect. Neither approach is justified merely because its estimate is numerically narrower than the no-assumption interval.

These are population bounds. In an application, \(p\), \(m_1\), and \(m_0\) are usually estimated from a sample, so substituting sample proportions produces estimated endpoints, not a confidence interval for the ATE. Sampling uncertainty requires separate inferential methods. A narrow confidence interval around estimated endpoints would not remedy a failure of the causal assumptions, and a broad identification interval is not the same thing as sampling imprecision.

Common Mistakes

  • Treating the observed difference as the causal effect. The reasoning “\(m_1-m_0\) is the ATE because those are the two observed means” ignores that the groups may differ before treatment. The correction is to distinguish factual outcomes from potential outcomes and use the stated bounds unless additional assumptions justify identification.
  • Using the same missing-outcome contribution for both potential outcomes. \(Y(1)\) is observed for treated individuals, while \(Y(0)\) is observed for untreated individuals. The unknown treated-group contribution to \(E[Y(0)]\) is weighted by \(p\); the unknown untreated-group contribution to \(E[Y(1)]\) is weighted by \(1-p\).
  • Calling the interval a confidence interval. The displayed interval describes values compatible with the population restrictions. Sample-based uncertainty around estimated endpoints is a separate issue.
  • Concluding that an interval containing zero proves no effect. It says zero is among the values not ruled out by these restrictions; it does not establish that the true effect equals zero.
Key takeaway: Binary outcomes restrict every missing potential-outcome mean to \([0,1]\). Combining those restrictions with treatment prevalence and the two observed group means yields an ATE interval. Without additional assumptions, the interval has width one and may include effects in both directions.

Check Your Understanding

Work each problem, then compare your result with the solutions below.

  1. A hypothetical population has \(p=0.50\), \(m_1=0.40\), and \(m_0=0.30\). Calculate the lower and upper ATE bounds.
  2. For \(p=0.25\), \(m_1=0.80\), and \(m_0=0.20\), calculate the bounds and state their units of interpretation.
  3. Why does an observed difference \(m_1-m_0=0\) not, by itself, imply an ATE of zero?
  4. If \(p=1\) and the observed treated mean is \(m_1=0.70\), what range for the ATE follows from the binary outcome restriction?

Solutions. (1) The lower endpoint is \(0.50(0.40)-[0.50(0.30)+0.50]=0.20-0.65=-0.45\). The upper endpoint is \([0.50(0.40)+0.50]-0.50(0.30)=0.70-0.15=0.55\). (2) The lower endpoint is \(0.25(0.80)-[0.75(0.20)+0.25]=0.20-0.40=-0.20\). The upper endpoint is \([0.25(0.80)+0.75]-0.75(0.20)=0.95-0.15=0.80\). These are risk differences, from a 20-percentage-point decrease to an 80-percentage-point increase. (3) The observed means describe different treatment groups. Without an assumption linking group membership to potential outcomes, they do not reveal the unobserved counterfactual outcomes needed for the ATE. (4) With everyone treated, \(E[Y(1)]=0.70\) and \(E[Y(0)]\in[0,1]\), so the ATE lies in \([0.70-1,0.70-0]=[-0.30,0.70]\).

Clinical Biostats