Tutorials › Response-Adaptive Randomization › Martingale Difference Sequences in Adaptive Trials

Asymptotic Theory for Adaptive Allocation · Tutorial 193 of 501

Martingale Difference Sequences in Adaptive Trials

You will derive the martingale difference conditions for response residuals and allocation innovations, then check them in Bernoulli and bounded-outcome adaptive trials.

Advanced 12 min read

What You'll Learn

  • Define a martingale difference sequence relative to a trial filtration.
  • Center observed responses using their conditional, arm-specific means.
  • Center allocation indicators using their probabilities given the past.
  • Verify conditional-mean and integrability conditions in Bernoulli examples.
  • Check bounded residuals and weighted response innovations.
  • Recognize when pooled centering or conditioning on the wrong information fails.

Why Centering Matters in an Adaptive Trial

In an adaptive trial, the probability of assigning the next participant to an arm may depend on allocations and responses already observed. The next allocation is therefore not generally independent of the trial’s past. Yet two useful kinds of random variation still have a precise mean-zero property once that past is taken into account: the difference between an allocation indicator and its conditional probability, and the difference between a response and its conditional mean. These are martingale differences.

The point is not to claim that successive allocations or responses are independent. Instead, each increment has conditional expectation zero given the information available just before it occurs. That property makes martingale difference sequences a basic way to represent uncertainty under adaptive sampling. Here the focus is on establishing the property and checking its conditions, not on drawing asymptotic conclusions from it.

Trial Information and the Definition

Suppose participants enter in order \(i=1,2,\ldots\), and there are \(K\) treatment arms. Let \(\mathcal{F}_{i-1}\) denote the information available immediately before participant \(i\) is assigned, and let \(\mathcal{F}_i\) include that participant’s assignment and observed response as well as the earlier information. Thus the information sets grow over time: \(\mathcal{F}_{i-1}\subseteq\mathcal{F}_i\). This is the timing convention needed to say what is known before each new random increment.

Write \(A_i\) for participant \(i\)’s assigned arm and \(I(A_i=k)\) for the indicator that the assignment is to arm \(k\). The design’s conditional allocation probability is \(\pi_{i,k}=P(A_i=k\mid\mathcal{F}_{i-1})\). It is determined by the observed past before the assignment is drawn, although its value may change from one participant to the next.

Definition: A sequence \(D_1,D_2,\ldots\) is a martingale difference sequence with respect to \(\mathcal{F}_0,\mathcal{F}_1,\ldots\) if each \(D_i\) is measurable with respect to \(\mathcal{F}_i\), \(E(|D_i|)<\infty\), and \(E(D_i\mid\mathcal{F}_{i-1})=0\) with probability one. Measurability means the increment is known once the information at time \(i\) is available; integrability ensures its expectation is well-defined.

The condition is conditional, not merely unconditional. It implies \(E(D_i)=0\) by taking expectations, but the reverse implication is not enough: an increment can have overall mean zero while its conditional mean varies with the history.

Allocation Innovations

For each arm \(k\), define the allocation innovation \(D^A_{i,k}=I(A_i=k)-\pi_{i,k}\). The allocation rule makes \(\pi_{i,k}\) the conditional probability of the indicator being one, given the pre-assignment information. Therefore \(E(I(A_i=k)\mid\mathcal{F}_{i-1})=\pi_{i,k}\), and subtracting \(\pi_{i,k}\) gives a conditional mean of zero. The innovation is known by time \(i\), and its absolute value is at most one, so it is integrable.

For a fixed participant, these arm-specific innovations are linked: summing them over all arms gives zero, because exactly one arm is assigned and the allocation probabilities sum to one. They are not independent across arms. This linkage does not prevent each component from being a martingale difference.

Worked Example: Allocation Innovations in a Two-Arm Trial

Suppose that, at a particular history, the design assigns the next participant to arm 1 with probability \(0.35\) and to arm 2 with probability \(0.65\). Let \(D^A_{i,1}=I(A_i=1)-0.35\). The conditional distribution given \(\mathcal{F}_{i-1}\) is \(D^A_{i,1}=0.65\) with probability \(0.35\), and \(D^A_{i,1}=-0.35\) with probability \(0.65\). Hence

$$ E(D^A_{i,1}\mid\mathcal{F}_{i-1}) =0.35(0.65)+0.65(-0.35) =0. $$

If the realized assignment is arm 2, the observed innovation is \(-0.35\). That realized value is not zero; the martingale difference condition concerns its conditional expectation before the draw, not the value of every increment. The innovation is bounded in absolute value by one, verifying integrability.

The conditional variance can also be checked directly: \(0.35(0.65)^2+0.65(0.35)^2=0.2275\). This is a property of the conditional distribution at this history; it does not change the zero conditional mean.

Response Residuals Centered at Conditional Means

