Tutorials › Group Sequential Methods and Sample Size Re-Estimation › Sequential Bonferroni Tests For Multiple Endpoints

Multiple Endpoints Across Interim Looks · Tutorial 602 of 699

Sequential Bonferroni Tests For Multiple Endpoints

You will derive and apply a simple Bonferroni procedure that controls false rejections across endpoints and interim looks.

Advanced 10 min read

What You'll Learn

  • Define strong familywise error control across endpoints and interim looks.
  • Allocate a global significance level among endpoint-specific hypotheses.
  • Combine endpoint allocations with valid group sequential tests.
  • Verify control under arbitrary dependence among endpoint statistics.
  • Recognize when conservative allocations or invalid interim analyses undermine the procedure.
  • Compare fixed Bonferroni allocations with closure-based testing.

Why Allocate Error by Endpoint?

A trial may monitor several endpoint hypotheses at interim analyses. If each endpoint is tested at the full significance level, the chance of at least one false rejection can exceed the intended familywise level. Interim monitoring adds another opportunity for a false rejection unless the test for each endpoint already accounts for its repeated looks. A direct solution is to divide the total error budget among endpoints, then use a valid group sequential test at each endpoint’s assigned level.

The procedure is deliberately simple: allocate the global level in advance, do not transfer unused level between endpoints, and let each endpoint’s sequential test account for all of its planned looks. The endpoint tests can be dependent; the error-control argument does not require independence.

Definition: For \(m\) endpoint hypotheses \(H_1,\ldots,H_m\), the familywise error rate is the probability of rejecting at least one hypothesis that is true. Strong control at level \(\alpha\) means this probability is at most \(\alpha\) for every configuration of true and false endpoint hypotheses, over all endpoints and all interim looks.

Fixed Endpoint Allocations and the Control Result

Choose nonnegative endpoint levels \(\alpha_1,\ldots,\alpha_m\) before examining the interim results, with \(\sum_{j=1}^m\alpha_j\leq\alpha\). For endpoint \(j\), use a group sequential test whose probability of ever rejecting \(H_j\), when \(H_j\) is true, is at most \(\alpha_j\). “Ever” includes rejection at any planned interim or final look. Each endpoint test must be calibrated for its own repeated looks; applying an ordinary single-look test repeatedly at level \(\alpha_j\) does not meet this condition in general.

Let \(R_j\) be the event that the procedure rejects \(H_j\) at any look, and let \(I_0\) be the set of indices of the true null hypotheses. A familywise error occurs when at least one of the events \(R_j\), \(j\in I_0\), occurs. The union bound gives the central justification:

$$ P\left(\bigcup_{j\in I_0}R_j\right) \leq \sum_{j\in I_0}P(R_j) \leq \sum_{j\in I_0}\alpha_j \leq \sum_{j=1}^{m}\alpha_j \leq \alpha. $$

The second inequality uses the validity of each endpoint’s group sequential test under its own null. The bound does not depend on the dependence among endpoints, nor does it require all null hypotheses to be true. It therefore establishes strong, rather than only global-null or weak, familywise error control.

A common allocation is equal Bonferroni allocation, \(\alpha_j=\alpha/m\). More generally, fixed nonnegative weights \(w_j\) with \(\sum_jw_j\leq1\) give \(\alpha_j=\alpha w_j\). These weights set priorities before the data are examined; this procedure does not recycle an endpoint’s unused level after that endpoint is rejected or reaches its final look.

Worked Example: Equal Allocation Across Three Endpoints

Consider a hypothetical trial with \(m=3\) endpoint hypotheses and a familywise target of \(\alpha=0.05\). Assign each endpoint the exact level \(\alpha_j=0.05/3\). Numerically, this is approximately \(0.01667\), but the allocation is the exact fraction, not the rounded display value.

Check the budget. The total assigned level is \(3(0.05/3)=0.05\), exactly. If the three group sequential tests each have probability at most \(0.05/3\) of ever rejecting a true endpoint null, then for any true-null set \(I_0\), \(P(\text{at least one true null rejected})\leq |I_0|(0.05/3)\leq0.05\).

For instance, if only endpoints 1 and 3 have true null hypotheses, the bound is \(2(0.05/3)=0.03333\), which is below the familywise target. If all three nulls are true, the bound is \(0.05\). These are upper bounds, not claims that the actual error probabilities equal those values.

