Tutorials › Response-Adaptive Randomization › Comparing Continuous-Outcome Power

Comparison With Fixed Randomization · Tutorial 82 of 501

Comparing Continuous-Outcome Power

Learn to calculate and compare power for a continuous endpoint under equal and unequal fixed allocation, and interpret the result under a stated model.

Intermediate 10 min read

What You'll Learn

  • State the common-variance normal model and the assumptions for the two-sample z test.
  • Calculate the standard error of a difference in means for fixed arm sizes.
  • Use the noncentrality parameter and critical value to calculate two-sided power.
  • Compare power for equal and unequal allocation at the same total sample size.
  • Explain how effect size, outcome variability, significance level, and allocation affect power.
  • Recognize when the common-variance known-variance calculation is not appropriate.

Why Calculate Continuous-Outcome Power?

For a continuous endpoint, the treatment effect is often expressed as a difference between two population means: for example, a change of 5 points on a symptom scale or 2 mmHg in mean blood pressure. A power calculation asks how often a planned test would reject the null hypothesis if a specified difference were truly present. Comparing equal and unequal fixed allocation makes the precision cost of an allocation imbalance explicit.

The comparison is meaningful only when the designs are evaluated under the same assumptions: the same total sample size, treatment effect, outcome variance, significance level, test, and analysis. Changing any of these while comparing allocation ratios can make it unclear which design feature caused a power difference.

Model and Power Formula

Consider a two-arm trial with \(n_1\) participants assigned to treatment 1 and \(n_0\) assigned to treatment 0, with fixed sample sizes \(n_1\) and \(n_0\). The total sample size is \(N=n_1+n_0\). Let \(X_{ij}\) be the continuous outcome for participant \(i\) in arm \(j\), where \(j\) is 0 or 1.

Model and assumptions: Assume independent outcomes within and between arms, with \(X_{ij}\sim N(\mu_j,\sigma^2)\) in each arm. Thus the two arms have a common population variance \(\sigma^2\). The variance is treated as known for the z-test calculation below, participants contribute independent observations, and the planned arm sizes are fixed. The treatment effect is \(\Delta=\mu_1-\mu_0\). The hypotheses for a two-sided test are \(H_0:\Delta=0\) and \(H_a:\Delta\ne 0\).

Under these assumptions, the difference in sample means, \(\bar{x}_1-\bar{x}_0\), has mean \(\Delta\) and standard error \(SE=\sigma\sqrt{1/n_1+1/n_0}\). For a two-sided test at significance level \(\alpha\), reject \(H_0\) when the absolute value of the standardized statistic exceeds \(z_{1-\alpha/2}\), the \(1-\alpha/2\) quantile of the standard normal distribution.

$$ Z=\frac{\bar{x}_1-\bar{x}_0}{\sigma\sqrt{1/n_1+1/n_0}}, \qquad c=z_{1-\alpha/2}. $$

Under a specified true effect \(\Delta\), this statistic has a normal distribution with variance 1 and mean \(\lambda=\Delta/SE\). The quantity \(\lambda\) is the signed noncentrality parameter: it expresses the true difference in standard-error units. Power is the probability, under that specified alternative, that the test statistic falls outside the null rejection region.

$$ \text{Power}(\Delta) =P(Z>c\mid\Delta)+P(Z<-c\mid\Delta) =1-\Phi(c-\lambda)+\Phi(-c-\lambda), \qquad \lambda=\frac{\Delta}{\sigma\sqrt{1/n_1+1/n_0}}. $$

Here \(\Phi\) is the standard normal cumulative distribution function. The formula includes both tails, including when \(\Delta\) is negative. For a positive effect large relative to its standard error, the lower-tail probability is often very small, but it is still part of the two-sided power.

Deriving the Standard Error and Power

The standard error follows from the variance of a difference of independent sample means. Under the model, each sample mean has variance \(\sigma^2/n_j\). Independence means the covariance between the sample means is zero, so their difference has variance equal to the sum of those variances. Dividing by its standard deviation gives the stated standardized statistic.

$$ \operatorname{Var}(\bar{x}_1-\bar{x}_0) =\frac{\sigma^2}{n_1}+\frac{\sigma^2}{n_0} =\sigma^2\left(\frac{1}{n_1}+\frac{1}{n_0}\right). $$

