Tutorials › Propensity Score Matching and Instrumental Variables › End-To-End Partial Identification Analysis

Partial Identification and Bounds · Tutorial 439 of 699

End-To-End Partial Identification Analysis

Learn to move from a causal estimand and explicit assumptions to baseline bounds, justified tightening restrictions, and a confidence set that reflects sampling uncertainty.

Advanced 12 min read

What You'll Learn

  • Define an average treatment effect when treatment is confounded and the outcome is bounded.
  • Derive support-based bounds from observed treated and untreated outcome means.
  • Use auxiliary restrictions to tighten bounds only when their assumptions are defensible.
  • Construct a conservative confidence set for the identified interval using bounded-outcome concentration bounds.
  • Diagnose incompatible restrictions, unsupported precision, and incorrect claims about what an interval establishes.
  • Distinguish an identified interval from a point estimate and from a confidence set.

From a Causal Question to a Reportable Interval

A partial-identification analysis is not complete when it produces two endpoints. The endpoints answer a causal question only if the estimand, outcome limits, restrictions, and uncertainty procedure all match the target population and the available data. This tutorial works through that chain for an average treatment effect when treatment may be confounded and the outcome is bounded.

Consider a hypothetical treatment and a bounded outcome, measured in points from 0 to 10. Let \(D\) indicate treatment, and let \(Y(1)\) and \(Y(0)\) denote the outcomes the same unit would have under treatment and no treatment. The observed outcome is \(Y=DY(1)+(1-D)Y(0)\). Consistency is required for this equality; no assumption that treatment is unconfounded is made.

Definition: The average treatment effect (ATE) in the target population is \(\tau=E[Y(1)-Y(0)]\). A parameter is partially identified when the stated data and assumptions determine a set of possible values rather than a single value. An identified interval is the interval between the smallest and largest values in that set.

The analysis below assumes the target population's treated share \(q=P(D=1)\) is known to be 0.60. The outcome satisfies \(0\leq Y(d)\leq 10\) for both treatment states and every unit. Both treatment groups occur in the target population. These support and population-share assumptions are substantive inputs, not consequences of the observed outcome means.

Deriving the Baseline Bounds

Write \(\mu_{11}=E[Y(1)\mid D=1]\) for the observed treated-group mean and \(\mu_{00}=E[Y(0)\mid D=0]\) for the observed untreated-group mean. The two counterfactual means, \(E[Y(0)\mid D=1]\) and \(E[Y(1)\mid D=0]\), are not observed. With no restrictions beyond the outcome support, each is between 0 and 10. Decomposing the ATE by observed treatment status gives

$$ \tau=q\{E[Y(1)\mid D=1]-E[Y(0)\mid D=1]\} +(1-q)\{E[Y(1)\mid D=0]-E[Y(0)\mid D=0]\}. $$

To make the ATE as small as possible, set each missing counterfactual mean to the value that lowers its group's treatment contrast: set \(E[Y(0)\mid D=1]=10\) and \(E[Y(1)\mid D=0]=0\). To make it as large as possible, set those means to 0 and 10, respectively. Thus the baseline identified interval is

$$ \left[ q(\mu_{11}-10)+(1-q)(0-\mu_{00}),\quad q(\mu_{11}-0)+(1-q)(10-\mu_{00}) \right]. $$

These endpoints are sharp under these assumptions: the missing potential outcomes can be set to their respective support limits without contradicting the observed means. This is a statement about what the assumptions permit, not a claim that the extreme values are plausible. The earlier material on linear programming for sharp causal bounds gives a general framework for more complicated constraints; the expression here is the simple two-group case.

Worked Example: Baseline Bounds From Outcome Limits

Suppose a hypothetical sample has 120 treated people and 80 untreated people, with sample outcome means of 7 points and 5 points, respectively. For this plug-in calculation, take the target treated share to be the known value \(q=0.60\), and use the support \(0\) to \(10\). The corresponding sample-based baseline bounds are

Lower endpoint. Assign the missing untreated potential outcome among treated people its maximum, 10, and the missing treated potential outcome among untreated people its minimum, 0:

$$ 0.60(7-10)+0.40(0-5) =-1.80-2.00=-3.80\text{ points}. $$

Upper endpoint. Reverse those assignments:

$$ 0.60(7-0)+0.40(10-5) =4.20+2.00=6.20\text{ points}. $$

The plug-in interval is \([-3.80,6.20]\) points. Its width reflects both the observed group contrast and the unknown outcomes under the treatments each group did not receive. Because the group means are estimated from a sample, this interval alone does not quantify sampling uncertainty.

Adding Restrictions That Have a Defensible Basis

