Tutorials › Pattern-Mixture and Selection Models › Worked Example: Treatment-Effect Conclusions Across Scenarios

Sensitivity Analysis for Treatment Effects · Tutorial 434 of 499

Worked Example: Treatment-Effect Conclusions Across Scenarios

You will calculate assumption-conditional treatment effects and confidence intervals for a hypothetical trial, then classify the evidence against a clinical-benefit threshold.

Advanced 11 min read

What You'll Learn

  • Define an arm-specific pattern-mixture shift and use it to adjust a trial’s estimated mean outcome.
  • Calculate a treatment-effect estimate and confidence interval under fixed sensitivity assumptions.
  • Distinguish statistical evidence of an effect from evidence that it exceeds a clinical-benefit threshold.
  • Compare conclusions across scenarios without treating one scenario as established fact.
  • Identify when a conditional confidence interval omits important sources of uncertainty.

From sensitivity assumptions to a clinical conclusion

A treatment-effect estimate is not a complete conclusion when some trial outcomes are missing. The estimated effect may depend on assumptions about the unobserved outcomes, so a useful analysis shows how the estimate and its uncertainty change under plausible alternatives. This tutorial works through a hypothetical two-arm trial and classifies the conclusions under several such assumptions.

The outcome is improvement in a symptom score at 12 weeks, measured in points; higher values mean greater improvement. The estimand is the intention-to-treat population mean improvement in the treatment arm minus the corresponding mean in the control arm. A positive contrast therefore favors treatment. Suppose clinicians have specified 0.5 points as the threshold for a clinically important benefit. This threshold is a benchmark for interpretation, not a value estimated from the trial.

Definition: In an arm-specific pattern-mixture sensitivity analysis, let \(\bar{Y}_{a,\mathrm{obs}}\) be the observed mean outcome in arm \(a\), and let \(\delta_a\) be the assumed difference between the mean among participants with missing outcomes and the observed mean in that same arm. Thus, the assumed missing-outcome mean is \(\bar{Y}_{a,\mathrm{obs}}+\delta_a\). If \(r_a\) is the fraction missing in arm \(a\), the arm’s adjusted mean is \(\bar{Y}_{a,\mathrm{obs}}+r_a\delta_a\). The treatment effect is the adjusted treatment mean minus the adjusted control mean.

The shift \(\delta_a\) is measured in outcome points. A negative shift says that, on average, participants with missing outcomes are assumed to have improved less than observed participants in that arm. The shifts are assumptions, not quantities identified by the observed outcomes alone. The introductory material on Pattern-Mixture Models provides the broad framework; here the emphasis is on carrying assumptions through to an effect estimate, interval, and clinical interpretation.

Trial data and a working interval calculation

Consider 100 randomized participants in each arm. In the treatment arm, 80 outcomes are observed, with mean improvement 8.0 points and sample standard deviation 5.0 points. In the control arm, 90 outcomes are observed, with mean improvement 6.0 points and sample standard deviation 4.5 points. Thus, the observed-data missing fractions are \(r_T=20/100=0.20\) and \(r_C=10/100=0.10\).

For the worked calculations, treat each scenario’s shifts and the missing fractions as fixed. Use the observed sample means and standard deviations to quantify sampling uncertainty, with independent arms and a large-sample normal approximation. Under this working approximation, the adjusted arm mean is the observed mean plus a fixed shift contribution, so its estimated standard error is the observed mean’s standard error. The standard error for the treatment-minus-control contrast is therefore

$$ SE(\hat{\Delta})= \sqrt{\frac{s_T^2}{n_{T,\mathrm{obs}}} +\frac{s_C^2}{n_{C,\mathrm{obs}}}}. $$

Here \(n_{T,\mathrm{obs}}=80\) and \(n_{C,\mathrm{obs}}=90\). This gives

$$ SE(\hat{\Delta})= \sqrt{\frac{5.0^2}{80}+\frac{4.5^2}{90}} =\sqrt{0.3125+0.2250} =\sqrt{0.5375} =0.7331\text{ points}. $$

We use the approximate 95% interval \(\hat{\Delta}\pm1.96SE(\hat{\Delta})\), with margin of error \(1.96(0.7331)=1.437\) points. The interval is conditional on the scenario shifts and this working variance calculation; it does not represent uncertainty about which shifts are appropriate.

