Tutorials › Multi-Arm Multi-Stage (MAMS) Designs › Averaging Expected Sample Size Across Scenarios

Sample Size and Expected Sample Size · Tutorial 214 of 499

Averaging Expected Sample Size Across Scenarios

Calculate and interpret a scenario-weighted expected enrollment while keeping the treatment-effect assumptions and scenario probabilities explicit.

Advanced 10 min read

What You'll Learn

  • Define scenario weights for a finite set of treatment-effect configurations.
  • Derive scenario-specific expected enrollment from stage continuation probabilities.
  • Calculate the weighted mean enrollment and check the arithmetic.
  • Distinguish a probability-weighted summary from an arbitrary average across scenarios.
  • Diagnose incoherent weights, incompatible scenarios, and implausible results.
  • Compare weighted expectations with scenario-specific and worst-case summaries.

Why Average Across Scenarios?

A MAMS design can have different expected enrollment under different treatment-effect configurations. For example, a design may tend to stop early when no treatment is active, but continue more often when several treatments are effective. Reporting expected enrollment under one configuration answers a useful question, but it does not summarize enrollment across a stated range of plausible configurations.

A scenario-weighted expected enrollment combines scenario-specific expectations using probabilities assigned to those configurations. The result depends on both ingredients: the expected enrollment within each scenario and the weights placed on the scenarios. It is therefore important to show those ingredients, rather than present the weighted mean as if it were a property of the design alone.

Definition: Let \(S\) denote a scenario label, with possible values \(1,\ldots,K\). Scenario \(s\) specifies a treatment-effect configuration, such as a global null, a global alternative, or partial treatment activity. Let \(w_s=P(S=s)\) be its stated scenario probability, where \(w_s\geq 0\) and \(\sum_{s=1}^K w_s=1\). If \(N\) is total enrollment, the scenario-weighted expected enrollment is $$ E[N]=\sum_{s=1}^K w_s E[N\mid S=s]. $$

The weights represent a probability distribution over the configurations for the purpose of the summary. They are not the probabilities of rejecting hypotheses, the probabilities that an individual treatment works, or the continuation probabilities at an interim analysis. Nor are scenario probabilities automatically supplied by a MAMS design. They must be chosen or specified for the planning question being answered.

From Continuation Probabilities to a Weighted Mean

Suppose a design has \(J\) stages, and \(d_j\) participants are enrolled at stage \(j\) if that stage is reached. Here \(d_j\) is an incremental stage enrollment, not the cumulative enrollment through stage \(j\). Let \(\pi_{j,s}\) be the probability, under scenario \(s\), that stage \(j\) is reached. The first stage is always reached, so \(\pi_{1,s}=1\). If \(I_j\) equals 1 when stage \(j\) is reached and 0 otherwise, then \(N=\sum_{j=1}^J d_jI_j\). Taking the expectation under scenario \(s\) gives:

$$ E[N\mid S=s] = \sum_{j=1}^J d_j E[I_j\mid S=s] = \sum_{j=1}^J d_j\pi_{j,s}. $$

The expected enrollment across scenarios follows by averaging these scenario-specific expectations with their weights. Equivalently, the probability of reaching stage \(j\) in the scenario mixture is \(\sum_s w_s\pi_{j,s}\), so:

$$ \sum_{s=1}^K w_s E[N\mid S=s] = \sum_{s=1}^K w_s\sum_{j=1}^J d_j\pi_{j,s} = \sum_{j=1}^J d_j\sum_{s=1}^K w_s\pi_{j,s}. $$

This identity is a direct application of the law of total expectation: the scenarios form a mutually exclusive and exhaustive partition, and the probability of each scenario is \(w_s\). It also shows two equivalent calculation routes. First calculate \(E[N\mid S=s]\) for each scenario and then weight those values, or first average the stage-reach probabilities and then calculate expected enrollment. Both routes must give the same result when they use the same stage sizes, scenarios, and weights.

