Tutorials › Biostatistics › Bayesian Borrowing of Historical Control Data

Bayesian Clinical Trial Design

Bayesian Borrowing of Historical Control Data

A practical and computational guide to incorporating historical control information into clinical trials using Bayesian methods, including prior construction, exchangeability, power priors, commensurate priors, robust borrowing, effective sample size, posterior inference, sensitivity analysis, and a complete worked example.

Advanced 18 min read

What You'll Learn

  • Why historical control data can improve precision while creating exchangeability concerns
  • How historical data are represented as Bayesian prior information
  • How power priors, commensurate priors, and robust priors implement borrowing
  • How effective sample size changes as historical information is discounted
  • How to calculate posterior probabilities using a Beta-Binomial model
  • How to perform sensitivity analyses and evaluate the operating characteristics of historical borrowing

Introduction

Clinical trials often contain less information than investigators would ideally like. This is particularly important when the control treatment has a well-established response rate from previous studies, registries, natural-history cohorts, or earlier trials.

Suppose a new Phase II study is evaluating an experimental treatment against a control treatment. A conventional analysis may use only the concurrently randomized control patients.

If the current trial contains relatively few control patients, however, a substantial amount of relevant information from previous studies may be left unused.

Bayesian borrowing provides a formal framework for incorporating historical control information into the current analysis.

Key idea: Historical control observations can be treated as prior information about the current control parameter. The current trial data then update that prior distribution to produce a posterior distribution.

The important qualification is that historical patients cannot automatically be assumed to be exchangeable with current patients.

Differences in patient population, disease severity, eligibility criteria, endpoint definitions, follow-up, standard of care, geography, calendar time, and trial conduct can make historical information less relevant to the current study.

Bayesian borrowing is therefore not simply a mechanism for increasing the sample size. It is a framework for quantifying how much historical information should influence the current inference.

The Basic Problem

Suppose a randomized clinical trial contains an experimental treatment arm and a concurrent control arm.

Let:

  • \(p_C\) = current control response probability
  • \(p_T\) = current treatment response probability

The treatment effect on the risk-difference scale is:

$$ \Delta=p_T-p_C $$

or on the risk-ratio scale:

$$ RR=\frac{p_T}{p_C} $$

Historical studies may contain information about:

$$ p_H $$

where \(p_H\) represents the historical control response probability.

The fundamental question is:

How similar is \(p_H\) to the control response probability that would be observed in the current trial?

If the historical and current populations are highly comparable, borrowing can be substantial.

If they differ materially, borrowing should be reduced or eliminated.

Why Borrow Historical Controls?

There are several reasons historical borrowing can be attractive.

  • Control patients may be difficult or expensive to recruit.
  • The disease may be rare.
  • Randomization may require a large number of concurrent controls.
  • Historical control information may be highly informative.
  • The control response rate may be relatively stable across studies.
  • Additional control enrollment may expose patients to an established treatment without directly increasing information about the experimental treatment.

For example, if 100 historical control patients have been studied under conditions very similar to the current trial, it may be inefficient to pretend that those observations contain no information whatsoever.

But historical data are not free information. Borrowing creates a bias-variance tradeoff. More borrowing can reduce posterior uncertainty when the historical information is relevant, but excessive borrowing can pull inference toward an inappropriate historical control rate when the current and historical populations differ.

Current Data Versus Historical Data

Suppose the current control arm contains:

$$ Y_C\mid p_C\sim\operatorname{Binomial}(n_C,p_C) $$

where \(Y_C\) is the number of responses among \(n_C\) concurrent control patients.

Suppose the historical dataset contains:

$$ Y_H\mid p_H\sim\operatorname{Binomial}(n_H,p_H) $$

The key Bayesian modeling question is whether:

$$ p_H=p_C $$

or whether they should instead be treated as related but potentially different parameters.

Complete Pooling: The Strongest Form of Borrowing

The simplest approach is to assume that historical and current controls have exactly the same response probability:

$$ p_H=p_C=p $$

Under this assumption, the historical and current observations can be pooled.

Suppose the historical study has:

$$ Y_H=60,\qquad n_H=100 $$

and the current trial has:

$$ Y_C=15,\qquad n_C=50 $$

Complete pooling would treat all 150 observations as coming from the same Bernoulli probability.

The combined observed response rate would be:

$$ \hat p = \frac{60+15}{100+50} = 0.50 $$

This can be statistically efficient if the exchangeability assumption is credible.

However, complete pooling can be dangerous if the historical study is not representative of the current trial.

Why Complete Pooling Can Fail

Imagine that the historical control response rate was 60%, while the current control response rate is actually closer to 30%.

If historical data are treated as fully exchangeable, the historical observations will pull the current estimate toward 60%.

That can create a seriously misleading analysis.

Historical borrowing should therefore be conditional on an explicit exchangeability assumption. The more questionable that assumption is, the more cautious the borrowing strategy should be.

Bayesian Prior Distributions

In Bayesian inference, uncertainty about an unknown parameter is represented with a probability distribution.

For a binary control endpoint, a natural prior is the Beta distribution:

$$ p_C\sim\operatorname{Beta}(a,b) $$

Its density is proportional to:

$$ \pi(p_C) \propto p_C^{a-1}(1-p_C)^{b-1} $$

The prior mean is:

$$ E(p_C)=\frac{a}{a+b} $$

and the prior effective sample size is commonly interpreted as:

$$ ESS_{\text{prior}}=a+b $$

This makes the Beta distribution especially convenient for explaining historical borrowing.

Historical Data as a Prior

Suppose historical controls produce:

$$ Y_H=y_H,\qquad n_H $$

Under a Beta-Binomial model, a simple historical-data prior is:

$$ p_C\sim \operatorname{Beta}(a+y_H,b+n_H-y_H) $$

when the baseline prior is:

$$ \operatorname{Beta}(a,b) $$

This is equivalent to treating the historical observations as prior evidence about the current control response probability.

However, using all historical observations with full weight is often too strong.

