Tutorials › Response-Adaptive Randomization › Why Plug-In Greedy Allocation Is Not Thompson Sampling

Thompson Sampling for Trials · Tutorial 258 of 501

Why Plug-In Greedy Allocation Is Not Thompson Sampling

Learn why choosing the arm with the largest posterior mean is not the same as sampling an arm according to its posterior probability of being best.

Advanced 11 min read

What You'll Learn

  • Define plug-in greedy allocation and Thompson sampling for binary response rates.
  • Calculate posterior best-arm probabilities from simple Beta posteriors.
  • Explain why posterior means alone do not determine Thompson sampling allocation.
  • Compare deterministic greedy decisions with uncertainty-sensitive sampling.
  • Identify assumptions and common implementation errors in the comparison.

Two Rules That Can Look Similar

Suppose a trial compares treatments with binary outcomes, such as response versus no response. After observing some outcomes, each treatment has a posterior distribution for its response probability. Two allocation rules may both use this posterior information, but they use it in different ways. Plug-in greedy allocation compares the posterior means and chooses the treatment with the largest mean. Thompson sampling accounts for the full posterior distributions: it draws a response probability for each treatment and allocates to the treatment whose draw is largest.

The distinction is not merely computational. A posterior mean summarizes the centre of a distribution, while a posterior probability of being best also reflects its uncertainty and shape. Plug-in greedy allocation turns a summary into a single choice. Thompson sampling randomizes in proportion to posterior uncertainty about which treatment is best. The allocation probabilities need not be equal to one another, and they are not generally determined by ranking the posterior means alone.

Definition: For treatment arm \(i\), let \(\theta_i\) be its binary response probability and let the current data and prior imply a posterior distribution for \(\theta_i\). Plug-in greedy allocation assigns to an arm with the largest posterior mean \(E(\theta_i\mid\text{data})\), using a specified tie rule if needed. Thompson sampling draws one value from each arm's posterior distribution and assigns to the arm with the largest draw. Thus, at the current posterior state, the probability that arm \(i\) is selected by Thompson sampling is the posterior probability that \(\theta_i\) is the largest response probability.

Model and Assumptions

For the examples, assume there are two arms, \(A\) and \(B\), and a patient's outcome on arm \(i\) is \(1\) for response and \(0\) for no response. Conditional on the arm's response probability \(\theta_i\), outcomes assigned to that arm are Bernoulli with success probability \(\theta_i\). Use independent Beta priors for the arms, and assume response probabilities remain stable during the period considered. Under these assumptions, if arm \(i\) has \(s_i\) responses and \(f_i\) nonresponses, a Beta prior with parameters \(a_i,b_i\) gives posterior \(\operatorname{Beta}(a_i+s_i,b_i+f_i)\). Its mean is \((a_i+s_i)/(a_i+b_i+s_i+f_i)\).

The posterior distributions describe uncertainty about the response probabilities, not variability in an individual patient's outcome. A posterior draw of \(\theta_i\) is therefore not a predicted response for the next patient. It is a draw representing uncertainty about the arm's underlying response rate. Thompson sampling compares these draws; plug-in greedy allocation compares posterior means.

The simple comparison assumes that a larger response probability is preferable, that the arms are evaluated on the same binary endpoint, and that the posterior distributions represent the intended model for the trial. If the clinical objective is not to maximize response probability, the relevant quantity for “best” would need to reflect that objective instead. These examples use independent posteriors so the probability calculations can be shown directly.

How the Two Rules Differ

Write the posterior means as \(\mu_A=E(\theta_A\mid\text{data})\) and \(\mu_B=E(\theta_B\mid\text{data})\). If \(\mu_A>\mu_B\), plug-in greedy allocation assigns to \(A\) (apart from any additional design restrictions). Thompson sampling instead selects \(A\) with probability \(P(\theta_A>\theta_B\mid\text{data})\), and selects \(B\) with probability \(P(\theta_B>\theta_A\mid\text{data})\). With continuous Beta posteriors, the probability of an exact tie is zero, so these two probabilities sum to one.

The central calculation for two arms is an average over the posterior distribution of one arm:

$$ P(\theta_A>\theta_B\mid\text{data}) =\int_0^1 f_A(x)\,F_B(x)\,dx, $$

where \(f_A\) is the posterior density of \(\theta_A\) and \(F_B(x)=P(\theta_B\leq x\mid\text{data})\) is the posterior cumulative distribution function for arm \(B\). To justify the expression, condition on \(\theta_A=x\). Given that value, the probability that \(B\)'s response probability is below it is \(F_B(x)\). Averaging that conditional probability over the posterior density \(f_A(x)\) yields the integral. This is a probability of being best, not a comparison of means.

