Tutorials › Bioavailability and Bioequivalence › Blinded Sample Size Re-Estimation from Pooled Variability

Sample Size and Power · Tutorial 140 of 399

Blinded Sample Size Re-Estimation from Pooled Variability

Use a prespecified pooled-variability estimate to adjust planned enrollment while keeping treatment identities and comparative results concealed.

Advanced 11 min read

What You'll Learn

  • Define the variability needed for sample-size re-estimation in a two-period crossover study
  • Relate residual variance to the variance of the within-subject treatment contrast
  • Recalculate enrollment using a blinded variance estimate and a prespecified rule
  • Distinguish a variance-only update from disclosure of comparative treatment results
  • Recognize when pooled variability is contaminated or too imprecise to guide enrollment
  • Compare blinded re-estimation with fixed enrollment and unblinded re-estimation

Why Re-Estimate Sample Size While Blinded?

Sample-size planning for bioequivalence depends strongly on variability. A study planned with an underestimated within-subject standard deviation may have less power than intended; one planned with an overestimated value may enroll more participants than necessary. Blinded sample-size re-estimation addresses this uncertainty by using accumulating variability data to update enrollment while keeping treatment identities concealed and withholding the estimated treatment difference from the decision process.

The method does not replace an initial design calculation. It adds a prespecified interim calculation: use a blinded estimate of the variance relevant to the primary endpoint, recalculate the required number of evaluable participants under the original design assumptions, and apply a stated minimum and maximum. The update is about variability, not whether the treatments appear equivalent.

Definition: Blinded sample-size re-estimation is a planned interim procedure that updates the target sample size using an estimate of variability obtained without disclosing treatment identities or using the estimated treatment contrast to decide enrollment. The data, estimator, timing, calculation rule, and enrollment limits should be specified before the study begins.

Which Variance Is Needed?

For a two-period crossover study, suppose the primary endpoint is analyzed on the log scale. Let \(Y_{iT}\) and \(Y_{iR}\) be participant \(i\)'s log endpoint under test and reference, and let \(D_i=Y_{iT}-Y_{iR}\) be the within-participant contrast. The sample-size calculation requires the standard deviation of these contrasts, denoted here by \(\sigma_D\). It is not the standard deviation of all individual log measurements pooled together.

A common model for a two-period, two-sequence crossover has fixed effects for sequence, period, and anonymously coded treatment, a participant effect nested within sequence, and a residual error. If residual errors in the two periods are independent with common variance \(\sigma_e^2\), the participant effect cancels from the within-participant contrast, and \(D_i\) has residual variance \(2\sigma_e^2\). Thus \(\sigma_D=\sqrt{2\sigma_e^2}\). A blinded analyst can estimate the residual variance from the prespecified model while the treatment labels remain masked, and provide the variance estimate without releasing a treatment contrast.

These relationships assume the model is appropriate: the endpoint is analyzed on the specified scale; participants are independent; residual variance is sufficiently stable across sequences and periods; and the residual errors have the structure used in the calculation. If the design or variance structure differs, the estimator must match that design rather than borrowing this conversion automatically.

Key relationship: In the equal-variance two-period crossover model, if \(MS_{\mathrm{res}}\) estimates the residual variance, then \(\hat{\sigma}_e^2=MS_{\mathrm{res}}\) and \(\hat{\sigma}_D=\sqrt{2MS_{\mathrm{res}}}\). The variance relevant to enrollment is the variance of the treatment contrast, not necessarily the residual variance itself.

Worked Example: Converting Residual Variance to Contrast Variability

Suppose an interim blinded analysis of a hypothetical two-period crossover study gives \(MS_{\mathrm{res}}=0.0648\) on the log scale. Assume the prespecified model has equal, independent residual errors across the two periods. Find the estimated standard deviation of the within-participant contrast.

The residual variance estimate is \(\hat{\sigma}_e^2=0.0648\). The contrast variance is twice this amount:

$$ \hat{\sigma}_D^2=2MS_{\mathrm{res}}=2(0.0648)=0.1296. $$

Taking the square root gives

$$ \hat{\sigma}_D=\sqrt{0.1296}=0.3600 $$

log units. This is the estimate to use in the variance-based update for the contrast-based design, subject to the original power calculation and any prespecified enrollment limits. It is not an estimate of the treatment difference.

How the Enrollment Update Works

