Tutorials › Bayesian Clinical Trials › Power Priors for Single-Arm Historical Benchmarks

Historical Controls and Power Priors · Tutorial 321 of 696

Power Priors for Single-Arm Historical Benchmarks

Learn to discount historical control evidence and quantify the posterior probability that a single-arm trial’s response rate exceeds its historical benchmark.

Advanced 11 min read

What You'll Learn

  • Derive the beta posterior for a historical benchmark under a fixed power-prior discount.
  • Update the current single-arm response rate with a binomial likelihood.
  • Calculate a posterior probability comparing uncertain current and historical rates.
  • Distinguish an uncertain historical benchmark from a fixed threshold based on its posterior mean.
  • Assess assumptions and limitations when interpreting a single-arm historical comparison.

Why Compare a Single-Arm Trial With a Historical Benchmark?

A single-arm trial estimates the response probability among participants receiving the study treatment, but it has no concurrent control group. A historical control rate can provide context: the question may be whether the current response probability is greater than a benchmark informed by earlier control data. A power prior offers one way to discount that historical evidence rather than treating it as fully equivalent to current information.

The resulting comparison is a posterior probability about two rates under a specified model. It is not, by itself, proof that treatment caused an improvement: differences in eligibility, care, outcome assessment, or calendar time can separate the historical benchmark from the rate that would have occurred in the current trial without treatment. The analysis therefore depends on both the probability model and the clinical credibility of the benchmark.

Definition: Let \(X_H\) be the number of responses among \(n_H\) historical control participants, and let \(p_H\) denote the historical control response probability. Let \(Y\) be the number of responses among \(n\) participants in a current single-arm trial, and let \(p_T\) denote the current treatment-arm response probability. A fixed power-prior discount \(a\), with \(0 \le a \le 1\), raises the historical likelihood to the power \(a\). The posterior comparison of interest can be \(P(p_T > p_H \mid X_H,Y)\), the probability that the current response rate exceeds the uncertain historical rate given both data sets and the specified priors.

Model and Assumptions

Assume the historical outcomes are conditionally independent Bernoulli observations with common response probability \(p_H\), so \(X_H \mid p_H \sim \operatorname{Binomial}(n_H,p_H)\). Assign \(p_H\) a beta prior, \(p_H \sim \operatorname{Beta}(\alpha_H,\beta_H)\), with positive shape parameters. For a fixed discount \(a\), the historical contribution to the posterior is proportional to

$$ \left[p_H^{X_H}(1-p_H)^{n_H-X_H}\right]^a $$

The current outcomes are assumed conditionally independent Bernoulli observations with common response probability \(p_T\), so \(Y \mid p_T \sim \operatorname{Binomial}(n,p_T)\). Assign \(p_T \sim \operatorname{Beta}(\alpha_T,\beta_T)\). For the simple comparison developed here, the two prior distributions are independent and the two data sets have independent likelihoods conditional on their respective response probabilities.

These assumptions require more than correct arithmetic. The endpoint and response definition must be meaningful across data sets; the binomial model treats participants within each group as having the same response probability; and the fixed discount must be chosen or justified for the analysis. Most importantly, borrowing historical information is credible only to the extent that the historical rate is relevant to the current clinical setting. Discounting reduces its mathematical influence, but cannot repair a mismatched endpoint or a materially different patient population.

Updating the Historical Benchmark and Current Rate

Conjugacy makes both updates explicit. Multiplying the historical beta prior by the powered likelihood adds \(aX_H\) to the first shape parameter and \(a(n_H-X_H)\) to the second. The current binomial likelihood updates its beta prior in the usual way. Under the stated independence assumptions, the resulting posteriors are independent:

$$ p_H \mid X_H \sim \operatorname{Beta}\!\left(\alpha_H+aX_H,\ \beta_H+a(n_H-X_H)\right), \qquad p_T \mid Y \sim \operatorname{Beta}\!\left(\alpha_T+Y,\ \beta_T+n-Y\right). $$

This calculation uses a fixed discount. When \(a=1\), the historical likelihood is used without discounting; when \(a=0\), it makes no contribution to the historical parameter’s posterior. For intermediate values, the historical successes and failures each contribute fractionally. This is the likelihood-weighting interpretation of the power prior; it does not make the historical participants actual participants in the current trial.

A derivation follows directly from the beta density. Ignoring factors that do not depend on \(p_H\), the posterior kernel is

