Tutorials › Response-Adaptive Randomization › Normalizing Multi-Arm Allocation Weights

Doubly Adaptive Biased Coin Designs · Tutorial 154 of 501

Normalizing Multi-Arm Allocation Weights

Learn how to calculate and verify normalized multi-arm DBCD probabilities from target shares, current allocation shares, and the tuning parameter.

Advanced 10 min read

What You'll Learn

  • Define target and current allocation shares for each treatment arm.
  • Derive the normalized multi-arm DBCD assignment probability.
  • Verify that the probabilities are nonnegative and sum to one.
  • Calculate allocation probabilities for three-arm examples.
  • Recognize zero-share and normalization errors that make a calculation invalid.
  • Compare ratio-corrected probabilities with direct randomization at the target shares.

Why Normalization Matters

With multiple treatment arms, the allocation rule must turn information about all arms into a single set of probabilities for the next participant. A target share says what long-run allocation the design aims to track; the current shares describe the allocation so far. The DBCD rule uses their ratios to adjust the next assignment probabilities. Because every arm’s probability depends on the other arms’ weights, those adjusted weights must be normalized jointly.

This tutorial treats the target allocation vector as given at the current assignment. It focuses on the calculation that converts that vector and the current allocation shares into probabilities. It does not require a particular method for constructing the target from response estimates.

Shares, Ratios, and Weights

Suppose there are \(K\) treatment arms. Before participant \(n+1\) is assigned, let \(N_k(n)\) be the number assigned to arm \(k\), and let \(n=\sum_{k=1}^K N_k(n)\) be the number already assigned. The current allocation share of arm \(k\) is \(x_k(n)=N_k(n)/n\). Let \(\pi_k\) be its target share, with \(\sum_{k=1}^K \pi_k=1\). The tuning parameter is \(\gamma\geq 0\).

Definition: For positive target shares and positive current allocation shares, the target-to-current ratio for arm \(k\) is \(\pi_k/x_k(n)\). Its DBCD allocation weight is \(\pi_k\{\pi_k/x_k(n)\}^{\gamma}\). The next-participant assignment probability for arm \(k\) is that weight divided by the sum of the weights over all \(K\) arms.

Writing \(p_k(n+1)\) for the probability of assigning participant \(n+1\) to arm \(k\), the rule is

$$ p_k(n+1) = \frac{\pi_k\left\{\pi_k/x_k(n)\right\}^{\gamma}} {\displaystyle\sum_{j=1}^{K}\pi_j\left\{\pi_j/x_j(n)\right\}^{\gamma}}, \qquad k=1,\ldots,K. $$

Equivalently, define the unnormalized weight \(w_k(n)=\pi_k\{\pi_k/x_k(n)\}^{\gamma}\), and let \(W(n)=\sum_{j=1}^K w_j(n)\). Then \(p_k(n+1)=w_k(n)/W(n)\). The ratios compare shares and are dimensionless. If a target share exceeds the current share, its ratio exceeds one; if it falls short, its ratio is below one. The power \(\gamma\) controls how strongly that ratio changes the weight, while the factor \(\pi_k\) retains the target share as the baseline weight.

Why These Weights Are Normalized This Way

The multi-arm extension needs both a correction and a baseline. The correction is the powered ratio \(\{\pi_k/x_k(n)\}^{\gamma}\); the baseline is \(\pi_k\). A natural construction is to multiply these arm-specific positive quantities to obtain weights, then divide each weight by their common total. The resulting probabilities are nonnegative and collectively describe one assignment.

$$ \sum_{k=1}^{K}p_k(n+1) = \sum_{k=1}^{K}\frac{w_k(n)}{W(n)} = \frac{\sum_{k=1}^{K}w_k(n)}{W(n)} =1. $$

This is the normalization check: it follows directly from using the same total weight in every denominator. Under the stated conditions, every \(w_k(n)>0\) and \(W(n)>0\), so each probability is positive as well. At \(\gamma=0\), the powered ratios equal one and the probabilities reduce to the target shares. If current shares equal target shares, every ratio is one, so the probabilities also equal the target shares for any \(\gamma\).

