Tutorials › Group Sequential Methods and Sample Size Re-Estimation › Sensitivity Analysis for Assumed Effects

Conditional and Predictive Power · Tutorial 354 of 699

Sensitivity Analysis for Assumed Effects

You will calculate conditional-power curves across plausible treatment effects and assess whether a decision is robust to the assumed effect.

Advanced 10 min read

What You'll Learn

  • Define conditional power as a function of the assumed future treatment effect.
  • Calculate plug-in and design-effect conditional power under a canonical normal model.
  • Interpret how quickly conditional power changes as the assumed effect varies.
  • Identify assumptions and diagnostics that can make a sensitivity analysis unreliable.
  • Decide whether conclusions hold across a plausible range of effects.

Why Vary the Assumed Effect?

Conditional power (CP) asks how likely the planned remaining analysis is to meet its success criterion, given the interim data and an assumption about the treatment effect from now on. That assumption can matter greatly. A plug-in calculation may use the treatment effect estimated at the interim look; a design-effect calculation may instead use the effect specified when planning the trial. If those choices lead to different assessments, a single CP value conceals an important source of uncertainty.

Sensitivity analysis makes that dependence explicit. Rather than asking only for CP at one assumed effect, calculate CP across a scientifically plausible range. The resulting curve helps show whether a conclusion is stable over that range or relies heavily on choosing an optimistic or pessimistic effect. It does not determine which effect is correct, and it does not turn CP into a probability distribution over effects.

Definition: A conditional-power sensitivity analysis evaluates conditional power as a function of the assumed treatment effect, while holding the interim data, remaining design, and success criterion fixed. A conclusion is sensitive when it changes materially across plausible assumed effects; it is robust over a stated range when it remains substantively the same throughout that range.

Model and Conditional-Power Curve

For a transparent calculation, consider a one-sided test of \(H_0:\theta=0\) against \(H_a:\theta>0\), with a fixed final critical value \(c\). Use the canonical normal model: at information \(I\), the standardized score has variance 1 and mean \(\theta\sqrt{I}\). Here \(\theta\) is the treatment effect in response units, and information is in reciprocal squared response units, so that \(\theta\sqrt{I}\) is dimensionless. The model assumes known or adequately approximated variance, independent information increments, and a correctly specified information schedule.

Let \(I_t\) be the information available at the interim analysis, \(I_T\) the planned final information, and \(Z_t=z_t\) the observed interim statistic. Assume \(0<I_t<I_T\). Conditional on the interim data and a fixed assumed effect \(\theta\), the final statistic has a normal distribution with mean and variance

$$ E_\theta(Z_T\mid Z_t=z_t) = \sqrt{\frac{I_t}{I_T}}\,z_t + \theta\frac{I_T-I_t}{\sqrt{I_T}}, \qquad \operatorname{Var}(Z_T\mid Z_t=z_t) = 1-\frac{I_t}{I_T}. $$

The first term in the mean carries forward the observed interim score; the second is the assumed mean contribution from the remaining information. For the one-sided final success rule \(Z_T\ge c\), the conditional power curve is therefore

$$ CP(\theta) = 1-\Phi\!\left( \frac{c-\sqrt{I_t/I_T}\,z_t-\theta(I_T-I_t)/\sqrt{I_T}} {\sqrt{1-I_t/I_T}} \right), $$

where \(\Phi\) is the standard normal cumulative distribution function. In a group sequential design with additional future looks, the success event is generally crossing one or more prespecified future boundaries. The same sensitivity principle applies, but CP must then be calculated for that full future crossing event; substituting a single final critical value would not represent the planned design.

Why the Curve Increases with the Assumed Effect

For a positive-direction test, the conditional standard deviation does not depend on \(\theta\), while the conditional mean increases with \(\theta\) whenever information remains. This gives a simple justification for the curve's direction. Define \(a=(I_T-I_t)/\sqrt{I_T}\) and \(s=\sqrt{1-I_t/I_T}\). Then \(a>0\), \(s>0\), and \(CP(\theta)=1-\Phi((c-m_0-a\theta)/s)\), where \(m_0=\sqrt{I_t/I_T}\,z_t\). Differentiating yields

