Tutorials › Response-Adaptive Randomization › Bias in the Difference of Adaptive Arm Means

Operating Characteristics and Controversies · Tutorial 321 of 501

Bias in the Difference of Adaptive Arm Means

Derive and interpret the finite-sample bias of a difference between adaptively sampled arm means, separating null-case cancellation from bias under treatment-effect alternatives.

Advanced 10 min read

What You'll Learn

  • Define the finite-sample bias of the difference between adaptive arm means relative to the treatment-effect estimand.
  • Identify how outcome-dependent future assignments can affect arm means through their random denominators.
  • Derive an exact bias formula for a three-patient Bernoulli example.
  • Distinguish equal-mean null hypotheses from stronger distributional symmetry conditions.
  • Compare adaptive arm means with an allocation-independent benchmark and diagnose when bias may arise.

Why a difference can be biased even when each arm is randomized

A difference between treatment-arm means looks like a familiar estimator of a treatment effect. Under adaptive assignment, however, the number assigned to each arm can depend on earlier outcomes. Those outcomes also appear in the arm means. Consequently, the final denominators can be related to the observations being averaged, and the expected difference need not equal the treatment-effect estimand in a finite trial.

The point is not that every adaptive difference is biased. Its bias depends on the allocation rule, the outcome distributions, and the estimand. In particular, a design that treats the arms symmetrically can have zero expected contrast under a suitably symmetric null, while the same design can have nonzero bias under an alternative.

Definition: Let \(A_t\in\{0,1\}\) be the treatment assigned to participant \(t\), and let \(Y_t(a)\) be that participant’s outcome under arm \(a\). Suppose \(N_a=\sum_{t=1}^n I(A_t=a)\) is positive for both arms. The adaptive arm mean is \(\bar{Y}_a=\sum_{t=1}^n I(A_t=a)Y_t(a)/N_a\). For the mean-difference estimand \(\Delta=\mu_1-\mu_0\), where \(\mu_a=E[Y_t(a)]\), the finite-sample bias of the difference in arm means is \(E[\bar{Y}_1-\bar{Y}_0]-\Delta\), with expectation over outcomes and assignments generated by the design.

For clarity, examples below use independent participants whose potential outcomes have the same distribution within each arm, and an allocation probability at each time that may depend on observed history but is set before the current participant’s outcome is observed. Both arms are kept nonempty. These conditions define the examples; they do not imply that the arm means are unbiased. If later assignments respond to earlier outcomes, the final \(N_a\) can depend indirectly on outcomes already included in \(\bar{Y}_a\).

Where the finite-sample bias enters

Write \(\varepsilon_t(a)=Y_t(a)-\mu_a\). Subtracting the estimand from the observed contrast gives the exact identity

$$ E[\bar{Y}_1-\bar{Y}_0]-\Delta = E\left[\sum_{t=1}^n \frac{I(A_t=1)\varepsilon_t(1)}{N_1}\right] - E\left[\sum_{t=1}^n \frac{I(A_t=0)\varepsilon_t(0)}{N_0}\right]. $$

The identity separates the two arm-specific contributions, but it does not make either expectation zero. Although \(A_t\) is decided before \(Y_t(A_t)\) is observed, \(N_a\) includes assignments made later. Those later assignments may depend on \(Y_t(A_t)\). Thus the ratio involving the final denominator can be associated with the current centered outcome.

With fixed arm counts, or with assignments independent of all outcomes, the denominator does not select observations according to their outcomes. Under the independent-outcome model above, each arm mean then has expectation \(\mu_a\), so their difference has expectation \(\Delta\). Adaptation changes this argument: it can make inclusion in a particular arm’s final average depend on previous responses. Even when that changes both arm means, their biases may cancel in the contrast—or may not.

A small Bernoulli design that exposes the bias

Consider a hypothetical three-participant trial. The first participant is assigned to arm 1 and the second to arm 0. Their outcomes are \(X\sim\mathrm{Bernoulli}(p_1)\) and \(Z\sim\mathrm{Bernoulli}(p_0)\), independently. The third participant is assigned to arm 1 with probability \(q(X)\), where \(q(1)=q_1\) and \(q(0)=q_0\). Otherwise the third participant is assigned to arm 0. Assume \(0<q_0,q_1<1\), and that the third outcome is an independent Bernoulli outcome from its assigned arm. Both arm means are defined for every assignment because the first two allocations put one participant in each arm.