Bayesian borrowing methods therefore introduce mechanisms for discounting historical information.

The Power Prior

One of the most intuitive approaches is the power prior.

Let the historical likelihood be:

$$ L_H(p_C) $$

The power prior raises the historical likelihood to a power \(\alpha\):

$$ \pi(p_C\mid\alpha,D_H) \propto L_H(p_C)^\alpha\pi_0(p_C) $$

where:

  • \(D_H\) = historical data
  • \(\alpha\) = borrowing weight
  • \(\pi_0(p_C)\) = baseline prior

The borrowing parameter satisfies:

$$ 0\le\alpha\le1 $$

Interpretation is straightforward:

Borrowing Weight Interpretation
\(\alpha=0\) No historical borrowing
\(\alpha=0.25\) Historical likelihood contributes at 25% of its original information
\(\alpha=0.50\) Approximately half-weight borrowing
\(\alpha=0.75\) Strong historical borrowing
\(\alpha=1\) Full historical borrowing

Power Prior for a Beta-Binomial Endpoint

Suppose:

$$ Y_H\sim\operatorname{Binomial}(n_H,p_C) $$

and:

$$ p_C\sim\operatorname{Beta}(a_0,b_0) $$

The historical likelihood is proportional to:

$$ L_H(p_C) \propto p_C^{y_H}(1-p_C)^{n_H-y_H} $$

After applying the power parameter:

$$ L_H(p_C)^\alpha \propto p_C^{\alpha y_H} (1-p_C)^{\alpha(n_H-y_H)} $$

Therefore the resulting prior is:

$$ p_C \sim \operatorname{Beta} \left( a_0+\alpha y_H,\; b_0+\alpha(n_H-y_H) \right) $$

This makes the interpretation particularly transparent.

Effective historical sample size: For the simple Beta-Binomial power prior, the historical contribution has an effective sample size of approximately \(\alpha n_H\).

A Worked Example

Consider a randomized Phase II study evaluating a new therapy for a rare disease.

The primary endpoint is binary response at a prespecified assessment time.

The current trial contains 50 concurrent control patients.

Historical studies contain 100 control patients who received the same nominal standard of care under broadly similar eligibility criteria.

The historical data contain:

$$ y_H=60,\qquad n_H=100 $$

so the historical response rate is:

$$ \hat p_H=\frac{60}{100}=0.60 $$

The current trial observes:

$$ y_C=15,\qquad n_C=50 $$

giving a current observed response rate of:

$$ \hat p_C=\frac{15}{50}=0.30 $$

The historical and current estimates are noticeably different.

This makes the example particularly useful because it illustrates why discounting historical information can matter.

Step 1: Specify a Weak Baseline Prior

Start with a weak Beta prior:

$$ p_C\sim\operatorname{Beta}(1,1) $$

This is the uniform distribution on the interval from 0 to 1.

Its effective sample size is:

$$ ESS_0=1+1=2 $$

The baseline prior therefore contributes relatively little information compared with the current and historical datasets.

Step 2: Apply a 50% Historical Borrowing Weight

Suppose the investigators decide that a power-prior weight of:

$$ \alpha=0.50 $$
will be used for this illustration.

The historical contribution becomes:

$$ \alpha y_H = 0.50(60) = 30 $$

and:

$$ \alpha(n_H-y_H) = 0.50(40) = 20 $$

Therefore the historical-data prior is:

$$ p_C \sim \operatorname{Beta}(1+30,1+20) $$

or:

$$ \boxed{ p_C\sim\operatorname{Beta}(31,21) } $$

The prior mean is:

$$ E(p_C) = \frac{31}{31+21} = 0.5962 $$

Thus the historical information centers the prior around approximately 59.6%.

Step 3: Add the Current Control Data

The current likelihood contributes:

$$ y_C=15 $$

responses and:

$$ n_C-y_C=35 $$

nonresponses.

Because the Beta prior is conjugate to the Binomial likelihood, the posterior is also Beta:

$$ p_C\mid D_H,D_C \sim \operatorname{Beta} (31+15,\;21+35) $$

Therefore:

$$ \boxed{ p_C\mid D_H,D_C \sim \operatorname{Beta}(46,56) } $$

Step 4: Calculate the Posterior Mean

The posterior mean is:

$$ E(p_C\mid D) = \frac{46}{46+56} $$
$$ \boxed{ E(p_C\mid D)\approx0.451 } $$

The current observed response rate was 30%, while the historical response rate was 60%.

The Bayesian posterior lies between these values because both sources of information contribute to the inference.

Important: The posterior is not simply an arithmetic average of 30% and 60%. The amount of influence each dataset has depends on its information content and the borrowing specification.

Step 5: Understand the Effective Sample Size

The historical effective sample size under the 50% power prior is:

$$ ESS_H=\alpha n_H $$

so:

$$ ESS_H=0.50(100)=50 $$

The current trial contributes 50 actual patients.

The approximate total information is therefore comparable to:

$$ 50+50=100 $$

current-equivalent observations, apart from the small baseline-prior contribution.

This illustrates why historical borrowing can materially increase precision.

Compare Different Borrowing Weights

One of the most useful ways to understand Bayesian borrowing is to repeat the analysis under different values of \(\alpha\).

\(\alpha\) Historical ESS Prior Mean Interpretation
0 0 50% No historical information
0.25 25 Approximately 59% Light borrowing
0.50 50 Approximately 60% Moderate borrowing
0.75 75 Approximately 60% Strong borrowing
1.00 100 Approximately 60% Full borrowing

The exact posterior distribution should be recalculated for each borrowing weight rather than interpreting the weight itself as a posterior probability.

No Borrowing: The Reference Analysis

A useful sensitivity analysis begins with:

$$ \alpha=0 $$

This means that the historical dataset contributes no information.

With the baseline prior:

$$ p_C\sim\operatorname{Beta}(1,1) $$

and current data of 15 responses among 50 patients:

$$ p_C\mid D_C \sim \operatorname{Beta}(16,36) $$

