Why Adaptation Changes the Error-Control Question
A sample-size re-estimation may change only the planned number of subjects, yet it can change the sampling distribution of the final test statistic. The key question is therefore not simply whether the adaptation used a blinded variance estimate. It is whether the complete procedure—including the rule for choosing the sample size and the final analysis—controls the chance of declaring equivalence when equivalence is false.
For average bioequivalence, let \(\theta\) be the population mean difference between test and reference on the log scale, and let \(L\) and \(U\) be the lower and upper equivalence limits, with \(L<U\). The alternative is that the mean difference lies inside the limits. In the usual two one-sided tests procedure (TOST), equivalence is declared only if both component null hypotheses are rejected.
The component hypotheses are \(H_{01}:\theta\leq L\) versus \(H_{a1}:\theta>L\), and \(H_{02}:\theta\geq U\) versus \(H_{a2}:\theta<U\). With a correctly calibrated fixed-sample test, each one-sided test is conducted at level \(\alpha\). Since equivalence requires rejecting both, the probability of declaring equivalence under either part of the null is no greater than the rejection probability of the corresponding component test, and thus no greater than \(\alpha\).
What Must Hold After Sample-Size Re-Estimation
Let \(A\) denote the outcome of the adaptation rule: for example, the selected sample size or a category of interim variance estimates. After adaptation, the final rejection rule may have a different conditional distribution for each value of \(A\). A direct sufficient condition for type I error control is that, for every possible adaptation outcome, the conditional probability of declaring equivalence is at most \(\alpha\) throughout each component of the null.
If the conditional probability inside the integral is at most \(\alpha\) for all \(a\) and every \(\theta\) in the null region, its average is also at most \(\alpha\). A less restrictive sufficient route is to establish conditional size at the relevant boundary and also establish monotonicity: for every adaptation outcome, the lower component test’s conditional rejection probability must be nondecreasing as \(\theta\) moves from the lower null region toward its alternative; for the upper component, it must be nonincreasing as \(\theta\) moves from the upper null region toward its alternative. Then the largest conditional rejection probability over each null component occurs at its boundary. Boundary control alone, without this monotonicity or control throughout the component, is not enough.
These are conditions on the conditional test, not merely on the fixed-sample test used after selecting a sample size. Conditioning on an adaptation outcome can alter the distribution of the statistic. A nominal critical value that is valid for a fixed sample size need not remain valid conditional on how that sample size was chosen.
Why Blinding Is Not a Complete Guarantee
A blinded variance-based rule avoids using the observed treatment difference as an explicit input to the sample-size decision. That can reduce the risk of directly selecting a larger study when interim results look favorable. It does not, by itself, prove type I error control. The relevant questions include whether the variance statistic is independent of the treatment-effect statistic under the model, whether the final statistic has its assumed distribution after adaptation, and whether the analysis accounts for any dependence created by using interim data both to adapt and to test.
For example, an independently conducted variance-only pilot can select the size of a subsequent confirmatory study. If the pilot outcome is independent of all confirmatory treatment-effect data, then conditional on the selected size the confirmatory data still follow the planned fixed-sample model. A different situation arises if the interim subjects contribute to the final estimate and their data also determine the sample size. In that case, validity needs to be established for the combined adaptive procedure rather than inferred from blinding.
Worked Example: Independent Variance Pilot and Conditional Control
Hypothetical setup. An independent, blinded pilot selects either \(n=24\) or \(n=40\) confirmatory subjects. The confirmatory estimate \(\bar{x}\) is normally distributed with known standard deviation \(\sigma=0.12\) log units, independently of the pilot. To examine the lower boundary, suppose the true value is \(\theta=L\), and use the one-sided z-test that rejects \(H_{01}\) when \(\bar{x}>L+1.645\sigma/\sqrt{n}\). The upper component is tested separately at the same one-sided level.
Conditional on either pilot outcome, \(n\) is fixed and the confirmatory data retain the specified normal distribution. At the lower boundary,
For any \(\theta<L\), the conditional rejection probability is
because \(L-\theta>0\), so the argument of \(\Phi\) is larger than 1.645. Thus the lower component test is conditionally controlled throughout its null component. Since declaring equivalence also requires rejecting the upper component, the probability of declaring equivalence when \(\theta\leq L\) cannot exceed this lower-test probability. The analogous argument applies when \(\theta\geq U\). The independent pilot’s choice between the two sizes therefore does not increase type I error under these stated assumptions.
This justification depends on independence and on the specified confirmatory analysis. It does not show that every blinded re-estimation using interim study data has the same property.
A Boundary Check Is Not Always Enough
The monotonicity requirement matters because an adaptation outcome can select data with unusual properties. The following hypothetical calculation shows the logical gap; it is not a model of a particular bioequivalence test.
Worked Example: Boundary Control Without Interior Control
Suppose an adaptive procedure has two equally likely adaptation outcomes at a particular null parameter value. At the lower boundary \(\theta=L\), its conditional probability of declaring equivalence is 0.05 for each outcome. The unconditional boundary probability is therefore
Now suppose that at a more distant point in the lower null component, the conditional probabilities are 0.08 and 0.02. Their average is
In this particular example, the unconditional probability is still 0.05, but one adaptation outcome has conditional rejection probability 0.08. If the outcomes were not equally likely at that parameter value—for instance, if their probabilities were 0.8 and 0.2—the unconditional probability would be \(0.8(0.08)+0.2(0.02)=0.068\), exceeding 0.05. Thus boundary control alone does not establish control over the entire null component. One needs conditional control throughout the component, or a justified monotonicity argument connecting the boundary to all other null values.
Methods for Protecting Type I Error
Several approaches can make an adaptive analysis defensible. The appropriate choice depends on what data inform the adaptation, what data enter the final estimate, and whether the planned test’s conditional distribution can be derived.
- Use an independent nuisance-parameter estimate. A separate pilot or genuinely independent variance source can select the confirmatory sample size while leaving the confirmatory test distribution unchanged, conditional on the selected size.
- Derive and use a conditionally calibrated test. The critical value or p-value is chosen to control the conditional error given the interim information. This requires a valid model for the conditional distribution, not simply reusing the fixed-sample critical value.
- Use a prespecified combination or sequential procedure. When interim data contribute to the final test, a method designed for the stages and decision rule can account for the information used at adaptation. Its error control must match the equivalence hypotheses and the actual adaptation rule.
- Keep the original fixed-sample design. This avoids adaptation-related distribution changes but may leave the study underpowered if the planning variance was inaccurate.
These approaches are not interchangeable shortcuts. A combination method or conditional-error method must be specified for the relevant test, adaptation timing, and data structure. An adjustment designed for a single one-sided test should not be assumed to control a two-component equivalence procedure without checking both components.
Worked Example: Correcting a Fixed-Sample Equivalence Interval
Consider hypothetical log-scale data with \(\bar{x}=0.010\) log units, \(s=0.120\) log units, and \(n=24\). Under the usual one-sample t model, suppose the fixed-sample TOST uses \(df=n-1=23\) and the 90% interval, with \(t_{0.95,23}\approx1.714\). The standard error is
The interval margin is
Therefore the 90% interval is
For prespecified log limits \(L\) and \(U\), this fixed-sample procedure declares equivalence when the entire interval lies within \((L,U)\). This calculation illustrates the fixed-sample rule only. If \(n\) had been selected using these same data, the interval’s nominal coverage and its associated test error would require separate justification; the arithmetic alone would not establish adaptive error control.
Diagnostics, Failure Modes, and Limitations
Before adopting a variance-based re-estimation, specify the adaptation rule, the maximum and minimum sample sizes, the final analysis, and the intended type I error argument. Then check whether the assumptions behind that argument match the actual design. In particular, distinguish an independent variance pilot from an interim estimate calculated using subjects whose treatment-effect information will later be analyzed.
Simulation under null parameter values is a useful diagnostic when analytic conditional distributions are difficult. Simulate the entire process—including the adaptation rule, sample-size caps, any rounding or allocation constraints, missing data handling, and final test—at the lower boundary, the upper boundary, and representative points farther into both null components. Estimate the probability of declaring equivalence at each value. Simulation can reveal inflation or implementation errors, but a finite simulation is not proof of control: its Monte Carlo uncertainty must be considered, and it cannot establish behavior for unexamined parameter values without additional reasoning.
Common failure modes include using an unadjusted fixed-sample test after an effect-informed sample-size decision; checking error only at the two boundaries without establishing conditional monotonicity; and treating blinding as a mathematical guarantee of independence. Another failure is to simulate only the boundary values while the procedure has a non-monotone rejection probability in the null interior. If these checks do not support a valid argument, the design or analysis needs adjustment rather than an assertion that the nominal level is preserved.
Check Your Understanding
Answer each question, then compare your reasoning with the solutions.
- What is the type I error event in an average bioequivalence test?
- Why does conditional error control for every adaptation outcome imply unconditional control?
- Why is conditional size at the lower boundary alone insufficient to establish control over \(\theta\leq L\)?
- A blinded variance pilot is independent of all confirmatory treatment-effect data. What additional feature of the final analysis supports the error-control argument?
- What should a null simulation include when checking an adaptive sample-size procedure?
Solutions. (1) It is declaring equivalence when \(\theta\leq L\) or \(\theta\geq U\). (2) The unconditional rejection probability is the average of the conditional probabilities over adaptation outcomes; an average of values no greater than \(\alpha\) is no greater than \(\alpha\). (3) Conditioning can change the statistic’s distribution, and its conditional rejection probability need not be monotone over the null component; an interior value could exceed the boundary value. (4) Conditional on the selected size, the confirmatory data must still follow the model used to calibrate the fixed-sample test, with the specified degrees of freedom and analysis. (5) It should reproduce the full prespecified procedure, including adaptation, sample-size limits, data reuse, allocation, missingness rules, and final analysis, and should examine boundary and interior null values.