Let \(Y_i\) be the response observed for participant \(i\). For each arm \(k\), suppose the conditional response mean is \(m_{i,k}=E(Y_i\mid\mathcal{F}_{i-1},A_i=k)\), whenever arm \(k\) has positive conditional assignment probability. This formulation allows the mean to depend on the information already observed. A common special case is a fixed arm mean, \(m_{i,k}=\mu_k\).

The observed response residual is \(D^Y_i=Y_i-m_{i,A_i}\). If response distributions are integrable and the conditional means are finite, then the residual is integrable. To verify its conditional mean, condition first on the possible assignments. For every arm that can be assigned, its response residual has conditional mean zero by the definition of \(m_{i,k}\). Averaging those zero means over the assignment probabilities also gives zero.

$$ \begin{aligned} E(D^Y_i\mid\mathcal{F}_{i-1}) &=\sum_{k=1}^K \pi_{i,k} E(Y_i-m_{i,k}\mid\mathcal{F}_{i-1},A_i=k)\\ &=\sum_{k=1}^K \pi_{i,k}(m_{i,k}-m_{i,k})\\ &=0. \end{aligned} $$

The same argument gives a useful arm-specific version: \(D^Y_{i,k}=I(A_i=k)(Y_i-m_{i,k})\). Its conditional expectation is the assignment probability multiplied by the within-arm conditional residual mean, so it is zero as well. If the response is integrable and the conditional mean is finite, the weighted residual is integrable. These statements require conditional centering at the appropriate arm-specific mean; they do not require identical responses across participants.

Condition to check: Before calling a response residual a martingale difference, specify its conditioning information and its centering mean. The response model must give \(E(Y_i\mid\mathcal{F}_{i-1},A_i=k)=m_{i,k}\), and the residual must have finite absolute expectation. For bounded responses these integrability conditions hold automatically when the means are finite.

Worked Examples: Response Conditions

Worked Example: Bernoulli Responses on Two Arms

Consider a hypothetical two-arm trial. At the current history, the assignment probabilities are \(0.60\) for arm 1 and \(0.40\) for arm 2. Conditional on assignment to arm 1, let the binary response be Bernoulli with success probability \(0.70\); on arm 2, let it be Bernoulli with success probability \(0.40\). Assume these response probabilities hold conditional on \(\mathcal{F}_{i-1}\) and the assigned arm. The conditional means are therefore \(m_{i,1}=0.70\) and \(m_{i,2}=0.40\).

Center the observed response by its assigned-arm mean: \(D^Y_i=Y_i-m_{i,A_i}\). Given assignment to arm 1, its mean is \(0.70-0.70=0\); given assignment to arm 2, its mean is \(0.40-0.40=0\). Averaging over assignment gives

$$ E(D^Y_i\mid\mathcal{F}_{i-1}) =0.60(0.70-0.70)+0.40(0.40-0.40)=0. $$

A Bernoulli response is either zero or one, and each centering mean lies between zero and one. Thus \(|D^Y_i|\leq 1\), so \(E(|D^Y_i|)<\infty\). If the observed assignment is arm 1 and the response is zero, the residual is \(-0.70\); if it is one, the residual is \(0.30\). Neither realized residual has to equal zero for the sequence to satisfy the martingale difference condition.

Worked Example: Bounded Responses with Conditional Means

In a hypothetical trial, suppose responses are measured on a scale from 0 to 10. Given the current history, assignment to arm 1 produces a response of 4 or 8 with equal probabilities, while assignment to arm 2 produces a response of 1 or 5 with equal probabilities. Let the assignment probabilities be \(0.25\) and \(0.75\), respectively. The conditional arm means are \(m_{i,1}=(4+8)/2=6\) and \(m_{i,2}=(1+5)/2=3\).

For the observed-response residual, the possible centered values are \(-2\) and \(2\) on either arm. Consequently, its conditional mean is

$$ E(D^Y_i\mid\mathcal{F}_{i-1}) =0.25\left[\tfrac12(4-6)+\tfrac12(8-6)\right] +0.75\left[\tfrac12(1-3)+\tfrac12(5-3)\right] =0. $$

The arithmetic verifies centering within each arm before averaging over assignments. Since the residual is always either \(-2\) or \(2\), its absolute value is at most 2 and it is integrable. For the arm-1 weighted residual, the conditional expectation is \(0.25[\tfrac12(-2)+\tfrac12(2)]=0\); it equals zero when arm 2 is assigned and has a finite absolute expectation. Thus both the observed residual and the arm-specific weighted residual meet the conditional-mean and integrability requirements.

Why the Conditional-Mean Calculation Works

The response calculation uses the law of iterated expectations: condition on the assignment, whose probabilities are known given \(\mathcal{F}_{i-1}\), and then average over those assignments. This short derivation is important because adaptive allocation changes the weights \(\pi_{i,k}\) as the trial proceeds. It does not invalidate the calculation: at each step, the probabilities used are the probabilities conditional on the actual past.