$$ \frac{d\,CP(\theta)}{d\theta} = \phi\!\left(\frac{c-m_0-a\theta}{s}\right)\frac{a}{s} >0, $$

where \(\phi\) is the standard normal density. Thus CP is increasing in the assumed beneficial effect under this model. This does not mean every group sequential success probability can be summarized by the final-look formula; it means the effect assumption shifts the distribution of future evidence in the beneficial direction while the specified design is held fixed.

Worked Example: Calculating CP at Several Assumed Effects

Hypothetical data. Suppose a one-sided design uses final critical value \(c=1.96\), with \(I_T=1.00\) and \(I_t=0.60\). The observed interim statistic is \(z_t=1.50\). Calculate CP for assumed effects \(\theta=0\), \(1\), and \(2\), in the response units used to define \(\theta\).

The assumptions for this calculation are the canonical normal model, independent increments, fixed planned information, and the one-sided final rule \(Z_T\ge1.96\). The interim information is less than the final information, so \(0.40\) information units remain. The conditional mean and standard deviation are

$$ m_\theta=\sqrt{0.60}(1.50)+0.40\theta =1.161895+0.40\theta, \qquad s=\sqrt{1-0.60}=\sqrt{0.40}=0.632456. $$

For \(\theta=0\), the standardized distance to the boundary is \( (1.96-1.161895)/0.632456=1.261915 \). Therefore \(CP(0)=1-\Phi(1.261915)\approx0.1035\).

For \(\theta=1\), the mean is \(m_1=1.561895\), so the distance is \( (1.96-1.561895)/0.632456=0.62946 \). Thus \(CP(1)=1-\Phi(0.62946)\approx0.2645\).

For \(\theta=2\), the mean is \(m_2=1.961895\), and the distance is \( (1.96-1.961895)/0.632456=-0.002996 \). Hence \(CP(2)=1-\Phi(-0.002996)\approx0.5012\).

The assumed effect changes CP from about 10.35% to about 50.12% across these values. The calculation does not say which effect is plausible. It shows that the CP conclusion is strongly dependent on that choice.

Plug-In and Design-Effect Values on the Same Curve

A plug-in CP calculation inserts an estimated treatment effect, often the interim estimate, into the CP formula as though it were the effect applying to the remaining data. A design-effect calculation inserts the prespecified effect used for planning. These are two evaluations of the same conditional-power function, not two different probability models. They answer different conditional questions because they use different effect assumptions.

In the known-variance canonical model above, a simple interim estimate is \(\hat{\theta}=z_t/\sqrt{I_t}\). This estimate can be noisy, especially when interim information is limited. Plug-in CP treats its value as fixed for the calculation; it does not average over uncertainty about \(\theta\). A predictive probability, by contrast, integrates future success probabilities over a distribution for the unknown effect. That is a distinct quantity.

Worked Example: Comparing Plug-In and Design-Effect CP

Hypothetical continuation of the design. Use the data and one-sided rule in the previous example. Suppose the design effect is \(\theta_D=1.50\). Calculate the plug-in and design-effect CP values, and compare them with the curve rather than treating either as the only possible answer.

The interim plug-in estimate is \( \hat{\theta}=z_t/\sqrt{I_t}=1.50/\sqrt{0.60}=1.936492 \). Using this assumed effect gives conditional mean \(1.161895+0.40(1.936492)=1.936492\). The standardized distance is \( (1.96-1.936492)/0.632456=0.03717 \), so \(CP(\hat{\theta})=1-\Phi(0.03717)\approx0.4852\).

For the design effect \(\theta_D=1.50\), the conditional mean is \(1.161895+0.40(1.50)=1.761895\). The standardized distance is \( (1.96-1.761895)/0.632456=0.31323 \), giving \(CP(\theta_D)=1-\Phi(0.31323)\approx0.3770\).

The plug-in value is higher because the interim estimate exceeds the design effect. Neither value is an unconditional chance of eventual success, nor does the higher plug-in value establish that the larger effect is credible. The difference is a reason to report the effect assumption and inspect the curve across a defensible range.

Assessing Whether a Conclusion Is Robust