To obtain power, standardize the difference in means using this standard error. Under the alternative its expectation is \(\Delta/SE=\lambda\), and its variance is 1. Thus \(Z\) has distribution \(N(\lambda,1)\) under the alternative. The probability of exceeding \(c\) is \(1-\Phi(c-\lambda)\); the probability of falling below \(-c\) is \(\Phi(-c-\lambda)\). Adding these disjoint rejection-tail probabilities gives the power formula.

This calculation makes the inputs visible. A larger absolute effect increases \(|\lambda|\); a larger \(\sigma\) or larger standard error reduces it. Increasing the total sample size generally reduces the standard error. A smaller \(\alpha\) makes the critical value \(c\) larger, so rejection is harder and power is lower when the other quantities stay fixed.

Worked Example: Equal Allocation With a Specified Effect

Suppose a hypothetical trial plans \(N=100\) participants, allocated equally, and the continuous outcome has common standard deviation \(\sigma=10\) points. The design targets a true mean difference of \(\Delta=5\) points and uses a two-sided test with \(\alpha=0.05\). Calculate power under the stated model.

The arm sizes are \(n_1=n_0=50\). The standard error and noncentrality parameter are

$$ SE=10\sqrt{\frac{1}{50}+\frac{1}{50}} =10\sqrt{0.04}=2.000\text{ points}, \qquad \lambda=\frac{5}{2}=2.500. $$

For \(\alpha=0.05\), \(c=z_{0.975}\approx1.960\). Therefore,

$$ \text{Power} =1-\Phi(1.960-2.500)+\Phi(-1.960-2.500) =1-\Phi(-0.540)+\Phi(-4.460) \approx0.705. $$

Under this model and if the true difference is 5 points, the planned two-sided test has about 70.5% power. This is not the probability that the treatment works, nor the probability that the null hypothesis is true; it is the rejection probability calculated under the specified alternative.

Comparing Equal and Unequal Allocation

At fixed total sample size, allocation determines the standard error through \(1/n_1+1/n_0\). To see why equal allocation minimizes that quantity, use \(n_1+n_0=N\):

$$ \frac{1}{n_1}+\frac{1}{n_0} =\frac{n_1+n_0}{n_1n_0} =\frac{N}{n_1n_0}. $$

For a fixed sum \(N\), the product \(n_1n_0\) is largest when the two arm sizes are equal (or as close to equal as integers allow). Consequently, equal allocation gives the smallest standard error in this common-variance model when the total sample size is held fixed. With a nonzero effect and the same two-sided test, a smaller standard error generally gives greater power. Unequal allocation can still be chosen for reasons beyond this particular power comparison, but its power should be calculated rather than assumed to match equal allocation.

Worked Example: Power With a 70:30 Allocation

Keep the hypothetical trial conditions from the previous example: \(N=100\), \(\sigma=10\) points, \(\Delta=5\) points, and a two-sided test with \(\alpha=0.05\). Now assign \(n_1=70\) and \(n_0=30\). Calculate power and compare it with equal allocation.

The standard error is

$$ SE=10\sqrt{\frac{1}{70}+\frac{1}{30}} =10\sqrt{0.047619}=2.182\text{ points}. $$

The noncentrality parameter is \(\lambda=5/2.182=2.291\). Using \(c=1.960\),

$$ \text{Power} =1-\Phi(1.960-2.291)+\Phi(-1.960-2.291) =1-\Phi(-0.331)+\Phi(-4.251) \approx0.630. $$

The calculated power is about 63.0%, compared with about 70.5% under equal allocation with the same \(N\), effect, variance, and test. The larger standard error in the 70:30 design explains the lower power under the assumed model. The figures do not imply that unequal allocation is always inappropriate; they quantify the power difference for these particular conditions.

Worked Example: A Different Effect and Allocation Ratio

A hypothetical trial plans \(N=120\) participants, with 80 assigned to treatment 1 and 40 to treatment 0. Suppose the common standard deviation is \(\sigma=12\) mmHg, the mean difference of interest is \(\Delta=4\) mmHg, and the planned test is two-sided with \(\alpha=0.05\). Calculate the power.

First calculate the standard error in the endpoint's units:

$$ SE=12\sqrt{\frac{1}{80}+\frac{1}{40}} =12\sqrt{0.0375} =2.324\text{ mmHg}. $$

Then \(\lambda=4/2.324=1.721\). With \(c=1.960\), the power is