$$ p_H^{\alpha_H-1}(1-p_H)^{\beta_H-1} p_H^{aX_H}(1-p_H)^{a(n_H-X_H)} = p_H^{\alpha_H+aX_H-1} (1-p_H)^{\beta_H+a(n_H-X_H)-1}. $$

This is the kernel of the beta distribution with the updated shape parameters shown above. Applying the same multiplication to the current prior and likelihood gives the stated posterior for \(p_T\). With independent priors and separate likelihoods, the joint posterior density factors into the product of these two beta densities.

Worked Example: Discounting a Historical Control Rate

Suppose a hypothetical historical control group has 12 responses among 40 participants. Use a \(\operatorname{Beta}(1,1)\) prior for \(p_H\) and a fixed power-prior discount \(a=0.5\). The binomial model treats all 40 outcomes as independent observations of the same historical response probability.

The discounted historical successes contribute \(0.5(12)=6\), and the discounted failures contribute \(0.5(40-12)=14\). Thus

$$ p_H \mid X_H \sim \operatorname{Beta}(1+6,\ 1+14) =\operatorname{Beta}(7,15). $$

The posterior mean benchmark is \(7/(7+15)=7/22=0.3182\), or about 31.8%. The historical data have added 20 discounted observations to the two prior shape parameters, rather than 40 full observations. This posterior mean is a summary of an uncertain benchmark, not a known control rate.

Comparing the Two Posterior Rates

A direct Bayesian comparison averages over uncertainty in both rates. If \(f_T\) is the posterior density of \(p_T\), and \(F_H\) is the posterior cumulative distribution function of \(p_H\), independence gives

$$ P(p_T>p_H \mid X_H,Y) = \int_0^1 f_T(t)F_H(t)\,dt = \int_0^1 f_H(h)\left[1-F_T(h)\right]\,dh. $$

For each possible current rate \(t\), \(F_H(t)\) is the posterior probability that the historical rate is below \(t\); averaging that quantity over the posterior for \(p_T\) yields the comparison probability. In general, beta distributions with different shape parameters do not reduce this comparison to a simple beta posterior. Numerical integration or draws from the two posterior distributions can evaluate it. In a simulation calculation, draw \(p_T^{(m)}\) and \(p_H^{(m)}\) independently from their respective posteriors, then estimate the probability by the proportion of draws for which \(p_T^{(m)}>p_H^{(m)}\).

Worked Example: Comparing Two Uncertain Rates

Consider hypothetical historical data with 18 responses among 50 controls, a \(\operatorname{Beta}(1,1)\) historical prior, and fixed discount \(a=0.4\). In a separate hypothetical single-arm trial, 30 of 60 participants respond, with an independent \(\operatorname{Beta}(1,1)\) prior for \(p_T\). Assume the binomial sampling and independence conditions stated above.

The historical posterior is

$$ p_H \mid X_H \sim \operatorname{Beta}(1+0.4(18),\ 1+0.4(32)) =\operatorname{Beta}(8.2,13.8). $$

Its mean is \(8.2/22=0.3727\). The current posterior is

$$ p_T \mid Y \sim \operatorname{Beta}(1+30,\ 1+30) =\operatorname{Beta}(31,31), $$

with mean \(31/62=0.5000\). To illustrate the scale of the comparison, approximate each beta posterior by a normal distribution with the same mean and variance. For a \(\operatorname{Beta}(u,v)\) distribution, the variance is \(uv/[(u+v)^2(u+v+1)]\). The historical variance is \(8.2(13.8)/[22^2(23)]=0.01016\), and the current variance is \(31(31)/[62^2(63)]=0.003968\). The difference \(p_T-p_H\) is therefore approximated by a normal distribution with mean \(0.5000-0.3727=0.1273\) and standard deviation \(\sqrt{0.01016+0.003968}=0.1189\).

Under this moment-matched normal approximation,

$$ P(p_T>p_H \mid X_H,Y) \approx \Phi\!\left(\frac{0.1273}{0.1189}\right) = \Phi(1.071) \approx 0.858. $$

Thus, under the specified model, the approximate posterior probability that the current response rate exceeds the uncertain historical rate is 85.8%. This is an approximation, not an exact beta comparison; numerical integration or posterior simulation can evaluate the integral directly. It is also not an 85.8% probability that treatment caused an improvement.

A Fixed Benchmark Is a Different Comparison

Sometimes the decision question is framed against a fixed clinical threshold \(p_0\), rather than an uncertain historical response probability. The posterior probability \(P(p_T>p_0 \mid Y)\) is then a one-posterior tail probability. It is not the same as \(P(p_T>p_H \mid X_H,Y)\): replacing \(p_H\) by its posterior mean ignores uncertainty in the benchmark and will generally change the result.