How Sequential Testing Fits the Allocation

The outer allocation handles multiplicity across endpoints. Inside each endpoint, the group sequential test handles repeated looks. The local test must have its type I error controlled over the entire sequence of looks under the endpoint null. It may use boundaries calibrated for a specified spending function or another valid group sequential design; the particular boundary construction is separate from the endpoint allocation.

One transparent, generally conservative construction is to combine Bonferroni adjustment across looks with the endpoint allocation. Suppose endpoint \(j\) has \(K_j\) planned looks and a valid single-look p-value is available at each look. Test at look \(k\) only if its p-value is at most \(\alpha_j/K_j\). Under a true \(H_j\), each fixed-look p-value satisfies \(P(p_{jk}\leq t)\leq t\). Therefore:

$$ P(\text{reject }H_j\text{ at any look}) \leq \sum_{k=1}^{K_j}P(p_{jk}\leq\alpha_j/K_j) \leq K_j(\alpha_j/K_j)=\alpha_j. $$

This remains valid if the statistics at different looks are correlated, because the union bound does not assume independence. A group sequential design calibrated directly to \(\alpha_j\) can account for the joint distribution across looks and may be less conservative. Either way, the local requirement is the same: control the probability of any false rejection for endpoint \(j\) over its full monitoring sequence.

Worked Example: Three Endpoints Monitored at Two Looks

Suppose a hypothetical trial has three endpoints, each observed at two planned looks, and uses the exact equal allocation \(\alpha_j=0.05/3\). For a simple conservative procedure, apply the single-look threshold \(\alpha_j/2=(0.05/3)/2=0.05/6=0.0083333\) at each look. Assume the p-values at each planned look are valid under their respective null hypotheses.

The hypothetical p-values are: endpoint 1, \(0.012\) at look 1 and \(0.006\) at look 2; endpoint 2, \(0.004\) at look 1; endpoint 3, \(0.010\) at look 1 and \(0.009\) at look 2. Stop testing an endpoint after its first rejection.

Endpoint 1: \(0.012>0.0083333\), so there is no rejection at look 1. At look 2, \(0.006\leq0.0083333\), so reject \(H_1\).
Endpoint 2: \(0.004\leq0.0083333\), so reject \(H_2\) at look 1.
Endpoint 3: both \(0.010\) and \(0.009\) exceed \(0.0083333\), so do not reject \(H_3\) at either look.

For each endpoint the probability of a false rejection over its two looks is at most \(2(0.0083333\ldots)=0.0166667\ldots=0.05/3\). Across the three endpoints, the probability of any false endpoint rejection is at most \(3(0.05/3)=0.05\). This calculation does not require independence across endpoints or looks. The stated rejection decisions describe the hypothetical data; they do not imply that the error bound is attained.

Unequal Allocations and Sequential P-Values

Equal allocation is not required. Investigators may assign more of the total level to endpoints considered more important, provided the levels are fixed in advance, nonnegative, and sum to no more than \(\alpha\). Allocation affects the local thresholds and hence the opportunity to reject each endpoint; it does not change the familywise guarantee if every local sequential test is valid at its assigned level.

A sequentially adjusted p-value for an endpoint summarizes the evidence along its monitoring path while accounting for the interim design. If \(q_j\) denotes such a p-value and is valid in the sense that \(P(q_j\leq t\mid H_j\text{ true})\leq t\), the endpoint can be rejected when \(q_j\leq\alpha_j\). It is not appropriate to substitute an unadjusted final-look p-value for \(q_j\) after repeatedly examining interim results.

Worked Example: Fixed Unequal Allocation

Suppose a hypothetical trial has three endpoint hypotheses, a global target \(\alpha=0.05\), and prespecified weights \(w_1=0.50\), \(w_2=0.30\), and \(w_3=0.20\). The weights sum to \(1\), so the endpoint levels are \(\alpha_1=0.05(0.50)=0.025\), \(\alpha_2=0.05(0.30)=0.015\), and \(\alpha_3=0.05(0.20)=0.010\). Their sum is \(0.025+0.015+0.010=0.050\).