Worked Example: The Reference Scenario

First use the reference assumption \(\delta_T=\delta_C=0\). Under this assumption, the mean among participants with missing outcomes is taken to equal the observed mean within each arm.

The adjusted treatment mean is \(8.0+0.20(0)=8.0\) points. The adjusted control mean is \(6.0+0.10(0)=6.0\) points. Thus the estimated treatment effect is \(8.0-6.0=2.0\) points.

The approximate 95% confidence interval is \[ 2.0\pm1.96(0.7331)=2.0\pm1.437=(0.563,\,3.437)\text{ points}. \] Rounded to three significant figures, this is \((0.563,\,3.44)\) points.

The interval excludes zero, so this scenario gives evidence of a positive mean effect under the stated assumptions and approximation. Its lower endpoint, 0.563 points, also exceeds the prespecified clinical threshold of 0.5 points. In this scenario the interval supports a benefit greater than that threshold; this conclusion remains conditional on the missing-outcome assumption.

How the assumptions change the effect

The adjustment follows directly from the arm-level mixture. If a fraction \(1-r_a\) has observed outcomes and a fraction \(r_a\) has missing outcomes, the assumed full-arm mean is

$$ (1-r_a)\bar{Y}_{a,\mathrm{obs}} +r_a(\bar{Y}_{a,\mathrm{obs}}+\delta_a) =\bar{Y}_{a,\mathrm{obs}}+r_a\delta_a. $$

Subtracting the adjusted control mean from the adjusted treatment mean gives

$$ \hat{\Delta} =(\bar{Y}_{T,\mathrm{obs}}-\bar{Y}_{C,\mathrm{obs}}) +r_T\delta_T-r_C\delta_C. $$

This derivation shows why the same shift need not have the same effect on the treatment contrast in both arms: the adjustment is multiplied by that arm’s missing fraction. It also makes the direction clear. A more unfavorable assumed shift in the treatment arm reduces the contrast, while a more unfavorable shift in the control arm increases it.

Worked Example: A Moderately Unfavorable Scenario

Suppose the missing treatment-arm outcomes are assumed to average 4 points below the observed treatment-arm mean, so \(\delta_T=-4\). Suppose missing control-arm outcomes are assumed to average 1 point below their observed mean, so \(\delta_C=-1\).

The adjusted treatment mean is \(8.0+0.20(-4)=8.0-0.8=7.2\) points. The adjusted control mean is \(6.0+0.10(-1)=6.0-0.1=5.9\) points. The treatment effect is \(7.2-5.9=1.3\) points. Equivalently, the observed contrast of 2.0 points changes by \(0.20(-4)-0.10(-1)=-0.7\) points, giving \(2.0-0.7=1.3\).

Using the same conditional standard error and margin of error, the approximate 95% confidence interval is \(1.3\pm1.437=(-0.137,\,2.737)\) points. The interval includes zero and the clinical threshold of 0.5 points. The point estimate exceeds 0.5, but the interval does not establish a positive effect or an effect larger than the clinical threshold. Nor does it rule out such an effect.

Worked Example: A Severe Treatment-Arm Shift

Now suppose the missing treatment-arm outcomes are assumed to average 8 points below the observed treatment-arm mean, while missing control outcomes are assumed to equal the observed control mean: \(\delta_T=-8\) and \(\delta_C=0\).

The adjusted treatment mean is \(8.0+0.20(-8)=6.4\) points. The adjusted control mean is \(6.0+0.10(0)=6.0\) points. The effect is \(6.4-6.0=0.4\) points. Its approximate 95% confidence interval is \(0.4\pm1.437=(-1.037,\,1.837)\) points.

This interval includes zero and the 0.5-point threshold. The data under this assumption do not establish a positive mean effect or clinically important benefit. However, the upper endpoint is above 0.5, so the interval also does not rule out a clinically important benefit. It is not correct to conclude that this scenario demonstrates no benefit merely because its point estimate is below the threshold.

Compare the conclusions, not just the estimates