Worked Example: Comparing With a Plug-In Historical Benchmark

Suppose hypothetical historical data show 14 responses among 50 controls. Use a \(\operatorname{Beta}(1,1)\) prior and \(a=0.4\). The historical posterior is \(\operatorname{Beta}(6.6,15.4)\), with mean \(6.6/22=0.3000\). For illustration only, treat this posterior mean as a fixed threshold. In a hypothetical current trial, 3 of 4 participants respond, with a \(\operatorname{Beta}(1,1)\) prior. The current posterior is \(\operatorname{Beta}(4,2)\).

For \(p_T \sim \operatorname{Beta}(4,2)\), the density is \(20t^3(1-t)\), so its cumulative distribution function is \(F_T(t)=5t^4-4t^5\). The posterior probability of exceeding the fixed 0.3000 threshold is

$$ P(p_T>0.3000 \mid Y) =1-F_T(0.3000) =1-\left[5(0.3000)^4-4(0.3000)^5\right] =1-0.03078 =0.9692. $$

This calculation gives a 96.92% posterior probability that \(p_T\) exceeds the fixed plug-in threshold. It does not give the probability that \(p_T\) exceeds the uncertain \(p_H\). With only four current participants, the current posterior remains broad; the high probability reflects the particular threshold and prior as well as the observed responses.

Interpretation, Diagnostics, and Limitations

The comparison probability is conditional on the likelihoods, priors, discount, and independence assumptions. A posterior probability such as 0.858 can be compared with a prespecified decision threshold, but that threshold is a decision rule, not the frequentist significance level \(\alpha\). The clinical consequences of false-positive and false-negative decisions, and the costs of treatment or further study, should inform the decision threshold.

Useful checks include verifying that the response counts lie between zero and their sample sizes, checking the discount and prior values used in both the reported posterior and any software, and confirming that the comparison integrates uncertainty in both rates when that is the intended estimand. Repeating the comparison over scientifically plausible discounts and prior choices can reveal whether the conclusion depends heavily on borrowing. A striking shift in the posterior probability as \(a\) changes is a sensitivity warning, not evidence that one discount is automatically correct.

The central limitation is that a historical comparison cannot rule out systematic differences between historical controls and current patients or care. A power prior changes the weight assigned to historical evidence; it does not model every source of between-study variation or remove confounding. Nor should the historical benchmark be confused with a directly observed concurrent control rate. Randomized comparisons address a different design question, while a fixed-rate analysis answers a different benchmark question.

Common pitfall: “The posterior mean historical rate is 30%, so the posterior probability that the current rate exceeds the historical rate is just \(P(p_T>0.30)\).” This substitutes a point estimate for an uncertain \(p_H\). The correction is to calculate \(P(p_T>p_H \mid X_H,Y)\) from the joint posterior when the target is comparison with an uncertain historical rate; use a fixed threshold only when that is explicitly the target.

Check Your Understanding

Answer each question using the model and conventions defined above. Solutions are included for checking your reasoning.

  1. With a \(\operatorname{Beta}(1,1)\) historical prior, \(X_H=20\), \(n_H=50\), and \(a=0.5\), what is the historical posterior and its mean?
  2. With a \(\operatorname{Beta}(1,1)\) current prior, \(Y=8\), and \(n=20\), what is the current posterior and its mean?
  3. Write an integral for \(P(p_T>p_H \mid X_H,Y)\) using the historical posterior density and the current posterior cumulative distribution function.
  4. Why is \(P(p_T>p_H \mid X_H,Y)\) not automatically the probability that treatment caused a higher response rate?

Solutions. (1) Discounted successes are \(0.5(20)=10\), and discounted failures are \(0.5(30)=15\). Thus \(p_H \mid X_H \sim \operatorname{Beta}(11,16)\), with mean \(11/27=0.4074\). (2) The current posterior is \(\operatorname{Beta}(1+8,1+12)=\operatorname{Beta}(9,13)\), with mean \(9/22=0.4091\). (3) If \(f_H\) is the historical posterior density and \(F_T\) the current posterior cumulative distribution function, the probability is \(\int_0^1 f_H(h)[1-F_T(h)]\,dh\). (4) The posterior comparison is conditional on a model and an external historical benchmark. Differences between historical and current populations, care, or measurement could explain a rate difference; a single-arm comparison does not by itself identify a causal treatment effect.

Clinical Biostats