Worked Example: A Larger Mean Does Not Mean Certain Selection

Suppose each arm begins with a \(\operatorname{Beta}(1,1)\) prior. In hypothetical trial data, \(A\) has one response and no nonresponses, while \(B\) has no responses and one nonresponse. The posterior distributions are \(A:\operatorname{Beta}(2,1)\) and \(B:\operatorname{Beta}(1,2)\).

The posterior means are \(E(\theta_A\mid\text{data})=2/(2+1)=2/3\) and \(E(\theta_B\mid\text{data})=1/(1+2)=1/3\). Plug-in greedy allocation therefore assigns to \(A\) with probability one, if it is implemented as a deterministic rule.

For Thompson sampling, the density of \(\theta_A\) is \(f_A(x)=2x\), and the cumulative distribution function for \(\theta_B\) is \(F_B(x)=1-(1-x)^2=2x-x^2\), for \(0\leq x\leq1\). Thus:

$$ P(\theta_A>\theta_B\mid\text{data}) =\int_0^1 2x(2x-x^2)\,dx =\left[\frac{4x^3}{3}-\frac{x^4}{2}\right]_0^1 =\frac{5}{6}\approx0.8333. $$

So Thompson sampling selects \(A\) with probability about \(0.8333\) and \(B\) with probability \(1-0.8333=0.1667\). Both rules favour \(A\), but they do not make the same allocation: the greedy rule always chooses \(A\), whereas Thompson sampling sometimes chooses \(B\) because its posterior still assigns \(B\) a chance of having the larger response rate.

Posterior Means and Best-Arm Probabilities

The calculation above also shows why replacing a posterior distribution with its mean changes the rule. A posterior mean is a single number for each arm. Once those numbers are compared, deterministic plug-in greedy allocation discards the rest of the distribution. In contrast, a best-arm probability depends on how much posterior mass each arm places above the other arm's possible values.

The posterior mean and the posterior probability of being best answer different questions. The mean describes the average response probability under the posterior. The best-arm probability describes how often the arm would have the highest response probability across draws from the joint posterior. A treatment can have a high posterior mean without having a correspondingly high probability of being best, particularly when its posterior is uncertain or asymmetric.

Worked Example: Identical Posteriors and a Greedy Tie

Suppose hypothetical data and priors lead to the same posterior for both arms: \(\theta_A,\theta_B\sim\operatorname{Beta}(2,2)\), independently. Each posterior mean is \(2/(2+2)=0.5\). A plug-in greedy rule now needs a tie-breaking convention. If its convention is “choose \(A\),” it assigns to \(A\) every time; if it breaks ties uniformly at random, it assigns to each arm with probability \(0.5\).

For Thompson sampling, the two independent posterior draws have identical continuous distributions. By symmetry, the event \(\theta_A>\theta_B\) and the event \(\theta_B>\theta_A\) have equal probabilities. Since an exact tie has probability zero, the two probabilities sum to one, so each is \(0.5\). Thus, with a uniform random tie rule, the two procedures happen to have the same marginal allocation probabilities in this special case. That numerical agreement does not make the rules equivalent: the greedy rule's tie convention is not generally determined by posterior uncertainty.

A Small Posterior Update Can Change the Sampling Probability

Greedy allocation can be insensitive to new evidence about how far ahead the leading arm is, as long as its posterior mean remains larger. Thompson sampling can respond to such evidence through a change in the probability that the leading arm is best. The next calculation illustrates this using the same pair of arms as before.

Worked Example: One More Response Changes the Best-Arm Probability

Start with independent \(\operatorname{Beta}(1,1)\) priors. Suppose, hypothetically, that arm \(A\) has two responses and no nonresponses, while \(B\) has no responses and one nonresponse. The posteriors are \(A:\operatorname{Beta}(3,1)\) and \(B:\operatorname{Beta}(1,2)\). Their means are \(3/4=0.75\) and \(1/3\), respectively, so plug-in greedy allocation selects \(A\).

The density for \(\theta_A\) is \(3x^2\), and \(F_B(x)=2x-x^2\) as before. Therefore:

$$ P(\theta_A>\theta_B\mid\text{data}) =\int_0^1 3x^2(2x-x^2)\,dx =\left[\frac{3x^4}{2}-\frac{3x^5}{5}\right]_0^1 =\frac{3}{2}-\frac{3}{5} =0.9000. $$