Assumptions for this calculation: The scenario labels are mutually exclusive and cover the configurations being averaged; the weights are nonnegative and sum to 1; each scenario-specific expectation uses the same defined enrollment measure and design; and the continuation probabilities are calculated under the treatment-effect configuration assigned to that scenario. If scenario-specific values are estimates from simulation, their numerical uncertainty should also be considered when interpreting the weighted result.

Worked Example: Weighting Three Treatment-Effect Configurations

Consider a hypothetical three-stage MAMS design. If reached, stages 1, 2, and 3 enroll 60, 80, and 100 participants, respectively. The scenario-specific probabilities of reaching stages 2 and 3 are as follows.

ScenarioWeight \(w_s\)\(\pi_{2,s}\)\(\pi_{3,s}\)
Global null0.200.250.08
Global alternative0.500.900.78
Partial treatment activity0.300.650.40

The probabilities of reaching stage 3 do not exceed those of reaching stage 2 in any row, as required: reaching stage 3 entails reaching stage 2. The weights are nonnegative and sum to \(0.20+0.50+0.30=1.00\).

Under the global null, expected enrollment is:

$$ E[N\mid S=\text{null}] = 60+80(0.25)+100(0.08) = 60+20+8 = 88\text{ participants}. $$

Under the global alternative, it is \(60+80(0.90)+100(0.78)=60+72+78=210\) participants. Under partial treatment activity, it is \(60+80(0.65)+100(0.40)=60+52+40=152\) participants. Thus the scenario-weighted expectation is:

$$ E[N] = 0.20(88)+0.50(210)+0.30(152) = 17.6+105+45.6 = 168.2\text{ participants}. $$

Under the stated scenario probabilities, the design’s average expected enrollment is 168.2 participants. This does not mean that an individual trial enrolls 168.2 people; actual enrollment is an integer-valued random outcome. Nor does the weighted mean say which scenario will occur. The scenario-specific values remain informative: they range from 88 to 210 participants.

Direct Weighting When Scenario Expectations Are Available

Sometimes the scenario-specific expected enrollments have already been calculated, for example from the continuation probabilities under each treatment-effect configuration. In that case, the averaging step is a weighted sum. Check the scenario definitions and probabilities before doing the arithmetic: a correct multiplication cannot repair weights that do not represent the intended scenario mixture.

Worked Example: Four Stages and Three Scenarios

For a hypothetical design, stage increments are 40, 60, 80, and 100 participants. The weights for three effect configurations are 0.20, 0.50, and 0.30. Suppose the respective stage-reach probabilities for stages 2 through 4 are \((0.30,0.12,0.05)\), \((0.70,0.55,0.40)\), and \((0.95,0.90,0.82)\). These probabilities are coherent because the chance of reaching a later stage is no greater than the chance of reaching an earlier one.

The scenario-specific expected enrollments are:

$$ \begin{aligned} E[N\mid S=1]&=40+60(0.30)+80(0.12)+100(0.05)=72.6,\\ E[N\mid S=2]&=40+60(0.70)+80(0.55)+100(0.40)=166,\\ E[N\mid S=3]&=40+60(0.95)+80(0.90)+100(0.82)=251. \end{aligned} $$

The weighted result is \(0.20(72.6)+0.50(166)+0.30(251)=14.52+83+75.30=172.82\) participants. This is below the maximum possible enrollment of \(40+60+80+100=280\) participants and above the minimum of 40 participants, as it should be.

Worked Example: Why the Weights Change the Answer

Suppose three hypothetical treatment-effect configurations have expected enrollments of 96, 154, and 208 participants, with stated probabilities 0.15, 0.35, and 0.50. The probabilities sum to 1, so the weighted expectation is:

$$ 0.15(96)+0.35(154)+0.50(208) = 14.4+53.9+104 = 172.3\text{ participants}. $$