These checks are useful for distinguishing the rule from a ratio-only formula. Normalizing \(\{\pi_k/x_k(n)\}^{\gamma}\) without the baseline factor \(\pi_k\) would give equal probabilities when \(\gamma=0\), not the target allocation in general. That is a different rule. The DBCD weights here preserve the target allocation as the \(\gamma=0\) case.

Worked Example: Three Arms With Unequal Current Shares

Suppose 20 participants have already been assigned: 8 to arm 1, 8 to arm 2, and 4 to arm 3. The target shares are \((0.50,0.30,0.20)\), and \(\gamma=1\). All current counts and target shares are positive, and the current counts sum to 20, so the ratios are defined.

The current shares are \(x_1=8/20=0.40\), \(x_2=8/20=0.40\), and \(x_3=4/20=0.20\). The weights are

$$ \begin{aligned} w_1&=0.50(0.50/0.40)=0.625,\\ w_2&=0.30(0.30/0.40)=0.225,\\ w_3&=0.20(0.20/0.20)=0.200. \end{aligned} $$

Their total is \(W=0.625+0.225+0.200=1.050\). Therefore,

$$ (p_1,p_2,p_3) = (0.625/1.050,\;0.225/1.050,\;0.200/1.050) \approx(0.5952,\;0.2143,\;0.1905). $$

The probabilities sum to \(0.5952+0.2143+0.1905=1.0000\), allowing for rounding. The first arm is below its 0.50 target share and receives the largest next-assignment probability. These are probabilities for the next participant, not a guarantee of the realized assignment or of the eventual allocation.

Derivation and Checks in Practice

The normalization step can be justified without assuming that the raw weights already sum to one. Start with positive weights \(w_1,\ldots,w_K\). A probability vector proportional to these weights has the form \(p_k=cw_k\), where \(c\) is a common positive constant. Requiring the probabilities to sum to one gives \(c\sum_j w_j=1\), hence \(c=1/\sum_j w_j\). Substitution gives \(p_k=w_k/\sum_jw_j\). Thus the common denominator is required by the probability constraint; it is not an optional adjustment.

Before using the calculation, check the inputs and output. The current counts should sum to \(n\); the target shares should sum to one; and each share used in a denominator should be positive. After computing the weights, their total should be positive. Finally, calculate all probabilities with the same total weight and verify that they sum to one before rounding. A sum slightly different from one after rounding can be harmless; a material discrepancy usually signals an arithmetic or transcription error.

Worked Example: Strong Correction When One Arm Is Underrepresented

Suppose \(n=30\), with counts \((6,15,9)\), so the current shares are \((0.20,0.50,0.30)\). Let the target shares be \((0.50,0.30,0.20)\), and use \(\gamma=2\). The counts sum to 30, and all current and target shares are positive.

The three ratios are \(0.50/0.20=2.5\), \(0.30/0.50=0.6\), and \(0.20/0.30=2/3\). Squaring each ratio and multiplying by its target share gives

$$ \begin{aligned} w_1&=0.50(2.5)^2=3.125,\\ w_2&=0.30(0.6)^2=0.108,\\ w_3&=0.20(2/3)^2\approx0.08889. \end{aligned} $$

Thus \(W\approx3.32189\), and

$$ (p_1,p_2,p_3) \approx (3.125/3.32189,\;0.108/3.32189,\;0.08889/3.32189) \approx(0.9407,\;0.03251,\;0.02676). $$

The probabilities sum to approximately \(1.0000\). The result assigns a high probability to arm 1 because its current share is much lower than its target share and \(\gamma=2\) magnifies that ratio. The calculation is mathematically valid, but it also illustrates why the computed probabilities should be inspected: a valid normalized vector can be highly concentrated.

Worked Example: At the Target Shares

Suppose 100 participants have been assigned, with counts \((50,30,20)\), and the target shares are \((0.50,0.30,0.20)\). Take \(\gamma=3\). The current shares equal the target shares, so every target-to-current ratio is one. The weights are therefore

