Tutorials › Response-Adaptive Randomization › Sequentially Updating Response Estimates

Doubly Adaptive Biased Coin Designs · Tutorial 150 of 501

Sequentially Updating Response Estimates

Follow the full sequential update from an observed response to revised estimates, target shares, and the next assignment probability, while keeping estimates and allocation counts distinct.

Advanced 11 min read

What You'll Learn

  • Distinguish the assignment counts used by the DBCD correction from the observations used to update response estimates.
  • Update a running treatment-specific mean without recomputing it from all prior responses.
  • Recalculate target shares and next-patient probabilities after each response is observed.
  • Identify which information must be available before an assignment probability is calculated.
  • Diagnose errors caused by confusing target shares with assignment probabilities or updating on pending outcomes.

From One Observed Response to the Next Assignment

A response-adaptive design changes its allocation as information accumulates. The operational unit is therefore not just an estimate or an allocation rule, but a sequence: record an assignment, observe its response, update the treatment-specific estimates, recalculate the target shares, and then calculate the probability for the next assignment. Keeping those steps in the correct order prevents information from a future response from influencing a probability that should already have been set.

This tutorial works with two treatment arms. The treatment-specific response estimator and the target map must be specified by the design; the update mechanics do not determine which response measure is scientifically appropriate. The examples use running means of numerical response scores and a simple normalized-estimate target solely to make each calculation explicit. They are hypothetical illustrations, not a recommendation to use that target map in a trial.

Definition: Immediately before the \((n+1)\)st assignment, let \(N_i(n)\) be the number assigned to arm \(i\) among the first \(n\) assignments, and let \(\hat{\theta}_i(n)\) be the current estimate of that arm’s response parameter, for \(i=1,2\). A prespecified target map converts the estimates into target shares \(\rho_1(n)\) and \(\rho_2(n)\), where \(\rho_1(n)+\rho_2(n)=1\). The DBCD assignment probability \(q_i(n+1)\) is the probability of assigning the next participant to arm \(i\), conditional on the information available before that assignment.

For two arms, the target shares and next-patient probabilities have different roles. The target share \(\rho_i(n)\) is the allocation proportion the design is currently tracking. The probability \(q_i(n+1)\) is the chance of the next assignment, adjusted using the current allocation counts. With tuning parameter \(\gamma\geq 0\), the two-arm DBCD rule is

$$ q_i(n+1)= \frac{\rho_i(n)\left\{\rho_i(n)/N_i(n)\right\}^{\gamma}} {\displaystyle\sum_{j=1}^{2}\rho_j(n)\left\{\rho_j(n)/N_j(n)\right\}^{\gamma}}, \qquad i=1,2. $$

The rule requires positive allocation counts and target shares so its ratios and weights are defined. Initialization addresses the early stage before estimates and counts are usable. At each later step, the calculation uses the current state; it does not assume the estimates or target shares remain fixed after a response.

The Update Order and the Information Set

The probability \(q_i(n+1)\) is calculated before the \((n+1)\)st assignment. Let \(A_{n+1}\) denote the arm assigned and \(Y_{n+1}\) its response once observed. The assignment updates the allocation count immediately: \(N_{A_{n+1}}(n+1)=N_{A_{n+1}}(n)+1\). The response estimator is updated when the response becomes available. Then the target map is applied to the updated estimates, and the DBCD formula gives the probability for assignment \(n+2\).

This order is a form of sequential bookkeeping. A probability must depend only on information available before its assignment. Thus, the response to assignment \(n+1\) cannot be used in \(q_i(n+1)\), but it can be used in \(q_i(n+2)\) if observed in time. If outcomes are delayed, the allocation count may already include a participant whose response has not yet updated the estimate. The design must state how such pending outcomes are handled; a response that has not been observed is not silently treated as a known response.

1
Start from the current state.
Record each arm’s allocation count and the response-estimator information available before the next assignment.
2
Calculate the next assignment probability.
Apply the current target shares and counts in the DBCD rule, then randomize the next participant using those probabilities.
3
Update after assignment and response.
Increase the assigned arm’s count. When its response is observed, update that arm’s estimator using the specified estimator rule.
4
Recompute for the following assignment.
Apply the target map to the updated estimates, then recalculate the DBCD probabilities using the updated counts.

Updating a Running Mean

One common estimator update for a numerical response is the running sample mean. If arm \(i\) has \(m_i\) observed responses with mean \(\bar{x}_i\), and its next observed response is \(y\), the updated mean is

$$ \bar{x}_{i,\mathrm{new}} =\frac{m_i\bar{x}_i+y}{m_i+1} =\bar{x}_i+\frac{y-\bar{x}_i}{m_i+1}. $$

