Tutorials › Response-Adaptive Randomization › Sensitivity to Unequal Outcome Variances

Operating Characteristics and Controversies · Tutorial 334 of 501

Sensitivity to Unequal Outcome Variances

See when a response-rate target aligns with the allocation that minimizes the variance of a treatment contrast—and when unequal arm variances reduce precision and power.

Advanced 11 min read

What You'll Learn

  • Distinguish the response rate that drives allocation from the variance of the endpoint used for inference.
  • Derive the allocation fraction that minimizes the variance of a two-arm mean difference.
  • Compare response-rate allocation with equal allocation and variance-optimal allocation.
  • Calculate how unequal variances change a contrast’s standard error and approximate power.
  • Recognize when variance affects allocation directly, indirectly, or not at all.
  • Diagnose variance-allocation mismatches without confusing them with response-rate estimation noise.

Why Unequal Variances Matter

A response-adaptive rule may assign more patients to the arm with the higher response rate. That can be attractive when the goal is to give more patients access to an arm that appears more beneficial. But the response-rate target does not necessarily minimize the uncertainty of an estimated treatment contrast. Precision also depends on how variable the outcomes are within each arm. If the more variable arm receives too few patients, its contribution to the variance of the contrast can be large.

To isolate this issue, first consider a setting where the rule targets arm-specific response probabilities, while the primary endpoint is a quantitative outcome. The response probabilities determine the intended allocation; the quantitative endpoint has arm-specific means and variances. These quantities need not be the same. A second example considers a binary endpoint that is itself the response used by the allocation rule, where response probabilities and outcome variances are linked.

Definition: Let \(N\) be the total sample size, with \(N_1\) patients assigned to arm 1 and \(N_0=N-N_1\) assigned to arm 0. The realized allocation fraction to arm 1 is \(q=N_1/N\). Let \(p_a\) be the response probability used by the allocation rule in arm \(a\), and let the inferential endpoint \(Y\) have mean \(\mu_a\) and variance \(\sigma_a^2\) in arm \(a\). A response-rate-proportional target is \(q_{\text{resp}}=p_1/(p_0+p_1)\). This target is defined when \(p_0+p_1>0\).

A realized adaptive allocation is random, so \(q\) need not equal its target exactly. The comparisons below first treat the arm counts as fixed at their planned or realized values. They describe the variance conditional on those counts, under independent observations within and across arms, finite arm-specific variances, and an estimator for which the usual arm-mean variance calculation is appropriate. For an adaptive trial, valid inference must also account for the assignment process where required; these calculations isolate the effect of unequal outcome variances rather than replacing the earlier material on adaptive inference.

Variance of the Treatment Contrast

For a difference in sample means, \(\bar{Y}_1-\bar{Y}_0\), independent observations give the conditional variance

$$ \operatorname{Var}(\bar{Y}_1-\bar{Y}_0\mid N_1,N_0) = \frac{\sigma_1^2}{N_1}+\frac{\sigma_0^2}{N_0} = \frac{1}{N}\left(\frac{\sigma_1^2}{q}+\frac{\sigma_0^2}{1-q}\right). $$

Thus, allocation balance alone does not determine precision. A smaller arm count is more costly when that arm also has a large outcome variance. With equal variances, equal allocation minimizes the variance. With unequal variances, the minimizing allocation generally gives more patients to the arm with the larger standard deviation.

Result: For fixed total sample size and positive \(\sigma_0,\sigma_1\), the allocation fraction minimizing the variance of \(\bar{Y}_1-\bar{Y}_0\) is \(q_{\text{var}}=\sigma_1/(\sigma_0+\sigma_1)\). It allocates in proportion to the arm-specific standard deviations, not the variances.

The result assumes that both arms receive a positive fraction of the sample. In practice, allocation constraints may prevent attaining this fraction exactly. It is a precision benchmark, not a claim that variance minimization should override patient benefit, ethical constraints, or the design’s scientific goals.

Why the Variance-Optimal Fraction Has This Form

