Tutorials › Causal Inference › Sensitivity Analysis as a Causal Robustness Problem

Sensitivity to Unmeasured Confounding · Tutorial 875 of 972

Sensitivity Analysis as a Causal Robustness Problem

Separate the assumptions used to identify a causal effect from the sensitivity parameters used to examine how robust an estimate is to unmeasured confounding.

Advanced 10 min read

What You'll Learn

  • Define a target causal estimand and an observed-data estimate
  • State the assumptions under which measured-covariate adjustment identifies an average treatment effect
  • Define an additive sensitivity parameter for residual confounding
  • Translate specified bias values into adjusted estimates and intervals
  • Distinguish sampling uncertainty from uncertainty about unmeasured confounding

Why Sensitivity Analysis Is a Causal Robustness Problem

An adjusted estimate from an observational study can be precise and still depend on an untestable claim: that the measured variables suffice to control confounding. Sensitivity analysis makes the consequences of relaxing that claim explicit. It asks how the causal conclusion would change if residual confounding produced a bias of a specified size, or within a specified range.

The first task is to keep three quantities separate: the target estimand (the causal effect we want to learn), the observed estimate (what the analysis calculates from the data under its identifying assumptions), and the sensitivity parameter (a value describing how far the target may differ from the quantity estimated under those assumptions). Sensitivity analysis does not make an untestable assumption disappear. It shows which conclusions depend on it.

Target Estimand and Identifying Assumptions

Definition: Let \(A\) be a binary treatment, with \(A=1\) for treatment and \(A=0\) for control. Let \(Y(1)\) and \(Y(0)\) be the potential outcomes for an individual under treatment and control. For a specified target population, the average treatment effect is \(\tau=E[Y(1)-Y(0)]\). The target population, outcome, treatment contrast, and time horizon are part of the estimand, not details to be chosen after seeing the result.

For example, an average effect on days alive and out of hospital over 30 days is not the same estimand as an effect on symptom score at six months. The units and direction of the outcome also matter: if a higher outcome is better, a positive \(\tau\) favors treatment; if a higher outcome is worse, it favors control.

In an observational study, one common identification strategy adjusts for measured pre-treatment covariates \(X\). The following assumptions state when that strategy identifies the average treatment effect:

  • Consistency: for a person who receives treatment \(A=a\), the observed outcome equals that person's potential outcome \(Y(a)\). The treatment must be sufficiently well-defined for this correspondence to make sense.
  • Conditional exchangeability: \(Y(a)\mathrel{\perp\!\!\!\perp} A\mid X\) for each \(a\in\{0,1\}\). Given the measured covariates, treatment assignment carries no further information about the potential outcomes.
  • Positivity: within covariate values represented in the target population, each treatment option has a positive probability: \(0<P(A=a\mid X=x)\) for \(a=0,1\).
  • No interference: one person's treatment does not change another person's potential outcome, unless the estimand and potential outcomes explicitly account for such effects.

With these assumptions, the causal mean under treatment level \(a\) is identified by averaging the observed conditional outcome mean over the target covariate distribution:

$$ E[Y(a)] = \sum_x E[Y\mid A=a,X=x]P_T(X=x), \qquad \tau = \sum_x \{E[Y\mid A=1,X=x]-E[Y\mid A=0,X=x]\}P_T(X=x). $$

Here \(P_T(X=x)\) denotes the distribution of covariates in the target population. For continuous covariates, the sum is replaced by an integral. In practice, an analysis estimates the conditional means and averages the resulting contrasts over the target population. Correct implementation also requires that the estimator adequately estimate those conditional means and that the data support the target covariate distribution. These are estimation and data-support concerns, distinct from the causal identification assumptions above.

Worked Example: State the Target Before Reading the Estimate

Suppose a hypothetical cohort study compares a new rehabilitation program with usual care. The outcome is the number of days alive and out of hospital during the 30 days after enrollment. Investigators standardize the two groups' outcomes over the covariate distribution of all eligible patients.

Target estimand. The target is the average difference in 30-day days alive and out of hospital if all eligible patients received the program rather than usual care: \(\tau=E_T[Y(1)-Y(0)]\). It is measured in days, and a positive value favors the program.

Suppose the standardized observed means are 24.6 days under the program and 22.1 days under usual care. The observed estimate is

$$ \hat{\theta}_{\mathrm{ID}}=24.6-22.1=2.5\text{ days}. $$

This calculation is an estimate of the target effect only if the identification assumptions hold, along with adequate estimation and support for the target population. The data alone do not establish conditional exchangeability. Thus 2.5 days is the observed-data estimate under the stated identifying strategy, not proof that the causal effect is 2.5 days.

Representing Residual Confounding With a Sensitivity Parameter

