Why Stress-Test a Re-Estimation Rule?
A sample-size re-estimation rule is a prespecified mapping from interim information to a final sample size. Its operating characteristics depend not only on the rule itself, but also on the assumptions used to generate the data. If the true outcome variance is larger than anticipated, events accrue more slowly than expected, or follow-up is shorter, a rule that performs well under its design assumptions may miss its power or feasibility targets.
Stress-testing uses targeted simulation to examine those vulnerabilities. The goal is not to search for one “correct” scenario. It is to learn which assumptions drive the rule’s performance, how severe plausible departures are, and what observable signs would indicate that the planned adaptation is inadequate. This complements the earlier material on simulating sample-size re-estimation rules and evaluating their operating characteristics: here the central task is choosing informative departures from the assumed data-generating conditions.
For each scenario, apply the complete planned rule to simulated trials, including the interim analysis, sample-size update, caps or floors, and final analysis. Summarize outcomes relevant to the design: probability of meeting the power objective under a stated treatment effect, probability of meeting an event target, achieved sample size, calendar duration, and other prespecified operating characteristics. Under a null treatment effect, rejection probability can also be recorded as a safety diagnostic; assessing or calibrating it is a separate task.
Designing Informative Simulation Scenarios
Begin with a reference scenario that represents the assumptions used to design the trial. Then create a small set of departures chosen because they are scientifically or operationally plausible. Change one factor at a time to isolate its influence, and include selected combined departures to assess whether risks compound. A stress test is not useful merely because it is extreme; the scenario should be interpretable and relevant to decisions about the rule.
For every scenario, specify the data-generating model and the rule being evaluated. For example, a continuous endpoint might be generated from independent normal outcomes in two randomized groups with common variance. A survival endpoint might be generated under proportional hazards, with an explicit event-time distribution, censoring mechanism, accrual pattern, and administrative end date. If these assumptions are not credible, the simulation may accurately evaluate the wrong model.
Use enough replications to make Monte Carlo uncertainty acceptably small. If \(R\) simulated trials produce an estimated probability \(\hat{q}\) of an outcome, its Monte Carlo standard error is \(SE_{\mathrm{MC}}=\sqrt{\hat{q}(1-\hat{q})/R}\). This quantifies simulation noise, not uncertainty about which scenario is realistic. At \(R=20{,}000\) and \(\hat{q}=0.80\), for example, the Monte Carlo standard error is about \(0.00283\), or 0.28 percentage points.
Worked Example: Stress-Testing a Variance-Based Rule
Suppose a hypothetical two-arm trial uses a continuous outcome, equal allocation, a two-sided 5% test, and 80% target power for a mean difference of 5 outcome units. Outcomes are independent and normally distributed with common standard deviation \(\sigma\). The rule estimates the variance at interim and sets the final sample size per group to the usual normal-approximation value, subject to a maximum of 140 per group. Assume the interim estimate is sufficiently precise for this illustration.
Under the planning value \(\sigma=12\), the required sample size per group is
so the rounded requirement is 91 participants per group. To stress-test variance misspecification, keep the treatment difference at 5 units but generate outcomes with \(\sigma=15\). If the interim estimate reflects that variance, the unrounded updated requirement is
The cap therefore limits enrollment to 140 per group. At that sample size, the standard error of the estimated difference is \(\,15\sqrt{2/140}=1.793\) units, and the noncentrality parameter is \(5/1.793=2.788\). The approximate two-sided power is \(\Phi(2.788-1.960)+\Phi(-2.788-1.960)\approx0.796\).
For a hypothetical simulation of 20,000 trials, suppose the estimated power is 0.795 with Monte Carlo standard error \(\sqrt{0.795(0.205)/20{,}000}=0.00285\). This is consistent with the calculation: the cap leaves power just below 80%. The key diagnostic is not simply that the rule increased sample size, but that a plausible variance departure pushes the required sample above the cap.
Events, Censoring, and Accrual Are Different Stressors
For an event-driven survival design, event count is often the information target. With equal allocation, proportional hazards, a log-rank test, a two-sided significance level \(\alpha\), and target power \(1-\beta\), a common approximation for required events is
where \(HR\) is the assumed hazard ratio. Translating \(D\) into enrollment requires assumptions about the probability of an observed event during the relevant follow-up. A lower observed event rate can arise from a lower underlying event incidence, more censoring, shorter follow-up, or some combination. Those mechanisms can have different implications for interpretation and duration, so they should not be collapsed into one unexplained “low event rate” scenario.
Worked Example: Lower Event Probability and a Sample-Size Cap
Consider a hypothetical event-driven trial with equal allocation, proportional hazards, a two-sided 5% test, 80% power, and a target \(HR=0.70\). Assume independent event times and non-informative censoring for the standard event-count approximation. The required number of events is
or 247 events. If the assumed probability of an observed event per enrolled participant is 0.40, approximately \(247/0.40=617.5\), or 618 participants, are needed. Now stress-test a lower event probability of 0.25, holding the hazard ratio and other assumptions fixed. The uncapped enrollment requirement becomes \(247/0.25=988\). With an 800-participant cap, the expected event count is only \(800(0.25)=200\).
Using the same approximation, the noncentrality parameter at 200 events is \(|\log(0.70)|\sqrt{200}/2=2.522\). Approximate power is therefore \(\Phi(2.522-1.960)+\Phi(-2.522-1.960)\approx0.713\), below the 80% target. In a hypothetical 20,000-trial simulation, suppose 71.4% of trials reject under this alternative; the Monte Carlo standard error is about 0.0032. That result is consistent with the approximation and indicates that the cap is consequential under this scenario.
Next, simulate the mechanisms separately. A lower underlying event incidence may reduce events even with complete follow-up. More censoring or shorter follow-up can also reduce observed events, but may additionally delay reaching the target. The simulation should represent the actual event-time and censoring processes; the event probability alone does not determine calendar duration.
Accrual is another distinct dimension. If a design requires 618 participants, uniform accrual at 40 participants per month implies about \(618/40=15.45\) months of enrollment. At 25 per month, enrollment takes about \(618/25=24.72\) months. If every participant then receives 12 months of follow-up, the corresponding end-to-end durations are about 27.45 and 36.72 months. These calculations assume steady accrual and complete follow-up; a simulation should include pauses, varying recruitment rates, or an administrative study end when those are plausible. With complete follow-up and the same event process per participant, slower accrual may lengthen duration without changing eventual statistical power. A fixed calendar deadline can change both by preventing the planned follow-up or event target from being reached.
Treatment-Effect Scenarios and Interpretation
A sample-size rule that updates a nuisance parameter, such as variance, does not thereby learn the true treatment effect. Stress-testing should therefore distinguish uncertainty in nuisance parameters from uncertainty about the effect used to define the design objective. Under a smaller true effect, a rule that retains its original target effect may select the same sample size while the probability of success falls. This is a sensitivity analysis of the design’s power, not a reason to change the primary treatment question after seeing interim results.
Worked Example: A Smaller True Treatment Effect
Return to the continuous-outcome rule with planning difference \(\delta=5\), standard deviation 12, and no binding cap at the planned value. The rule selects 91 participants per group when the interim variance estimate is near 144. Now simulate a scenario in which the true difference is 4 units, while the variance remains 144. Assume the variance estimator is based on within-group variation and that the rule does not re-estimate the treatment effect.
The standard error at 91 per group is \(12\sqrt{2/91}=1.779\) units. Under the true difference of 4 units, the noncentrality parameter is \(4/1.779=2.249\), giving approximate two-sided power \(\Phi(2.249-1.960)+\Phi(-2.249-1.960)\approx0.614\). The rule still selects about 91 per group because the variance and its assumed target difference have not changed. A hypothetical simulation estimate of 0.615 power across 20,000 trials has Monte Carlo standard error about 0.0034 and agrees with the approximation.
This scenario reveals a limitation, not a computational defect: the adaptation addresses variance uncertainty, while power is sensitive to the true effect. Report results across a scientifically meaningful range of effects rather than treating the design effect as known with certainty.
Diagnostics, Comparisons, and Limitations
Useful stress-test summaries include the distribution of final sample size, the fraction of trials hitting the cap or floor, achieved power, event shortfall, and calendar duration. Inspect more than averages: a mean sample size can conceal a substantial fraction of trials at a cap, and an average duration can conceal a long tail. For each simulated performance probability, report Monte Carlo uncertainty; for duration or sample size, examine quantiles as well as the mean.
A warning sign is a performance change that cannot be explained by the scenario’s mechanism. For example, a variance-only change should not alter the simulated treatment effect, while an accrual-only change should not alter eventual power if every participant receives the same complete follow-up and the event process is unchanged. If it does, inspect the simulation code for unintended changes such as a calendar cutoff or altered censoring. Also check that each simulated trial applies the interim rule at the correct time and respects its sample-size bounds.
One-factor-at-a-time scenarios are easy to interpret but can miss interactions; combined scenarios reveal whether, for example, slow accrual and heavier censoring jointly make an event target infeasible. Analytic power and sample-size calculations are useful checks for simple models, but they cannot represent complex interim rules, caps, variable accrual, and censoring as flexibly as simulation. Conversely, simulation is only as credible as its data-generating assumptions. It does not establish that a scenario is plausible, and it cannot replace clinical judgment about meaningful effects or operational constraints.
A common mistake is to conclude that a rule is robust because it meets its target in the reference scenario. That reasoning tests only the assumptions under which it was designed. The correction is to prespecify plausible departures, evaluate each with the full rule, and report where performance degrades or feasibility fails. Another mistake is to read a difference between two simulated success rates as meaningful without considering Monte Carlo error; increase the number of replications or report the uncertainty before interpreting small differences.
Check Your Understanding
Answer each question, then compare your reasoning with the solutions.
- Why should a stress test change one assumption at a time before examining combined departures?
- A simulated success probability is 0.70 from 10,000 trials. Calculate its Monte Carlo standard error.
- In an event-driven trial, why can lower incidence and heavier censoring produce different operational consequences even if both reduce the observed event count?
- A variance-based rule selects the same sample size when the true treatment effect is smaller than planned. Does that show the rule is working as intended, or that the simulation failed?
- What should be examined in addition to the average final sample size when a rule has a maximum enrollment cap?
Solutions. (1) Isolating factors helps attribute changes in performance to a particular assumption; combined scenarios can then assess interactions. (2) \(SE_{\mathrm{MC}}=\sqrt{0.70(0.30)/10{,}000}=0.00458\), or about 0.46 percentage points. (3) Censoring and follow-up affect when events become observable and may delay reaching the target, while lower incidence changes event occurrence; the time course and mechanism matter. (4) It can be an expected result if the rule updates variance but not treatment effect; power can fall because the true effect is smaller. (5) Examine the fraction of trials hitting the cap, the distribution and upper quantiles of sample size, and achieved power under the stress scenarios.