Sensitivity analysis is most useful when it is tied to a decision or interpretation. For example, an analysis may label CP below 0.30 as low. That threshold is a decision convention, not a universal statistical boundary. To assess robustness, find whether plausible assumed effects lie on both sides of the threshold, and report the range over which the classification changes. The choice of plausible range should have a scientific rationale; selecting it after seeing the curve can make the result misleading.

Worked Example: Locating a CP Threshold

Hypothetical decision rule. For the same interim data and design, suppose the team uses \(CP=0.30\) as a descriptive threshold and considers \(\theta\) values from \(0.80\) to \(1.40\) scientifically plausible. Find approximately where the CP curve reaches 0.30 and determine whether the classification is robust over that range.

Using \(m_0=1.161895\), \(a=0.40\), and \(s=0.632456\), solve \(1-\Phi((1.96-m_0-0.40\theta)/s)=0.30\). Since \(\Phi^{-1}(0.70)\approx0.5244\), the threshold occurs when

$$ \frac{1.96-1.161895-0.40\theta}{0.632456}=0.5244, \qquad \theta = \frac{1.96-1.161895-0.632456(0.5244)}{0.40} \approx1.166. $$

Because CP increases with \(\theta\), values below about \(1.166\) give CP below 0.30, and values above it give CP above 0.30. The plausible range \(0.80\) to \(1.40\) crosses that point, so the threshold classification is not robust over the stated range. This conclusion is about sensitivity to the assumption, not evidence that either side of the threshold is the right decision.

Diagnostics, Alternatives, and Limitations

Before interpreting a curve, check that it represents the planned analysis. Verify the interim information fraction, remaining information, observed statistic, direction of the test, and the applicable boundary or boundaries. Under the simple final-look model, the curve should increase with \(\theta\) for a beneficial effect and should remain between 0 and 1. A decreasing curve, values outside that range, or a discontinuity in a smooth one-look calculation signals a sign, scaling, or implementation error. With multiple future looks, the curve may need numerical evaluation of correlated future statistics, but the success event must still match the design.

The canonical calculation also relies on its variance and information assumptions. If the variance is estimated with substantial uncertainty, if information is mismeasured, or if the assumed effect is not an appropriate model for future data, the displayed probabilities can be poorly calibrated. Changing the assumed effect alone does not address those problems. Nor should one silently replace the original group sequential boundaries with a nominal final-look cutoff; that changes the design's rejection rule and can change its error properties.

A sensitivity curve complements, rather than replaces, the main alternatives. A single plug-in value is compact but can be dominated by a noisy interim estimate. A single design-effect value is tied to the planning assumption but does not show how conclusions shift away from it. A predictive probability represents uncertainty about the effect by averaging conditional success probabilities over an effect distribution; it answers a different question and depends on that distribution. Sensitivity analysis instead displays conditional power at fixed effect values, without assigning probabilities to those values.

The essential limitation is that a curve does not choose the scientifically credible effect range or a stopping threshold. Those choices require substantive justification and should be stated. The useful output is not simply a graph: it is a contextual statement identifying the assumptions, range examined, and whether the decision-relevant conclusion changes across that range.

Key takeaway: Report conditional power as a function of the assumed treatment effect when that assumption could alter the interpretation. A plug-in or design-effect value is one point on the curve; robustness can only be assessed relative to a justified effect range and a clearly stated decision criterion.

Check Your Understanding

Answer each question, then compare with the brief solutions below.

  1. In the canonical model, what happens to one-sided conditional power as the assumed beneficial effect increases, and why?
  2. Why is plug-in CP not the same as predictive probability?
  3. If every plausible effect value gives CP above a stated descriptive threshold, what can be concluded about robustness, and what cannot be concluded?
  4. Name two checks that could reveal a mis-specified conditional-power calculation in a group sequential design.

Solutions. (1) CP increases because the assumed effect raises the conditional mean of future evidence while the conditional standard deviation stays fixed. (2) Plug-in CP treats one estimated effect as fixed; predictive probability averages CP over a distribution for the unknown effect. (3) The classification is robust over the specified range and under the stated model, but this does not prove that the effect range is correct or that the decision threshold is universally appropriate. (4) Check that the information fraction and remaining information match the design, and that the success event uses the planned boundaries; also check the test direction and that calculated probabilities lie between 0 and 1.

Clinical Biostats