Tutorials › Group Sequential Methods and Sample Size Re-Estimation › Nuisance-Driven Re-Estimation and Error Control

Type I Error Control in Re-Estimation · Tutorial 452 of 699

Nuisance-Driven Re-Estimation and Error Control

Learn to separate nuisance-driven planning from treatment-effect-driven adaptation and identify the conditional validity requirements that protect type I error.

Advanced 12 min read

What You'll Learn

  • Distinguish nuisance parameters from treatment-effect estimates in sample size re-estimation.
  • Calculate a second-stage sample size from an estimated variance or event rate.
  • State the conditional validity conditions needed for nuisance-based adaptation.
  • Explain why blinding alone does not guarantee type I error control.
  • Recognize when conditioning or an independent nuisance estimate supports valid combination testing.
  • Diagnose outcome-dependent sample size choices that can invalidate a nominal cutoff.

Why Nuisance-Driven Re-Estimation Needs Its Own Validity Check

A second-stage sample size may be changed because the observed variability is larger than expected, or because the event rate suggests that the original plan will yield less information than anticipated. These are nuisance-driven changes: they use quantities that affect precision or information, rather than the estimated treatment effect itself. That distinction is useful, but it does not by itself guarantee type I error control. The key question is whether, after the nuisance information is observed and the sample size is chosen, the test still has a valid null distribution.

The relevant standard is conditional validity. If the adaptation information is denoted by \(A\), then for each possible value \(a\), the testing procedure should control its conditional probability of rejection under \(H_0\). Averaging valid conditional error rates over possible values of \(A\) then controls the unconditional type I error. If the nuisance estimate is related to the treatment contrast, however, conditioning on it can change the distribution of the test statistic. Blinding may conceal the treatment labels from the person making the choice, but concealment alone does not establish independence or conditional validity.

Definition: A nuisance-driven sample size re-estimation chooses the second-stage sample size using an estimate of a parameter that affects the variability, event frequency, or information of the test, rather than using the estimated treatment effect as the target for that choice. Conditional validity—calibration under the null conditional on the information used for adaptation—is a sufficient guarantee of unconditional type I error control, but it is not necessary. Unconditional type I error control requires only that the rejection probability, averaged over the adaptation information, is at most the nominal level.

Sample Size Calculations Use Nuisance Quantities

For planning, nuisance estimates commonly enter a standard-error or information calculation. For example, for two independent groups with normally distributed outcomes, equal group sizes, common known standard deviation \(\sigma\), and a one-sided test at level \(\alpha\), the approximate per-group sample size for power \(1-\beta\) at a target difference \(\delta>0\) is

$$ n \approx \frac{2\sigma^2(z_{1-\alpha}+z_{1-\beta})^2}{\delta^2}. $$

Here \(z_q\) is the \(q\)-quantile of the standard normal distribution. The formula makes clear why an estimated \(\sigma^2\) may change the planned sample size: larger variance requires more observations to achieve the same target precision and power. But this planning calculation is not itself a type I error argument. That argument concerns how the final rejection rule behaves under \(H_0\), after the sample size has been selected.

Worked Example: Re-Estimating Sample Size From an External Variance Estimate

Suppose a hypothetical two-group study plans to detect a mean difference of \(4\) units, using a one-sided significance level of \(0.025\) and power \(0.90\). An independent external variance study suggests a common standard deviation of \(5\) units. Under the normal planning approximation, calculate the required number per group.

The assumptions for this calculation are independent observations, equal group sizes, a common standard deviation, and a normal-approximation planning formula. Use \(z_{0.975}=1.960\) and \(z_{0.90}=1.2816\):

$$ n \approx \frac{2(5^2)(1.960+1.2816)^2}{4^2} = \frac{50(3.2416)^2}{16} = \frac{50(10.50797)}{16} \approx 32.84. $$

Rounding up gives \(33\) participants per group. This is a planning estimate, not a claim that a test using \(33\) per group automatically has exact type I error control. In this example, the variance estimate comes from an independent external source, so it does not select the sample size according to the study's observed treatment contrast. The final test must still use a properly calibrated rejection rule.

The Conditional Validity Argument

A convenient way to see the protection is through a one-sided inverse-normal combination test. Let \(P_1\) and \(P_2\) be stagewise p-values, and define \(Z_i=\Phi^{-1}(1-P_i)\), where \(\Phi\) is the standard normal cumulative distribution function. Suppose the adaptation information \(A\) determines the second-stage sample size and therefore determines weights \(w_1(A)\) and \(w_2(A)\). A combination statistic is

$$ C=w_1(A)Z_1+w_2(A)Z_2, \qquad w_1(A)^2+w_2(A)^2=1. $$

