Tutorials › Bayesian Clinical Trials › Sequential Beta-Binomial Updating At Interim Looks

Posterior Computation · Tutorial 108 of 696

Sequential Beta-Binomial Updating At Interim Looks

Learn to update a response-rate posterior at each interim look and understand when successive updates produce the same result as one pooled update.

Intermediate 10 min read

What You'll Learn

  • State the assumptions needed for sequential beta-binomial updating.
  • Update beta parameters after each batch of trial outcomes.
  • Derive why the final posterior depends only on total successes and failures.
  • Calculate the posterior predictive probability of a later batch.
  • Identify data and design features that can invalidate the simple update.
  • Distinguish sequential Bayesian updating from frequentist sequential error control.

Why Update at Interim Looks?

In a clinical trial, outcomes may become available in batches rather than all at once. Investigators may want a current summary of the response probability after each batch has been observed. For a binary endpoint, Bayesian updating makes this possible: the posterior after one look becomes the prior information carried forward to the next look. The procedure is sequential in its bookkeeping, even though—under the model below—the final posterior is the same as if all observed outcomes had been analyzed together.

This equality is useful for checking computations and communicating how evidence accumulates. It does not mean that interim data can be reused as if they were new observations, nor that the model is automatically appropriate for every trial with repeated looks.

Model and Assumptions

Let \(p\) be the probability that a participant has the specified binary response by the specified assessment time. For each observed participant, let the response be 1 for success and 0 for failure. Assume outcomes are conditionally independent given \(p\), and that every participant in this analysis has the same response probability \(p\). Begin with a proper beta prior with positive parameters \(a\) and \(b\).

Definition: If \(p\sim\operatorname{Beta}(a,b)\), its density on \(0<p<1\) is proportional to \(p^{a-1}(1-p)^{b-1}\), where \(a>0\) and \(b>0\). For \(m\) observed participants with \(y\) responses, the beta-binomial update is \(p\mid y,m\sim\operatorname{Beta}(a+y,b+m-y)\).

For sequential updating, divide the observed participants into \(K\) non-overlapping batches. At look \(k\), let \(m_k\) be the number of newly observed outcomes and \(y_k\) the number of successes in that batch. Each participant must appear in exactly one batch, and the recorded outcome must follow the same definition and assessment window throughout. The simple model also assumes that response probabilities do not change across batches after conditioning on \(p\). For example, a major change in eligibility, treatment, or response ascertainment could make one common \(p\) inappropriate.

The Recursive Update

Write the beta parameters available just before look \(k\) as \(a_{k-1}\) and \(b_{k-1}\). The initial values are \(a_0=a\) and \(b_0=b\). The batch contributes \(y_k\) successes and \(m_k-y_k\) failures, so the recursive update is

$$ a_k=a_{k-1}+y_k,\qquad b_k=b_{k-1}+(m_k-y_k). $$

Thus after each interim look, the success parameter increases by the number of new successes, and the failure parameter increases by the number of new failures. These updated parameters summarize the prior information together with all outcomes observed so far; they are not a replacement prior to be combined again with those same outcomes.

Worked Example: Updating at Two Interim Looks

Suppose the response probability has prior \(p\sim\operatorname{Beta}(2,3)\). At the first look, 4 responses are observed among 10 participants. At the second look, a new batch has 3 responses among 5 participants. Assume the binary endpoint and assessment time are unchanged, and the participants in the two batches are distinct.

At the first look, there are \(y_1=4\) successes and \(m_1-y_1=10-4=6\) failures. Therefore,

$$ p\mid\text{first look}\sim \operatorname{Beta}(2+4,\;3+6) =\operatorname{Beta}(6,9). $$

The second batch adds 3 successes and \(5-3=2\) failures. Using the first-look posterior as the current distribution gives

$$ p\mid\text{both looks}\sim \operatorname{Beta}(6+3,\;9+2) =\operatorname{Beta}(9,11). $$