$$ (w_1,w_2,w_3) = (0.50(1)^3,\;0.30(1)^3,\;0.20(1)^3) =(0.50,0.30,0.20). $$

The total weight is one, and the next-assignment probabilities are \((0.50,0.30,0.20)\), which sum to one. This fixed-point check holds for any nonnegative \(\gamma\): when current shares match the target, the ratio correction leaves the target probabilities unchanged.

Failure Modes and Alternatives

The formula as stated requires positive current and target shares. If an arm has \(N_k(n)=0\), its current share is zero and \(\pi_k/x_k(n)\) is undefined. A design needs a suitable initialization or another explicitly specified boundary rule before applying this formula; silently substituting zero, an arbitrary small value, or an infinite weight does not follow from the rule. A zero target share also falls outside the positive-target condition used here and may require a separate design definition. These cases should be handled deliberately, not by pretending the ordinary calculation still applies.

The principal alternative at a given state is to randomize directly with probabilities equal to the target shares, \(p_k(n+1)=\pi_k\), without correcting for current allocation shares. That is simple and remains well-defined for positive target shares, but it does not use the current imbalance in this assignment step. The normalized DBCD rule incorporates that imbalance through the ratios. Neither calculation changes the target vector itself; the key distinction is whether current shares are used to adjust the next assignment probabilities.

Another common error is to compute weights and report them as probabilities. Weights such as \((0.625,0.225,0.200)\) are not probabilities unless they sum to one. Divide every weight by the same total. Also avoid rounding ratios or weights too early: small rounding errors across several arms can obscure whether normalization is correct.

Common mistake: “The target-to-current ratios are the allocation probabilities, so I can use them directly.” Ratios need not sum to one, and they omit the target-share baseline. Use \(w_k=\pi_k(\pi_k/x_k)^\gamma\) for each arm, then set \(p_k=w_k/\sum_jw_j\). Check the resulting probabilities sum to one.

Normalization guarantees a valid probability vector under the stated conditions; it does not guarantee that the allocation will equal the target at each step, nor does it remove uncertainty in the realized sequence of assignments. The formula is also sensitive to its inputs: small current shares can produce large ratios, especially when \(\gamma\) is large. Thus, verify the computed vector and its assumptions rather than treating successful normalization as evidence that every design choice is appropriate.

1
Check the state.
Confirm that the current counts sum to \(n\), the target shares sum to one, and the shares required in denominators are positive.
2
Compute each weight.
For every arm, calculate \(w_k=\pi_k(\pi_k/x_k(n))^\gamma\), using the same \(\gamma\).
3
Normalize and verify.
Divide every weight by \(W=\sum_jw_j\), then check that the probabilities are nonnegative and sum to one, allowing only for rounding.

Check Your Understanding

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

  1. With \(\gamma=0\), what are the normalized assignment probabilities if the target shares are \((0.6,0.3,0.1)\)?
  2. For target shares \((0.4,0.4,0.2)\), current shares \((0.2,0.5,0.3)\), and \(\gamma=1\), calculate the weights and probabilities.
  3. Why can’t the ordinary formula be evaluated when an arm has current share zero?
  4. A calculation gives positive weights summing to 2.4. What should be done before reporting assignment probabilities?

Solutions. (1) At \(\gamma=0\), each powered ratio equals one, so the weights are the target shares and sum to one. The probabilities are \((0.6,0.3,0.1)\). (2) The ratios are \(2\), \(0.8\), and \(2/3\). The weights are \(0.4(2)=0.8\), \(0.4(0.8)=0.32\), and \(0.2(2/3)\approx0.1333\). Their total is approximately \(1.2533\), giving probabilities approximately \((0.6383,0.2553,0.1064)\), which sum to one. (3) The target-to-current ratio divides by the current share, so a zero current share makes that ratio undefined. Initialization or a separately defined boundary rule is needed. (4) Divide every weight by 2.4. The normalized probabilities will sum to one, subject to rounding.

Clinical Biostats