Bounds can become more informative when credible information constrains the missing potential outcomes. The restriction must be justified for the same outcome definition and target population. It cannot be selected merely because it produces a narrower interval. For example, a carefully matched auxiliary data source might support bounds on the mean untreated outcome among treated people and the mean treated outcome among untreated people. Such bounds are moment restrictions: they constrain group averages, not each person's unobserved outcome.

For illustration, suppose external evidence that is applicable to this hypothetical target population supports \(4\leq E[Y(0)\mid D=1]\leq 6\) and \(6\leq E[Y(1)\mid D=0]\leq 8\). Assume these restrictions are valid and known for the target population. The lower endpoint uses the largest permitted missing \(Y(0)\) mean in the treated group and the smallest permitted missing \(Y(1)\) mean in the untreated group. The upper endpoint uses the opposite ends. The support-only bounds still apply, but these restrictions tighten them.

Worked Example: Tightening Bounds With Auxiliary Restrictions

Using the same hypothetical sample means, \(q=0.60\), and the restrictions above, calculate the plug-in interval. For the lower endpoint, use 6 for the missing untreated mean among treated people and 6 for the missing treated mean among untreated people:

$$ 0.60(7-6)+0.40(6-5) =0.60+0.40=1.00\text{ point}. $$

For the upper endpoint, use 4 for the missing untreated mean among treated people and 8 for the missing treated mean among untreated people:

$$ 0.60(7-4)+0.40(8-5) =1.80+1.20=3.00\text{ points}. $$

The restricted plug-in identified interval is \([1.00,3.00]\) points. It lies inside the baseline interval \([-3.80,6.20]\), as it must if the added restrictions are compatible and genuinely narrow the set of admissible potential-outcome distributions. The calculation does not establish that the effect is positive in the population: sampling uncertainty in the observed means remains to be addressed.

In a larger problem, restrictions can be combined with observed-data constraints and optimized over the feasible set, for example by linear programming. Before interpreting the resulting endpoints, check that the restrictions are mutually compatible with the outcome support and observed data. If they imply an empty feasible set, the problem is not evidence of an impossible causal effect; it signals that at least one input, definition, or restriction needs review.

Quantifying Sampling Uncertainty

An identified interval describes what the population assumptions permit; a confidence set accounts for the fact that the observed group means are estimated. To show the calculation transparently, suppose the 120 treated and 80 untreated observations are independent random samples within their respective treatment groups, and each outcome lies in \([0,10]\). The target share \(q=0.60\) and the auxiliary restrictions are treated as known. If the share or restrictions are estimated, their uncertainty must also be included.

For a mean from \(n_g\) independent observations bounded in an interval of width \(R\), Hoeffding's inequality gives \(P(|\bar{Y}-\mu|\geq r)\leq 2\exp(-2n_g r^2/R^2)\). Set \(\alpha=0.05\), \(R=10\), and use the conservative radius \(r_g=R\sqrt{\log(4/\alpha)/n_g}\). For either group the probability of missing its interval is at most \(2\exp\{-2\log(4/\alpha)\}=\alpha^2/8\). By the union bound, the probability that either group mean falls outside its interval is at most \(\alpha^2/4=0.000625\), which is below 0.05. Thus the two intervals cover both population means simultaneously with probability at least 95% (in fact, this bound guarantees at least 99.9375%). This deliberately conservative choice makes the coverage guarantee clear.

Worked Example: A Confidence Set for the Restricted Bounds

For the treated mean, the radius is \(r_1=10\sqrt{\log(80)/120}\approx1.910934\). For the untreated mean, it is \(r_0=10\sqrt{\log(80)/80}\approx2.3404\). Use the unrounded radii in endpoint calculations. The corresponding mean intervals are approximately \([5.0891,8.9109]\) and \([2.6596,7.3404]\). If displaying two decimal places, round these outward, for example to \([5.08,8.92]\) and \([2.65,7.35]\); rounding inward could invalidate the stated coverage guarantee.

For the baseline bounds, substitute the lower treated mean and upper untreated mean into the lower endpoint, and the upper treated mean and lower untreated mean into the upper endpoint:

$$ L_{\mathrm{base}}=0.60(5.089066-10)+0.40(0-7.3404) \approx-5.88272, $$
$$ U_{\mathrm{base}}=0.60(8.910934-0)+0.40(10-2.6596) \approx8.28272. $$

The confidence set for the baseline identified interval is therefore reported conservatively as \([-5.883,8.283]\) points. In particular, the upper endpoint is rounded upward, not down to 8.282.