The posterior mean is:

$$ \frac{16}{52} \approx0.3077 $$

which is close to the observed current response rate of 30%.

Full Borrowing

At the other extreme:

$$ \alpha=1 $$

the historical observations are given full weight.

The prior becomes:

$$ p_C \sim \operatorname{Beta}(61,41) $$

After adding the current data:

$$ p_C\mid D \sim \operatorname{Beta}(76,76) $$

with posterior mean:

$$ E(p_C\mid D) = \frac{76}{152} = 0.50 $$

Notice what happened.

The current data suggested a control response rate of 30%, while the historical data suggested 60%.

Full borrowing produces a posterior centered at 50%.

Lesson: When historical and current data disagree, full borrowing can substantially move the current inference toward the historical estimate. Whether this is appropriate depends on the scientific plausibility of exchangeability.

Dynamic Borrowing

A major development beyond fixed power priors is dynamic borrowing.

Instead of deciding in advance that:

$$ \alpha=0.50 $$

the model allows the amount of borrowing to depend on how compatible the historical and current data appear to be.

Conceptually:

$$ \text{similar data} \Rightarrow \text{more borrowing} $$

while:

$$ \text{inconsistent data} \Rightarrow \text{less borrowing} $$

This is particularly attractive because the investigator does not need to pretend that historical information is either perfectly exchangeable or completely irrelevant.

Commensurate Priors

The commensurate prior approach explicitly models the similarity between the historical and current parameters.

Suppose:

$$ p_H $$

represents the historical control rate and:

$$ p_C $$

represents the current control rate.

A commensurate model might specify:

$$ p_C\mid p_H,\tau \sim \text{distribution centered around }p_H $$

where \(\tau\) controls how closely the two parameters are expected to agree.

A high degree of commensurability means:

$$ p_C\approx p_H $$

while weak commensurability allows:

$$ p_C $$

to differ substantially from:

$$ p_H $$

The model can therefore learn how compatible the datasets are.

The Concept of a Commensurability Parameter

A useful conceptual model is:

$$ \theta_C\mid\theta_H,\tau \sim N(\theta_H,\tau^{-1}) $$

where:

  • \(\theta_H\) = historical parameter
  • \(\theta_C\) = current parameter
  • \(\tau\) = precision describing how closely they are expected to agree

Large \(\tau\) implies strong similarity.

Small \(\tau\) implies greater uncertainty about whether the historical and current parameters are comparable.

Intuition: A commensurate prior does not simply ask whether historical data should be used. It asks how strongly the historical parameter should inform the current parameter.

Robust Bayesian Borrowing

Another important strategy is the robust mixture prior.

Instead of assuming that the historical prior is always appropriate, construct a mixture of an informative historical prior and a weakly informative or vague component.

For example:

$$ \pi(p_C) = w\pi_H(p_C) + (1-w)\pi_V(p_C) $$

where:

  • \(\pi_H\) = historical-data prior
  • \(\pi_V\) = weakly informative prior
  • \(w\) = prior probability assigned to the historical component

This creates protection against complete dependence on historical information.

Why Robust Priors Are Useful

Suppose historical data strongly suggest a 60% response rate.

But there is a possibility that the historical population is not comparable to the current population.

A robust prior can represent both possibilities:

1
Commensurate component: the historical study is relevant and should be borrowed.
2
Weak component: the historical study is not sufficiently comparable, so the current data should dominate.

The posterior probability assigned to the components can then shift as current data accumulate.

Historical Borrowing Does Not Mean Replacing Randomization

Historical borrowing is often discussed in the context of randomized trials, but it should not be confused with abandoning concurrent controls.

A randomized trial may still contain a concurrent control arm while borrowing additional information about the expected control response.

The current randomized comparison remains the primary source of protection against many forms of bias.

Important distinction: Bayesian borrowing can augment the information available about a control parameter. It does not automatically eliminate the need for a concurrent control group when a concurrent control is scientifically and operationally appropriate.

Borrowing and Treatment Effect Estimation

Suppose the treatment arm contains:

Y_T\sim\operatorname{Binomial}(n_T,p_T) $$

and the current control arm contains:

$$ Y_C\sim\operatorname{Binomial}(n_C,p_C) $$

The treatment effect is:

$$ \Delta=p_T-p_C $$

Bayesian inference produces a posterior distribution for both parameters.

The clinically important quantity may therefore be:

$$ P(\Delta>0\mid D) $$

or, more usefully, the probability that the treatment exceeds a clinically meaningful effect:

$$ P(\Delta>\delta\mid D) $$

where \(\delta\) is the minimum clinically important treatment effect.

Example: Probability of a Clinically Meaningful Benefit

Suppose the clinically meaningful improvement is:

$$ \delta=0.10 $$

The Bayesian decision criterion could be:

$$ P(p_T-p_C>0.10\mid D)>0.95 $$

This is fundamentally different from asking whether a conventional p-value is less than 0.05.

The Bayesian analysis directly describes the posterior probability that the treatment effect exceeds a clinically meaningful threshold.

Historical Borrowing and Effective Sample Size

Effective sample size is one of the most intuitive ways to communicate the strength of historical borrowing.

For a power prior:

$$ ESS_H\approx\alpha n_H $$

For example, if:

$$ n_H=200 $$

and:

$$ \alpha=0.25 $$

then:

$$ ESS_H\approx50 $$

The historical dataset contains 200 patients, but under the selected borrowing scheme it contributes information approximately equivalent to 50 fully exchangeable observations.

Caution: Effective sample size is model-dependent. It should not be interpreted as a literal replacement of 50 new patients in every inferential respect.

Historical Data Quality Matters

Before borrowing historical information, investigators should examine the historical dataset itself.

Relevant characteristics include:

  • Eligibility criteria
  • Baseline disease severity
  • Demographic characteristics
  • Prior treatment history
  • Concomitant medications
  • Endpoint definition
  • Assessment schedule
  • Follow-up duration
  • Protocol requirements
  • Standard of care
  • Geographic region
  • Calendar period
  • Study conduct