The calculations can be summarized with separate statements about direction, statistical evidence, and the clinical benchmark. “Not established” does not mean “absent”: an interval may be too wide to distinguish an important benefit from no benefit or harm.

ScenarioEffect (95% CI), pointsConclusion relative to zeroConclusion relative to 0.5 points
Both shifts zero2.00 (0.563, 3.44)Interval supports a positive effectInterval supports an effect above the threshold
\(\delta_T=-4,\ \delta_C=-1\)1.30 (-0.137, 2.74)Positive effect not establishedBenefit above threshold not established or ruled out
\(\delta_T=-8,\ \delta_C=0\)0.400 (-1.04, 1.84)Positive effect not establishedBenefit above threshold not established or ruled out

The first scenario supports both a positive effect and an effect exceeding the clinical threshold, conditional on its assumptions. The other two scenarios weaken those conclusions, but neither interval shows that clinically important benefit is impossible. The result is therefore not uniformly robust across these three assumptions. A scenario table should be interpreted as a set of conditional results, not as evidence that one unobserved-outcome assumption is correct.

Assumptions, diagnostics, and limitations

The calculations rely on several assumptions and approximations. Randomization supports comparison of treatment assignments, but it does not determine the outcomes that would have been observed for participants who are missing. Each scenario must therefore state its shifts, including their units, direction, and clinical rationale. The observed means and standard deviations must also be meaningful summaries for the analyzed outcome, and the large-sample normal approximation should be adequate for the observed sample sizes and outcome distributions.

The interval calculation above treats missing fractions and shifts as fixed and uses observed-outcome sampling variability. It does not add uncertainty from estimating shifts, uncertainty about the missingness process, or extra imputation variability. It is a transparent working calculation for the example, not a general claim that this variance formula is appropriate for every fitted sensitivity model. In an actual analysis, the interval procedure must match the estimator and account for relevant sources of uncertainty.

Useful checks include verifying that missing fractions are calculated from the randomized arm totals, that the shift signs match the stated clinical assumptions, and that all adjusted means use the same outcome scale. Recompute the treatment contrast from the two adjusted arm means as well as from the shift formula. If these disagree, the arithmetic or sign convention is wrong. Also check whether the intervals are so wide that the scenario labels (“moderate” or “severe”) create more confidence than the data warrant.

The reference scenario is often described as a missing-at-random-style reference, but it should not be treated as verified simply because it is familiar. Complete-case analysis, model-based imputation, and selection models make different modeling choices; none can identify the missing outcomes without assumptions. Pattern-mixture shifts make one set of those assumptions explicit and allow their consequences to be compared. They do not prove that the chosen scenarios span all plausible departures.

Common pitfall: “The estimate is above the clinical threshold, so the trial establishes clinically important benefit.” A point estimate alone does not express its uncertainty. Compare the confidence interval with both zero and the clinical threshold. Conversely, an interval that crosses the threshold does not prove the true effect is below it.

A sound report states the scenario assumptions, the effect and interval in outcome units, and what each interval does and does not support. The scenarios should be selected and summarized according to a prespecified plan where possible; highlighting only the most favorable result can make conditional evidence appear unconditional. Statistical significance and clinical importance are distinct, and conclusions can change when assumptions change even if the observed trial data stay fixed.

Check Your Understanding

Answer each question using the trial and conventions above. Solutions are included after the questions.

  1. If \(\delta_T=-2\) and \(\delta_C=0\), what is the adjusted treatment effect?
  2. Under the moderately unfavorable scenario, does the confidence interval establish an effect above 0.5 points? Explain.
  3. Why does the severe-shift interval not establish that clinically important benefit is absent?
  4. Name one source of uncertainty not represented by the working interval calculation.

Solutions. (1) The effect is \(2.0+0.20(-2)-0.10(0)=1.6\) points. (2) No. Its interval, \((-0.137,\,2.737)\) points, includes 0.5, so an effect above the threshold is not established; the interval also does not rule one out. (3) Its interval extends to 1.837 points, which is above 0.5, so values representing clinically important benefit remain compatible with the interval. (4) For example, uncertainty about the shift parameters or additional imputation variability is omitted because shifts are treated as fixed and the variance uses observed-outcome sampling variability.

Clinical Biostats