If the third participant goes to arm 0, the expected contrast conditional on \(X,Z\) is \(X-p_0\): the control mean averages \(Z\) and a fresh outcome with mean \(p_0\). If instead the third participant goes to arm 1, the conditional expected contrast is \((X+p_1)/2-Z\). The difference between these two conditional expectations is \((p_1-X)/2+(p_0-Z)/2\). Since assignment probability depends on \(X\), but not on \(Z\), averaging over \(Z\) removes its centered contribution. Therefore

$$ E[\bar{Y}_1-\bar{Y}_0] = p_1-p_0+\frac{1}{2}E[q(X)(p_1-X)] = \Delta+\frac{p_1(1-p_1)}{2}(q_0-q_1). $$

The final equality follows by splitting over \(X=1\) and \(X=0\): \(E[q(X)(p_1-X)]=p_1q_1(p_1-1)+(1-p_1)q_0p_1 =p_1(1-p_1)(q_0-q_1)\). This formula is exact for this three-participant design, not an approximation.

Worked Example: Bias Under Equal Bernoulli Means

Suppose both arms have response probability \(p_1=p_0=0.5\). The third participant has probability \(q_1=0.8\) of assignment to arm 1 after an initial arm-1 success, and probability \(q_0=0.2\) after an initial arm-1 failure. The rule is randomized at both histories, and the true mean difference is \(\Delta=0\).

Using the derived expression, the bias is

$$ \operatorname{Bias} = \frac{p_1(1-p_1)}{2}(q_0-q_1) = \frac{0.5(0.5)}{2}(0.2-0.8) = 0.125(-0.6) = -0.075. $$

Thus \(E[\bar{Y}_1-\bar{Y}_0]=-0.075\), in response-probability units, despite equal arm means. The design tends to send the third participant to arm 1 after a success and to control after a failure. In this setup that outcome-dependent rule creates unequal finite-sample selection effects.

Worked Example: Bias Under a Treatment-Effect Alternative

Keep the same allocation rule, but suppose \(p_1=0.6\) and \(p_0=0.4\), so the true difference is \(\Delta=0.2\). The design still uses \(q_1=0.8\) and \(q_0=0.2\).

Substitution gives

$$ \operatorname{Bias} = \frac{0.6(1-0.6)}{2}(0.2-0.8) = \frac{0.24}{2}(-0.6) = -0.072. $$

The expected arm-mean difference is therefore \(0.2-0.072=0.128\). This is below the true difference of \(0.2\) by \(0.072\). It does not mean that every realized trial will underestimate the effect; bias describes the average over repetitions of this specified design and outcome model.

What null symmetry does—and does not—guarantee

A null hypothesis of equal means, \(H_0:\mu_1=\mu_0\), is not by itself enough to ensure that an adaptive difference has expectation zero. The first worked example has equal Bernoulli means, and therefore identical outcome distributions, but the allocation rule favors arm 1 after a success in that arm and control after a failure. The rule is not unchanged by swapping the arm labels, so equal distributions do not force a zero expected contrast.

A stronger design symmetry can give a useful null result. Suppose the two arms have identical outcome distributions, both remain represented, and the allocation rule and estimator are invariant to swapping arm labels. Swapping the labels changes the sign of \(\bar{Y}_1-\bar{Y}_0\) but leaves its distribution unchanged. Hence its expectation equals its own negative and must be zero. This is a symmetry argument about the complete design, not a consequence of equal means alone. In particular, if the outcome distributions differ despite having equal means, identical-distribution symmetry is unavailable.

Under an alternative, even label-symmetric allocation does not generally protect the difference from bias. Treatment and control outcomes now come from different distributions, so swapping labels is not a symmetry of the data-generating process. A small symmetric-rule calculation makes the distinction concrete.

Worked Example: A Symmetric Rule Under an Alternative

Use the same initial assignments, but let the third participant’s probability of assignment to arm 1 be \(0.8\) when \((X,Z)=(1,0)\), \(0.2\) when \((X,Z)=(0,1)\), and \(0.5\) when \(X=Z\). This rule swaps its arm-1 probability with its arm-0 probability when the arm labels and their outcomes are exchanged. For this rule, the conditional contrast if the third participant is assigned to arm 0 is \(X-(Z+p_0)/2\); assigning to arm 1 changes that conditional expectation by \((p_1+p_0-X-Z)/2\).