As a check, the pooled data contain \(4+3=7\) responses among \(10+5=15\) participants, hence \(15-7=8\) failures. Updating the original prior once gives \(\operatorname{Beta}(2+7,3+8)=\operatorname{Beta}(9,11)\), the same result. The posterior mean is \(9/(9+11)=0.45\), a posterior summary of the response probability, not the observed response proportion, which is \(7/15\approx0.4667\).

Why the Order of Batches Does Not Matter

The path-independence result follows directly from the likelihood. For a batch with \(m_k\) participants and \(y_k\) successes, the binomial likelihood as a function of \(p\) is proportional to \(p^{y_k}(1-p)^{m_k-y_k}\). The binomial coefficient does not depend on \(p\), so it does not affect the posterior density. Multiplying the likelihood contributions from all batches gives the justification below.

Result: Under the common-\(p\), conditionally independent Bernoulli model, after \(K\) disjoint batches the posterior is \(p\mid\text{data}\sim\operatorname{Beta}\left(a+\sum_{k=1}^{K}y_k,\;b+\sum_{k=1}^{K}(m_k-y_k)\right)\). Consequently, the final posterior depends on the total numbers of successes and failures, not on the order in which the batches arrived.

To verify the result, multiply the prior density by the likelihood for each batch. Let \(Y=\sum_{k=1}^{K}y_k\) and \(F=\sum_{k=1}^{K}(m_k-y_k)\). Then the product is proportional to

$$ p^{a-1}(1-p)^{b-1} \prod_{k=1}^{K}p^{y_k}(1-p)^{m_k-y_k} = p^{a+Y-1}(1-p)^{b+F-1}. $$

This is the kernel of a \(\operatorname{Beta}(a+Y,b+F)\) density. The same result follows from the recursion because adding the batch success counts in any order gives \(a+\sum_k y_k\), and adding the failure counts gives \(b+\sum_k(m_k-y_k)\). This derivation requires that each batch contribute only its new outcomes. Counting earlier participants again would multiply in their likelihoods more than once and would produce incorrect parameters.

Worked Example: Reversing the Batch Order

Let the prior be \(\operatorname{Beta}(1.5,2.5)\). Two distinct batches contain 2 responses among 6 participants and 5 responses among 8 participants. Compare the posterior obtained in each order.

If the 2-of-6 batch arrives first, it adds 2 successes and 4 failures: \(\operatorname{Beta}(1.5+2,2.5+4)=\operatorname{Beta}(3.5,6.5)\). The 5-of-8 batch then adds 5 successes and 3 failures, giving \(\operatorname{Beta}(3.5+5,6.5+3)=\operatorname{Beta}(8.5,9.5)\).

If the 5-of-8 batch arrives first, the interim posterior is \(\operatorname{Beta}(1.5+5,2.5+3)=\operatorname{Beta}(6.5,5.5)\). Adding the 2 successes and 4 failures from the other batch gives \(\operatorname{Beta}(6.5+2,5.5+4)=\operatorname{Beta}(8.5,9.5)\).

Both orders yield 7 successes and \(4+3=7\) failures in total. The matching posterior is therefore also obtained directly as \(\operatorname{Beta}(1.5+7,2.5+7)\). The interim posteriors differ because different data are available at those interim times, but the final posterior agrees.

Predicting a Later Batch

An interim posterior can also describe what outcomes are plausible in a future batch before those outcomes are observed. If the current posterior is \(\operatorname{Beta}(A,B)\), the probability of exactly \(z\) successes in a new batch of \(r\) participants is obtained by averaging the binomial probability over uncertainty about \(p\):

$$ P(Z=z\mid\text{current data}) = \binom{r}{z} \frac{B(A+z,B+r-z)}{B(A,B)},\qquad z=0,\ldots,r, $$

where \(B(\cdot,\cdot)\) is the beta function. This is a beta-binomial predictive probability. It is not generally the same as substituting the posterior mean of \(p\) into a binomial model, because it averages over the full posterior distribution.

