Docs/Group Sequential Design

Group Sequential Design (GSD)

Technical documentation for interim monitoring with early stopping rules. This page covers the error spending framework, conditional power, futility monitoring, DMC considerations, regulatory alignment, spending function selection, and validation benchmarks against industry-standard software.

1. Theoretical Foundation

Group Sequential Designs (GSD) allow for the interim monitoring of a clinical trial to permit early stopping for efficacy (success) or futility (failure). Zetyra utilizes the α\alpha-spending function approach (Lan & DeMets, 1983), which maintains the overall family-wise error rate (FWER) while allowing for flexible timing of interim analyses.

The Error Spending Framework

To prevent the inflation of the Type I Error rate, Zetyra spends a portion of the total significance level (α\alpha) at each look. The boundaries are computed by Zetyra's own Python implementation of the Lan-DeMets error-spending framework, which is validated against the peer-reviewed gsDesign R package (Anderson, 2023), widely used in regulatory submissions. gsDesign is the external reference the engine is checked against, not a dependency it calls: see the validation appendix below for the agreement achieved on each design.

Alpha-Spending Function Definition

For information fraction t[0,1]t \in [0, 1], the spending functionα(t)\alpha^*(t) defines cumulative Type I error spent:

α(t)={22Φ(z1α/2/t)O’Brien-Fleming (Lan-DeMets)αlog ⁣[1+(e1)t]Pocock-like (Lan-DeMets)α1eγt1eγHwang-Shih-DeCaniαtLinear / uniform spending\alpha^*(t) = \begin{cases} 2 - 2\Phi(z_{1-\alpha/2}/\sqrt{t}) & \text{O'Brien-Fleming (Lan-DeMets)} \\ \alpha \cdot \log\!\left[1 + (e - 1)t\right] & \text{Pocock-like (Lan-DeMets)} \\ \alpha \cdot \frac{1 - e^{-\gamma t}}{1 - e^{-\gamma}} & \text{Hwang-Shih-DeCani} \\ \alpha \cdot t & \text{Linear / uniform spending} \end{cases}

Key principle: By carefully allocating how much of the total α=0.05\alpha = 0.05 is “spent” at each interim analysis, we can perform multiple hypothesis tests while maintaining the overall Type I error rate. The Lan-DeMets framework allows flexible timing—actual analysis times need not match planned times.

2. Conditional Power & Futility Monitoring

Conditional Power Definition

At interim analysis kk, conditional power (CP) is the probability of achieving statistical significance at the final analysis, given the data observed so far.

CPk(θ)=P(ZK>z1αKZk=zk,θ)CP_k(\theta) = P\left(Z_K > z_{1-\alpha_K} \mid Z_k = z_k, \theta\right)

Where ZkZ_k is the observed test statistic at look kk, and θ\theta is the assumed treatment effect for the remaining patients.

Effect Size Assumptions for CP

Under Current Trend

θ=θ^k\theta = \hat{\theta}_k (MLE from interim data). Most commonly used for futility assessment.

Under Design Alternative

θ=θ1\theta = \theta_1 (original design assumption). Conservative for futility decisions.

Under Null

θ=0\theta = 0. Used to verify Type I error control; rarely used for decision-making.

Futility Boundary Types

Non-Binding Futility (Default)

The trial may continue even if the futility boundary is crossed. Type I error is calculated assuming the trial always continues.

  • Advantage: Preserves nominal Type I error regardless of decision
  • Advantage: Provides flexibility for DMC judgment
  • Disadvantage: Slightly larger sample size than binding

Binding Futility

The trial must stop if the futility boundary is crossed. Type I error accounts for mandatory stopping.

  • Advantage: Smaller sample size (uses alpha “saved” from futility)
  • Disadvantage: Less flexibility—must stop even if external evidence changes
  • Disadvantage: Rarely recommended by FDA for pivotal trials

DMC Decision Point: When to Stop for Futility

Common thresholds in practice:

  • CP < 5% under current trend: Strong evidence of futility
  • CP < 10% under design alternative: Futility if even optimistic assumptions fail
  • CP < 20%: Consider stopping if operational costs are high

Note: These are guidelines, not rules. DMC should consider clinical context, safety data, and external evidence.

Sample Size Re-estimation (SSR) Considerations

Zetyra Does NOT Support Adaptive SSR

Zetyra's GSD calculator provides classical group sequential designs with fixed maximum sample sizes. Blinded or unblinded sample size re-estimation (SSR) based on interim variance or effect estimates requires additional statistical methodology:

  • Blinded SSR: Re-estimate pooled variance without unblinding—generally acceptable but requires pre-specification
  • Unblinded SSR: Requires combination tests (e.g., inverse normal method) or conditional error functions to maintain Type I error
  • Promising Zone designs: Proschan & Hunsberger (1995) methodology not currently implemented

Reference: ICH E20 Adaptive Designs for Clinical Trials (2025).

3. Statistical Assumptions & Requirements

Prespecification (Mandatory)

Per ICH E9 (1998) and ICH E20 (2025), the following must be documented in the Statistical Analysis Plan (SAP) before unblinding:

  • • Number of interim analyses
  • • Planned timing (information fractions or calendar time)
  • • Alpha-spending function with all parameters
  • • Efficacy and futility boundaries (binding vs. non-binding)
  • • Decision rules and DMC charter references

Independent Increments

Test statistics at successive looks must follow a multivariate normal distribution with independent increments:

(Z1,Z2,,ZK)N(θ(I1,I2,,IK),Σ)(Z_1, Z_2, \ldots, Z_K) \sim N\left(\theta \cdot (\sqrt{I_1}, \sqrt{I_2}, \ldots, \sqrt{I_K}), \Sigma\right)

where Σjk=Ij/Ik\Sigma_{jk} = \sqrt{I_j / I_k} for jkj \leq k.

Practical Implication: This holds when subjects are randomized independently and outcomes are measured without bias. Violations can occur with:

  • • Time-varying treatment effects
  • • Informative censoring (survival endpoints)
  • • Cluster randomization without proper adjustment

Information Time vs. Calendar Time

Analyses should be timed by information fraction, not calendar time:

tk=IkIK=nkN (for continuous endpoints)t_k = \frac{I_k}{I_K} = \frac{n_k}{N} \text{ (for continuous endpoints)}

For survival endpoints, information fraction is based on number of events:tk=dk/Dt_k = d_k / D where dkd_k is observed events and DD is target events.

4. Spending Function Selection Guide

The choice of spending function affects how aggressively alpha is allocated to early looks. This decision should balance statistical power, expected sample size savings, and regulatory acceptability.

FunctionEarly StoppingPower ImpactWhen to Choose

O'Brien-Fleming

α(t)=22Φ(z1α/2/t)\alpha^*(t) = 2 - 2\Phi(z_{1-\alpha/2}/\sqrt{t})

Conservative (Hard)

<0.1% alpha at 50% info

Minimal (~1-2%)Default for Phase III. Preserves power for final analysis while allowing early stopping only for overwhelming efficacy. FDA's preferred choice.

Pocock

α(t)=αt\alpha^*(t) = \alpha \cdot t

Aggressive (Easy)

~2.5% alpha at 50% info

LowerTime-critical trials. When early termination is high priority (e.g., pandemic response, diseases with rapid progression). Requires larger max N.

Hwang-Shih-DeCani

γ=4\gamma = -4 (OBF-like)

ConservativeMinimalApproximately O'Brien-Fleming behavior with closed-form formula.

Hwang-Shih-DeCani

γ=1\gamma = 1 (Pocock-like)

AggressiveModerateApproximately Pocock behavior with closed-form formula.

Hwang-Shih-DeCani

γ=0\gamma = 0

LinearModerateCompromise option. Linear spending between OBF and Pocock. Good when moderate early stopping is desired.

“Lower” power means Pocock requires a larger maximum sample size than O'Brien-Fleming to achieve the same power. At the same max N, Pocock has lower power.

Recommendation: For most Phase III confirmatory trials, O'Brien-Fleming is the standard choice. It preserves statistical power while allowing early stopping only when treatment effects are substantially larger than planned. Pocock should be reserved for situations where early stopping has exceptional operational or ethical value.

Limitations & When Not to Use GSD

Group sequential designs add operational complexity. Consider whether the benefits outweigh the costs in your specific trial context.

Small Trials (N < 100)

GSD overhead (DMC meetings, unblinding logistics) may exceed sample size savings. With small N, the expected number of subjects saved is often <10.

Very Short Trials (<6 months)

By the time you organize a DMC meeting and analyze interim data, enrollment may be nearly complete. The operational delay can negate early stopping benefits.

Long-Term Endpoints (e.g., 5-year survival)

Early looks occur when most patients haven't reached the endpoint. Information fraction lags calendar time significantly, limiting utility of early interims.

Regulatory Complexity

Some regulatory pathways (especially first-in-class drugs) may scrutinize GSD more heavily. If your trial already faces regulatory hurdles, a fixed design may be simpler to defend.

Rule of thumb: GSD is most beneficial when (1) the trial is large enough that even a 10-20% reduction in expected N is meaningful, (2) enrollment is slow enough that interim analyses can influence the trial, and (3) there's genuine uncertainty about treatment effect magnitude.

5. Handling Timing Deviations

One key advantage of the Lan-DeMets error spending approach is its flexibility in handling deviations from planned interim analysis timing. Unlike the original Pocock and O'Brien-Fleming group sequential boundaries (which require equal spacing), the spending function approach maintains Type I error control even when actual information fractions differ from planned values.

Recalculating Boundaries at Actual Information Times

When the actual information fraction tkactualt_k^{actual} differs from planned tkplannedt_k^{planned}:

  1. 1. Calculate α(tkactual)\alpha^*(t_k^{actual}) using the pre-specified spending function
  2. 2. Compute incremental alpha: αk=α(tkactual)α(tk1actual)\alpha_k = \alpha^*(t_k^{actual}) - \alpha^*(t_{k-1}^{actual})
  3. 3. Derive boundary zkz_k to spend exactly αk\alpha_k given correlation structure

Acceptable Deviations

  • • Interim at 48% vs. planned 50% information
  • • Slight delays due to enrollment variability
  • • Event-driven analyses arriving early/late

The spending function automatically adjusts boundaries to maintain Type I error control.

Problematic Deviations

  • • Skipping a planned interim entirely
  • • Adding unplanned interim analyses
  • • Changing the spending function after unblinding

These require protocol amendments and may raise regulatory concerns.

Documentation Requirement

When actual timing differs from planned, document in the interim analysis report: (1) the reason for the deviation, (2) the actual information fraction achieved, (3) the recalculated boundary using the pre-specified spending function, and (4) confirmation that the spending function was not modified.

6. DMC Perspective

The Data Monitoring Committee (DMC) reviews unblinded interim data and makes recommendations about trial continuation. While the statistical boundaries provide quantitative guidance, DMC decisions involve broader considerations.

What the DMC Sees at Each Interim Analysis

Efficacy Assessment

  • • Test statistic (Z-score) vs. efficacy boundary
  • • Point estimate and confidence interval
  • • Conditional power under current trend and design alternative
  • • Predicted probability of success at final analysis

Futility Assessment

  • • Test statistic vs. futility boundary (if specified)
  • • Conditional power (if <10-20%, strong futility signal)
  • • Trend in treatment effect across time
  • • Comparison to external trials or historical data

Possible DMC Recommendations

Continue as Planned

Data are within expected range; no boundary crossed. This is the most common outcome at early interims.

Stop for Efficacy

Efficacy boundary crossed with strong treatment effect. Consider external validity, subgroup consistency, and safety profile before recommending.

Stop for Futility

Conditional power extremely low (<5-10%). Continuing would expose patients to experimental treatment with minimal scientific benefit.

Stop for Safety

Unacceptable adverse event profile. This decision is independent of efficacy boundaries and takes precedence.

Beyond the Boundaries: Clinical Judgment

Statistical boundaries are guidelines, not mandates. The DMC may consider:

  • External evidence: Results from competing trials or emerging safety signals
  • Subgroup heterogeneity: Is the effect driven by a specific population?
  • Clinical meaningfulness: Is the effect size clinically relevant, even if statistically significant?
  • Regulatory context: Would early stopping be accepted by the target regulatory agency?

7. Practical Considerations

Enrollment Overrun

In fast-enrolling trials, patients may be randomized between the data cutoff for an interim analysis and the DMC meeting/decision. This “pipeline” enrollment must be considered.

Handling Overrun

  • Option 1: Include in analysis — Analyze all randomized patients. Information fraction may exceed planned interim (e.g., 55% vs. 50%). The spending function handles this automatically.
  • Option 2: Exclude pipeline patients — Analyze only patients enrolled before cutoff. Cleaner information fraction but may raise questions about selection.
  • Recommendation: Pre-specify the approach in the DMC charter. Option 1 is generally preferred for intent-to-treat integrity.

Boundary Conversion: Z-scores, P-values, and Confidence Intervals

Boundaries can be expressed in multiple equivalent forms:

Z-score Boundary

Z = 2.963

Raw test statistic threshold

Nominal P-value

p = 0.0030

p=2(1Φ(z))p = 2(1 - \Phi(z)) (two-sided)

Repeated CI

99.70% CI

(1p)×100%(1 - p) \times 100\% coverage

Important: The “nominal p-value” at an interim is NOT directly comparable to the conventional 0.05 threshold. Always compare to the boundary p-value for that specific interim analysis.

Alpha Allocation for Efficacy vs. Futility

When designing asymmetric boundaries with different spending functions for efficacy and futility:

  • Efficacy: Typically O'Brien-Fleming to preserve power
  • Futility: Often more aggressive than the efficacy rule, to enable earlier stopping when the trial is unlikely to succeed. In this calculator the aggressive choice is Pocock beta spending (the default when efficacy uses O'Brien-Fleming); HSD is available at a user-selected γ\gamma, and at its defaultγ=4\gamma = -4 it is close to O'Brien-Fleming in shape and therefore stops for futility far less readily
  • Non-binding futility: Does not “spend” alpha; Type I error calculated assuming trial continues regardless of futility boundary

Number of Interim Analyses

# LooksTypical UseInflation, efficacy-only (OBF)Considerations
2Small trials, limited budget for DMC~1%Simple; limited opportunity for early stopping
3Standard Phase III~2%Balance of flexibility and simplicity
4-5Large, long trials; safety monitoring~3%More operational burden; diminishing returns
ContinuousAlpha spending only (no fixed schedule)~4%Maximum flexibility but complex operations

These figures are for efficacy-only monitoring. They are the cost of spending α across looks and nothing else. Adding a futility boundary costs considerably more, because a trial that stops for futility forfeits its chance to reject and that power has to be bought back with a larger maximum N. For the calculator's default O'Brien–Fleming efficacy with Pocock futility (α = 0.025 one-sided, 90% power) the inflation is:

  • 2 looks — +15.4% (vs ~1% efficacy-only)
  • 3 looks — +22.3% (vs ~2%)
  • 4 looks — +26.3% (vs ~3%)
  • 5 looks — +28.9% (vs ~4%)

This matches gsDesign test.type=4 to within ~0.002. The trade is deliberate: you pay for the maximum N up front and recover it on average, since expected N under H0 falls to roughly 48% of the maximum for a 3-look design.

8. Regulatory Documentation Checklist

Per ICH E20 Adaptive Designs for Clinical Trials (2025), the following must be pre-specified in your SAP:

Number and timing of analyses

Fixed or information-based.

Example SAP text: "Three interim analyses at 33%, 67%, and 100% of target events (200, 400, and 600 events respectively)."

Alpha-spending function

Including all parameters.

Example SAP text: "O'Brien-Fleming spending function for efficacy (α = 0.025 one-sided). Pocock spending for futility (non-binding)."

Stopping boundaries

Both efficacy and futility, if applicable.

Example SAP text: "Efficacy boundaries: Z = 3.710, 2.509, 1.994 at looks 1-3. Futility boundary: Stop if conditional power < 10% under current trend at any interim."

Binding vs. non-binding futility

Explicitly state commitment level.

Example SAP text: "Futility boundaries are non-binding. The Type I error rate is calculated assuming the trial continues regardless of futility boundary crossing."

Decision rules

Criteria for each possible outcome.

Example SAP text: "If efficacy boundary crossed: DMC recommends early stopping. If futility boundary crossed: DMC may recommend stopping but is not required."

DMC Charter reference

Who makes interim decisions and with what information.

Example SAP text: "Per DMC Charter v2.0 (dated XX/XX/XXXX), the DMC will review unblinded efficacy and safety data. Sponsor will receive blinded recommendations only."

Handling of timing deviations

Procedure if actual timing differs from planned.

Example SAP text: "If actual information fraction differs from planned by more than 10%, boundaries will be recalculated using the pre-specified spending function evaluated at actual information times."

Final analysis adjustment

Specification for p-values and confidence intervals.

Example SAP text: "Final p-values will be calculated using the stagewise ordering. Confidence intervals will use the repeated confidence interval approach of Jennison & Turnbull."

9. Validation Appendix: GSD Benchmarking

Boundaries and sample sizes are validated against gsDesign (R), the open-source reference implementation. Each comparison is like-for-like: the same boundary construction, the same tailwise α\alpha, the same information timing, and the same futility regime. Where a construction has no direct counterpart it is reported separately rather than as a deviation. Cases 1–2 validate the UPPER boundary of the two-sided α=0.05\alpha=0.05 design against test.type=1 at α=0.025\alpha=0.025; symmetry supplies the lower boundary. Both tools compute max N via MVN power integration (bisection search for the smallest N where power meets the target). Zetyra's classical Pocock construction is compared with gsDesign(sfu="Pocock"); sfLDPocock is shown separately as a different design, not scored as a deviation. After matching constructions, residual differences are assessed against prespecified tolerances without assuming a conservative direction.

Case 1: O'Brien-Fleming (3-Look, Efficacy Only)

Inputs: α=0.05\alpha=0.05 (two-sided), Power=0.80, 1:1 Allocation, Continuous Endpoint, δ=0.25\delta=0.25, σ=1.0\sigma=1.0.

Zetyra offers two O'Brien-Fleming algorithms. The Lan-DeMets variant (default) uses the OBF spending function with per-look root-finding, matching gsDesign sfLDOF. The Classical variant uses the parametric form bk=c/tkb_k = c/\sqrt{t_k} from O'Brien & Fleming (1979).

Lan-DeMets OBF Spending (default)

LookInfo FracEfficacy ZNominal pCumulative α\alpha
133%3.7100.00020.0002
267%2.5090.01210.0121
3100%1.9940.04610.0500

Classical Parametric OBF

LookInfo FracEfficacy ZNominal pCrossing Prob
133%3.4690.00050.0005
267%2.4530.01420.0147
3100%2.0030.04520.0500

“Crossing Prob” is the cumulative boundary-crossing probability under H0, computed from the actual c/tc/\sqrt{t} boundaries. This differs from the Lan-DeMets spending schedule shown in the table above because the classical form does not follow a spending function — it controls total α\alpha directly.

Cross-Software Comparison (Max N Total):

AlgorithmgsDesign (R)ZetyraDeviation
Lan-DeMets sfLDOF5095090 (exact)
Classical (c/tc/\sqrt{t})512

gsDesign does not implement the classical parametric OBF form. Both tools use MVN power integration for sample size, giving near-identical results for shared spending functions.

Case 2: Pocock (3-Look, Efficacy Only)

Inputs: Same as Case 1. Classical constant-boundary Pocock (Pocock 1977), not Lan–DeMets Pocock spending.

LookInfo FractionEfficacy ZNominal pCumulative α\alpha
133%2.28950.0220500.022050
267%2.28950.0220500.037936
3100%2.28950.0220500.049997

Cross-Software Comparison

ConstructionBoundariesMax N
Zetyra — classical constant-boundary Pocock2.2895 (constant)586
gsDesign sfu="Pocock" — like-for-like2.289478 (constant)586
gsDesign sfLDPocock — different construction, context only2.279428, 2.294910, 2.295939588

Primary result: ΔN = 0 subjects; ΔZ < 0.0001 against the like-for-like specification. Fixed-design N = 503.

Zetyra's Pocock efficacy option uses the classical constant-boundary construction. The like-for-like gsDesign specification, sfu="Pocock", gives the same rounded boundary and the same integer maximum sample size: 586. sfLDPocock is a different, Pocock-like Lan–DeMets spending construction; it produces nonconstant boundaries and max N 588. That result is shown separately and is not treated as a numerical deviation.

Case 3: Hwang-Shih-DeCani with Futility (4-Look, One-Sided)

Inputs: α=0.025\alpha=0.025 (one-sided — the same upper-tail efficacy error as 0.05 two-sided, though the design itself is one-sided), Power=0.90, 1:1 Allocation, Continuous Endpoint, δ=0.25\delta=0.25, σ=1.0\sigma=1.0. Efficacy and futility both HSD γ=4\gamma=-4 (non-binding), which is the default; γ\gamma is now a user control, so this case pins it explicitly for reproducibility. This case was previously documented as two-sided with a separate gamma = -2 futility function, which the calculator cannot produce: a two-sided test already rejects in the lower tail and so carries no futility boundary. Restated one-sided so the numbers are reproducible.

LookInfo FractionEfficacy ZFutility ZCum. α\alpha (Eff.)Cum. β\beta (Fut.)
125%3.1554-1.06880.0008010.003206
250%2.81830.01690.0029800.011920
375%2.43911.00410.0089020.035609
4100%2.01362.01360.0250000.100000

Cross-Software Comparison (Max N Total):

gsDesign (R)

703

Zetyra

703

Deviation

0 (exact)

Zetyra and gsDesign agree here on the integer maximum sample size and on every boundary to under 0.002 Z units, within the prespecified validation tolerance.

Note: Futility boundaries are shown as Z-scores and are derived by spending the β-spending function under the alternative, so for typical designs they are positive at interim looks (e.g. Z = 0.377, 1.279 for a 3-look OBF/Pocock design): the trial stops for futility when the observed effect is present but too small to be on track. Only designs with many looks push the first bound slightly negative (a 5-look design starts at −0.24), where so little information has accrued that almost no result is discouraging enough to stop. At the final look the futility and efficacy bounds meet — reaching the end without crossing efficacy isthe futility outcome.

Interpretation of Cross-Software Comparisons

  • Construction matching (Case 2): Zetyra's Pocock is the classical constant-boundary form, so it is compared with gsDesign(sfu="Pocock") — the same construction. sfLDPocock is a different, Pocock-like spending design; it is reported separately for context and is not scored as a deviation.
  • Numerical tolerance (all cases): multivariate-normal probabilities are computed by adaptive integration, so results can differ in the last digit between implementations. Any residual difference reported above is within the prespecified validation tolerance; no direction is claimed for it.

Residual differences are numerical, not methodological. They are reported against the prespecified tolerance rather than characterised as safe in a particular direction; a design whose numbers are treated as conservative by assumption is one whose discrepancies stop being investigated.

11. API Quick Reference

POST /api/v1/calculators/group-sequential

Key Parameters

ParameterTypeDescription
kintNumber of analyses/looks (default: 3)
spending_functionstring"OBrienFleming" | "Pocock" | "HwangShihDecani"
efficacy_hsd_gammafloatHSD shape parameter for efficacy, range [−10, 10] (default: −4). Consulted only when spending_function is "HwangShihDecani". Negative spends conservatively early, 0 spends linearly, positive spends aggressively early.
futility_hsd_gammafloat | nullHSD shape parameter for futility. null (default) = inherit efficacy_hsd_gamma. An explicit value, including 0, overrides the inherited one. Because the default inherits, changing the efficacy γ also changes the futility rule unless this is set explicitly.
beta_spending_functionstring | nullFutility spending. null = auto (Pocock when efficacy uses OBF). "none" = efficacy-only design: no futility boundary, and none of the sample-size inflation that paying for futility stopping costs (~22% for a 3-look OBF design).
effect_sizefloatStandardized effect size (default: 0.3)
alpha, powerfloatDefaults: 0.025, 0.90
test_typestring"one_sided" | "two_sided"

Key Response Fields

  • max_sample_size — Maximum sample size needed
  • expected_sample_size — Expected sample under H₁
  • boundaries.efficacy — Z-score boundaries for efficacy stopping
  • boundaries.futility — Z-score boundaries for futility stopping
  • stopping_probabilities — Cumulative stopping probabilities
View full API documentation →

12. References

  1. Pocock SJ. Group sequential methods in the design and analysis of clinical trials. Biometrika. 1977;64(2):191-199.
  2. O'Brien PC, Fleming TR. A multiple testing procedure for clinical trials. Biometrics. 1979;35(3):549-556.
  3. Lan KKG, DeMets DL. Discrete sequential boundaries for clinical trials. Biometrika. 1983;70(3):659-663.
  4. Kim K, DeMets DL. Design and analysis of group sequential tests based on the type I error spending rate function. Biometrika. 1987;74(1):149-154.
  5. Hwang IK, Shih WJ, De Cani JS. Group sequential designs using a family of type I error probability spending functions. Statistics in Medicine. 1990;9(12):1439-1445.
  6. Jennison C, Turnbull BW. Group Sequential Methods with Applications to Clinical Trials. Chapman & Hall/CRC Press; 2000.
  7. Anderson K. gsDesign: Group Sequential Design. R package version 3.6.4; 2024.
  8. Wassmer G, Pahlke F. rpact: Confirmatory Adaptive Clinical Trial Design and Analysis. R package version 4.0; 2024.
  9. U.S. Food and Drug Administration. Adaptive Designs for Clinical Trials of Drugs and Biologics: Guidance for Industry. November 2019.
  10. International Council for Harmonisation (ICH). E9 Statistical Principles for Clinical Trials. February 1998.

Last updated: July 2026

Ready to design?

Use our Group Sequential Design Calculator for interim analyses with O'Brien-Fleming, Pocock, and HSD spending functions, sample-size inflation factors, and SAP-ready boundary tables.

Open Group Sequential Design Calculator