A historical response rate should not be treated as a universal biological constant.

Calendar-Time Drift

One particularly important source of historical-data incompatibility is calendar time.

Standard-of-care therapy can change substantially over several years.

Diagnostic criteria can change.

Supportive care can improve.

Patient referral patterns can change.

The apparent response rate can therefore change even if the underlying treatment has not changed.

Historical borrowing should account for temporal drift. Data from a highly similar study conducted recently may be much more relevant than data from an apparently similar study conducted a decade earlier.

Endpoint Compatibility

Historical borrowing is especially sensitive to differences in endpoint definitions.

For example, "response at 12 weeks" is not necessarily interchangeable with "best overall response during treatment."

Similarly:

  • Independent central review may differ from investigator assessment.
  • Different response criteria may produce different response rates.
  • Different assessment windows can change the observed event probability.
  • Different missing-data rules can alter the estimand.

Historical data should therefore be harmonized to the extent possible before they are incorporated into the model.

Exchangeability

The statistical concept underlying historical borrowing is exchangeability.

Informally, exchangeability means that after accounting for relevant characteristics, the historical and current observations can reasonably be treated as arising from related populations.

It does not necessarily mean that every patient has exactly the same risk.

Rather, it represents a modeling assumption about the relationship between datasets.

A useful conceptual hierarchy is:

Assumption Borrowing Behavior
Strong exchangeability Substantial borrowing
Moderate exchangeability Discounted borrowing
Weak exchangeability Limited borrowing
No plausible exchangeability No borrowing

Historical Data Can Be Highly Informative—and Highly Misleading

This is one of the most important concepts in Bayesian historical borrowing.

Suppose the historical sample size is enormous.

One might assume that this automatically makes the historical information valuable.

But if the historical study is systematically different from the current study, the large sample size can actually make the resulting inference more strongly misleading.

More information is not necessarily better information. The value of historical data depends on both its statistical precision and its relevance to the current estimand.

Prior-Data Conflict

A useful Bayesian diagnostic is prior-data conflict.

Suppose the historical data strongly favor:

$$ p_C\approx0.60 $$

while current data favor:

$$ p_C\approx0.30 $$

This discrepancy should not simply be hidden by the model.

Instead, investigators should investigate whether it reflects:

  • Random variation
  • Population differences
  • Changes in standard of care
  • Different endpoint definitions
  • Different follow-up
  • Changes in trial conduct
  • True temporal changes in the disease or treatment effect

Sequential Updating

One advantage of the Bayesian framework is that historical information can be combined with current information sequentially.

Conceptually:

$$ \text{Historical data} \rightarrow \text{Prior} \rightarrow \text{Current data} \rightarrow \text{Posterior} $$

As additional current patients are enrolled, the posterior can be updated.

The historical information does not have to be treated as if it were newly collected current data.

Power Prior Versus Commensurate Prior

Feature Power Prior Commensurate Prior
Historical information Discounted by a weight Linked through a similarity model
Borrowing mechanism Explicit power parameter Precision/commensurability parameter
Fixed borrowing possible? Yes Yes, depending on specification
Dynamic borrowing Possible Natural feature
Interpretation Fractional historical likelihood Similarity between historical and current parameters

Robust Mixture Priors

Feature Standard Historical Prior Robust Mixture Prior
Historical component Yes Yes
Weak component Not necessarily Yes
Protection against incompatibility Limited Explicitly modeled
Borrowing Potentially strong Can adapt through mixture weighting

Bayesian Borrowing Is a Modeling Decision

It is tempting to think of historical borrowing as a computational trick for reducing sample size.

That is not the correct perspective.

The borrowing model encodes a scientific belief:

\[ \text{How much should what we learned previously influence what we believe now?} \]

That belief must be justified using clinical and statistical knowledge.

Type I Error Under Bayesian Decision Rules

Bayesian analyses do not automatically have a frequentist type I error of 5%.

For example, suppose a decision rule is:

$$ P(\Delta>0\mid D)>0.95 $$

The posterior probability threshold of 0.95 is not automatically equivalent to a frequentist one-sided type I error of 5%.

The frequentist operating characteristics must be evaluated through simulation or exact calculation under a specified data-generating process.

Critical distinction: A Bayesian posterior probability threshold and a frequentist type I error threshold are different quantities. If a trial requires frequentist operating characteristics, they should be evaluated explicitly under the proposed Bayesian design.

Simulation-Based Operating Characteristics

Suppose the trial declares success when:

$$ P(\Delta>\delta\mid D)>0.95 $$

To evaluate frequentist type I error, simulate trials under a null scenario in which:

$$ \Delta\le\delta $$

For each simulated trial:

  1. Generate historical data according to the assumed historical scenario.
  2. Generate current control data.
  3. Generate current treatment data.
  4. Fit the Bayesian model.
  5. Calculate the posterior probability.
  6. Record whether the decision criterion is crossed.

After many simulated trials:

$$ \widehat{\alpha} = \frac{\text{number of false declarations}} {\text{number of simulated trials}} $$

provides an empirical estimate of the type I error.

Power Simulation

Similarly, power can be estimated by simulating under a clinically meaningful alternative.

For example:

$$ \Delta=\delta_{\text{target}} $$

Then:

$$ \widehat{\text{Power}} = \frac{\text{number of successful simulated trials}} {\text{number of simulated trials}} $$

This is often essential for Bayesian clinical trial designs involving nonconjugate models, dynamic borrowing, multiple endpoints, or hierarchical structures.

R: Basic Beta-Binomial Historical Borrowing

The simple conjugate example can be implemented directly in R.

y_hist <- 60
n_hist <- 100

y_curr <- 15
n_curr <- 50

alpha_borrow <- 0.50

a0 <- 1
b0 <- 1

a_prior <- a0 + alpha_borrow * y_hist
b_prior <- b0 + alpha_borrow * (n_hist - y_hist)