Assume valid group sequential tests at these levels produce sequentially adjusted p-values \(q_1=0.021\), \(q_2=0.017\), and \(q_3=0.008\). Compare each value with its own allocation: \(0.021\leq0.025\), so reject \(H_1\); \(0.017>0.015\), so do not reject \(H_2\); and \(0.008\leq0.010\), so reject \(H_3\).

If \(H_1\) and \(H_3\) are in fact true, the probability of falsely rejecting at least one of them is bounded by \(0.025+0.010=0.035\). If all three are true, the bound is \(0.050\). Either way, strong familywise control holds. The weights are not estimated from these p-values: choosing or changing them after seeing the results would invalidate the stated allocation argument.

Diagnostics, Alternatives, and Limitations

Before applying the procedure, verify both levels of calibration. Across endpoints, check that the exact assigned levels sum to no more than \(\alpha\). Within each endpoint, confirm that its test controls rejection over all planned looks at its assigned level, under the endpoint null and the actual monitoring plan. The endpoint test’s validity must reflect relevant design features, such as the analysis schedule and any adaptation that affects the test statistic or its distribution.

Warning signs include comparing every interim p-value with the full endpoint level, using a final-look p-value as though no interim analyses occurred, or rounding levels upward so their sum exceeds the global target. Another failure is using a local sequential procedure whose calibration assumes a different number or timing of looks. Bonferroni allocation cannot repair a locally invalid sequential test.

The method is easy to audit and remains valid under arbitrary dependence among endpoints, but it can be conservative. The union bound ignores the extent to which rejection events overlap, and a fixed allocation cannot make use of evidence from one endpoint to alter the level for another. The closure principle provides a more general framework for multiplicity control by testing intersection hypotheses; this simple allocation is attractive when a transparent, non-recycling rule is sufficient. A single global test may answer a different scientific question and, by itself, does not necessarily provide controlled endpoint-specific claims.

Common pitfall: “Each of three endpoints gets level 0.01667, so the total is 0.05.” If \(0.01667\) is treated as the literal allocation, then \(3(0.01667)=0.05001\), which exceeds \(0.05\). Use the exact allocation \(0.05/3\) and display \(0.01667\) only as a rounded approximation, or choose levels whose exact sum is at most \(0.05\).

The procedure’s scope is also important. It controls the probability of one or more false endpoint rejections, not the probability that every endpoint decision is correct, and it does not guarantee adequate power for any particular endpoint. If endpoint hypotheses, weights, or the monitoring plan are changed in response to accumulating data, the original simple proof applies only if the revised procedure is itself shown to preserve the local error bounds and the total allocation.

Check Your Understanding

Answer each question using the definitions and bounds above. Solutions follow the questions.

  1. Four endpoint hypotheses are assigned \(\alpha_j=0.0125\) each. What familywise level is guaranteed if each local group sequential test is valid at its assigned level?
  2. Two endpoints are monitored at three looks each. With endpoint levels \(\alpha_1=0.03\) and \(\alpha_2=0.02\), what single-look threshold results from equal Bonferroni allocation over looks within each endpoint?
  3. Why does the proof still work when endpoint statistics are correlated?
  4. With a global level of \(0.05\), three literal allocations of \(0.01667\) each are proposed. Do they meet the budget? Give a corrected allocation.
  5. One endpoint’s group sequential test repeatedly compares ordinary single-look p-values with its full allocated level. Does the endpoint allocation alone guarantee control?

Solutions. (1) The sum is \(4(0.0125)=0.05\), so the procedure guarantees strong familywise control at level \(0.05\). (2) For endpoint 1, \(0.03/3=0.01\) per look; for endpoint 2, \(0.02/3=0.006667\) per look. Their local error bounds are \(0.03\) and \(0.02\), respectively, and sum to \(0.05\). (3) The union bound \(P(\cup_jR_j)\leq\sum_jP(R_j)\) does not require independent events. (4) No: \(3(0.01667)=0.05001>0.05\). Use the exact levels \(0.05/3\) each, displaying each as approximately \(0.01667\), or choose rounded-down levels such as \(0.01666\) each. (5) No. Repeated single-look testing at the full local level may exceed that level over the sequence. Use a valid group sequential test at the local level or a procedure that adjusts for the looks.

Clinical Biostats