$$ \text{Power} =1-\Phi(1.960-1.721)+\Phi(-1.960-1.721) =1-\Phi(0.239)+\Phi(-3.681) \approx0.406. $$

The design has about 40.6% power for a true 4 mmHg difference under the model. The standardized effect is \(\Delta/\sigma=4/12=0.333\), but the final power calculation also depends on the sample size and allocation through the standard error. Reporting the effect in mmHg and the standard deviation in mmHg keeps the assumed effect interpretable in context.

Checking a Power Comparison

Before interpreting the result, check that the calculation matches the planned design and analysis. In particular, confirm the following:

  • Endpoint and effect: State the outcome scale and the target mean difference \(\Delta\) in its units. A power value without a specified effect is incomplete.
  • Variance: Specify the assumed common standard deviation \(\sigma\), not just a variance estimate from an unrelated outcome or scale.
  • Allocation and total: Verify \(n_1+n_0=N\) and use the actual arm sizes in \(SE\). Do not substitute \(N\) into an equal-allocation formula when the planned allocation is unequal.
  • Test and significance level: Match the critical value to the planned one-sided or two-sided test and to \(\alpha\). This tutorial's main formula is two-sided.
  • Model: Consider whether independence, approximate normality of the sample means, and the common-variance assumption are defensible for the planned analysis.

The calculations above use a z test with known \(\sigma\). In practice, \(\sigma\) is commonly estimated and the analysis may use a t test; that procedure has a different reference distribution, especially when samples are small. The next tutorial on small-sample power comparisons addresses that setting. If arm variances differ, the common-variance standard error used here is not the appropriate general calculation. Likewise, if allocation is changed in response to accumulating outcomes, fixed arm sizes do not describe the design, so this fixed-allocation calculation alone does not characterize its power.

Common Mistakes and Design Interpretation

Common mistake: “The unequal-allocation trial has the same total sample size, so it must have the same power.” Equal total sample sizes do not guarantee equal standard errors. The correction is to use \(SE=\sigma\sqrt{1/n_1+1/n_0}\) with the actual arm sizes, then calculate power using the specified effect and test.

A second error is to report a power number without its assumptions. For example, “power is 80%” is not a complete design statement unless the assumed effect, variance, allocation, total sample size, and test are identified. Power is a property of the specified model and testing procedure under a specified alternative, not a permanent property of a treatment or endpoint.

The main comparison here is between fixed equal and fixed unequal allocation at matched conditions. Equal allocation minimizes the standard error for a difference in means under the common-variance model at fixed \(N\). Unequal allocation raises that standard error, with the amount depending on the ratio. A fair comparison preserves the other inputs and reports both the power values and the assumptions that produced them.

Check Your Understanding

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

  1. A hypothetical equal-allocation trial has \(N=64\), common standard deviation \(\sigma=8\) points, target difference \(\Delta=4\) points, and a two-sided \(\alpha=0.05\) test. Calculate its power.
  2. For \(n_1=48\), \(n_0=32\), \(\sigma=6\) units, and \(\Delta=3\) units, calculate the standard error and two-sided power at \(\alpha=0.05\).
  3. Why is it incorrect to use the common-variance formula without qualification when the treatment and control outcome variances are substantially different?
  4. With the effect, variance, allocation, and total sample size held fixed, what happens to two-sided power when \(\alpha\) decreases from 0.05 to 0.01? Explain using the rejection threshold.

Solutions. (1) The arms each have 32 participants. Thus \(SE=8\sqrt{1/32+1/32}=2.000\) points and \(\lambda=4/2=2.000\). With \(c=1.960\), power is \(1-\Phi(-0.040)+\Phi(-3.960)\approx0.516\), or 51.6%. (2) \(SE=6\sqrt{1/48+1/32}=6\sqrt{0.052083}=1.369\) units. Then \(\lambda=3/1.369=2.191\), and power is \(1-\Phi(1.960-2.191)+\Phi(-1.960-2.191)\approx0.591\), or 59.1%. (3) The formula assumes a common variance, so it does not represent the sampling variance of the difference under substantially unequal arm variances. The analysis and power calculation must instead match a method and variance model that account for that difference. (4) Power decreases: reducing \(\alpha\) raises the two-sided critical value from about 1.960 to 2.576, making the rejection region narrower under the alternative.

Clinical Biostats