Why Baseline Correlation Matters for Enrollment
A baseline measurement can improve the precision of a continuous-outcome comparison when it predicts the outcome at follow-up. In an analysis of covariance (ANCOVA), the outcome is adjusted for the baseline value; the resulting residual variance can be smaller than the variance used in an unadjusted comparison. That reduction can translate into a lower planned enrollment for the same power, provided the design and analysis assumptions are appropriate.
In a blinded internal pilot, treatment labels are unavailable for the sample-size update. Baseline and outcome measurements are still paired for each participant, so their covariance can be estimated without revealing assignments. But the raw pooled outcome variance includes variation due to the treatment effect. A raw pooled baseline-outcome correlation can therefore understate the within-treatment correlation relevant to ANCOVA. The estimate must use a treatment-effect correction consistent with the design assumptions.
Model and Assumptions
Consider a randomized two-group trial. Suppose the baseline and outcome follow a common-slope model:
Here \(g\) identifies the treatment group, \(\mu_g\) is its outcome mean at the baseline reference value, \(\beta\) is the common baseline slope, and \(\mu_X\) is the baseline mean. Assume that the baseline distribution is the same across groups in the target population, randomization makes treatment assignment independent of baseline, and the residual variance is common. For the variance reduction below, the linear relationship should be an adequate description over the observed baseline range, with residuals uncorrelated with \(X\). Baseline must be measured before treatment, and missingness or selection must not create a different baseline-outcome relationship in the analyzed data.
Within a treatment group, write the marginal outcome variance as \(\sigma_Y^2\). The model gives \(\sigma_Y^2=\beta^2\sigma_X^2+\sigma^2\), where \(\sigma_X^2=\operatorname{Var}(X)\). The covariance is \(\operatorname{Cov}(X,Y)=\beta\sigma_X^2\), so the correlation satisfies \(\rho^2=\beta^2\sigma_X^2/\sigma_Y^2\). These identities connect the correlation to the variance remaining after linear adjustment.
Deriving the ANCOVA Variance Factor
The adjusted outcome removes the part of \(Y\) linearly predicted by \(X\). Under the common-slope model, that predicted component has variance \(\beta^2\sigma_X^2\). The residual variance after adjustment is therefore
Thus, relative to an unadjusted comparison using the within-treatment outcome variance \(\sigma_Y^2\), ANCOVA uses the variance fraction \(1-\rho^2\). In a standard two-group sample-size calculation with the same target treatment effect, significance level, power, and allocation ratio, sample size is proportional to the outcome variance. The adjusted-design target is consequently approximately the unadjusted target multiplied by \(1-\rho^2\). This is a planning relationship, not a guarantee that a realized trial will attain that exact reduction.
This approximation presumes the baseline covariate is available for the participants included in the analysis, and that the analysis uses the corresponding ANCOVA adjustment. Finite-sample degrees of freedom, rounding, missing baseline data, and other design constraints can affect the final enrollment target. The correlation-based reduction should not be applied to a design that will not actually adjust for baseline.
Estimating the Correlation While Treatment Is Masked
Let \(\pi\) be the planned fraction allocated to one group, and let \(\Delta_0\) be the treatment effect assumed for the blinded variance correction. If the treatment effect is constant and assignment is independent of baseline, the population pooled outcome variance is the within-treatment outcome variance plus the between-treatment contribution \(\pi(1-\pi)\Delta_0^2\). The pooled baseline-outcome covariance has no corresponding treatment-effect contribution in expectation, because baseline is independent of randomized assignment.
Using masked sample summaries, a practical moment-based estimate is therefore
Here \(s_X^2\) is the pooled baseline variance, \(s_{Y,\mathrm{mask}}^2\) is the pooled outcome variance with treatment labels masked, and \(s_{XY}\) is the pooled baseline-outcome covariance computed on the paired observations. The corrected variance \(\widehat{\sigma}_{Y,w}^{\,2}\) estimates the within-treatment marginal outcome variance, not the ANCOVA residual variance. The covariance and baseline variance then determine the squared correlation. Small-sample implementations may need design-specific refinements, but the distinction between the raw pooled variance and the within-treatment variance remains essential.
The correction relies on the assumed treatment effect and allocation fraction. In a finite randomized sample, chance baseline imbalance can contribute to the pooled covariance, and unequal allocation changes the contamination term. If the effect varies with baseline, treatment groups have different baseline distributions, or outcome variance differs by group, this simple pooled-moment calculation may not estimate a single common within-group correlation reliably.
Worked Example: Estimating Correlation From Masked Summaries
Suppose a hypothetical internal pilot has masked summaries \(s_X^2=100\) baseline units squared, \(s_{XY}=30\) baseline units times outcome points, and \(s_{Y,\mathrm{mask}}^2=169\) outcome points squared. The planned allocation is equal, so \(\pi=0.5\), and the assumed treatment effect for the correction is \(\Delta_0=6\) outcome points. Assume randomization, a constant effect, and a common baseline-outcome slope and residual variance.
First remove the assumed between-treatment contribution from the pooled outcome variance:
Then estimate the correlation:
The estimated squared correlation is \(\hat{\rho}^{\,2}=0.05625\). The corresponding estimated variance factor is \(1-0.05625=0.94375\). In this hypothetical example, adjustment for baseline is estimated to reduce the within-treatment outcome variance by about 5.63%, relative to using that variance without adjustment. This is a planning estimate, not evidence that the treatment effect is 6 points.
From Variance Reduction to Enrollment
When all other sample-size inputs are held fixed, multiply the unadjusted total enrollment target by the estimated variance factor. The resulting number is a starting target for an adjusted design; apply the trial's allocation, minimum-enrollment, and rounding rules afterward. If those rules require an even total for equal allocation, round up to an even number so the planned target is not undershot.
Worked Example: Updating the Planned Total
Suppose a hypothetical unadjusted two-group design requires 240 participants in total. Use the correlation estimate from the preceding example, and suppose the final analysis will use ANCOVA with the same baseline measurement. Assume no additional minimum or attrition adjustment is needed for this calculation.
The unrounded adjusted target is
For equal allocation, the planned total must be an even integer and must not be below 226.5. Rounding up gives 228 participants, or 114 per group. This calculation assumes the estimated variance factor is suitable for planning and that no other sample-size constraint overrides the result.
Worked Example: Seeing the Cost of Using Raw Pooled Correlation
Consider another hypothetical pilot with \(s_X^2=25\) baseline units squared, \(s_{XY}=20\) baseline units times outcome points, and masked pooled outcome variance \(s_{Y,\mathrm{mask}}^2=125\) outcome points squared. Allocation is equal and the assumed treatment effect is 10 outcome points. Assume the same common-slope and constant-effect conditions as above.
The treatment-effect contribution is \((0.5)(0.5)(10^2)=25\) outcome points squared, so the corrected within-treatment outcome variance is \(125-25=100\) points squared. The corrected correlation estimate is
If someone instead uses the raw pooled variance in the denominator, the result is \(20/\sqrt{(25)(125)}\approx0.3578\), with square about \(0.1280\). That would imply a variance factor of about \(0.8720\), rather than the corrected factor \(0.8400\). For an unadjusted total target of 300, the corrected calculation gives \(300(0.8400)=252\) participants; the raw-correlation calculation gives about \(300(0.8720)=261.6\), or 262 after rounding up. The raw calculation has diluted the apparent correlation by retaining treatment-effect variation in the outcome denominator.
Diagnostics, Alternatives, and Limitations
Check that the corrected within-treatment outcome variance is positive and that the estimated correlation is in the range \([-1,1]\). Also check the implied ANCOVA residual variance:
A negative value, or \(|\hat{\rho}|>1\), is not a meaningful variance reduction. It indicates that the observed masked moments and the assumed correction are incompatible, or that sampling variability is substantial. Investigate the summaries, assumed effect, allocation, outliers, data pairing, and missingness rather than silently forcing the estimate into range. A value near zero can also be unstable when the internal pilot is small; reporting it to many decimal places does not make it precise.
The principal alternatives answer different planning questions. An unadjusted endpoint design avoids relying on an estimated baseline correlation, but does not claim the ANCOVA efficiency gain. An ANCOVA design using external information about correlation avoids a noisy internal-pilot estimate, but depends on whether that external population and measurement process are relevant. A raw masked correlation is easy to calculate but is not a suitable substitute when treatment-effect variation inflates the pooled outcome variance. The masked correction is useful when its assumptions are credible, yet it inherits uncertainty both from the covariance estimate and from the assumed effect used in the variance correction.
The simple factor \(1-\rho^2\) is limited to the common linear adjustment and variance structure described here. A nonlinear baseline-outcome relationship, treatment-by-baseline interaction, markedly different baseline distributions, heteroscedasticity, or selective loss to follow-up can invalidate the single-correlation summary. The estimate also does not establish that the treatment effect is absent or present; it is a nuisance-parameter estimate for planning. In an interim review, document the masking, correction assumptions, estimated factor, and any decision limits applied to the recalculated enrollment.
Check Your Understanding
Answer each question, then compare your reasoning with the solutions.
- In one sentence, explain why a raw pooled masked outcome variance may be too large for estimating the within-treatment baseline-outcome correlation.
- A hypothetical pilot has \(s_X^2=16\), \(s_{XY}=12\), \(s_{Y,\mathrm{mask}}^2=100\), equal allocation, and assumed effect \(\Delta_0=8\). Find the corrected variance, \(\hat{\rho}\), and the ANCOVA variance factor.
- If an estimated squared correlation is 0.25 and the unadjusted target is 400 participants, what is the unrounded adjusted target, before design-specific rounding?
- What should be concluded if the correction produces \(\widehat{\sigma}_{Y,w}^{\,2}\leq 0\) or \(|\hat{\rho}|>1\)?
Solutions. (1) The pooled variance can include between-treatment variation from the treatment effect, while the relevant correlation uses within-treatment outcome variance. (2) The contamination is \((0.5)(0.5)(8^2)=16\), so the corrected variance is \(100-16=84\). Then \(\hat{\rho}=12/\sqrt{(16)(84)}=12/36.66\approx0.3273\), \(\hat{\rho}^{\,2}\approx0.1071\), and the variance factor is about \(0.8929\). (3) \(400(1-0.25)=300\) participants. (4) The moments and assumptions do not support a valid correlation or residual-variance estimate; investigate the data and correction rather than clipping the estimate without a pre-specified rule.