The exact argument below requires that, under \(H_0\), conditional on \(A\), the stagewise p-values are independent and uniform on \((0,1)\), and that the weights are fixed once \(A\) is given. In that case, conditional on \(A\), \(Z_1\) and \(Z_2\) are independent standard normal variables. Thus \(C\), conditional on \(A\), is standard normal. For a level-\(\alpha\) upper-tail test, reject when \(C>z_{1-\alpha}\). The conditional rejection probability is then \(\alpha\) for every value of \(A\). Averaging over \(A\) gives the unconditional result:

$$ P_{H_0}(C>z_{1-\alpha}) = E_{H_0}\!\left[ P_{H_0}(C>z_{1-\alpha}\mid A) \right] = E_{H_0}[\alpha] = \alpha. $$

This is a specific sufficient condition, not a blanket guarantee for every nuisance-based rule. In particular, if the first-stage p-value is not uniform conditional on the information used to select the weights or sample size, the normal combination argument does not follow. For discrete tests, p-values may instead be super-uniform; the test then needs a calibration that accounts for that fact rather than assuming continuous uniformity without justification.

Worked Example: A Nuisance-Based Choice of Combination Weights

Consider a hypothetical two-stage test with a one-sided level of \(0.025\). An independent nuisance assessment determines the second-stage size and, according to a prespecified rule, fixes \(w_1=0.6\) and \(w_2=0.8\). Suppose the observed stagewise p-values are \(P_1=0.0800\) and \(P_2=0.0300\). For this illustration, assume the stagewise p-values are independent and conditionally uniform under \(H_0\) given the information used to choose the weights.

The weights satisfy the required normalization:

$$ 0.6^2+0.8^2=0.36+0.64=1. $$

Transform the p-values to normal scores. Using \(\Phi^{-1}(0.9200)\approx1.405\) and \(\Phi^{-1}(0.9700)\approx1.881\), the combination statistic is

$$ C=0.6(1.405)+0.8(1.881) =0.8430+1.5048 =2.3478. $$

The one-sided p-value is approximately \(1-\Phi(2.3478)=0.00945\). Since \(0.00945<0.025\), this hypothetical result rejects \(H_0\) at the specified level. The type I error justification is not that the weights happened to be based on a nuisance quantity; it is that the weights were fixed conditional on that information and the stagewise p-values satisfy the stated conditional null assumptions.

Event Rates and Blinded Information

Binary outcomes introduce a related issue. An event-rate estimate can affect the variance of a difference in proportions, and hence the planned information. A blinded pooled event rate may be attractive because it does not identify which treatment group has the higher observed rate. Nevertheless, pooled information is still computed from outcomes. It cannot be declared independent of the treatment statistic merely because treatment labels were hidden.

One way to make the reasoning explicit is to use an independent nuisance source, such as a prespecified external estimate or a separate blinded run-in sample, for the sample-size choice. Another is to use a test and conditioning argument that establishes conditional null validity given the information used in adaptation. For example, conditional tests for binary outcomes can condition on appropriate event totals; discreteness then often makes the p-values conservative rather than exactly uniform. The particular test and conditioning set must match the adaptation rule.

Worked Example: Event-Rate Planning and the Limits of Blinding

Suppose a hypothetical trial plans to compare two proportions, targets an absolute difference of \(0.10\), and uses a one-sided significance level of \(0.025\) with power \(0.90\). An independent pilot suggests a common event rate of \(q=0.20\) for planning. Using the common-rate normal approximation, the per-group sample size is

$$ n \approx \frac{2q(1-q)(z_{0.975}+z_{0.90})^2}{(0.10)^2} = \frac{2(0.20)(0.80)(3.2416)^2}{0.01} = \frac{0.32(10.50797)}{0.01} \approx 336.26. $$

Rounding up gives \(337\) participants per group under this approximation. It assumes independent binary outcomes and uses the same event rate in both groups for the variance calculation; it is a planning approximation, not an exact sample-size formula for every alternative.

Because the pilot is independent, this calculation does not use the trial's treatment contrast to select the trial size. By contrast, suppose the same calculation is repeated using the trial's pooled event rate at an interim look. Hiding the group labels prevents a direct treatment-effect estimate from being used, but it does not prove that the first-stage test statistic is conditionally valid given that pooled rate. Before using a nominal combination cutoff, the analyst must show that the test's null calibration holds given the pooled information, or use a procedure that preserves the appropriate conditional error.

Why Outcome-Dependent Choices Can Break the Argument

The danger is not limited to deliberately choosing a sample size from a treatment-effect estimate. It arises whenever adaptation information carries information about the statistic used for testing and the final procedure treats that information as if it were irrelevant. Suppose the first-stage contrast tends to be positive under \(H_0\) simply through random variation, and a rule responds to that contrast by selecting a second-stage size or combination weight. Conditional on that choice, the distribution of the first-stage evidence is no longer the unconditional null distribution used by a fixed-design cutoff.