Now apply the auxiliary restrictions. The smallest restricted ATE uses the lower treated mean, the upper permitted missing \(Y(0)\) mean of 6, the lower permitted missing \(Y(1)\) mean of 6, and the upper untreated mean:

$$ L_{\mathrm{restricted}} =0.60(5.089066-6)+0.40(6-7.3404) \approx-1.08272. $$

The largest uses the upper treated mean, the lower permitted missing \(Y(0)\) mean of 4, the upper permitted missing \(Y(1)\) mean of 8, and the lower untreated mean:

$$ U_{\mathrm{restricted}} =0.60(8.910934-4)+0.40(8-2.6596) \approx5.08272. $$

A conservative report is \([-1.083,5.083]\) points. Under the stated sampling and restriction assumptions, this confidence set contains the population identified interval with at least 95% repeated-sampling coverage. It is not a 95% probability statement about a fixed effect, and it does not show that the population effect is positive. Its width reflects both sampling uncertainty and the remaining partial identification.

Diagnostics, Alternatives, and Limits

Before reporting an interval, audit the analysis in the order the quantities enter it. Confirm that the outcome support is clinically or substantively possible for every relevant potential outcome, that the target population and \(q\) are specified, and that the observed group means use the same outcome scale. For every added restriction, record its source, target population, and whether it is treated as exact or estimated. Check that the restricted interval is contained within the support-only interval. A violation often indicates a coding error, mismatched population, or incompatible restriction.

The main alternatives answer different questions. A point estimate requires assumptions strong enough to identify one effect, such as assumptions underlying an adjustment or matching analysis; it is not automatically preferable to a wide but honest interval. Baseline bounds require fewer restrictions but can be uninformative. Auxiliary constraints can sharpen bounds, but only to the extent their validity is credible. A confidence set addresses sampling error; it does not repair confounding assumptions or make weak restrictions persuasive. The introductory material on reporting identified intervals emphasizes the same distinction: report what is identified, not a single preferred value from inside the interval.

A frequent mistake is to read the sample-based interval \([1,3]\) as proof that the true ATE is positive, or to describe \([-1.083,5.083]\) as a range in which the effect has a 95% probability of lying. The first confuses a plug-in identified set with population inference; the second assigns a probability interpretation that a frequentist confidence procedure does not provide. The correction is to report the confidence set and its assumptions, then state that it contains the population identified set with the specified repeated-sampling coverage. If the interval includes zero, these data and assumptions do not rule out a zero ATE.

Finally, the numerical coverage guarantee depends on the stated sampling design, bounded outcomes, and fixed external inputs. If observations are clustered, sampled with unequal probabilities, or dependent, the displayed Hoeffding calculation may not apply as written. If auxiliary restrictions are estimated, their uncertainty must be incorporated rather than silently treating them as known. A narrow interval is not a substitute for these checks.

1
Define the estimand and target.
State the treatment, outcome scale, target population, and the ATE being bounded.
2
Write the assumptions.
Record consistency, outcome limits, treatment share, sampling conditions, and any auxiliary restrictions.
3
Compute and audit the bounds.
Calculate support-only endpoints first; add only justified, compatible restrictions.
4
Account for estimation and report carefully.
Construct a confidence set, preserve coverage when rounding, and distinguish it from the identified interval.

Check Your Understanding

Answer each question, then compare with the solutions below.

  1. With \(q=0.50\), treated mean 6, untreated mean 4, and outcome support \([0,10]\), what are the baseline plug-in bounds?
  2. For that same setting, suppose the missing mean \(E[Y(0)\mid D=1]\) is restricted to \([3,5]\), and \(E[Y(1)\mid D=0]\) to \([5,7]\). What are the restricted plug-in bounds?
  3. Why is it not valid to treat a restriction estimated from an auxiliary sample as known without qualification?
  4. What does a 95% confidence set for an identified interval mean under repeated sampling?

Solutions. 1. The lower endpoint is \(0.50(6-10)+0.50(0-4)=-2-2=-4\) points; the upper is \(0.50(6-0)+0.50(10-4)=3+3=6\) points. 2. The lower endpoint uses missing means 5 and 5: \(0.50(6-5)+0.50(5-4)=1\) point. The upper uses 3 and 7: \(0.50(6-3)+0.50(7-4)=3\) points. 3. The auxiliary restriction has sampling uncertainty; treating it as fixed can produce a confidence set with too little coverage. Its uncertainty must be propagated or a justified deterministic restriction must be established. 4. Across repetitions under the sampling assumptions, the procedure produces a set that contains the population identified interval at least 95% of the time. This is not a 95% probability that a fixed parameter lies in the realized set.

Clinical Biostats