The first expression adds the new response to the previous total and divides by the new number of observed responses. The second shows the direction and size of the change: the estimate moves toward \(y\), with a smaller change when more responses have already contributed to the mean. This update is valid for a running mean of the observed numerical responses; other estimators may require different stored information and update equations.

Worked Example: Updating Estimates and the Next Probability

Suppose eight hypothetical participants have been assigned, with four in each arm. The observed response scores have means \(\bar{x}_1=6.0\) points and \(\bar{x}_2=8.0\) points. For illustration, let the target shares be the estimates normalized to sum to one, and take \(\gamma=1\). Calculate the next assignment probability, then update after arm 2 receives a response of 10 points.

Before the next assignment. The target shares are \(\rho_1=6/(6+8)=0.4286\) and \(\rho_2=8/(6+8)=0.5714\). With \(N_1=N_2=4\), the unnormalized DBCD weights are

$$ w_1=\frac{\rho_1^2}{N_1} =\frac{0.4286^2}{4}=0.04592, \qquad w_2=\frac{\rho_2^2}{N_2} =\frac{0.5714^2}{4}=0.08163. $$

Therefore \(q_1=0.04592/(0.04592+0.08163)=0.3600\) and \(q_2=0.6400\). The probability of assigning the next participant to arm 1 is 0.3600; it is not the same as arm 1’s target share of 0.4286.

Suppose the next participant is assigned to arm 2 and has a response of 10 points. Arm 2’s updated mean is \((4(8.0)+10)/5=42/5=8.4\) points. The allocation counts are now \(N_1=4\) and \(N_2=5\). The new target shares are \(\rho_1=6.0/(6.0+8.4)=0.4167\) and \(\rho_2=0.5833\). For the following assignment,

$$ w_1=\frac{0.4167^2}{4}=0.04340, \qquad w_2=\frac{0.5833^2}{5}=0.06806, \qquad q_1=\frac{0.04340}{0.04340+0.06806}=0.3894. $$

Thus, the next probability for arm 1 rises from 0.3600 to 0.3894. Both the target share and the allocation-count correction contribute to the change.

A Second Sequential Update

The calculation repeats after every observed response. Counts change because of assignments; estimates change because of observed outcomes; target shares are recalculated from those estimates. These state changes need not move the target share and the next-patient probability in the same direction.

Worked Example: Recalculating After an Arm 1 Response

In a separate hypothetical state, suppose \(N_1=N_2=5\), with running means of 7 points for arm 1 and 5 points for arm 2. Use the same illustrative normalized-estimate target and \(\gamma=1\). Find the probability for the next assignment. Then suppose arm 1 is assigned and its response is 4 points; calculate the probability for the assignment after that response.

Initially, \(\rho_1=7/12=0.5833\) and \(\rho_2=5/12=0.4167\). The weights are \(w_1=0.5833^2/5=0.06806\) and \(w_2=0.4167^2/5=0.03472\). Hence \(q_1=0.06806/(0.06806+0.03472)=0.6622\).

After assignment to arm 1 and observation of 4 points, its mean becomes \((5(7)+4)/6=39/6=6.5\) points. The counts are now \(N_1=6\) and \(N_2=5\), and the revised shares are \(\rho_1=6.5/(6.5+5)=0.5652\) and \(\rho_2=0.4348\). The next weights are

$$ w_1=\frac{0.5652^2}{6}=0.05324, \qquad w_2=\frac{0.4348^2}{5}=0.03781. $$

Therefore the probability of arm 1 for the next assignment is \(q_1=0.05324/(0.05324+0.03781)=0.5848\). The target share for arm 1 declined from 0.5833 to 0.5652, while its next-patient probability declined from 0.6622 to 0.5848. The initial probability was above the target share because arm 1’s allocation count was not yet in the target proportion.

Pending Responses and Sequential Validity

When responses are not immediate, the design needs a precise rule for distinguishing participants assigned from responses observed. For example, suppose an arm has five observed responses with mean 6 points and a sixth participant has been assigned to that arm but has not yet responded. The count used to describe assignments is six; the running mean still summarizes five observed responses. When the sixth response arrives, the estimator can be updated using all six observed values.

A response-adaptive procedure should specify whether it waits for an outcome before making another assignment, permits multiple pending outcomes, or uses a prespecified estimator that accommodates incomplete follow-up. The relevant requirement for the sequential probability is that it be determined from the information actually available at that point. Updating an estimate with an unobserved value, or retrospectively changing a probability using a response that arrived later, breaks that ordering.

Worked Example: A Response Arrives After Its Assignment

Suppose five responses in arm 1 have mean 6.0 points, four participants have been assigned to arm 2, and the current estimates are 6.0 and 8.0 points. A sixth participant is assigned to arm 1, but that participant’s response is pending. Under the illustrative normalized target and with \(\gamma=1\), determine the target and next probability before the pending response arrives. Then update when the response is 9 points.