a_post <- a_prior + y_curr
b_post <- b_prior + (n_curr - y_curr)

a_prior
b_prior
a_post
b_post

For the example:

a_prior
# 31

b_prior
# 21

a_post
# 46

b_post
# 56

Posterior Mean

posterior_mean <- a_post /
  (a_post + b_post)

posterior_mean
# approximately 0.451

The posterior mean is approximately 45.1%.

Posterior Probability That the Control Rate Exceeds 40%

A Bayesian analysis can answer direct probability questions.

For example:

$$ P(p_C>0.40\mid D) $$

In R:

prob_above_40 <- 1 - pbeta(
  0.40,
  a_post,
  b_post
)

prob_above_40

This is a posterior probability, not a p-value.

Posterior Credible Interval

A 95% equal-tail credible interval can be calculated with:

qbeta(
  c(0.025, 0.975),
  a_post,
  b_post
)

This interval has a Bayesian interpretation: conditional on the model and prior assumptions, it contains 95% of the posterior probability for the parameter.

That interpretation differs from the frequentist interpretation of a 95% confidence interval.

Compare Borrowing Weights in R

posterior_summary <- function(
  borrow,
  y_hist,
  n_hist,
  y_curr,
  n_curr,
  a0 = 1,
  b0 = 1
) {

  a_prior <- a0 + borrow * y_hist
  b_prior <- b0 + borrow * (n_hist - y_hist)

  a_post <- a_prior + y_curr
  b_post <- b_prior + n_curr - y_curr

  c(
    borrowing = borrow,
    prior_ess = a_prior + b_prior,
    posterior_mean =
      a_post / (a_post + b_post),
    lower =
      qbeta(0.025, a_post, b_post),
    upper =
      qbeta(0.975, a_post, b_post)
  )
}

sapply(
  c(0, 0.25, 0.50, 0.75, 1),
  posterior_summary,
  y_hist = y_hist,
  n_hist = n_hist,
  y_curr = y_curr,
  n_curr = n_curr
)

This provides a simple sensitivity analysis showing how the posterior changes as historical borrowing increases.

Plotting the Prior and Posterior

p <- seq(0, 1, length.out = 1000)

prior_density <- dbeta(
  p,
  a_prior,
  b_prior
)

posterior_density <- dbeta(
  p,
  a_post,
  b_post
)

plot(
  p,
  prior_density,
  type = "l",
  lwd = 2,
  xlab = "Control response probability",
  ylab = "Density"
)

lines(
  p,
  posterior_density,
  lwd = 2,
  lty = 2
)

legend(
  "topright",
  legend = c("Historical-data prior", "Posterior"),
  lty = c(1, 2),
  lwd = 2
)

The prior represents what the historical information contributes before the current data are observed.

The posterior represents the updated inference after incorporating the current control data.

Dynamic Borrowing With a Grid Over the Power Parameter

A simple way to explore dynamic borrowing is to evaluate several possible borrowing weights.

weights <- seq(0, 1, by = 0.05)

posterior_means <- sapply(
  weights,
  function(w) {

    a_prior <- 1 + w * y_hist
    b_prior <- 1 + w * (n_hist - y_hist)

    a_post <- a_prior + y_curr
    b_post <- b_prior + n_curr - y_curr

    a_post / (a_post + b_post)
  }
)

data.frame(
  borrowing = weights,
  posterior_mean = posterior_means
)

This creates a borrowing sensitivity curve.

Historical Borrowing in a Randomized Trial

Now suppose the experimental treatment arm contains:

$$ Y_T=25,\qquad n_T=50 $$

The observed treatment response rate is therefore:

$$ \hat p_T=\frac{25}{50}=0.50 $$

The current control arm contains:

$$ Y_C=15,\qquad n_C=50 $$

with an observed response rate of 30%.

The naive observed risk difference is:

$$ \hat\Delta = 0.50-0.30 = 0.20 $$

Bayesian borrowing changes the uncertainty around the control parameter and therefore changes the posterior distribution of the treatment effect.

Bayesian Treatment Effect Model

One simple model is:

$$ p_T\sim\operatorname{Beta}(a_T,b_T) $$

and:

$$ p_C\sim\pi_C(p_C) $$

where \(\pi_C\) is constructed using historical borrowing.

After observing current data:

$$ p_T\mid D_T \sim \operatorname{Beta}(a_T+y_T,b_T+n_T-y_T) $$

while the control posterior incorporates both historical and current control information.

Posterior draws can then be used to calculate:

$$ \Delta^{(s)} = p_T^{(s)}-p_C^{(s)} $$

for posterior draw \(s\).

Monte Carlo Calculation of Treatment Benefit

set.seed(123)

S <- 100000

p_treatment <- rbeta(
  S,
  1 + y_curr_treatment,
  1 + n_curr_treatment - y_curr_treatment
)

p_control <- rbeta(
  S,
  a_post,
  b_post
)

delta <- p_treatment - p_control

mean(delta > 0.10)

quantile(
  delta,
  c(0.025, 0.50, 0.975)
)

The first result estimates:

$$ P(p_T-p_C>0.10\mid D) $$

This is often an intuitive quantity for clinical decision-making.

Borrowing Can Increase Precision

Suppose the historical and current control populations are sufficiently comparable.

Borrowing then effectively adds information about the control response rate.

The posterior distribution of \(p_C\) can become narrower than the posterior based only on the current control patients.

This can increase precision for the treatment effect:

$$ \operatorname{Var}(\Delta\mid D) $$

may decrease because:

$$ \operatorname{Var}(p_C\mid D) $$

is reduced.

But Precision Is Not the Same as Accuracy

This distinction is fundamental.

A narrower posterior distribution is not necessarily better if it is centered around the wrong value.

If historical controls are systematically different from current controls, aggressive borrowing can create a highly precise but biased inference.

The objective is not maximum borrowing. The objective is an appropriate balance between information and robustness to historical-data incompatibility.

Conflict Scenario

Consider two hypothetical situations.