Let \(n\) denote the number of evaluable participants, \(\alpha\) the one-sided significance level for each equivalence test, and \(1-\beta\) the target power. Let \(L\) and \(U\) be the lower and upper equivalence limits on the log scale, and let \(\theta_0\) be the assumed true mean contrast for planning. A normal-approximation calculation illustrates the dependence of sample size on variance:

$$ n \approx \left[ \frac{(z_{1-\alpha}+z_{1-\beta})\sigma_D} {\min(U-\theta_0,\ \theta_0-L)} \right]^2. $$

Here \(z_q\) is the \(q\) quantile of the standard normal distribution. The expression uses the distance from the assumed true contrast to the nearer equivalence boundary; it is a planning approximation, not a substitute for an exact power calculation where that is required. With the margins, \(\alpha\), target power, and \(\theta_0\) held fixed, it shows why the required sample size scales with \(\sigma_D^2\). In practice, the planned calculation should use the same power method before and after re-estimation.

Re-estimation rule: If the planning calculation is variance-scaled and all other inputs remain fixed, an interim estimate \(\hat{\sigma}_{D,\mathrm{new}}\) updates the initial target \(n_{\mathrm{old}}\) approximately as \(n_{\mathrm{new}}=n_{\mathrm{old}}(\hat{\sigma}_{D,\mathrm{new}}/\hat{\sigma}_{D,\mathrm{old}})^2\). Round according to the prespecified rule, respect the minimum and maximum, and use the design’s allocation requirements. For an exact power calculation, recompute power with the new variance and find the smallest eligible sample size meeting the target.

Worked Example: Updating the Target After a Variance Increase

A hypothetical crossover study initially plans for 24 evaluable participants using \(\hat{\sigma}_{D,\mathrm{old}}=0.300\) log units. At the prespecified interim point, a blinded estimate is \(0.360\) log units. For this illustration, use the variance-scaled approximation, keep all other planning inputs fixed, round up to a whole participant, and impose a maximum of 40 evaluable participants.

The variance ratio is

$$ \left(\frac{0.360}{0.300}\right)^2 =(1.200)^2=1.440. $$

The unrounded updated target is

$$ n_{\mathrm{new}}=24(1.440)=34.56. $$

Rounding up gives 35 evaluable participants, below the maximum of 40. The variance increase therefore raises the target by 11 evaluable participants. This calculation does not show that the treatments differ; it only shows that the interim variance estimate is larger than the value used for initial planning.

If the study separately expects 10% of randomized participants not to be evaluable, an operational target intended to yield 35 evaluable participants is \(\lceil 35/0.90\rceil=\lceil38.89\rceil=39\) randomized participants. That dropout allowance is a separate planning step: the re-estimation target above is in evaluable participants.

Derivation of the Variance-Scaling Rule

The scaling rule follows directly from the normal-approximation formula when all quantities except \(\sigma_D\) are fixed. Write \(c=(z_{1-\alpha}+z_{1-\beta})/\min(U-\theta_0,\theta_0-L)\). Then \(n\approx c^2\sigma_D^2\). Applying the same expression to the old and new standard deviations gives

$$ \frac{n_{\mathrm{new}}}{n_{\mathrm{old}}} \approx \frac{c^2\hat{\sigma}_{D,\mathrm{new}}^2} {c^2\hat{\sigma}_{D,\mathrm{old}}^2} = \left(\frac{\hat{\sigma}_{D,\mathrm{new}}} {\hat{\sigma}_{D,\mathrm{old}}}\right)^2. $$

The constant \(c^2\) cancels. This derivation explains both the convenience and the limitation of the shortcut: it is valid for the stated approximation with unchanged design inputs. Small-sample power calculations, integer sequence balance, or a changed target power can prevent this ratio from giving the exact revised sample size. In those settings, recompute power using the original method rather than assuming the scaling identity is exact.

Keeping the Update Blinded

Blinding is a data-handling procedure, not a property guaranteed by the word “pooled.” A defensible process keeps treatment identities masked, uses a prespecified model that accounts for the design’s treatment and period structure, and returns only the variance information needed for the size decision. The people deciding whether to expand enrollment should not receive treatment-specific summaries or the estimated treatment contrast.

The decision rule also needs clear operational details. Specify when the interim data are available, how many participants or complete profiles are required, how the variance estimator is calculated, the target power and other inputs that remain fixed, how rounding and sequence balance are handled, and the minimum and maximum enrollment. A rule such as “increase if the interim result looks unfavorable” is not blinded variance re-estimation: it is an undefined response to accumulating results.

Worked Example: Why Pooling Raw Measurements Can Mislead