Before observing the pending response, the estimates remain 6.0 and 8.0 points, so the target shares remain \(3/7=0.4286\) and \(4/7=0.5714\). The assignment counts are \(N_1=6\) and \(N_2=4\). Thus \(w_1=0.4286^2/6=0.03061\), while \(w_2=0.5714^2/4=0.08163\), giving \(q_1=0.03061/(0.03061+0.08163)=0.2727\).

When the pending response of 9 points is observed, arm 1’s running mean becomes \((5(6)+9)/6=39/6=6.5\) points. The allocation counts do not change at the time the response arrives: they were already \(N_1=6\) and \(N_2=4\). The revised target shares are \(6.5/14.5=0.4483\) and \(8/14.5=0.5517\). The new weights are \(0.4483^2/6=0.03349\) and \(0.5517^2/4=0.07610\), so the following arm 1 probability is \(0.03349/(0.03349+0.07610)=0.3056\).

The estimate and target changed when the response was observed; the allocation count did not. This distinction is essential whenever assignment and outcome observation occur at different times.

Why the Update Rule Is Coherent

The DBCD formula produces valid probabilities at each step when its inputs are valid. To see this, define \(w_i(n)=\rho_i(n)\{\rho_i(n)/N_i(n)\}^{\gamma}\). If \(N_i(n)>0\), \(\rho_i(n)>0\), and \(\gamma\geq0\), then \(w_i(n)>0\). Dividing each weight by the sum of the two positive weights gives positive probabilities whose sum is one:

$$ q_1(n+1)+q_2(n+1) =\frac{w_1(n)+w_2(n)}{w_1(n)+w_2(n)}=1. $$

The update sequence also preserves the intended timing: the current state determines the next randomization, and newly observed data can affect later randomizations. This is not a claim that the point estimates are unbiased, stable, or scientifically adequate. Those properties depend on the estimator, the response process, and the target map. The mechanics ensure only that the specified estimator and target are carried forward in the correct sequence.

Diagnostics, Alternatives, and Limitations

A practical implementation should retain an audit trail for each assignment: the counts, observed-response information, estimates, target shares, tuning parameter, computed probabilities, randomization, and later outcome update. Recomputing a few steps independently is a useful check. At every step, verify that the target shares sum to one, the assignment probabilities sum to one, and the probability used was calculated before the corresponding assignment.

  • Using a target share as the randomization probability. The target is not generally the next-patient probability. The count correction in the DBCD rule changes the probability when current allocations differ from the target.
  • Updating an estimate at assignment rather than observation. An assignment supplies no response value. Keep the estimate unchanged until the response is observed or use a clearly specified method for pending data.
  • Mixing counts and estimator sample sizes. Allocation counts and numbers of observed responses may differ with delayed outcomes. Use each quantity for its intended role.
  • Changing the target map midstream without specification. The response estimates alone do not define target shares. A map must be prespecified or any adaptation to it must itself be explicitly defined.

A fixed allocation probability does not respond to the estimates. Directly randomizing according to the target shares also differs from DBCD because it omits the correction for current allocation imbalance. Batch updating can be appropriate in designs that deliberately randomize groups of participants, but it is not the same sequence as recalculating after each available response. Each alternative changes the design’s allocation behavior and should not be substituted silently.

Finally, an online update is not a guarantee of good performance. Early estimates can be unstable; an extreme estimate can shift a target sharply if the target map is sensitive; and delayed or missing responses can make the current estimates less current than the assignment counts. The method described here specifies how to propagate estimates through a DBCD rule. It does not by itself solve estimator selection, uncertainty quantification, missing outcomes, or the scientific choice of target.

Check Your Understanding

Work through the questions, then compare with the solutions.

  1. With \(N_1=N_2=3\), estimates 4 and 6, the normalized-estimate target, and \(\gamma=1\), what is the next probability of arm 1?
  2. An arm has four observed responses with mean 5 points. A new observed response is 9 points. What is the updated running mean?
  3. A participant has been assigned, but the response is pending. Should the estimate be updated using that response? What may already have changed?
  4. Why can a DBCD next-patient probability differ from the corresponding target share?

Solutions. (1) The target shares are \(4/10=0.4\) and \(6/10=0.6\). The weights are \(0.4^2/3=0.05333\) and \(0.6^2/3=0.1200\), so \(q_1=0.05333/(0.05333+0.1200)=0.3077\). (2) The updated mean is \((4(5)+9)/5=29/5=5.8\) points. (3) No. The assignment count has increased, but the response estimate remains based on observed responses until the pending value is available. (4) The probability includes the ratio correction using current allocation counts, whereas the target share describes the allocation proportion being tracked.

Clinical Biostats