To examine a possible failure of conditional exchangeability, define a discrepancy between the causal target and the quantity identified by the measured-covariate analysis. Let \(\theta_{\mathrm{ID}}\) be the observed-data estimand identified by that analysis if its identifying assumptions hold. Define the additive sensitivity parameter

Definition: The residual-confounding bias parameter is \(b=\theta_{\mathrm{ID}}-\tau\). A positive \(b\) means the measured-covariate analysis overstates the target effect; a negative \(b\) means it understates it. The corresponding sensitivity-adjusted estimate at a specified value of \(b\) is \(\hat{\tau}(b)=\hat{\theta}_{\mathrm{ID}}-b\).

The parameter \(b\) is expressed in the outcome's units. It is not generally identifiable from the observed data: different levels of unmeasured confounding can be compatible with the same observed distribution. A sensitivity analysis therefore specifies values or ranges for \(b\), explains why they are relevant, and reports the resulting causal estimates. The parameter is a transparent way to represent a departure from the identifying assumptions; it is not a correction automatically learned from the data.

This additive representation is useful when the effect and possible bias can be compared on a meaningful scale. It is a general sensitivity framework, not a claim that every form of confounding produces a constant additive shift in every population. More elaborate sensitivity models may define the departure differently. Whatever the parameterization, the analysis must say what the parameter means, its units or scale, and how it relates to the target estimand.

Worked Example: Translate a Specified Bias Into an Effect

Return to the hypothetical estimate of 2.5 additional days. Suppose investigators examine a scenario in which unmeasured confounding makes the measured-covariate estimate too high by \(b=0.8\) days. Under that scenario, the adjusted estimate is

$$ \hat{\tau}(0.8)=\hat{\theta}_{\mathrm{ID}}-b =2.5-0.8=1.7\text{ days}. $$

The interpretation is conditional: if the residual bias is 0.8 days in this direction, the implied average effect is 1.7 additional days alive and out of hospital. The calculation does not show that the bias is 0.8 days, and it does not establish that the adjusted value is the true effect. It makes the consequence of that assumption explicit.

Sampling Uncertainty and Sensitivity Uncertainty

An estimate varies across samples even when the identifying assumptions hold. A confidence interval describes sampling uncertainty under the statistical model and assumptions used to construct it. Sensitivity analysis addresses a different uncertainty: whether the identifying assumptions are adequate, and how large the resulting discrepancy \(b\) might be. A narrow confidence interval does not, by itself, rule out substantial unmeasured confounding.

If \(b\) is treated as a fixed specified value, subtracting it from an estimate shifts the entire interval by the same amount. If the original interval is \([L,U]\), the shifted interval is \([L-b,U-b]\). This interval still reflects sampling uncertainty conditional on that value of \(b\); it does not incorporate uncertainty about whether that value is plausible. If \(b\) is allowed to vary over a range, reporting the results across that range shows how the conclusion changes, but the range itself needs justification.

Worked Example: Separate a Significance Threshold From a Null-Effect Threshold

Suppose a hypothetical analysis estimates a treatment effect of 2.0 points on a quality-of-life scale, with a 95% confidence interval of 0.4 to 3.6 points under its identifying assumptions. Assume higher scores are better. Consider a fixed positive bias \(b\), meaning the measured-covariate analysis overstates the causal effect by \(b\) points.

For \(b=0.7\), the adjusted estimate is \(2.0-0.7=1.3\) points. The interval shifts to \([0.4-0.7,\;3.6-0.7]=[-0.3,2.9]\) points. Conditional on this fixed bias value, the shifted interval includes zero; the data no longer exclude a null effect at the corresponding 95% confidence level.

The adjusted point estimate reaches zero when

$$ 2.0-b=0 \quad\Longrightarrow\quad b=2.0\text{ points}. $$

That is the bias required to move the point estimate to the null. It is not the bias required for the confidence interval first to include zero. The lower endpoint reaches zero when \(0.4-b=0\), or \(b=0.4\) points. These are different questions: the first concerns the point estimate, while the second concerns whether the interval excludes zero. Neither calculation establishes how plausible a bias of that size is.

Identification Assumptions Versus Robustness Assumptions

An identifying assumption is part of the argument that observed data determine a causal estimand. For the measured-covariate strategy above, conditional exchangeability is such an assumption. A sensitivity assumption instead specifies how the target might depart from the estimate if that identifying assumption fails. The two play different roles: the first supplies identification; the second evaluates dependence on that identification claim.

QuestionIdentifying assumptionSensitivity assumption
What is being asserted?Measured covariates suffice to control confounding, conditional on them.Residual confounding would shift the effect by a specified amount or within a stated range.
What role does it play?It connects the observed-data functional to the causal target.It describes how the target could differ if the identifying claim is imperfect.
Is it learned from the observed outcome and treatment data alone?Usually not fully testable from those data.Not generally; values or ranges require explicit substantive justification.