For fixed \(N\), minimizing the contrast variance is equivalent to minimizing \(g(q)=\sigma_1^2/q+\sigma_0^2/(1-q)\) over \(0<q<1\). Differentiating gives \(g'(q)=-\sigma_1^2/q^2+\sigma_0^2/(1-q)^2\). At an interior minimum, setting the derivative to zero and using positive standard deviations yields \(\sigma_1/q=\sigma_0/(1-q)\). Rearranging gives \(q=\sigma_1/(\sigma_0+\sigma_1)\). The second derivative is \(g''(q)=2\sigma_1^2/q^3+2\sigma_0^2/(1-q)^3>0\), so this stationary point is the unique minimum.

Consequently, a response-rate-proportional target is variance-optimal for this contrast only when

$$ \frac{p_1}{p_0+p_1} = \frac{\sigma_1}{\sigma_0+\sigma_1}, \qquad\text{equivalently}\qquad \frac{p_1}{p_0}=\frac{\sigma_1}{\sigma_0}, $$

assuming positive response probabilities and standard deviations. If those ratios differ, the response-rate target and the precision target differ. This is not automatically a design error: they optimize different objectives.

Worked Example: When the Two Targets Agree

Suppose a rule targets response probabilities \(p_0=0.25\) and \(p_1=0.50\). The quantitative endpoint has standard deviations \(\sigma_0=10\) points and \(\sigma_1=20\) points. With \(N=120\), compare the response-rate target with equal allocation and calculate the contrast standard error.

The response-rate target is \(q_{\text{resp}}=0.50/(0.25+0.50)=2/3\). The variance-optimal fraction is \(q_{\text{var}}=20/(10+20)=2/3\), so the targets agree. The corresponding counts are \(N_1=80\) and \(N_0=40\).

The conditional variance of the difference is \(\frac{20^2}{80}+\frac{10^2}{40}=\frac{400}{80}+\frac{100}{40}=5+2.5=7.5\) points squared. Thus \(SE=\sqrt{7.5}=2.739\) points. Under equal allocation, with 60 patients per arm, the variance is \(400/60+100/60=8.333\) points squared and the standard error is \(2.887\) points.

Here the higher-variance arm also has the higher response-rate target, and the response-based allocation improves precision relative to equal allocation. The agreement is contingent on these particular ratios; it is not a general property of response-adaptive randomization.

When the Targets Do Not Agree

A mismatch can leave the high-variance arm relatively undersampled. The resulting increase in standard error reduces the separation between the expected test statistic under a nonzero effect and its rejection threshold. Under a normal approximation for a two-sided level-\(\alpha\) test, the noncentrality parameter is approximately \(\lvert\mu_1-\mu_0\rvert/SE\). A larger standard error therefore generally means lower power at the same effect size and sample size. This is an approximation; finite-sample power depends on the test and its reference distribution.

Worked Example: A Response Target That Undersamples the Noisier Arm

Suppose \(p_0=0.40\) and \(p_1=0.60\), so a response-rate-proportional rule targets \(q_{\text{resp}}=0.60\). For a continuous endpoint, let \(\sigma_0=1\) point and \(\sigma_1=3\) points, with \(N=100\). Compare response-rate allocation, equal allocation, and the variance-optimal allocation.

At the response-rate target, \(N_1=60\) and \(N_0=40\). The contrast variance is \(3^2/60+1^2/40=9/60+1/40=0.175\) points squared, giving \(SE=\sqrt{0.175}=0.4183\) points.

With equal allocation, the variance is \(9/50+1/50=0.2000\) points squared, so \(SE=0.4472\) points. The response-rate target is more precise than equal allocation here, but it is not variance-optimal. The variance-optimal fraction is \(q_{\text{var}}=3/(3+1)=0.75\). At \(N_1=75\), \(N_0=25\), the variance is \(9/75+1/25=0.1600\) points squared, giving \(SE=0.4000\) points.

The response target places too few patients on arm 1 relative to the variance-optimal benchmark, but not so few that it performs worse than equal allocation. For a fixed effect size, its power would lie between the power under equal and variance-optimal allocation, under the stated approximate testing conditions.

Binary Outcomes: Variance Can Be Tied to Response

If the endpoint is a binary response \(Y\), then within arm \(a\), \(\mu_a=p_a\) and \(\sigma_a^2=p_a(1-p_a)\). The response rate and variance are no longer separate inputs. However, targeting allocation in proportion to \(p_a\) still does not generally minimize the variance of the difference in response proportions. The variance-optimal fraction is proportional to the standard deviations \(\sqrt{p_a(1-p_a)}\), not to \(p_a\). For example, the response rate can be high while its Bernoulli variance is modest; variance is largest at a response probability of one half.

Worked Example: A Binary Response Target Is Not Variance-Optimal

Suppose the binary endpoint is also the response used by the allocation rule, with \(p_0=0.60\) and \(p_1=0.40\), and \(N=100\). The response-rate target assigns \(q_{\text{resp}}=0.40/(0.60+0.40)=0.40\) to arm 1: 40 patients to arm 1 and 60 to arm 0.

The Bernoulli variances are \(\sigma_1^2=0.40(0.60)=0.24\) and \(\sigma_0^2=0.60(0.40)=0.24\). The contrast variance under the response target is \(0.24/40+0.24/60=0.0100\), so \(SE=0.1000\) on the proportion scale. Since the variances are equal, equal allocation is variance-optimal: its variance is \(0.24/50+0.24/50=0.009600\), giving \(SE=0.0980\).

The target favors arm 0 because its response probability is higher, although both arms have the same outcome variance. For this endpoint’s contrast, that imbalance slightly reduces precision. With a difference of 0.20, the approximate noncentrality is \(0.20/0.1000=2.000\) under the response target and \(0.20/0.0980=2.041\) under equal allocation. For a two-sided 5% normal test, these correspond to approximate powers of 0.516 and 0.532, respectively. These are illustrative normal-approximation calculations, not exact small-sample power guarantees.

What to Compare and What to Check

Unequal variances should be evaluated alongside allocation, not treated as a standalone property. A useful analysis reports the arm counts, the estimated arm-specific variances, the resulting standard error, and the intended allocation objective. For a difference in means, a heteroscedastic standard error uses the separate arm variances, \(s_1^2/N_1+s_0^2/N_0\); it does not pool the variances as if they were equal. A Welch-type test can use this standard error with an approximate degrees-of-freedom adjustment. In an adaptive trial, the design and analysis must still address the adaptive assignment mechanism.

Allocation approachWhat determines the targetImplication when variances differ
Equal allocationSame planned fraction in each armVariance-optimal only when the arm standard deviations are equal; remains a useful benchmark.
Response-rate-proportional allocationArm response probabilitiesCan improve patient-benefit allocation, but does not generally minimize contrast variance.
Variance-optimal allocationArm standard deviations for the specified contrastMinimizes the fixed-sample variance under the stated model, but may conflict with patient-benefit or other design goals.

Several diagnostics help identify a consequential mismatch:

  • Compare \(q_{\text{resp}}\) with \(q_{\text{var}}\), rather than judging imbalance from the arm counts alone.
  • Calculate the contrast variance at the response target and at relevant alternatives, keeping the total sample size fixed.
  • Check whether a small arm is also the arm with the larger endpoint variance; this combination can inflate uncertainty.
  • Assess whether variance estimates are stable enough to support the comparison. Highly uncertain estimates make a plug-in variance-optimal target unreliable.
  • When actual adaptive counts vary, examine operating characteristics across possible allocation paths rather than relying only on one expected allocation fraction.

A common mistake is to say, “The response rate is higher in arm 1, so allocating more patients there must give the most powerful test.” This confuses a potential patient-benefit objective with a precision objective. The correction is to calculate the contrast variance using both arm-specific variances and counts, then compare the resulting power under the scientific assumptions of interest.

Another mistake is to assume that unequal variances automatically invalidate a comparison. They do not. They change the appropriate standard error and can make a particular allocation inefficient; inference can still be appropriate when the analysis accommodates heteroscedasticity and the adaptive design. Conversely, merely reporting a variance-robust test does not remove the precision cost of an allocation that leaves a noisy arm small.

Key takeaway: A response-rate target allocates according to response, whereas the variance-optimal allocation for a mean contrast assigns patients in proportion to endpoint standard deviations. These targets coincide only under a specific relationship between response probabilities and outcome variability. Evaluate both precision and the design’s patient-benefit objective.

Check Your Understanding

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

  1. For a mean contrast, let \(\sigma_0=2\) and \(\sigma_1=6\). What fraction of a fixed sample minimizes the variance, and why?
  2. A response-rate target assigns \(q=0.50\), while the endpoint standard deviations are \(\sigma_0=1\) and \(\sigma_1=2\). Is the response target variance-optimal?
  3. For a binary endpoint with \(p_0=0.20\) and \(p_1=0.80\), calculate the variance-optimal allocation fraction to arm 1. Compare it with allocation proportional to response rates.
  4. Why can a larger standard error reduce power even if the estimated treatment effect is unchanged?

Solutions. (1) \(q_{\text{var}}=6/(2+6)=0.75\), so 75% goes to arm 1, the arm with the larger standard deviation. (2) No. \(q_{\text{var}}=2/(1+2)=2/3\), not 0.50; the response target undersamples the more variable arm relative to the variance-optimal benchmark. (3) Both Bernoulli variances are \(0.20(0.80)=0.16\), so the standard deviations are equal and \(q_{\text{var}}=0.50\). Proportional-to-response allocation gives \(q_{\text{resp}}=0.80/(0.20+0.80)=0.80\), which is not variance-optimal. (4) The test statistic’s expected separation from its null rejection threshold is smaller when its standard error is larger: for effect \(\Delta\), the approximate noncentrality is \(\lvert\Delta\rvert/SE\). Lower separation generally means a lower probability of rejection under the alternative.

Clinical Biostats