A common error is to say, “The second-stage sample size depends only on variance, so type I error is unchanged.” The missing step is a dependence check. If the variance estimate is statistically independent of the first-stage test statistic under the null, and the stagewise testing construction remains valid conditional on that estimate, the argument can work. If the estimate shares data with the test statistic, this independence or conditional calibration must be established—not assumed. Conversely, an outcome-dependent choice is not automatically invalid if the final procedure explicitly preserves conditional error; it is invalid to use a fixed-design calibration without showing that the calibration still applies.

1
Specify the adaptation information.
Write down exactly which data or estimate determines the second-stage sample size, including whether it uses trial outcomes, external data, or both.
2
Check its relation to the test statistic.
Under \(H_0\), establish independence or derive the conditional null distribution given that information. Blinding alone is not such a derivation.
3
Check the final calibration.
Confirm that stagewise p-values and any data-dependent weights meet the conditions of the combination procedure conditional on the adaptation information.

Diagnostics, Alternatives, and Limitations

Practical diagnostics should focus on the joint behaviour of the adaptation statistic and the testing statistic under the null model. Simulation can be useful: generate data under representative null parameter values, repeat the full sample-size selection and testing algorithm, and estimate the rejection frequency. This is a check on the complete design, not a substitute for a proof when a formal guarantee is required. Look for unexpected association between the nuisance estimate and the stage-one signed statistic, changes in the conditional distribution of stagewise p-values across adaptation values, and differences between simulated and claimed conditional rejection rates.

The main alternatives are to use an independent nuisance estimate, to use a test that is valid conditional on the observed nuisance information, or to recalibrate the final test through conditional-error methods. The first can simplify the independence argument but may be less precise or less current. Conditioning can give rigorous control but may cost power, especially with discrete outcomes. Conditional-error recalibration can accommodate broader adaptations, but it requires computing and preserving the relevant conditional rejection probability. These approaches are not interchangeable without checking their assumptions.

Nuisance-based re-estimation also has limits beyond type I error. A variance estimate from a small interim sample can be unstable; a second-stage size based on it may be too small or unnecessarily large. A sample-size cap may prevent the planned power from being attained, even when type I error remains controlled. Misspecified variance or event-rate models can make the information calculation inaccurate. Therefore report the estimate used, the rule mapping it to sample size, any bounds or rounding, and the final testing calibration. Sample-size adaptivity can improve planning, but it cannot guarantee power under every nuisance value or repair a misspecified outcome model.

Common Pitfall

  • Wrong reasoning: “The sample size is based on a blinded pooled variance or event rate, not the estimated treatment effect, so the nominal type I error is automatically preserved.”
  • Correction: Show that the test remains calibrated conditional on the information used for adaptation, or use a procedure that explicitly preserves the relevant conditional error. Blinding may help prevent a direct effect-driven choice, but it does not itself prove conditional independence or conditional validity.
Key takeaway: A nuisance estimate may determine the second-stage sample size without inflating type I error, but only when the final test is valid under the null after accounting for that choice. The operative question is conditional validity, not whether the adaptation variable has been labelled “nuisance” or kept blinded.

Check Your Understanding

Answer each question, then compare with the solutions below.

  1. A second-stage size is selected from an independent external variance estimate. What additional condition about the final test is still needed?
  2. Why does blinding treatment labels not, on its own, prove that a pooled event-rate estimate is safe to use for re-estimation?
  3. In an inverse-normal combination test, weights are chosen from adaptation information \(A\). State the conditional assumptions that yield a standard normal combination statistic.
  4. What is the per-group sample size from the common-rate approximation if \(q=0.25\), the target difference is \(0.10\), and \(z_{1-\alpha}+z_{1-\beta}=3.0\)?
  5. A nuisance estimate is associated with the first-stage signed test statistic under \(H_0\). Is a fixed-design cutoff justified solely because the estimate measures variance? Explain.

Solutions. (1) The final test must remain calibrated under \(H_0\); for the conditional combination argument, its stagewise p-values must meet the required null assumptions and the rejection rule must account for the selected design. (2) The pooled estimate is still a function of outcomes, and hiding labels does not establish independence or conditional validity relative to the test statistic. (3) Conditional on \(A\), the stagewise p-values must be independent uniform null variables, and the weights must be fixed with \(w_1(A)^2+w_2(A)^2=1\). (4) \(n\approx 2(0.25)(0.75)(3.0)^2/(0.10)^2=3.375/0.01=337.5\) participants per group. (5) No. The association means the test statistic's conditional null distribution may differ from the fixed-design one; establish conditional calibration or use an error-controlled recalibration.

Clinical Biostats