Why Information Scaling Matters
An O’Brien-Fleming boundary is high at early analyses and approaches its terminal reference threshold as information accumulates. The shape has a direct explanation: the score is standardized by the square root of information, so a fixed threshold on the accumulated score becomes a larger threshold on the standardized scale when information is small. This tutorial derives that relationship and clarifies exactly what it does—and does not—say about a calibrated group sequential design.
The argument uses the information-standardized score and the information fraction. The earlier material on the O’Brien-Fleming boundary shape establishes the qualitative pattern; here the aim is to explain its approximate scaling algebraically. The scaling is a useful description of the boundary’s shape, not a substitute for calibrating the planned analyses to control the design’s overall Type I error.
For clarity, take positive values of \(\theta\) to favor the treatment and consider an upper efficacy boundary. A two-sided design applies the corresponding argument to \(|Z_k|\); a lower one-sided boundary follows by reversing the sign. Write \(c_*\) for the selected terminal reference critical value. For example, the single-analysis two-sided normal reference at significance level \(\alpha=0.05\) is approximately \(1.96\). In a group sequential design, however, the actual terminal critical value may be chosen as part of a joint calibration and need not equal that single-analysis reference.
Deriving the Inverse-Square-Root Shape
The O’Brien-Fleming scaling approximation sets the standardized threshold at look \(k\) to \(c_*/\sqrt{t_k}\). Substituting \(t_k=I_k/I_K\) and using \(Z_k=U_k/\sqrt{I_k}\) shows what this approximation means on the unstandardized score scale.
The last expression is the key justification: under this shape approximation, the efficacy threshold on the cumulative score scale is constant across looks. Dividing that fixed score threshold by \(\sqrt{I_k}\) to obtain a standardized threshold produces the factor \(1/\sqrt{I_k}\), or equivalently \(1/\sqrt{t_k}\) after expressing information relative to its final value. At the final analysis, \(t_K=1\), so the approximation gives \(c(1)=c_*\).
This derivation is algebraically exact once the approximate boundary \(c_*/\sqrt{t_k}\) has been posited. It does not prove that these values are the exact critical values that give a chosen familywise Type I error for every set of planned looks. The standardized scores across looks are correlated, and the probability of crossing at any look depends on their joint distribution. A design must calibrate its actual boundaries accordingly.
Worked Calculations
Worked Example: Calculating Thresholds From Information Fractions
Suppose a hypothetical two-sided design uses \(c_*=1.96\) as its terminal reference value. Calculate the approximate standardized thresholds at information fractions \(t=0.20\), \(0.50\), and \(1.00\).
Use \(c(t)=c_*/\sqrt{t}\). At \(t=0.20\), \(\sqrt{0.20}\approx0.4472\), so \(c(0.20)=1.96/0.4472\approx4.383\). At \(t=0.50\), \(\sqrt{0.50}\approx0.7071\), so \(c(0.50)=1.96/0.7071\approx2.772\). At the terminal look, \(c(1.00)=1.96/1=1.960\).
Thus the approximate two-sided crossing criteria are \(|Z|\geq4.383\), \(|Z|\geq2.772\), and \(|Z|\geq1.960\), respectively. These are shape calculations using the stated reference value. They are not, by themselves, proof that a multi-look design using these cutoffs has exactly a 5% overall Type I error.
Worked Example: Checking the Score-Scale Meaning
Consider a hypothetical study with planned final information \(I_K=400\) and an interim analysis at \(I_k=100\), in the study’s information units. Use \(c_*=1.96\). Find the approximate interim threshold on both the standardized-score and cumulative-score scales.
The information fraction is \(t_k=100/400=0.25\). Therefore, the approximate standardized threshold is \(c(0.25)=1.96/\sqrt{0.25}=1.96/0.5=3.920\). Since \(Z_k=U_k/\sqrt{100}=U_k/10\), crossing the upper boundary requires \(U_k/10\geq3.920\), or \(U_k\geq39.20\).
The fixed score threshold predicted by the derivation is \(c_*\sqrt{I_K}=1.96\sqrt{400}=1.96(20)=39.20\), which agrees. At the final information, the corresponding score threshold is \(1.96\sqrt{400}=39.20\) again, while the standardized threshold is \(39.20/20=1.960\). The interim threshold is higher only on the standardized scale because the same score is divided by a smaller square root of information.
Worked Example: Converting a Boundary to an Effect-Estimate Threshold
Suppose a hypothetical analysis has \(I_K=625\), interim information \(I_k=225\), and terminal reference value \(c_*=1.96\). Under the canonical score model, the effect estimate is \(\hat{\theta}_k=U_k/I_k\). Calculate the approximate interim boundary and the corresponding effect-estimate threshold.
First, \(t_k=225/625=0.36\), so \(\sqrt{t_k}=0.60\). The standardized boundary is \(c(0.36)=1.96/0.60=3.267\). The interim standard error on the information scale is \(1/\sqrt{225}=1/15\), and \(Z_k=\hat{\theta}_k\sqrt{I_k}\). Thus the upper boundary requires \(\hat{\theta}_k\geq3.267/15\approx0.2178\), in the units of \(\theta\).
At the final analysis, the corresponding estimate threshold is \(1.96/\sqrt{625}=1.96/25=0.07840\) in the same units. The interim estimate threshold is larger: \(0.2178/0.07840\approx2.778\), which equals \(1/t_k=1/0.36\) up to rounding. This is consistent with a fixed score threshold: an estimate based on less information must be larger to generate the same cumulative score.
Assumptions and What the Approximation Leaves Out
The calculation rests on a particular statistical structure. The \(U_k\) values must be cumulative scores for the same treatment-effect parameter, and \(I_k\) must represent their corresponding cumulative information. Under the canonical model, the score has mean approximately \(\theta I_k\) and variance approximately \(I_k\), so \(Z_k=U_k/\sqrt{I_k}\) has approximately unit variance under \(H_0\). The planned information fractions must be meaningful for the analysis sequence; in particular, the simple expression treats \(I_k\) and \(I_K\) as specified information quantities rather than arbitrary sample-size percentages.
For a standard group sequential calibration, the joint behavior across looks also matters. In the canonical model, earlier and later scores are dependent; under the null, their correlation is determined by the ratio of the earlier to later information. The probability of crossing at least one boundary is therefore not obtained by treating the analyses as independent or by assigning the full \(\alpha\) separately at every look. The inverse-square-root rule explains the O’Brien-Fleming shape, while a joint calibration determines the actual cutoffs needed for a specified overall error rate.
Nor does the algebra make an approximate boundary robust to a misspecified information model. If the actual variance structure differs from the one used to calculate information, then \(Z_k\) may not be standardized as assumed. Likewise, if the information fractions are inaccurately planned or estimated, substituting them into \(c_*/\sqrt{t_k}\) can yield a boundary that does not match the intended design. The formula is a guide to shape under the stated model, not a repair for invalid scores or information estimates.
Comparison, Diagnostics, and Pitfalls
A flat, Pocock-style boundary keeps the standardized critical value roughly constant across looks. O’Brien-Fleming scaling instead demands a much larger standardized score early and moves toward the terminal reference later. Both shapes require a design-level calibration when used for group sequential testing; neither comparison changes the need to account for repeated looks. A nominal single-look cutoff applied independently at every analysis is not generally a valid substitute for that calibration.
Useful checks follow directly from the scaling argument:
- At \(t=1\), the formula must return \(c_*\).
- For \(0<t<1\), \(c_*/\sqrt{t}\) must exceed \(c_*\) when \(c_*>0\).
- Halving the information fraction multiplies the approximate threshold by \(\sqrt{2}\), not by \(2\).
- Multiplying the standardized boundary by \(\sqrt{I_k}\) should recover the same score threshold \(c_*\sqrt{I_K}\) at each look.
A common mistake is to say that the required standardized threshold is inversely proportional to the information fraction itself. That would give \(c_*/t\), making the threshold too high. The score standardization uses \(\sqrt{I_k}\), so the correct relationship is \(c_*/\sqrt{t_k}\). A second mistake is to interpret this formula as an exact boundary that automatically spends the desired Type I error. The correction is to separate the shape approximation from joint calibration: use the scaling to understand the geometry, then use a valid group sequential calibration for the specified design and error rate.
The approximation is most useful for understanding why early O’Brien-Fleming thresholds are stringent and for checking whether a proposed boundary has the expected scale. It does not provide exact calibrated values for every finite set of looks, replace the assumed score model, or determine how a design should be adapted when information changes. Those questions require methods beyond this scaling derivation.
Check Your Understanding
Answer each question using the scaling argument. Brief solutions follow the questions.
- With \(c_*=2.00\) and \(t=0.25\), what is the approximate standardized threshold?
- If the information fraction increases from \(0.25\) to \(1.00\), by what factor does the approximate threshold change?
- For \(I_K=900\) and \(c_*=2.00\), what cumulative-score threshold does the derivation imply at every look?
- Why does an early analysis require a larger effect estimate than the final analysis to reach the approximate upper boundary?
- Does using \(c_*/\sqrt{t_k}\) alone guarantee a specified overall Type I error across multiple looks? Explain.
Solutions. (1) \(c(0.25)=2.00/\sqrt{0.25}=2.00/0.50=4.00\). (2) The threshold at \(t=1\) is \(2.00\), one-half the threshold at \(t=0.25\), so it decreases by a factor of \(2\). (3) The implied score threshold is \(c_*\sqrt{I_K}=2.00\sqrt{900}=2.00(30)=60.0\). (4) The effect estimate is \(U_k/I_k\); with less information, a larger estimate is needed to produce the same score threshold. (5) No. Overall Type I error depends on the joint behavior of the correlated statistics and requires joint calibration of the design’s boundaries.