Scenario Historical Rate Current Rate Implication
A 58% 57% Strong compatibility
B 60% 30% Potential prior-data conflict

The appropriate borrowing strategy may be very different between these situations.

Scenario A provides little evidence that the historical and current settings differ.

Scenario B requires investigation and sensitivity analysis.

Historical Borrowing Sensitivity Analysis

A robust analysis should generally examine how conclusions change under different assumptions about historical relevance.

For example:

  • No borrowing
  • Weak borrowing
  • Moderate borrowing
  • Strong borrowing
  • Robust dynamic borrowing

The investigator should then examine whether the clinical conclusion is stable.

A Useful Sensitivity Table

Analysis Historical Weight Purpose
Primary Prespecified dynamic borrowing Primary Bayesian analysis
Sensitivity 1 0 Current data only
Sensitivity 2 0.25 Limited historical borrowing
Sensitivity 3 0.50 Moderate borrowing
Sensitivity 4 1.00 Full historical borrowing

If the conclusion changes dramatically between these analyses, the dependence on the historical-data assumption should be explicitly discussed.

Historical Borrowing and Sample Size

Historical borrowing can also affect prospective trial design.

Suppose the current trial requires a target precision for the control response rate.

If sufficiently relevant historical information is available, the required number of concurrent control patients may be reduced.

However, sample-size calculations should be based on the actual Bayesian model and decision criterion.

It is not generally sufficient to say:

$$ N_{\text{effective}} = N_{\text{current}}+N_{\text{historical}} $$

because the historical observations may be discounted or dynamically weighted.

Bayesian Sample Size by Simulation

A typical simulation-based sample-size workflow is:

1
Specify the historical-data model.
2
Specify the current treatment and control models.
3
Specify the borrowing mechanism.
4
Define the Bayesian decision criterion.
5
Simulate trials under clinically relevant scenarios.
6
Estimate type I error and power.
7
Repeat for alternative sample sizes.
8
Select a design satisfying the prespecified operating characteristics.

What Scenarios Should Be Simulated?

A comprehensive simulation should not assume that the historical data are always perfectly compatible.

For example, consider:

Scenario Historical Parameter Current Control Parameter Purpose
1 0.30 0.30 Exact agreement
2 0.30 0.35 Small drift
3 0.30 0.45 Moderate drift
4 0.30 0.60 Strong disagreement
5 0.60 0.30 Historical rate substantially higher

These scenarios reveal whether the proposed borrowing method behaves appropriately when exchangeability is violated.

Bias Under Historical Mismatch

Suppose the true current control response rate is 30%, but the historical response rate is 60%.

If the historical sample is large and receives substantial weight, the posterior for the current control rate can be pulled upward.

This can make the experimental treatment appear less effective if the treatment response is compared against the borrowed control.

The opposite can also occur.

If the historical control rate is lower than the current control rate, excessive borrowing can make the treatment appear more effective than it actually is relative to the current control population.

Historical borrowing can introduce bias in either direction. The direction depends on the relationship between the historical and current control parameters.

Borrowing Multiple Historical Studies

Often there is not one historical study but several.

Suppose there are:

  • Historical Study 1
  • Historical Study 2
  • Historical Study 3
  • Historical Study 4

A naive approach would pool all observations.

A hierarchical Bayesian model can instead allow study-specific control rates:

$$ p_j\sim\mathcal{D}(\mu,\tau) $$

where \(j\) indexes historical studies.

The studies then share information through the population-level distribution rather than being assumed to have exactly identical response rates.

Hierarchical Historical Borrowing

A conceptual normal-scale model is:

$$ \theta_j\mid\mu,\tau \sim N(\mu,\tau^{-1}) $$

where:

  • \(\theta_j\) = study-specific treatment or control parameter
  • \(\mu\) = overall population mean
  • \(\tau\) = between-study precision

The amount of shrinkage depends on the between-study heterogeneity.

If studies are highly consistent, they borrow strongly from one another.

If studies are heterogeneous, the model reduces the amount of pooling.

Partial Pooling

This is known as partial pooling.

Compare three approaches:

Approach Assumption
No pooling Each study has a completely independent parameter
Complete pooling All studies share exactly the same parameter
Partial pooling Study parameters are related but allowed to differ

Bayesian hierarchical models are particularly useful when several historical datasets exist.

Historical Borrowing for Continuous Endpoints

The same principles apply to continuous endpoints.

Suppose the current control mean is:

$$ \mu_C $$

and the historical control mean is:

$$ \mu_H $$

A normal-normal model might specify:

$$ \mu_C\sim N(\mu_0,\sigma_0^2) $$

with historical information incorporated through the prior distribution.

The same concerns remain:

  • Exchangeability
  • Temporal drift
  • Population differences
  • Endpoint compatibility
  • Borrowing strength
  • Prior-data conflict

Historical Borrowing for Time-to-Event Endpoints

Historical borrowing can also be used for survival endpoints.

For example, historical information may be incorporated into a model for a baseline hazard or regression parameter.

A simplified proportional-hazards representation is:

$$ h(t\mid x) = h_0(t)\exp(x\beta) $$

Historical data might inform:

  • Baseline hazard parameters
  • Regression coefficients
  • Control-arm survival parameters
  • Frailty or random-effect distributions

The modeling is considerably more complex than the conjugate Beta-Binomial example, but the conceptual principle remains the same.

Historical Borrowing and Estimands

An important issue is that the historical and current studies must address compatible estimands.

Suppose the current trial targets:

$$ \text{Response probability at Week 12 under the current treatment strategy} $$

Historical data should provide information relevant to that same quantity.

If the historical study instead measures a different endpoint, population, or treatment strategy, borrowing may not be scientifically justified even if the numbers appear superficially similar.

Historical Borrowing and Missing Data

Historical datasets may have different patterns of missingness than the current trial.

For example, the historical response rate may have been calculated among evaluable patients while the current trial uses an intention-to-treat estimand.

Simply borrowing the historical response count without accounting for the different estimand can create incompatibility.