1
Specify the information available.
Use \(\mathcal{F}_{i-1}\) for the history before assignment and \(\mathcal{F}_i\) for the history after observing assignment and response.
2
State the conditional model.
Give \(\pi_{i,k}=P(A_i=k\mid\mathcal{F}_{i-1})\) and \(m_{i,k}=E(Y_i\mid\mathcal{F}_{i-1},A_i=k)\).
3
Center each increment.
Subtract \(\pi_{i,k}\) from the allocation indicator, or \(m_{i,A_i}\) from the observed response.
4
Check both requirements.
Calculate the conditional expectation given \(\mathcal{F}_{i-1}\), then establish finite absolute expectation.

This argument needs no independence between increments. It does require that the conditional means and assignment probabilities match the specified sequential model. For example, if the response mean depends on earlier observed history, a fixed arm mean is a valid centering value only if it still equals the conditional mean at each relevant history.

Comparison, Diagnostics, and Common Pitfalls

SequenceWhat is subtracted?Key check
Allocation innovationConditional assignment probability \(\pi_{i,k}\)\(E[I(A_i=k)\mid\mathcal{F}_{i-1}]=\pi_{i,k}\)
Observed-response residualAssigned-arm conditional mean \(m_{i,A_i}\)\(E[Y_i\mid\mathcal{F}_{i-1},A_i=k]=m_{i,k}\)
Arm-specific weighted residualArm-specific conditional mean, multiplied by \(I(A_i=k)\)Within-arm residual mean is zero and the increment is integrable

An ordinary independent, mean-zero sequence is one special case of a martingale difference sequence relative to a suitable past-information filtration. The martingale condition is more flexible: it allows the distribution of an increment to depend on the past, provided its conditional mean is zero. Conversely, an unconditional mean-zero calculation alone is weaker and does not establish the martingale difference property.

A frequent mistake is to center every response by one pooled mean. In the Bernoulli example above, the mean of the response given the history is \(0.60(0.70)+0.40(0.40)=0.58\). If the assignment probability changes to \(0.80\) for arm 1, the conditional response mean becomes \(0.80(0.70)+0.20(0.40)=0.64\). Subtracting the old pooled value \(0.58\) would then give conditional mean \(0.06\), not zero. Subtracting the assigned-arm mean instead gives zero within either arm, regardless of the current allocation probabilities.

Other useful diagnostics are direct checks of the timing and model assumptions. Is \(\pi_{i,k}\) calculated from information available before assignment, rather than from the observed assignment itself? Does the stated response mean condition on both the relevant history and the assigned arm? Is the response integrable? For bounded outcomes, integrability follows from the bound; for unbounded outcomes, it must be checked from the response model. A wrong filtration, a misspecified conditional mean, or a response distribution with no finite absolute expectation can invalidate the claimed sequence.

The martingale difference property does not assert that residuals are independent, that their variances are constant, or that an arm will be sampled often. Nor does it by itself establish a law of large numbers or a central limit theorem. Those require additional conditions. The practical value here is more basic: it separates variation that is unpredictable on average given the past from predictable changes in allocation or response means.

Key takeaway: Allocation indicators become martingale differences when centered by their conditional assignment probabilities. Observed responses become martingale differences when centered by their conditional means given the past and assigned arm. In both cases, verify the conditional expectation and integrability rather than relying on an unconditional mean-zero statement.

Check Your Understanding

For each question, identify the relevant conditioning information and verify the claim. Solutions follow the questions.

  1. At a history, an arm has assignment probability \(0.20\). What are the two possible values of its allocation innovation, and what is its conditional mean?
  2. A binary response has success probability \(0.30\) on the assigned arm. Give the two possible centered response residuals and explain why they are integrable.
  3. In a two-arm trial, conditional arm means are 2 and 5, and assignment probabilities are \(0.40\) and \(0.60\). Verify that the observed-response residual has conditional mean zero.
  4. Why does \(E(D_i)=0\) alone not prove that \(D_i\) is a martingale difference?
  5. Responses are bounded between 0 and 10, but the analyst centers them by a constant that differs from their conditional mean. Which condition may fail, and which condition is still assured by boundedness?

Solutions. (1) The innovation is \(1-0.20=0.80\) with probability \(0.20\), and \(0-0.20=-0.20\) with probability \(0.80\). Its conditional mean is \(0.20(0.80)+0.80(-0.20)=0\). (2) The residual is \(1-0.30=0.70\) with probability \(0.30\), or \(0-0.30=-0.30\) with probability \(0.70\). It is bounded in absolute value by 1, so it is integrable. (3) The conditional mean of the residual is \(0.40(2-2)+0.60(5-5)=0\), conditioning on the pre-assignment history and averaging over assignments. (4) Unconditional mean zero does not guarantee \(E(D_i\mid\mathcal{F}_{i-1})=0\); the conditional mean could vary with the past. (5) The conditional-mean-zero condition may fail because the centering constant is not the response’s conditional mean. Integrability is assured: bounded responses and a finite constant centering value give a bounded residual.

Clinical Biostats