Why Summarize the Whole Stopping Pattern?
A group sequential design is not described completely by its maximum planned information or its overall power. Under a particular data-generating scenario, some simulated trials may stop at the first look, others may continue, and others may reach the maximum number of looks. The stopping pattern shows where trial paths end; the rejection pattern shows at which analysis the design rejects \(H_0\); and the sample-size distribution describes how much information or how many participants the design uses across paths.
These summaries answer different planning questions. A high probability of stopping early may reduce average sample-size use, but it does not by itself indicate high power: early stopping may be for efficacy or futility. Likewise, a design's average sample size does not say how often it uses close to its maximum. Each summary must be tied to the simulated scenario and to the stopping and analysis rules that generated the paths.
Path-Level Quantities and Their Summaries
Suppose the design has \(K\) planned analyses. For one simulated trial path, let \(J\) be the look at which the trial stops, so \(J\in\{1,\ldots,K\}\). The final look is an administrative stopping point: a path that has not stopped earlier ends at look \(K\). Let \(N_j\) denote the sample size at look \(j\), and let \(N_J\) be the sample size at which that path stops. If the schedule is specified in information rather than participants, the same summaries can be formed using the information at each look.
The conditional stopping probability has a different denominator from the unconditional stopping probability. At look \(j\), only paths that reach that look can stop there. A rejection at look \(j\) is also an event on a path that stops at \(j\); a path cannot reject at multiple looks if the design stops upon rejection. Consequently, analysis-specific rejection events are disjoint, and overall rejection probability is their sum.
These summaries rely on a clearly specified simulation scenario and correctly implemented design. Each simulated path must be generated under the scenario's assumptions, evaluated using the planned looks and stopping rules, and followed through to its actual stopping point. For a simulated frequency to estimate a probability, the paths should be independent replicates from the same specified scenario. Counts should include every completed path exactly once. If several scenarios are simulated, report the summaries separately for each one.
Estimating Stopping and Rejection Probabilities
For \(R\) simulated paths, let \(C_j\) be the number that stop at look \(j\), and let \(A_j\) be the number that reject \(H_0\) at look \(j\). The estimated unconditional stopping probability is \(C_j/R\). The estimated conditional stopping probability is \(C_j/\sum_{\ell=j}^{K}C_\ell\), because the denominator counts paths that reached look \(j\). The estimated analysis-specific rejection probability is \(A_j/R\), not \(A_j/C_j\). The latter is the rejection fraction among paths stopping at that look, a useful conditional description but not the probability of rejection on a randomly selected trial path.
Since each path stops once, the stopping counts must satisfy \(\sum_{j=1}^{K}C_j=R\). Since rejection at look \(j\) requires stopping at look \(j\), \(0\le A_j\le C_j\). The simulated overall rejection probability, or power under an alternative scenario, is \(\sum_{j=1}^{K}A_j/R\). Under a null scenario, the corresponding total estimates the design's Type I error probability for the specified null data-generating conditions.
Worked Example: Estimating the Stopping Pattern
Suppose 10,000 hypothetical paths are simulated under one specified alternative scenario for a design with three looks. The paths stop as follows: 1,800 at look 1, 2,700 at look 2, and 5,500 at look 3.
The estimated unconditional stopping probabilities are \(1{,}800/10{,}000=0.18\), \(2{,}700/10{,}000=0.27\), and \(5{,}500/10{,}000=0.55\). They sum to \(1.00\), as required because every path stops once.
At look 1, all paths reached the look, so the conditional stopping probability is \(1{,}800/10{,}000=0.18\). At look 2, \(10{,}000-1{,}800=8{,}200\) paths reached it, giving \(2{,}700/8{,}200=0.3293\). At look 3, all 5,500 paths reaching it stop there, so the conditional probability is \(5{,}500/5{,}500=1.00\).
Thus, 55% of all simulated paths stopped at the final look, while about 32.9% of paths that reached look 2 stopped at that look. These answer different questions; the conditional percentage at look 2 is not 27%.
Worked Example: Rejection Probability at Each Analysis
For the same hypothetical set of 10,000 paths, suppose rejection occurred at look 1 on 900 paths, at look 2 on 600 paths, and at look 3 on 1,400 paths. These counts do not exceed the corresponding stopping counts.
The analysis-specific rejection probability estimates are \(900/10{,}000=0.09\) at look 1, \(600/10{,}000=0.06\) at look 2, and \(1{,}400/10{,}000=0.14\) at look 3. Adding these disjoint events gives an estimated overall rejection probability of \(0.09+0.06+0.14=0.29\), or 29% under this scenario.
Among paths stopping at each look, the rejection fractions are \(900/1{,}800=0.50\), \(600/2{,}700=0.2222\), and \(1{,}400/5{,}500=0.2545\). For example, the 50% at look 1 is conditional on stopping there. It is not the overall probability that a simulated trial rejects at look 1, which is 9%. The remaining paths stopping at each look did not reject there.
Summarizing the Sample-Size Distribution
If stopping at look \(j\) entails sample size \(N_j\), the simulated distribution of \(N_J\) is obtained by assigning each path the sample size at its stopping look. With a fixed schedule of increasing sample sizes, its mean can be calculated directly from the stopping probabilities. Its median and other quantiles can be found by accumulating probability across the ordered sample sizes. A complete summary often reports the look-specific stopping frequencies, sample sizes at those looks, the mean and median sample size, and the maximum planned sample size. The mean describes average use; the maximum describes the planned upper end, not typical use.
To justify the second expression, write the sample size on a path that stops at \(J\) as successive increments. The increment from \(N_{j-1}\) to \(N_j\) is included precisely when the path reaches look \(j\), or \(J\ge j\). Thus,
Taking expectations and using \(E[\mathbf{1}(J\ge j)]=P(J\ge j)\) yields the stated expression. This form emphasizes that expected sample-size savings depend on how often paths avoid the additional sample required to reach later looks. If sample size varies between paths even at the same look, summarize the actual path-level sample sizes instead of treating \(N_j\) as a fixed value.
Worked Example: Mean and Median Sample Size
Suppose the three looks in the previous examples occur at sample sizes of 100, 200, and 300 participants. Using stopping probabilities 0.18, 0.27, and 0.55, the expected stopping sample size is
The distribution assigns probability 0.18 to 100 participants and a cumulative probability of \(0.18+0.27=0.45\) to sample sizes no greater than 200. The cumulative probability reaches 1.00 at 300, so the median is 300 participants: at least half the probability is at or below 300, while less than half is at or below 200. The mean of 237 participants is below the maximum of 300, but it does not mean that a typical path stops at exactly 237 participants. In this discrete example, no path stops at 237.
Interpreting, Comparing, and Checking the Summaries
Stopping distributions are most informative when the reason for stopping is retained. If efficacy and futility rules can both stop a trial, report the mutually exclusive stopping category assigned to each path, such as efficacy, futility, or reaching the final analysis without an earlier stop. The category definitions must reflect the implemented rule, including what happens if more than one criterion appears to be met at a look. Counts for mutually exclusive categories should add to the number of paths. Do not combine different reasons into a single “early stop” frequency if the distinction matters for the planning question.
Analysis-specific rejection probabilities complement stopping frequencies. A stopping table can show that many paths end early, but only the rejection table indicates how much rejection probability is contributed at each analysis. Comparing scenarios can reveal a shift in the timing of rejection or in average sample-size use, even when their overall rejection probabilities are similar. Comparisons are meaningful only when the scenarios and design rules are stated; a simulated frequency is not a universal property of the design.
Several checks can expose errors in summaries. Confirm that stopping counts sum to \(R\), that each path has one stopping look, that rejection counts do not exceed stopping counts, and that the total rejection probability agrees with the sum of analysis-specific probabilities. Check that conditional stopping denominators count paths reaching the look rather than all paths. Verify that the sample size assigned to each path matches its actual stopping look, and that reported units are clear. A discrepancy may indicate duplicated or missing paths, overlapping stopping categories, an incorrect denominator, or a mismatch between the stopping rule and the summary code.
A common alternative is to report only the maximum sample size or only the expected sample size. The maximum is essential for describing the planned upper limit, but it cannot show how often that limit is reached. The mean gives a compact measure of average use, but conceals the shape of the distribution: two designs can have the same mean and quite different probabilities of early stopping or of using nearly the maximum. Look-specific frequencies and selected quantiles add that context. These summaries do not replace the overall power or Type I error estimate; they describe different features of the design.
All simulation summaries are estimates under the assumed data-generating model and implementation. They do not establish that the assumptions are appropriate for a future trial, and they can be misleading if paths are generated or classified incorrectly. Report the scenario, number of simulated paths, look schedule, stopping categories, and the definition of sample size used. Interpret small differences cautiously, especially when they are not accompanied by the uncertainty inherent in a finite simulation.
Check Your Understanding
Answer each question using the definitions and calculations in this tutorial. Solutions follow the questions.
- In 2,000 simulated paths, 300 stop at look 1 and 500 stop at look 2. What is the estimated conditional probability of stopping at look 2, given that a path reaches it?
- Of 5,000 paths, 150 reject at look 1 and 250 reject at look 2. Assuming no path rejects at both looks, what is the estimated overall rejection probability? What does the look-2 value \(250/5{,}000\) estimate?
- A design stops at sample sizes 80, 160, and 240 with probabilities 0.20, 0.30, and 0.50. Calculate expected sample size.
- Why is the fraction rejecting among paths that stop at look 1 not the same quantity as the overall probability of rejection at look 1?
Solutions. (1) The number reaching look 2 is \(2{,}000-300=1{,}700\). The conditional stopping probability is \(500/1{,}700=0.2941\), or about 29.4%. (2) The overall rejection estimate is \(150/5{,}000+250/5{,}000=0.03+0.05=0.08\), or 8%. The look-2 value estimates the unconditional probability of rejection at look 2 among all simulated paths. (3) \(80(0.20)+160(0.30)+240(0.50)=16+48+120=184\) participants. (4) The fraction among look-1 stoppers uses only paths that stopped there as its denominator; the overall probability uses all simulated paths. It is conditional rather than unconditional.