Consider a hypothetical equal-allocation parallel study with two masked groups. Suppose each group has within-group variance \(0.0400\) on the analysis scale, and the two group means differ by \(0.300\), although the identities and direction of the difference are concealed. What is the expected variance of all observations around the grand mean if the groups are simply pooled?

With equal group sizes, each group mean is \(0.150\) from the grand mean. The pooled variance contains within-group variance plus the between-group contribution:

$$ 0.0400+(0.150)^2=0.0400+0.0225=0.0625. $$

The resulting standard deviation is \(\sqrt{0.0625}=0.250\), rather than the within-group standard deviation \(\sqrt{0.0400}=0.200\). Thus, pooling raw observations can mix treatment separation with residual variability. It may produce a larger estimate even though no treatment identity or direction has been disclosed. The appropriate blinded estimator depends on the design and must separate the variance component relevant to the planned test from mean structure as specified in advance.

Diagnostics, Alternatives, and Limitations

Before using an interim estimate, check whether it is numerically and scientifically credible. Review the number of complete profiles contributing to it, model convergence, residual variation by sequence and period, influential observations, protocol deviations, and whether the endpoint transformation and analysis population match the plan. A striking estimate based on very few participants may reflect ordinary sampling noise; a variance estimate is itself uncertain. Data quality problems should be investigated rather than treated automatically as evidence that enrollment must increase.

A blinded update is not the only option. Fixed enrollment avoids interim estimation and its operational complexity, but cannot respond to a meaningful mismatch between planned and observed variability. Unblinded re-estimation may use more direct information about treatment-specific effects, but it reveals comparative information to those conducting the update and requires additional safeguards. Blinded re-estimation preserves concealment of treatment identities and the comparative result, but its validity still depends on the estimator, model, and prespecified decision rule; “blinded” does not by itself guarantee correct statistical behavior.

This procedure also has limits. It cannot repair a misspecified endpoint, biased measurements, incorrect equivalence limits, or a design with inadequate power for reasons other than variance. A hard maximum may leave the study below its target power if variability is higher than anticipated. A cap can be necessary for feasibility, but its consequences should be evaluated during planning. Nor should the updated estimate be interpreted as evidence for or against bioequivalence: the interim calculation concerns enrollment, while the final equivalence conclusion follows from the prespecified final analysis.

Common Mistake

  • Wrong reasoning: “The pooled data have a larger standard deviation, so the treatments must be different; I will use the observed difference to decide how many participants to add.”
  • Correction: A blinded sample-size update uses a prespecified estimate of the variance relevant to the design and does not use the treatment contrast to make the enrollment decision. Pooled raw variability can include mean differences, so use a design-appropriate estimator and keep comparative summaries concealed.
  • Wrong reasoning: “The standard deviation increased by 20%, so the sample size also increases by 20%.”
  • Correction: Under the variance-scaled approximation, sample size is proportional to the square of the standard deviation. A 20% increase gives a factor of \(1.20^2=1.44\), or approximately a 44% increase before rounding and applying limits.
Key takeaway: Blinded re-estimation can adapt enrollment to accumulating variability without disclosing comparative treatment results, but only when the variance estimator and update rule are matched to the design, prespecified, and separated from treatment-effect interpretation.

Check Your Understanding

Answer each question, then compare your reasoning with the solutions.

  1. A crossover model gives \(MS_{\mathrm{res}}=0.0500\). Under the equal, independent residual-error assumption, what is the estimated standard deviation of the within-participant contrast?
  2. An initial plan requires 30 evaluable participants at a contrast standard deviation of \(0.250\). The interim estimate is \(0.300\). Using variance scaling, what is the unrounded revised target?
  3. Why is the variance of all raw measurements not automatically a valid estimate of within-subject contrast variability?
  4. Name two safeguards that help keep a blinded re-estimation from revealing comparative treatment results.

Solutions. (1) \(\hat{\sigma}_D=\sqrt{2(0.0500)}=\sqrt{0.1000}=0.3162\) log units. (2) The variance ratio is \((0.300/0.250)^2=(1.200)^2=1.440\), so the unrounded target is \(30(1.440)=43.20\) evaluable participants; the prespecified rounding and maximum rules determine the operational target. (3) Raw pooled variability can include differences in group means or other modeled effects, so it need not estimate the variance component relevant to the contrast. (4) Examples include keeping treatment identities masked, using a prespecified design-appropriate variance model, returning only the variance estimate to the enrollment decision group, and withholding treatment-specific summaries and the estimated contrast.