An equal-weight average would instead be \((96+154+208)/3=152.7\) participants (to 4 significant figures). That is a valid descriptive average only if equal weights are intended; it is not the result implied by the stated probabilities. Here the highest-enrollment scenario receives the largest weight, so the probability-weighted result is higher. Changing the weights changes the summary even though the design and its three scenario-specific expectations have not changed.

Interpretation, Checks, and Alternatives

A scenario-weighted expectation is useful when a single planning summary is needed and the scenario probabilities have a defensible interpretation. For example, the weights might encode a specified distribution over plausible treatment-effect configurations. If the weights are instead chosen merely to summarize a scenario grid, the result should be labelled as a weighted summary under those chosen weights, not as a predictive expectation based on known probabilities.

A weighted mean does not replace the scenario-specific results. Two designs can have the same weighted expected enrollment but behave quite differently across scenarios: one may have similar expectations in every scenario, while another may have very low enrollment in one and very high enrollment in another. Reporting the weights and scenario-specific values makes that variation visible.

Other summaries answer different questions. A single expected enrollment under the global null or global alternative describes that particular configuration, not a mixture. An unweighted mean gives every listed scenario equal influence regardless of its stated probability. A maximum-enrollment figure describes the design’s cap, not its average under a scenario distribution. These summaries are not interchangeable; choose the one that matches the planning question, and show the scenario-specific results when the spread matters.

Diagnostic checks:
  • Confirm that every \(w_s\) is nonnegative and that the weights sum to 1, allowing only minor rounding discrepancy.
  • Check that scenarios do not overlap and that the set of scenarios matches the stated averaging purpose.
  • Use a consistent definition of enrollment and consistent stage increments in every scenario.
  • Verify that continuation probabilities are between 0 and 1 and do not increase across successive stages.
  • Check that the weighted result lies between the smallest and largest scenario-specific expected enrollments.

A frequent mistake is to list several plausible scenarios, assign equal weights without explanation, and call the arithmetic mean the expected enrollment. The correction is to state what the weights mean. Equal weights can be appropriate when the aim is an explicitly equal-weight summary of the listed scenarios; they do not become scenario probabilities merely because they sum to 1. Similarly, assigning probabilities to scenarios is a modelling or planning choice, not evidence that treatment effects are randomly drawn from those configurations in the actual trial.

The weighted mean also has limits. It compresses a range of possible scenario-specific behavior into one number, depends on the chosen scenarios and weights, and can conceal sensitivity to either choice. A useful report therefore gives the weighted expectation together with the scenario definitions, weights, and conditional expectations. If those inputs are uncertain, show how the summary changes under reasonable alternative weights rather than implying unwarranted precision.

Check Your Understanding

Use the stated weights as probabilities for the configurations in each question. Try each calculation before reading the solutions.

  1. Two scenarios have expected enrollments of 100 and 180 participants, with weights 0.75 and 0.25. What is the weighted expected enrollment?
  2. Three scenarios have expected enrollments 80, 140, and 220 participants. A colleague averages them equally. What additional information is needed to interpret that average as a scenario-weighted expectation?
  3. A two-stage design enrolls 50 participants at stage 1 and 70 at stage 2 if reached. Under one scenario, the probability of reaching stage 2 is 0.40. Find the expected enrollment in that scenario.
  4. Scenario-specific expected enrollments are 90 and 210 participants, weighted 0.60 and 0.40. A reported weighted expectation is 150. Is it consistent with the stated inputs?

Solutions. (1) \(0.75(100)+0.25(180)=75+45=120\) participants. (2) To interpret the equal-weight average as a probability-weighted expectation, the scenario probabilities must be justified as equal, and the scenarios must form a mutually exclusive, exhaustive set for the intended mixture. Otherwise it is only an equal-weight descriptive summary: its value is \((80+140+220)/3=146.7\) participants. (3) \(50+70(0.40)=78\) participants. (4) No. The stated inputs give \(0.60(90)+0.40(210)=54+84=138\) participants, not 150. Recheck the weights, scenario expectations, or reported arithmetic.

Clinical Biostats