Thompson sampling selects \(A\) with probability \(0.9000\) and \(B\) with probability \(0.1000\). Compared with the preceding example, the additional response on \(A\) raised its posterior probability of being best from \(5/6\) to \(0.9000\). The greedy rule still selects \(A\) with probability one in both cases. The numerical allocation probabilities in this example describe the current posterior comparison; they are not a claim that future outcomes or future posterior states are fixed.

Why the Distinction Matters

Thompson sampling does not allocate in proportion to posterior means. It allocates according to posterior best-arm probabilities. This gives some allocation to an arm with a lower posterior mean whenever that arm has a nonzero chance of being best, while favouring arms that are more likely to be best. A deterministic greedy rule, by contrast, can stop exploring an arm as soon as its posterior mean falls below another's, even if substantial uncertainty remains.

This distinction is useful when checking a design implementation. A program that calculates posterior means and always assigns to the largest one implements plug-in greedy allocation, not Thompson sampling. A Thompson sampling implementation must draw from the posterior distributions (or otherwise calculate the corresponding best-arm probabilities) and let the resulting posterior comparison determine the selected arm. The draw-based procedure approximates the best-arm probabilities over repeated draws; finite Monte Carlo computation can add numerical error, but that error is separate from the conceptual difference between the two rules.

FeaturePlug-in greedy allocationThompson sampling
Posterior information usedPosterior mean for each armFull posterior distributions
Selection ruleChoose an arm with the largest meanDraw arm values and choose the largest draw
Selection probability for a leading-mean armOne under deterministic tie-free selectionIts posterior probability of being best
Role of uncertaintyDiscarded after calculating the meansReflected in the probability of being best

Assumptions, Diagnostics, and Limitations

The comparison is meaningful only if the posterior distributions are appropriate for the trial question. For the Beta-Bernoulli calculations here, outcomes must be represented as binary responses, the likelihood must be Bernoulli conditional on the arm-specific response probability, and the Beta update must use the correct counts. Independent posteriors are an assumption of the worked calculations, not a requirement for the general distinction: with dependent posteriors, Thompson sampling still compares draws from the joint posterior, but the best-arm probabilities must respect that dependence.

A useful diagnostic is to inspect the procedure's actual decision rule, not just its reported posterior means. If allocation always goes to the largest mean, then the procedure is greedy, even if it was described as “Bayesian.” If the implementation uses posterior draws, check that it draws from each arm's current posterior and compares those draws, rather than comparing the draws' averages or substituting the posterior means afterward. Also check that the binary outcome coding and success/failure counts match the model.

Neither rule is automatically appropriate for every clinical objective. Plug-in greedy allocation may under-explore arms with uncertain response rates. Thompson sampling can allocate to an arm that is not best in a particular posterior draw, and it can be sensitive to the model and prior used to form the posterior. The examples assume stable response probabilities and do not establish that either design meets any particular trial's ethical, inferential, or operational requirements. Those judgments require considering the trial's goals and assumptions, not just the allocation formula.

Common Mistake

Common mistake: “Arm \(A\) has the larger posterior mean, so Thompson sampling assigns every patient to \(A\).” This is plug-in greedy reasoning. Thompson sampling assigns to \(A\) according to \(P(\theta_A>\theta_B\mid\text{data})\), which can be less than one even when \(A\)'s posterior mean is larger. The correction is to distinguish the posterior mean from the posterior probability of being best, then apply the allocation rule that the design actually specifies.

Check Your Understanding

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

  1. Arm \(A\) has posterior mean \(0.68\), and arm \(B\) has posterior mean \(0.52\). What does a deterministic plug-in greedy rule do? What additional quantity determines Thompson sampling's probability of choosing \(A\)?
  2. In the first worked example, what is the posterior probability that \(B\) is best?
  3. Why can two identical continuous posterior distributions give each arm a Thompson sampling best-arm probability of \(0.5\)?
  4. A programmer draws posterior values for the arms but then chooses the arm with the largest posterior mean. What allocation rule has the programmer implemented?

Solutions. (1) The deterministic greedy rule chooses \(A\), because \(0.68>0.52\). Thompson sampling's probability of choosing \(A\) is \(P(\theta_A>\theta_B\mid\text{data})\), which depends on both posterior distributions, not just their means. (2) It is \(1-5/6=1/6\approx0.1667\), since the continuous posteriors have zero probability of a tie. (3) By symmetry, the two arms have equal probabilities of being larger; those probabilities sum to one because a tie has probability zero, so each is \(0.5\). (4) Plug-in greedy allocation: the posterior draws are not determining the choice if the programmer discards them and selects by posterior means.

Clinical Biostats