Harmonize the estimand before harmonizing the data. The question being answered should be comparable before the historical observations are treated as information about the current parameter.

Historical Borrowing and Prior Predictive Checks

Before analyzing current data, investigators can examine what the historical prior predicts.

For a Beta prior:

$$ p_C\sim\operatorname{Beta}(a,b) $$

the prior predictive distribution for a future sample can be calculated.

This asks:

Before seeing the current trial data, what control response counts does the historical prior consider plausible?

If the prior predicts implausible current outcomes, the prior may be too strong or poorly constructed.

Prior Predictive Simulation in R

set.seed(123)

S <- 100000

p_prior <- rbeta(
  S,
  a_prior,
  b_prior
)

y_prior_predictive <- rbinom(
  S,
  size = n_curr,
  prob = p_prior
)

quantile(
  y_prior_predictive,
  c(0.025, 0.50, 0.975)
)

This provides a simple diagnostic for whether the historical-data prior implies reasonable numbers of responses in a future sample of the current trial's size.

Posterior Predictive Distribution

After observing the current data, the posterior can be used to predict future control outcomes.

For a future sample of size \(m\):

$$ Y_{\text{future}}\mid D \sim \operatorname{Beta\text{-}Binomial} (m,a_{\text{post}},b_{\text{post}}) $$

This allows questions such as:

$$ P(Y_{\text{future}}\ge k\mid D) $$

which can be useful for predictive trial planning and decision-making.

Borrowing and Predictive Probability

Historical information can also be incorporated into predictive probability calculations.

For example, suppose a trial is currently underpowered based only on concurrent controls.

A Bayesian model can calculate:

$$ P(\text{eventual success}\mid D_{\text{current}},D_{\text{historical}}) $$

This can support adaptive decisions such as:

  • Continue enrollment
  • Stop for futility
  • Modify the design
  • Declare success

provided those rules were appropriately prespecified and their operating characteristics evaluated.

Common Mistakes

  1. Treating historical data as automatically exchangeable. Similarity must be justified rather than assumed.
  2. Using all historical observations at full weight without sensitivity analysis. Historical mismatch can produce substantial bias.
  3. Choosing a borrowing weight after seeing the current outcome. This can introduce undesirable operating characteristics unless the adaptation is explicitly incorporated into the design.
  4. Calling a posterior probability a p-value. They have different interpretations.
  5. Assuming a 95% posterior probability implies 5% type I error. That equivalence does not automatically hold.
  6. Ignoring historical study differences. Endpoint definitions, eligibility, standard of care, and follow-up can all matter.
  7. Reporting the posterior without describing the prior. The posterior depends on both the prior and the likelihood.
  8. Reporting an effective sample size as if it were literal enrollment. ESS is a model-dependent measure of information.
  9. Failing to simulate prior-data conflict. Historical mismatch should be part of design evaluation.
  10. Optimizing only under the assumption of perfect historical compatibility. A borrowing method should also be examined under plausible departures from exchangeability.

A Practical Historical Borrowing Workflow

1
Define the current trial estimand.
2
Identify potentially relevant historical datasets.
3
Assess population, endpoint, treatment, and calendar-time compatibility.
4
Specify a baseline prior.
5
Select a historical borrowing framework.
6
Specify how borrowing responds to historical-current disagreement.
7
Perform prior predictive checks.
8
Simulate operating characteristics under compatible and incompatible scenarios.
9
Evaluate type I error, power, bias, coverage, and precision.
10
Perform sensitivity analyses across plausible borrowing assumptions.
11
Prespecify the final Bayesian model and decision rule.
12
Document the historical datasets and rationale for borrowing.

What Should Be Prespecified?

A Bayesian historical-borrowing analysis should be fully described before the current outcome data are available.

The statistical analysis plan should specify:

  • Historical datasets included in the analysis
  • Eligibility criteria for historical observations
  • Endpoint harmonization rules
  • Baseline prior
  • Borrowing model
  • Hyperpriors
  • Borrowing parameters
  • Decision criterion
  • Posterior probability threshold
  • Effective sample-size diagnostics
  • Prior predictive checks
  • Sensitivity analyses
  • Simulation scenarios
  • Handling of missing historical data
  • Handling of missing current data
  • Computational methods
  • Convergence diagnostics for MCMC models, when applicable

When MCMC Is Required

The simple Beta-Binomial power-prior example is conjugate and can be solved analytically.

More sophisticated borrowing models are often not.

MCMC may be required when the model contains:

  • Commensurate priors
  • Hierarchical study effects
  • Robust mixture priors
  • Multiple historical datasets
  • Covariate adjustment
  • Time-to-event models
  • Nonlinear treatment effects
  • Complex missing-data models

In those settings, the posterior is represented by a large collection of simulated draws rather than a closed-form distribution.

MCMC Diagnostics

When MCMC is used, investigators should examine:

  • Trace plots
  • Effective sample size
  • Potential scale reduction statistics
  • Autocorrelation
  • Monte Carlo error
  • Posterior predictive checks

A model that produces apparently plausible posterior estimates but has not converged adequately should not be considered reliable.

Borrowing and Regulatory Considerations

Historical borrowing can be particularly attractive in rare diseases, pediatric populations, and other settings where randomized control enrollment is difficult.

However, the use of historical information requires careful justification.

The statistical documentation should make it possible to understand:

  • Where the historical data came from
  • Why they are relevant
  • How they were incorporated
  • How much they were allowed to influence the current analysis
  • What happened when historical and current data disagreed
  • How operating characteristics were evaluated
  • How sensitive conclusions were to borrowing assumptions
Transparency matters: The strength of a Bayesian historical-borrowing analysis comes from making its assumptions explicit rather than hiding them inside a conventional analysis.

Bayesian Borrowing Versus Ignoring Historical Data

Feature No Historical Borrowing Bayesian Historical Borrowing
Current data used Yes Yes
Historical data used No Yes
Potential precision gain Limited to current data Potentially greater
Dependence on exchangeability Low Explicit
Risk from historical mismatch Low Must be assessed
Sensitivity analysis Conventional Essential for borrowing assumptions

