Tutorials › Group Sequential Methods and Sample Size Re-Estimation › The O’Brien-Fleming Boundary Shape

Classical Boundaries · Tutorial 96 of 699

The O'Brien-Fleming Boundary Shape

You will learn to calculate and interpret O’Brien-Fleming critical values across planned information fractions, and distinguish the boundary’s shape from its error-rate calibration.

Intermediate 10 min read

What You'll Learn

  • Express the O’Brien-Fleming boundary as a critical value for a standardized statistic at each information fraction.
  • Explain mathematically why the critical value is largest at the earliest look and declines toward the final analysis.
  • Calculate illustrative critical values and assess whether hypothetical statistics cross them.
  • Distinguish boundary shape from calibration of the familywise Type I error rate.
  • Recognize common scale, interpretation, and implementation errors.

Why the O’Brien-Fleming Boundary Has a Distinctive Shape

A group sequential design permits planned interim analyses before the target information is reached. At each look, an efficacy decision is based on a standardized statistic and a critical value. The O’Brien-Fleming boundary is distinctive because it demands unusually strong standardized evidence early, then relaxes as the trial approaches its final analysis. This shape allows interim stopping for efficacy while making early crossing uncommon under the null hypothesis.

The boundary is expressed on the standardized-statistic scale: it gives a critical value for \(Z\), not directly a treatment-effect difference in its original units. Its shape is determined by the information fraction at each look. The constant that sets the boundary’s overall height must be calibrated for the design; the shape formula alone does not establish a particular Type I error rate.

Definition: Let \(I_k\) be the planned statistical information at look \(k\), and let \(I_K\) be the planned information at the final look \(K\). The information fraction is \(t_k=I_k/I_K\), with \(0<t_1<\cdots<t_K=1\). For a one-sided upper efficacy boundary, the canonical O’Brien-Fleming-shaped critical value is \(b_k=c/\sqrt{t_k}\), where \(c>0\) is a design-specific constant. The trial crosses the boundary at look \(k\) if \(Z_k\ge b_k\). For a symmetric two-sided design, crossing is defined by \(|Z_k|\ge b_k\).

Assumptions and What the Formula Specifies

The usual group sequential framework uses a standardized statistic \(Z_k\) at each planned look. Under the null hypothesis, the canonical joint model treats the vector of statistics as multivariate normal, with each statistic having mean zero and variance one, and with correlation \(\operatorname{Corr}(Z_i,Z_j)=\sqrt{t_i/t_j}\) for \(i<j\). This model reflects the fact that analyses use accumulating data and are therefore correlated. In applications, the normal model may be an approximation, so its adequacy depends on the outcome, sample size, and chosen statistic.

The information fractions must describe the actual planned information at the looks, and the sequence must be increasing, ending at \(t_K=1\). The formula applies to a specified efficacy direction: an upper boundary tests for evidence in the positive direction; the symmetric two-sided version also allows evidence in the negative direction. A design with a lower efficacy direction can reverse the sign convention. The statistic and direction must be defined consistently.

The coefficient \(c\) is not universally fixed. To control the probability of at least one false efficacy crossing under \(H_0\) at the target familywise significance level \(\alpha\), the design calibrates \(c\) using the joint null distribution of the planned statistics and the complete set of information fractions. Therefore, two designs can share the O’Brien-Fleming shape but have different numerical critical values. A value of \(c\) used for illustrating the shape should not be described as a calibrated design unless that calibration has actually been performed.

Worked Example: Calculating the Shape Across Three Looks

Suppose a hypothetical three-look design has information fractions \(t_1=0.25\), \(t_2=0.60\), and \(t_3=1.00\). For illustration, suppose the design uses the stipulated coefficient \(c=2.02\). The calculations below show the shape; they do not, by themselves, verify the design’s familywise Type I error rate.

At the first look,

$$ b_1=\frac{2.02}{\sqrt{0.25}}=\frac{2.02}{0.50}=4.040. $$

At the second look,

$$ b_2=\frac{2.02}{\sqrt{0.60}}=\frac{2.02}{0.7745967}\approx 2.608. $$

At the final look,

$$ b_3=\frac{2.02}{\sqrt{1.00}}=2.020. $$

Thus the upper critical values are 4.040, 2.608, and 2.020 standardized-statistic units. They decline as information accumulates. The final value equals \(c\), because the final information fraction is one. These are cutoffs for \(Z\), not treatment effects measured in, for example, milligrams or points.

Why the Critical Values Decline

The shape follows directly from the dependence of the cutoff on \(t\). For \(c>0\), define \(b(t)=c t^{-1/2}\) for \(0<t\le 1\). Differentiation gives

$$ \frac{d b(t)}{dt}=-\frac{c}{2t^{3/2}}<0. $$

Consequently, \(b(t)\) is strictly decreasing as the information fraction increases. Also, \(b(1)=c\), while \(b(t)>c\) whenever \(t<1\). As \(t\) approaches one, the cutoff approaches \(c\). The formula therefore gives both defining features of the shape: a higher threshold at every interim information fraction and a decline toward the final critical value.

This mathematical description concerns the standardized-statistic boundary. If a statistic is formed by dividing an estimated treatment effect by its standard error, a critical value for \(Z\) can be translated into a corresponding threshold for the estimated effect only after specifying that standard error. The standardized cutoff alone does not state a fixed effect-size requirement in the outcome’s measurement units.