Worked Example: Predicting Three Responses in the Next Four Participants

Suppose the current posterior is \(\operatorname{Beta}(8,6)\), for example after a \(\operatorname{Beta}(2,2)\) prior and 6 responses among 10 participants. For a new batch of \(r=4\), calculate the predictive probability of exactly \(z=3\) responses.

Substitute \(A=8\), \(B=6\), \(r=4\), and \(z=3\):

$$ P(Z=3\mid\text{data}) =\binom{4}{3}\frac{B(11,7)}{B(8,6)} =4\frac{(8)(9)(10)(6)}{(14)(15)(16)(17)} =\frac{17280}{57120} \approx0.3025. $$

Thus the posterior predictive probability of exactly 3 responses in the next 4 participants is about 30.25%. If that batch is subsequently observed to contain 3 responses and 1 failure, the posterior becomes \(\operatorname{Beta}(8+3,6+1)=\operatorname{Beta}(11,7)\). The prediction describes possible future data; it does not change the posterior before that data are observed.

Interim Checks, Alternatives, and Limitations

At each look, check that the batch contains only new participants, that successes plus failures equal the batch size, and that the posterior parameters equal the original prior parameters plus the cumulative counts. The posterior parameters must not decrease as valid new observations are added. A mismatch often signals a duplicated batch, an incorrect denominator, or a success/failure coding error.

Also assess whether the assumptions make sense across enrollment time. Different response rates by site, calendar period, dose, or participant subgroup can signal that a single common response probability is inadequate. Dependence between participants or outcome-dependent missingness can likewise undermine the simple likelihood. Incomplete follow-up requires particular care: a participant whose response is not yet known should not automatically be counted as a failure. The analysis may need a model that accounts for the way observation or follow-up occurs.

The computational alternative is to pool all observed outcomes and update once. Under the stated model, pooling and recursive updating are exactly equivalent; sequential bookkeeping is convenient when interim summaries are needed, while pooled updating is a useful independent check. A different model may be needed if the response probability varies over time or across groups. Frequentist sequential procedures address a different issue: repeatedly examining data can affect frequentist error rates, whereas Bayesian posterior updating is a probability calculation conditional on the model and observed data. The equality of the two Bayesian calculations does not by itself establish frequentist operating characteristics or justify a trial decision rule.

The central limitation is therefore not the order of arithmetic but the assumptions behind the likelihood. If batches reflect different response mechanisms, or if participants enter the analysis selectively based on unobserved outcomes, totals alone may not contain all relevant information. A correct sequential update is only as credible as the model and data collection process it represents.

Check Your Understanding

Answer each question, then compare with the brief solutions.

  1. Starting from \(\operatorname{Beta}(3,2)\), a batch contains 5 responses among 12 participants. What is the posterior?
  2. Two successive batches have 1 response among 4 and 6 responses among 9. Explain why reversing their order does not change the final posterior for any fixed beta prior.
  3. A data analyst uses \(\operatorname{Beta}(4,7)\) after the first look, then combines that distribution with all outcomes from both the first and second looks. What is the error?
  4. Give one reason the common-\(p\) beta-binomial model might be inappropriate across interim batches.

Solutions. (1) There are \(12-5=7\) failures, so the posterior is \(\operatorname{Beta}(3+5,2+7)=\operatorname{Beta}(8,9)\). (2) The batches contribute 7 successes and \((4-1)+(9-6)=6\) failures in total. The final posterior is the original \(\operatorname{Beta}(a,b)\) prior updated to \(\operatorname{Beta}(a+7,b+6)\), regardless of the order of addition. (3) The first-look outcomes are counted twice. The second update must use only the new second-batch outcomes, or the analyst must update the original prior once with the pooled data. (4) For example, a change in treatment, eligibility, or response ascertainment could make the response probability differ by batch, violating the single common-\(p\) assumption.

Clinical Biostats