The Bias-Variance Tradeoff

Historical borrowing can be understood through the familiar bias-variance tradeoff.

Without borrowing:

\[ \text{less model dependence} \quad\text{but}\quad \text{more uncertainty} \]

With strong appropriate borrowing:

\[ \text{potentially less uncertainty} \quad\text{with acceptable bias} \]

With strong inappropriate borrowing:

\[ \text{less uncertainty} \quad\text{but potentially substantial bias} \]

The statistical objective is therefore not simply to minimize posterior variance.

A Conceptual Comparison of Borrowing Strategies

Method Basic Idea Adaptation to Conflict
Complete pooling Treat historical and current data as exchangeable None
Power prior Discount historical likelihood Depends on borrowing parameter
Commensurate prior Model similarity between historical and current parameters Yes
Robust mixture prior Mix historical and weak components Yes
Hierarchical model Allow multiple studies to share information Through between-study heterogeneity

What Makes a Good Historical Dataset?

Historical information is particularly attractive when:

  • The disease definition is the same.
  • The patient population is highly comparable.
  • The treatment regimen is the same.
  • The control treatment is the same.
  • The endpoint is defined identically.
  • The assessment schedule is comparable.
  • The standard of care has not materially changed.
  • The historical study was conducted recently enough to remain relevant.
  • The historical analysis population is compatible with the current estimand.

What Makes Historical Data Less Suitable?

Borrowing should be treated cautiously when:

  • Eligibility criteria differ substantially.
  • Endpoint definitions differ.
  • Control treatment changed.
  • Supportive care changed.
  • The disease population changed.
  • Historical follow-up differs substantially.
  • There is evidence of strong temporal drift.
  • The historical study had substantially different ascertainment procedures.
  • The historical population is much more or less severe than the current population.

The Most Important Diagnostic: Does the Conclusion Depend on Borrowing?

Suppose the primary Bayesian analysis concludes:

$$ P(\Delta>\delta\mid D)=0.97 $$

but the analysis without historical borrowing produces:

$$ P(\Delta>\delta\mid D)=0.82 $$

That difference is important.

It means the conclusion depends materially on the historical information.

That does not automatically invalidate the Bayesian analysis.

It does mean that the relevance of the historical data becomes especially important to the scientific interpretation.

A Robust Reporting Framework

A transparent report could present:

Analysis Historical Borrowing Posterior Probability of Benefit
Primary Dynamic borrowing Report value
Sensitivity No borrowing Report value
Sensitivity Moderate fixed borrowing Report value
Sensitivity Full borrowing Report value

This makes the contribution of the historical information transparent.

Worked Example Summary

Component Value
Endpoint Binary response
Historical responses 60 / 100
Historical response rate 60%
Current control responses 15 / 50
Current response rate 30%
Baseline prior \(\operatorname{Beta}(1,1)\)
Illustrative borrowing weight 0.50
Historical effective sample size 50
Historical-data prior \(\operatorname{Beta}(31,21)\)
Posterior \(\operatorname{Beta}(46,56)\)
Posterior mean Approximately 45.1%
No-borrowing posterior mean Approximately 30.8%
Full-borrowing posterior mean 50.0%

The Most Important Concept

Bayesian borrowing of historical control data is fundamentally about information relevance.

Historical observations can contain valuable information about what the control response rate is likely to be in the current trial.

But their value depends on how exchangeable they are with the current population.

A useful conceptual sequence is:

1
Historical data provide information about a control parameter.
2
A Bayesian prior represents that information.
3
A borrowing mechanism determines how strongly the historical data influence the current analysis.
4
Current trial data update the historical information.
5
The resulting posterior describes uncertainty about the current control parameter.
6
The posterior can then be propagated into the treatment-effect analysis.

Bottom Line

Bayesian borrowing of historical control data is a formal method for using information from previous studies to improve inference about the current control population. The simplest implementation uses a power prior that discounts the historical likelihood. More sophisticated approaches—including commensurate priors, robust mixture priors, and hierarchical models—allow the amount of borrowing to depend on the compatibility of historical and current data. The central statistical challenge is not computational. It is determining whether the historical data are sufficiently relevant to justify borrowing. When they are, borrowing can increase precision and reduce the amount of concurrent control information required. When they are not, excessive borrowing can introduce bias. For this reason, a credible Bayesian borrowing strategy should include explicit assumptions about exchangeability, prior predictive checks, sensitivity analyses, and operating-characteristic simulations under both compatible and incompatible historical-data scenarios.

References

Ibrahim, J.G. & Chen, M.H. (2000). Power prior distributions for regression models. Statistical Science, 15(1), 46–60.

Hobbs, B.P., Carlin, B.P., Mandrekar, S.J. & Sargent, D.J. (2011). Hierarchical commensurate and power prior models for adaptive incorporation of historical information in clinical trials. Biometrics, 67(3), 1047–1056.

Neuenschwander, B., Capkun-Niggli, G., Branson, M. & Spiegelhalter, D.J. (2010). Summarizing historical information on controls in clinical trials. Clinical Trials, 7(1), 5–18.

Viele, K., Berry, S., Neuenschwander, B., et al. (2014). Use of historical control data for assessing treatment effects in clinical trials. Pharmaceutical Statistics, 13(1), 41–54.

Schmidli, H., Gsteiger, S., Roychoudhury, S., O'Hagan, A., Spiegelhalter, D. & Neuenschwander, B. (2014). Robust meta-analytic-predictive priors in clinical trials with historical control information. Biometrics, 70(4), 1023–1032.

Spiegelhalter, D.J., Abrams, K.R. & Myles, J.P. (2004). Bayesian Approaches to Clinical Trials and Health-Care Evaluation. Wiley.

Berry, S.M., Carlin, B.P., Lee, J.J. & Muller, P. (2010). Bayesian Adaptive Methods for Clinical Trials. CRC Press.