Worked Example: Deciding Whether an Interim Statistic Crosses

Use the hypothetical fractions and coefficient from the preceding example. Suppose the observed upper-direction statistics are \(Z_1=3.20\) at \(t_1=0.25\) and \(Z_2=2.70\) at \(t_2=0.60\). Assume the first analysis has not crossed and that these statistics and information fractions are calculated as planned.

At the first look, \(3.20<4.040\), so the statistic does not cross the O’Brien-Fleming-shaped boundary. At the second look, compare the statistic with the second critical value:

$$ 2.70>2.608. $$

The statistic crosses the upper boundary at the second look. Under the design’s prespecified decision rule, this is an efficacy-boundary crossing; whether to stop or take another action is governed by that rule. The comparison does not establish that the illustrative coefficient controls a particular \(\alpha\), nor does it mean there is a 2.70-unit treatment effect. The observed number is a standardized statistic.

Reading the Boundary Without Confusing Its Scales

A critical value answers a specific question: how large must the standardized statistic be at this look to cross the boundary? It is not a treatment effect, a standard error, or a nominal \(p\)-value. At the same information fraction, a larger positive \(Z_k\) gives stronger evidence in the upper direction, but the decision rule compares it with the boundary assigned to that look.

For instance, with \(c=2.02\) and \(t=0.40\), the cutoff is

$$ b(0.40)=\frac{2.02}{\sqrt{0.40}} =\frac{2.02}{0.6324555} \approx 3.194. $$

An observed upper-direction statistic of 2.30 does not cross this interim boundary because \(2.30<3.194\). Comparing it instead with the final cutoff, 2.020, would answer the wrong question: it would ignore the planned interim boundary. A group sequential decision must use the critical value for the look at which the statistic was observed.

Comparison, Diagnostics, and Limitations

A useful contrast is Pocock’s constant critical-value boundary. In the standard Pocock shape, the standardized cutoff is the same at each planned look; in the O’Brien-Fleming shape, it starts higher and declines. Both require calibration under the joint distribution of the repeated statistics to meet a specified familywise error rate. The practical distinction is where the design places its stringency: O’Brien-Fleming demands more extreme standardized evidence at early looks, whereas a flat boundary does not increase its cutoff at those early looks in the same way.

A fixed-sample analysis uses a critical value for one planned analysis, rather than a sequence of look-specific group sequential cutoffs. Reusing that fixed-sample cutoff at every interim look is not equivalent to applying a calibrated group sequential design: repeated opportunities to cross can change the probability of a false rejection. The O’Brien-Fleming shape is therefore part of a planned sequential decision rule, not a justification for unplanned repeated testing.

Several checks help identify implementation problems:

  • Confirm that each \(t_k\) is calculated relative to the planned final information and that the final fraction is one.
  • Verify that the computed cutoffs decline as \(t_k\) increases, and that the final cutoff equals \(c\).
  • Check that the observed statistic, boundary direction, and one-sided or two-sided rule agree.
  • Use the cutoff for the current look rather than reusing an interim value or the final value at every analysis.
  • Ensure that \(c\) was calibrated for the actual number and timing of looks, target \(\alpha\), and specified design—not merely borrowed from an illustrative calculation.

The shape does not itself choose the number or timing of analyses, establish that the assumed joint model is adequate, or guarantee good operating characteristics under alternatives. Nor does it supply a treatment-effect threshold without the relevant standard error. The design must specify and validate those features separately. These limitations are especially important when actual information fractions depart from the planned ones: inserting an incorrect fraction into the formula can produce a numerically precise but inappropriate cutoff.

Key takeaway: On the standardized-statistic scale, the O’Brien-Fleming-shaped efficacy cutoff is \(b_k=c/\sqrt{t_k}\). Because this function decreases with information, early looks require larger \(Z\)-statistics, and the final cutoff is \(c\). The shape and the calibration of familywise error are related but distinct parts of a design.

Check Your Understanding

Answer each question using the boundary definition and the distinction between shape and calibration.

  1. For an O’Brien-Fleming-shaped boundary with \(c=2.10\), calculate the critical value at \(t=0.25\).
    Solution: \(b=2.10/\sqrt{0.25}=2.10/0.50=4.200\).
  2. Why must the cutoff at \(t=0.40\) exceed the final cutoff, when \(c>0\)?
    Solution: Since \(\sqrt{0.40}<1\), \(c/\sqrt{0.40}>c\). At the final look \(t=1\), the cutoff is \(c/\sqrt{1}=c\).
  3. With \(c=2.00\), \(t=0.64\), and an upper-direction statistic \(Z=2.40\), does the statistic cross the boundary?
    Solution: The cutoff is \(2.00/\sqrt{0.64}=2.00/0.80=2.50\). Since \(2.40<2.50\), it does not cross.
  4. Does specifying \(b_k=c/\sqrt{t_k}\) alone prove that a design controls familywise Type I error at a stated \(\alpha\)?
    Solution: No. The coefficient \(c\) must be calibrated for the design’s joint null distribution, information fractions, and target error rate.
  5. What is the main shape difference between O’Brien-Fleming and Pocock boundaries on the standardized-statistic scale?
    Solution: O’Brien-Fleming cutoffs decline from larger early values toward the final value; a standard Pocock boundary uses a constant critical value across looks.
Clinical Biostats