Take \(p_1=0.7\), \(p_0=0.4\), so \(\Delta=0.3\). The four initial outcome pairs have probabilities \(0.28\), \(0.18\), \(0.42\), and \(0.12\) for \((1,1),(0,0),(1,0),(0,1)\), respectively. The expected bias is

$$ \begin{aligned} \operatorname{Bias} &=0.28(0.5)(-0.45)+0.18(0.5)(0.55)\\ &\quad+0.42(0.8)(0.05)+0.12(0.2)(0.05)\\ &=-0.063+0.0495+0.0168+0.0012=0.0045. \end{aligned} $$

The expected difference is \(0.3+0.0045=0.3045\). The bias is small here, but it is not zero. Under identical Bernoulli means, the same label-symmetric rule would give zero expected difference by the symmetry argument; that null cancellation does not extend automatically to the alternative.

Comparison, diagnostics, and limitations

The ordinary difference in adaptive arm means is transparent and targets a familiar contrast, but its finite-sample expectation depends on the interaction between outcomes and future allocation. A fixed-count randomized design avoids that particular random-denominator mechanism under the independent-outcome model: the count in each arm is set without using responses, and the arm means are unbiased for their respective means. An outcome-independent randomized allocation can also leave counts random while retaining unbiased arm means, provided assignment and counts are independent of outcomes and both means are defined. These are benchmarks, not claims that fixed allocation is preferable for every clinical objective.

Inverse-probability weighting is a different estimator family intended to address outcome-dependent assignment probabilities; it is not the same quantity as the unweighted difference studied here. Its construction and properties require separate conditions. Likewise, a small bias in this contrast does not establish good uncertainty quantification or valid testing; those are distinct operating characteristics.

In practice, the source code or written allocation rule should be inspected for whether assignment probabilities depend on earlier outcomes, whether the rule treats arm labels symmetrically, and whether both arm counts are guaranteed positive. Repeated simulations under specified null and alternative outcome models can estimate the design-specific bias by averaging \(\bar{Y}_1-\bar{Y}_0-\Delta\) across trials. Such a simulation should preserve the full sequential assignment mechanism rather than condition on a final allocation pattern.

A finite-sample bias calculation is specific to the model and policy. The Bernoulli formula above cannot simply be carried over to continuous outcomes, different histories, or other sample sizes. It also says nothing by itself about the variance of the estimator or whether an inferential procedure accounts correctly for adaptation. The key conclusion is narrower: the unweighted arm-mean difference is not automatically unbiased merely because each participant’s current assignment is randomized.

Common pitfall: “The null says the two means are equal, so the expected difference between their sample means must be zero.” This reasoning ignores the allocation rule. Equal means alone do not impose arm-label symmetry or remove outcome-dependent denominators. Check the full design and outcome assumptions; then evaluate \(E[\bar{Y}_1-\bar{Y}_0]-\Delta\), rather than inferring it from \(H_0\) alone.

Check Your Understanding

Answer each question using the distinction between the treatment-effect estimand and the expected adaptive contrast.

  1. In the three-participant Bernoulli design, what feature of the allocation rule can make the final arm-1 denominator related to the first arm-1 outcome?
  2. If \(p_1=p_0=0.4\), \(q_1=0.7\), and \(q_0=0.3\) in that design, calculate the bias.
  3. Does equal mean response guarantee zero bias for every adaptive allocation rule? State the relevant distinction.
  4. Why can a label-symmetric design have zero expected contrast under an identical-distribution null but nonzero bias under an alternative?

Solutions. (1) The third assignment probability depends on the first arm-1 response, so that response helps determine whether arm 1 receives a second participant and hence determines its final denominator. (2) The bias is \(0.4(0.6)(0.3-0.7)/2=0.12(-0.4)=-0.048\). (3) No. Equal means do not ensure that the allocation rule is label-symmetric; nor do equal means alone ensure identical outcome distributions. (4) Under the identical-distribution null, swapping labels preserves the design and data distribution while reversing the contrast, forcing its expectation to zero. Under an alternative, the arm distributions differ, so that swap no longer preserves the data-generating process and cancellation need not hold.

Clinical Biostats