A useful sensitivity analysis keeps the estimand fixed while varying the assumed departure. Changing the target population, outcome, or treatment contrast changes the question, rather than merely testing robustness to confounding. It is also important to state which assumptions are held fixed across scenarios. For example, varying \(b\) does not automatically address positivity problems, an incorrectly defined treatment, or model misspecification.

Worked Example: Report a Range of Plausible Biases

Suppose a hypothetical study estimates 1.6 fewer symptom-days with treatment than control, using a scale on which a negative contrast favors treatment. For clarity, define the outcome as symptom-days avoided, so larger values are better; the observed estimate is then \(\hat{\theta}_{\mathrm{ID}}=1.6\) days avoided. Investigators consider a range \(0.3\leq b\leq 1.0\) days, where positive \(b\) means the observed estimate overstates benefit.

At the upper end of the proposed range, the adjusted estimate is \(1.6-1.0=0.6\) days avoided. At the lower end it is \(1.6-0.3=1.3\) days avoided. Since \(\hat{\tau}(b)=1.6-b\) decreases as \(b\) increases, the estimates over the full range are

$$ 1.6-1.0\leq \hat{\tau}(b)\leq 1.6-0.3, \qquad 0.6\leq \hat{\tau}(b)\leq 1.3\text{ days avoided}. $$

The estimated benefit remains positive throughout this specified range. The appropriate conclusion is that the point-estimate direction is robust to biases no larger than one day in the stated direction, under the chosen additive model. It is not a conclusion that all possible unmeasured confounding has been ruled out. If plausible bias could be larger, have the opposite direction, or vary in a way not represented by this model, this calculation would not cover that possibility.

A Short Derivation and a Practical Workflow

The additive adjustment follows directly from the definition \(b=\theta_{\mathrm{ID}}-\tau\). Rearranging gives \(\tau=\theta_{\mathrm{ID}}-b\). Replacing the unknown \(\theta_{\mathrm{ID}}\) with its observed estimate gives \(\hat{\tau}(b)=\hat{\theta}_{\mathrm{ID}}-b\). If a confidence interval for \(\theta_{\mathrm{ID}}\) is \([L,U]\) and \(b\) is fixed, subtracting the same constant from each possible endpoint gives \([L-b,U-b]\). If instead \(b\) is only known to lie between \(b_L\) and \(b_U\), the range of adjusted point estimates is \([\hat{\theta}_{\mathrm{ID}}-b_U,\hat{\theta}_{\mathrm{ID}}-b_L]\), because the adjusted estimate decreases as \(b\) increases.

1
Specify the target.
Name the treatment contrast, outcome, time horizon, target population, and effect measure.
2
State the identifying strategy.
List the causal assumptions and describe how the observed data are used to estimate the target under those assumptions.
3
Define the sensitivity parameter.
Give its meaning, direction, and units. State explicitly what kind of departure from the identifying assumptions it represents.
4
Evaluate and report scenarios.
Calculate adjusted estimates, and intervals when appropriate, over substantively justified values or ranges. Distinguish sensitivity to bias from sampling uncertainty.

Common Mistakes and Why This Matters

  • Calling the observed estimate the causal effect without qualification. The estimate has a causal interpretation only under the stated identification assumptions and adequate estimation. Label it as an observed-data estimate under a specified strategy.
  • Changing the estimand while claiming to test robustness. A different target population or outcome asks a different causal question. Keep the target fixed when varying a confounding sensitivity parameter.
  • Treating a sensitivity value as an estimated fact. The observed data generally do not reveal the amount of bias from an unmeasured variable. Describe the value as a scenario, and justify any range.
  • Confusing interval width with protection against confounding. A confidence interval can be narrow while the identifying assumptions remain questionable. Report sampling uncertainty and sensitivity results as distinct parts of the analysis.
  • Interpreting a favorable sensitivity range too broadly. Robustness over a chosen range means robust over that range and parameterization. It does not prove robustness to every unmeasured-confounding structure.

Sensitivity analysis matters because causal conclusions are often conditional on assumptions that cannot be verified from the observed data alone. Making the target, estimate, and sensitivity parameter explicit lets readers see both what the analysis identifies under its assumptions and how much a stated violation would change the result.

Check Your Understanding

Use the definitions and calculations in this tutorial to answer the following.

  1. In your own words, distinguish a target estimand from an observed estimate.
  2. For an average treatment effect, what does conditional exchangeability assert, and what role does it play?
  3. An estimate is 4.0 units and \(b=1.2\) units, with \(b\) defined as estimate minus target. What is the sensitivity-adjusted estimate?
  4. A confidence interval is \([0.5,5.5]\) and the specified bias is \(b=1.0\). What is the shifted interval, and what uncertainty does it represent?
  5. Why does robustness over a chosen range of \(b\) not establish robustness to all possible unmeasured confounding?
Clinical Biostats