Tutorials › Biostatistics › Prior Distributions: Choosing and Justifying Priors

Bayesian Statistics

Prior Distributions: Choosing and Justifying Priors

A practical and mathematical guide to selecting Bayesian prior distributions, including informative, weakly informative, skeptical, diffuse, and reference priors, prior elicitation, effective sample size, prior predictive checks, sensitivity analysis, clinical trial examples, and transparent prior justification.

Intermediate 20 min read

What You'll Learn

  • What a Bayesian prior distribution represents and how it enters Bayes' theorem
  • How to distinguish informative, weakly informative, skeptical, diffuse, and reference priors
  • How to choose priors using historical data, expert knowledge, and scientific constraints
  • How to interpret prior parameters and effective sample size
  • How to use prior predictive checks and sensitivity analyses to evaluate priors
  • How to justify and document priors in clinical trials and regulatory analyses

Introduction

A Bayesian analysis begins with a statistical model for the data and a probability distribution describing uncertainty about the unknown parameters before the current data are incorporated.

That probability distribution is the prior distribution.

For a parameter \(\theta\), the prior is commonly written:

$$ \pi(\theta) $$

After observing data \(D\), Bayesian inference combines the prior with the likelihood:

$$ \boxed{ \pi(\theta\mid D) \propto L(\theta;D)\pi(\theta) } $$

The result is the posterior distribution.

Choosing the prior is therefore not a cosmetic step in Bayesian analysis. It is part of the statistical model.

Key idea: A prior distribution should be scientifically defensible, mathematically appropriate for the parameter being modeled, transparent to readers, and evaluated for its potential influence on the posterior inference.

What Is a Prior Distribution?

Suppose \(\theta\) is an unknown treatment effect. Before observing the current trial, there is uncertainty about the value of \(\theta\). A Bayesian analysis represents that uncertainty with a probability distribution:

$$ \theta\sim\pi(\theta) $$

The prior distribution describes how plausible different values of \(\theta\) are before incorporating the current data.

This does not necessarily mean the prior represents personal belief in a subjective or informal sense.

It may instead be based on:

  • Previous clinical trials
  • Historical control data
  • Meta-analyses
  • External observational studies
  • Pharmacological knowledge
  • Biological constraints
  • Expert elicitation
  • Weakly informative regularization
  • Conservative or skeptical assumptions

The Prior Is Part of the Model

A common mistake is to think of the prior as an optional adjustment added after the statistical model has been constructed.

It is more accurate to think of a Bayesian model as having at least two important components:

$$ \text{Likelihood} \quad+\quad \text{Prior} $$

The likelihood describes how the observed data arise conditional on the parameters.

The prior describes uncertainty about the parameters before incorporating the current data.

Together they determine the posterior.

Important: The prior is not a substitute for the likelihood. It does not replace the current trial data. It supplies additional structure about parameter uncertainty that is then updated by the observed data.

Bayes' Theorem

For a parameter \(\theta\) and data \(D\):

$$ \pi(\theta\mid D) = \frac{ L(\theta;D)\pi(\theta) }{ \int L(\theta;D)\pi(\theta)\,d\theta } $$

The numerator combines:

  • The likelihood \(L(\theta;D)\)
  • The prior \(\pi(\theta)\)

The denominator normalizes the resulting distribution so that it integrates to one.

The central relationship is:

$$ \boxed{ \text{Posterior} \propto \text{Likelihood} \times \text{Prior} } $$

Why Prior Choice Matters

The effect of a prior depends heavily on how much information the current data contain.

If the dataset is extremely informative, reasonable priors may have little effect on the posterior.

If the dataset is small or sparse, the prior may have a much larger influence.

This can be understood conceptually as:

$$ \text{Posterior information} \approx \text{Prior information} + \text{Information from current data} $$

This is not a universal algebraic identity, but it is a useful way to understand why prior strength matters.

The smaller the information content of the current data, the more important it becomes to justify the prior carefully.

A Simple Beta-Binomial Example

Suppose a clinical trial studies a binary response probability \(p\). The likelihood is:

$$ X\sim\operatorname{Binomial}(n,p) $$

Because \(p\) lies between 0 and 1, a natural prior is a Beta distribution:

$$ p\sim\operatorname{Beta}(\alpha,\beta) $$

The Beta density is:

$$ \pi(p) \propto p^{\alpha-1}(1-p)^{\beta-1} $$

The prior mean is:

$$ E(p) = \frac{\alpha}{\alpha+\beta} $$

The prior variance is:

$$ \operatorname{Var}(p) = \frac{ \alpha\beta }{ (\alpha+\beta)^2(\alpha+\beta+1) } $$

The parameters \(\alpha\) and \(\beta\) therefore control both the location and concentration of the prior.

Prior Mean and Prior Strength

Suppose:

$$ p\sim\operatorname{Beta}(2,8) $$

The prior mean is:

$$ E(p)=\frac{2}{10}=0.20 $$

The quantity:

$$ \alpha+\beta $$

is often interpreted as the prior's equivalent sample size in the Beta-binomial conjugate setting.

Thus:

$$ 2+8=10 $$

can be viewed as approximately 10 prior observations worth of information, with 2 corresponding to prior successes and 8 to prior failures.

Important qualification: The "equivalent sample size" interpretation is especially intuitive for conjugate Beta-binomial models. For more complex priors and models, prior strength is not generally reducible to a literal number of observations.

Prior Effective Sample Size

Effective sample size is a useful way to communicate how strongly a prior can influence an analysis.

For the Beta prior:

$$ \operatorname{Beta}(\alpha,\beta) $$

a common conjugate interpretation is:

$$ ESS_{\text{prior}} = \alpha+\beta $$

For example:

Prior Prior Mean Approx. ESS
\(\operatorname{Beta}(1,1)\) 0.50 2
\(\operatorname{Beta}(2,8)\) 0.20 10
\(\operatorname{Beta}(10,40)\) 0.20 50
\(\operatorname{Beta}(100,400)\) 0.20 500

Notice that the last three priors have the same prior mean.

They differ substantially in how concentrated they are.

Therefore: Choosing a prior mean is not enough. The strength or precision of the prior must also be justified.

Same Mean, Different Prior Strength

Consider two priors:

$$ p\sim\operatorname{Beta}(2,8) $$

and:

$$ p\sim\operatorname{Beta}(20,80) $$

Both have mean:

$$ 0.20 $$

But their concentrations are very different.

The second prior is much more concentrated around 20%.

Consequently, it represents substantially stronger prior information.

Practical lesson: A statement such as "we chose a prior centered at 20%" is incomplete. You should also explain how uncertain the prior is around 20% and why that degree of uncertainty is scientifically reasonable.

Informative Priors

An informative prior contributes substantial information to the posterior.

For example:

$$ \theta\sim N(0.20,0.05^2) $$

expresses relatively strong prior concentration around an effect of 0.20.

Informative priors can be appropriate when high-quality external information exists.

Examples include:

  • A well-established historical control response rate
  • A validated pharmacological parameter
  • Multiple high-quality previous trials
  • Strong biological knowledge
  • Reliable external registry information

The stronger the prior, the stronger the justification should generally be.

Weakly Informative Priors

A weakly informative prior provides some structure without attempting to dominate the current data.

For example:

$$ \beta\sim N(0,2.5^2) $$

might be used for a standardized regression coefficient when very large effects are scientifically implausible but small and moderate effects remain quite plausible.

Weakly informative priors can be particularly useful for:

  • Regularizing unstable estimates
  • Preventing implausibly extreme parameter values
  • Improving computational stability
  • Handling sparse data
  • Reducing sensitivity to separation in logistic regression
Weakly informative does not mean completely noninformative. A weakly informative prior deliberately rules out or downweights parameter values that are inconsistent with reasonable scientific expectations.

Skeptical Priors

A skeptical prior deliberately places greater probability on small or negligible treatment effects.

Suppose \(\theta\) is the treatment effect and:

$$ \theta\sim N(0,\sigma^2) $$

The prior is centered on no treatment effect.

This does not mean the analysis assumes that the treatment has no effect.

Rather, it expresses the idea that large effects are considered less plausible before observing the current trial.

Skeptical priors can be useful in sensitivity analyses or when historical evidence suggests that large effects are uncommon.

Optimistic Priors

An optimistic prior places more prior probability on clinically favorable effects.

For example:

$$ \theta\sim N(0.15,0.10^2) $$

might represent an expectation that the treatment is likely to have a positive effect.

Optimistic priors can be scientifically appropriate when supported by substantial external evidence.

However, they should not be chosen merely because the investigator wants a positive result.

Prior choice must be separated from desired conclusions. A prior should represent defensible pre-data information, not the result the analyst hopes to obtain.

Diffuse and Vague Priors

A diffuse prior attempts to spread probability over a broad range of parameter values.

For example:

$$ \theta\sim N(0,100^2) $$

may appear extremely broad on a scale where plausible treatment effects are close to zero.

However, simply making a prior very broad does not automatically make it "objective" or harmless.

Diffuse priors can place substantial probability on scientifically implausible values.

They can also behave poorly after nonlinear transformations.

Do not confuse "very broad" with "uninformative." A prior can be mathematically broad while still having important consequences for the resulting posterior.

Noninformative and Reference Priors

Some Bayesian analyses use priors designed to have minimal influence under specific theoretical criteria.

These are sometimes called reference priors or objective priors, depending on the construction.

Examples include Jeffreys-type priors.

For a Bernoulli probability \(p\), the Jeffreys prior is:

$$ p\sim\operatorname{Beta}\left(\frac12,\frac12\right) $$

This prior is not literally uniform. The uniform prior is:

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

The two priors therefore represent different constructions.

Important: "Noninformative" is not a universal mathematical property. Whether a prior is weakly influential can depend on the parameterization, likelihood, sample size, and inferential target.

Uniform Priors Are Not Always Neutral

Suppose an investigator assigns:

$$ p\sim\operatorname{Uniform}(0,1) $$

This may look completely neutral.

But consider the transformed parameter:

$$ \operatorname{logit}(p) = \log\left(\frac{p}{1-p}\right) $$

A uniform prior on \(p\) does not correspond to a uniform prior on \(\operatorname{logit}(p)\).

Thus, "uniform" depends on the scale of the parameter.

Key lesson: There is generally no prior that is automatically neutral under every reparameterization.

Choose the Parameterization Carefully

The prior should be considered on the scale where the scientific interpretation is most meaningful.

For example, in logistic regression:

$$ \log\left( \frac{p_i}{1-p_i} \right) = \alpha+\beta X_i $$

A prior on \(\beta\) is a prior on the log-odds scale.

If:

$$ \beta\sim N(0,1) $$

then the implied prior distribution on the odds ratio:

$$ OR=e^\beta $$

is log-normal.

Therefore, prior evaluation should consider both the model parameterization and the scientific scale on which effects will ultimately be interpreted.

Prior Choice Should Begin With the Scientific Question

A useful prior-selection process starts before choosing a distribution.

First ask: What parameter are we trying to describe?

For example:

  • Response probability
  • Mean treatment difference
  • Log hazard ratio
  • Log odds ratio
  • Regression coefficient
  • Variance component
  • Between-study heterogeneity
  • Pharmacokinetic clearance

The scientific meaning of the parameter should guide the prior.

Do not begin with "Which prior distribution should I use?" Begin with "What values of this parameter are scientifically plausible, and how certain am I about that range?"

Step 1: Define a Plausible Range

Suppose \(\theta\) is a treatment effect measured as a difference in response rates. Before selecting a mathematical distribution, identify a plausible range.

For example:

$$ -0.20\le\theta\le0.30 $$

The range might be based on:

  • Historical clinical trials
  • Clinical expertise
  • Biological constraints
  • Regulatory precedent
  • Published literature
  • Mechanistic knowledge

The mathematical prior should then reflect this scientific information.

Step 2: Identify the Center of the Prior

The next question is: Where should the prior distribution be centered?

For example, suppose historical studies suggest an average response rate near 25%.

A prior might therefore be centered around:

$$ p=0.25 $$

But the center alone is insufficient. The analyst must also quantify uncertainty around that value.

Step 3: Quantify Prior Uncertainty

Suppose the prior mean is 0.25. We might believe that values between 0.10 and 0.40 are reasonably plausible.

The prior distribution should place most of its probability within a scientifically defensible range.

This is often more informative than choosing a standard deviation simply because it is mathematically convenient.

A good prior justification describes probability in terms of scientifically meaningful quantities. For example: "The prior places approximately 95% probability between X and Y" is often easier to defend than simply saying "we selected a standard deviation of 0.15."

Prior Elicitation

Prior elicitation is the process of translating expert knowledge into a probability distribution.

Suppose experts are asked about an unknown treatment effect \(\theta\). Instead of asking: "Choose a standard deviation for a normal distribution." it is often better to ask questions such as:

  • What value do you consider most plausible?
  • What values would be surprising?
  • What value would you regard as clinically important?
  • What is the probability that the effect is positive?
  • What is the probability that the effect exceeds a specified threshold?

These answers can then be used to construct or calibrate a prior distribution.

Quantile-Based Prior Elicitation

Suppose an expert believes:

$$ P(\theta<0.05)=0.025 $$

and:

$$ P(\theta<0.45)=0.975 $$

This suggests that approximately 95% of prior probability lies between 0.05 and 0.45.

A candidate prior distribution can then be calibrated to match those quantiles.

This is often more transparent than selecting a distribution first and attempting to invent a scientific explanation afterward.

Prior Predictive Distribution

One of the most important tools for evaluating a prior is the prior predictive distribution.

The prior predictive distribution is:

$$ p(D) = \int p(D\mid\theta)\pi(\theta)\,d\theta $$

It describes the kinds of datasets that the model and prior predict before observing the actual data.

This is extremely useful because a prior can look reasonable when viewed only on the parameter scale but produce absurd predictions on the observable-data scale.

Prior predictive question: "If I generated a study from this model and prior before seeing the data, would the resulting datasets look scientifically plausible?"

Prior Predictive Checking for a Binary Endpoint

Suppose:

$$ p\sim\operatorname{Beta}(2,8) $$

and:

$$ X\sim\operatorname{Binomial}(50,p) $$

The prior predictive distribution describes the number of responses we would expect to observe in 50 patients before seeing the actual trial.

If the prior implies that 0 responses and 50 responses are both highly plausible, despite strong clinical knowledge that neither scenario is realistic, the prior may be poorly calibrated.

Prior Predictive Checking in Clinical Trials

Suppose a historical response rate is approximately 20%. An analyst proposes a very broad prior:

$$ p\sim\operatorname{Beta}(0.5,0.5) $$

The prior itself places substantial probability near both 0 and 1.

That may be inconsistent with the historical clinical context.

A prior predictive simulation may reveal that the model frequently generates trials with extremely high or extremely low response counts.

This does not automatically prove the prior is wrong. It does, however, force the analyst to explain why those predictions are scientifically reasonable.

Prior Predictive Checks vs. Posterior Predictive Checks

Check When? Main Question
Prior predictive Before current data influence inference Does the prior/model generate plausible data?
Posterior predictive After observing data Does the fitted model reproduce important features of the observed data?

Both are valuable, but they answer different questions.

Prior Sensitivity Analysis

A responsible Bayesian analysis should determine whether its conclusions are sensitive to reasonable alternative priors.

Suppose the primary prior is:

$$ \pi_1(\theta) $$

The analyst might also examine:

$$ \pi_2(\theta) \qquad\text{and}\qquad \pi_3(\theta) $$

The resulting posterior distributions can then be compared.

The key question is: Would reasonable prior alternatives change the scientific conclusion?

What Counts as a Reasonable Alternative Prior?

Sensitivity analysis should not merely compare the primary prior with absurd alternatives.

Useful alternatives may include:

  • A skeptical prior
  • A neutral prior
  • A moderately informative prior
  • A weakly informative prior
  • A prior based on a different credible historical source

The purpose is to determine whether conclusions depend on plausible modeling choices.

Good sensitivity analysis asks whether the conclusion is robust to reasonable uncertainty about prior specification. It is not simply an exercise in finding a prior that produces the desired answer.

Worked Example: Response Rate

Suppose a Phase II study enrolls 20 patients and observes 6 responses.

$$ x=6,\qquad n=20 $$

The observed response rate is:

$$ \hat p=\frac{6}{20}=0.30 $$

Suppose the historical response rate is approximately 20%. A moderately informative prior could be:

$$ p\sim\operatorname{Beta}(4,16) $$

The prior mean is:

$$ E(p)=\frac{4}{20}=0.20 $$

The prior effective sample size is approximately:

$$ ESS_{\text{prior}}=20 $$

After observing 6 responses among 20 patients:

$$ p\mid D \sim \operatorname{Beta}(4+6,\;16+14) $$
$$ \boxed{ p\mid D \sim \operatorname{Beta}(10,30) } $$

The posterior mean is:

$$ E(p\mid D) = \frac{10}{40} = 0.25 $$

The observed proportion was 30%, but the posterior mean is 25% because the historical prior contributes information centered at 20%.

Changing the Prior Strength

Now consider a much weaker prior:

$$ p\sim\operatorname{Beta}(1,4) $$

This also has prior mean:

$$ E(p)=0.20 $$

but the effective sample size is only:

$$ ESS_{\text{prior}}=5 $$

The posterior becomes:

$$ p\mid D \sim \operatorname{Beta}(7,18) $$

with posterior mean:

$$ E(p\mid D) = \frac{7}{25} = 0.28 $$

The prior mean is identical to the previous example, but the posterior is closer to the observed 30% because the prior is weaker.

This example illustrates an essential principle: Two priors can encode the same prior expectation but produce different posterior conclusions because they assign different amounts of uncertainty around that expectation.

Historical Data as Prior Information

Historical clinical data are a common source of Bayesian prior information.

Suppose previous studies provide information about a control response rate. A prior can be constructed from those data.

However, historical borrowing should not be automatic. The analyst should consider whether the historical studies are sufficiently comparable with the current study.

Important dimensions include:

  • Patient population
  • Eligibility criteria
  • Disease severity
  • Endpoint definition
  • Assessment timing
  • Standard of care
  • Concomitant therapy
  • Study conduct
  • Geographic region
  • Calendar period
Historical data are informative only to the extent that they are relevant to the current target population and estimand. More historical observations do not automatically mean more valid prior information.

Historical Borrowing Can Be Dangerous

Suppose a historical control group has a response rate of 20%. A new trial has a response rate of 35%. If the historical population is systematically different from the new population, blindly borrowing historical information could bias the posterior toward an inappropriate value.

The problem is not that Bayesian methods "use prior information." The problem is that the prior information may not be transportable.

This is why historical-data priors should be accompanied by sensitivity analyses and, when appropriate, models that allow historical information to be discounted.

Exchangeability and Historical Borrowing

Bayesian hierarchical models can explicitly represent uncertainty about whether different studies are sufficiently similar.

For example:

$$ \theta_i\sim N(\mu,\tau^2) $$

where \(\theta_i\) represents the treatment or control parameter in study \(i\).

The heterogeneity parameter \(\tau\) determines how strongly studies are expected to resemble one another.

If between-study heterogeneity is large, the model naturally permits less information sharing.

This is often preferable to treating historical data as if they were literally observed in the current trial.

Commensurate and Dynamic Borrowing

More sophisticated Bayesian models can allow the amount of borrowing from historical data to depend on how compatible those data are with the current study.

This is sometimes described as dynamic borrowing.

The basic concept is:

$$ \text{More compatible historical data} \Rightarrow \text{more borrowing} $$

while:

$$ \text{Less compatible historical data} \Rightarrow \text{less borrowing} $$

Such models can be useful, but their operating characteristics should be carefully evaluated through simulation before being used in a clinical trial.

Prior Choice for Treatment Effects

Suppose \(\theta\) is a treatment difference. A common prior is:

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

Here:

  • \(\mu_0\) is the prior mean.
  • \(\sigma_0\) is the prior standard deviation.

For example:

$$ \theta\sim N(0,0.10^2) $$

places the prior center at no treatment difference while allowing clinically meaningful positive or negative effects.

If approximately 95% of the prior probability is desired between \(-0.20\) and \(0.20\), then a rough normal approximation is:

$$ \sigma_0\approx\frac{0.20}{1.96} \approx0.102 $$

This provides a transparent way to connect the standard deviation to a scientifically meaningful range.

Skeptical Priors for Clinical Effects

Suppose large treatment effects are considered unlikely. A skeptical prior might be:

$$ \theta\sim N(0,0.05^2) $$

Approximately 95% of the prior probability lies within:

$$ -0.098 \lesssim \theta \lesssim 0.098 $$

Thus, effects larger than approximately 10 percentage points in magnitude are given relatively low prior probability.

The appropriateness of such a prior depends entirely on the scientific context.

Prior Choice for Odds Ratios

Clinical investigators often interpret logistic regression results using odds ratios.

Because the odds ratio is positive:

$$ OR>0 $$

it is often convenient to specify a prior on the log odds-ratio:

$$ \log(OR)=\beta $$

with:

$$ \beta\sim N(0,\sigma^2) $$

This implies:

$$ OR\sim\operatorname{Lognormal}(0,\sigma^2) $$

For example, if:

$$ \beta\sim N(0,0.5^2) $$

then the prior places substantial probability around odds ratios close to 1 while still allowing meaningful treatment effects.

The prior should always be examined on the odds-ratio scale as well as the log-odds scale.

Prior Choice for Hazard Ratios

The same principle applies to survival analysis.

If:

$$ HR=e^\beta $$

then a prior such as:

$$ \beta\sim N(0,\sigma^2) $$

induces a prior distribution for the hazard ratio.

For example, if:

$$ \beta\sim N(0,0.25^2) $$

then approximately 95% of \(\beta\) lies between:

$$ -0.49 \quad\text{and}\quad 0.49 $$

which corresponds approximately to:

$$ 0.61 < HR < 1.63 $$

This makes the scientific meaning of the prior much easier to assess.

Always inspect priors on the clinically interpreted scale. A normal prior may look harmless on a regression coefficient scale while implying implausibly extreme odds ratios or hazard ratios after transformation.

Prior Choice for Variance Parameters

Variance parameters require special care because they must be positive.

For example:

$$ \sigma>0 $$

Possible priors include:

  • Half-normal distributions
  • Half-\(t\) distributions
  • Log-normal distributions
  • Other positive-support distributions

A half-normal prior can be written:

$$ \sigma\sim\operatorname{HalfNormal}(s) $$

The scale \(s\) should be chosen based on the expected magnitude of the standard deviation on the original measurement scale.

Variance priors should not be selected solely because they are computationally convenient.

Prior Choice for Regression Coefficients

In a linear regression:

$$ Y_i = \beta_0+\beta_1X_i+\epsilon_i $$

one might use:

$$ \beta_1\sim N(0,\sigma_\beta^2) $$

But the appropriate value of \(\sigma_\beta\) depends heavily on the scale of \(X\) and \(Y\).

A prior that is reasonable when \(X\) is measured in kilograms may be inappropriate when \(X\) is measured in grams.

Standardization matters. Prior scales should be interpreted in relation to the units and scaling of the model covariates.

Standardizing Predictors Can Simplify Prior Specification

Suppose a predictor \(X\) is standardized:

$$ X^* = \frac{X-\bar X}{SD(X)} $$

Then a coefficient prior can be specified on a scale where one unit corresponds to approximately one standard deviation of the predictor.

This can make the prior easier to interpret and compare across coefficients.

However, the transformation should be prespecified appropriately and the final results should still be translated back to clinically meaningful units.

Prior Predictive Checking for Regression

Suppose:

$$ \beta\sim N(0,10^2) $$

This may appear "noninformative."

But if \(X\) ranges from 0 to 100, the implied variation in:

$$ \beta X $$

could be enormous.

The resulting prior predictive distribution for \(Y\) might imply impossible values.

Prior predictive simulation would reveal this immediately.

Prior predictive checking is often more informative than looking at the prior density alone. Ask what the prior implies for observable clinical outcomes.

Prior Choice for Logistic Regression

Logistic regression can be especially sensitive to prior specification when data are sparse.

The model is:

$$ \operatorname{logit}(p_i) = \alpha+\beta_1X_{i1}+\cdots+\beta_kX_{ik} $$

With sparse data, maximum likelihood estimates can become extremely large or even fail to exist because of complete or quasi-complete separation.

A weakly informative prior on regression coefficients can regularize the analysis.

For example:

$$ \beta_j\sim N(0,2.5^2) $$

may place meaningful probability on moderate effects while strongly downweighting extreme coefficients.

The exact scale should depend on predictor coding and the scientific context.

Prior Choice in Hierarchical Models

Hierarchical models often contain parameters describing between-group heterogeneity.

For example:

$$ \theta_j\sim N(\mu,\tau^2) $$

where \(\tau\) describes between-group variability.

A prior for \(\tau\) should reflect how much heterogeneity is plausible.

A common approach is:

$$ \tau\sim\operatorname{HalfNormal}(s) $$

or a half-\(t\) distribution.

The choice can have substantial influence when there are few groups.

Hierarchical priors deserve special attention. With only a small number of groups, the data may provide limited information about between-group heterogeneity, making the prior for \(\tau\) influential.

Prior Choice in Subgroup Analysis

Suppose a trial estimates treatment effects in several subgroups:

$$ \theta_1,\theta_2,\ldots,\theta_K $$

A hierarchical prior might allow the subgroup effects to be related:

$$ \theta_k\sim N(\mu,\tau^2) $$

This can stabilize estimates for small subgroups through partial pooling.

However, if there is strong evidence that treatment effects differ radically across subgroups, aggressive pooling may be inappropriate.

The prior structure should therefore reflect plausible biological and clinical relationships among the populations.

Prior Choice Should Consider the Estimand

The same clinical question can have different estimands.

For example, a treatment effect might be:

  • A treatment-policy estimand
  • A hypothetical estimand
  • A while-on-treatment estimand
  • A difference in response probabilities
  • A log hazard ratio
  • A restricted mean survival-time difference

The prior should be specified for the parameter associated with the actual estimand.

Prior specification cannot be separated from estimand definition. A prior for a treatment effect under one estimand is not automatically appropriate for a different treatment-effect definition.

Prior Choice in a Confirmatory Clinical Trial

Suppose a Phase III trial evaluates a treatment effect \(\theta\). A Bayesian confirmatory analysis might specify:

$$ \theta\sim N(0,\sigma^2) $$

The choice of \(\sigma\) should be justified using relevant external evidence, clinical expertise, or a prespecified weakly informative rationale.

The statistical analysis plan should explain:

  • Why the prior was selected
  • What external information informed it
  • How prior strength was determined
  • What prior predictive checks were performed
  • What sensitivity analyses were prespecified
  • How the posterior will be summarized
  • What decision threshold will be used

Prior Choice and Regulatory Transparency

In a regulated clinical development program, prior selection should be particularly transparent.

A reviewer should be able to understand:

  • Where the prior information came from
  • Why it is relevant to the current population
  • How it was translated into a mathematical distribution
  • How strongly it influences the analysis
  • Whether alternative reasonable priors change the conclusion
  • How the operating characteristics of the design were evaluated
For confirmatory use, prior selection should be treated as part of the prespecified statistical methodology rather than an analysis decision made after seeing the results.

Prior Selection Before Looking at the Current Data

One of the most important principles is temporal separation.

The prior should be specified using information that is legitimately available before the current data are used for the analysis.

Using the current trial results to choose a prior and then presenting that prior as if it were genuinely pre-data information can distort the Bayesian analysis.

For example, an analyst should not observe a surprisingly large treatment effect and then select an optimistic prior specifically because it produces a desired posterior.

Prior information should be defined independently of the current outcome data whenever the analysis is intended to represent genuine pre-data uncertainty.

Data-Dependent Priors

There are specialized Bayesian methods in which prior information can depend on other data.

For example, empirical Bayes methods estimate hyperparameters from observed data.

These methods can be useful, but they should not be casually described as ordinary fixed priors chosen independently of the data.

The inferential properties and interpretation are different.

Robust Priors

Sometimes an analyst wants to use historical information without allowing a small amount of prior-data conflict to dominate the current trial.

A robust prior can combine an informative component with a more diffuse component:

$$ \pi(\theta) = w\pi_{\text{informative}}(\theta) + (1-w)\pi_{\text{robust}}(\theta) $$

where \(0\le w\le1\).

The robust component provides protection against prior-data conflict.

This approach can be especially useful when historical information is relevant but its transportability is uncertain.

Prior-Data Conflict

A prior-data conflict occurs when the observed data are substantially inconsistent with the prior.

For example, suppose the prior strongly centers a response probability near 20%, but the current trial produces overwhelming evidence for a response rate near 60%.

The analyst should not simply ignore the discrepancy.

Possible explanations include:

  • The historical population differs from the current population
  • The endpoint changed
  • The treatment context changed
  • The prior was too strong
  • The current trial has an unusual result
  • The assumed likelihood is misspecified
Prior-data conflict is information. It should prompt investigation rather than automatic deletion of either the prior or the current data.

Prior Predictive Probability of Extreme Results

A useful diagnostic is to calculate the prior predictive probability of observing results at least as extreme as the actual study.

Conceptually:

$$ P(T(D)\ge T(D_{\text{obs}})) $$

where \(T(D)\) is a suitable test statistic or discrepancy measure.

If this probability is extremely small under the prior predictive distribution, the prior/model may be in substantial conflict with the observed data.

Prior Choice Should Be Evaluated on Multiple Scales

A prior can be examined on at least three scales:

  1. The model parameter scale
  2. The clinically interpretable effect scale
  3. The observable-data scale

For example, a prior on a logistic regression coefficient should be examined both as a distribution for \(\beta\) and as a distribution for the odds ratio.

It should also be examined through its implications for predicted response probabilities.

A prior is easier to defend when its implications are understandable to both statisticians and clinical scientists.

Prior Calibration Using Clinically Important Differences

Suppose the minimum clinically important difference is:

$$ \delta=0.10 $$

A prior can be evaluated by asking:

$$ P(\theta>\delta) $$

and:

$$ P(\theta<-\delta) $$

before observing the current data.

These probabilities provide a useful language for discussing the prior with clinical experts.

For example: "Before the trial, we believe there is approximately a 20% probability that the true treatment benefit exceeds the clinically important threshold."

This is often much easier for collaborators to evaluate than an unexplained standard deviation.

Prior Probability of Benefit

Suppose the prior is:

$$ \theta\sim N(0.05,0.10^2) $$

If the clinically important threshold is:

$$ \delta=0.10 $$

then the prior probability of clinically important benefit is:

$$ P(\theta>0.10) $$

This probability can be calculated from the normal distribution.

The important point is that the prior can be communicated directly in terms of the clinical decision.

Prior Probability of Harm

The same prior can be used to quantify concern about harm.

For example:

$$ P(\theta<-0.10) $$

might represent the prior probability that the treatment is meaningfully harmful.

A clinically useful prior therefore does not need to focus exclusively on benefit. It should represent the full range of scientifically plausible outcomes.

Prior Distributions Should Not Be Chosen by Convention Alone

There is no universal rule such as: "Always use a normal prior." or: "Always use a weakly informative prior."

The correct prior depends on:

  • The parameter
  • The parameterization
  • The scientific context
  • The amount of external information
  • The current sample size
  • The likelihood
  • The intended decision

A conventional prior can be reasonable, but convention is not a substitute for justification.

Common Prior Distributions

Parameter Common Prior Choices Reason
Probability \(p\) Beta Support restricted to 0–1
Unrestricted treatment difference Normal Real-valued support
Positive parameter Log-normal, half-normal, gamma Positive support
Odds ratio Normal prior on log(OR) Convenient unconstrained scale
Hazard ratio Normal prior on log(HR) Convenient unconstrained scale
Standard deviation Half-normal, half-\(t\) Positive support
Count rate Gamma or log-normal Positive support
Correlation Specialized correlation priors Restricted parameter space

Choosing a Distribution Family

The distribution family should respect the parameter's support.

For example, a normal prior is generally inappropriate for a probability because a normal distribution allows values below 0 and above 1.

Likewise, a standard normal prior is not directly appropriate for a variance because variance parameters must be positive.

This sounds obvious, but prior support becomes particularly important in complex hierarchical models.

Support Is Only the First Requirement

A prior can have the correct mathematical support and still be scientifically poor.

For example, both:

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

and:

$$ p\sim\operatorname{Beta}(100,100) $$

respect the support \(0

But the second prior is far more concentrated around 0.50.

Thus prior selection requires both:

  • Correct support
  • Scientifically reasonable concentration

Prior Choice for Rare Events

Rare-event endpoints require particular care.

Suppose a safety event probability is expected to be below 1%. A prior centered around 50% would be obviously inappropriate.

Instead, the prior should be calibrated to plausible rare-event probabilities. For example:

$$ p\sim\operatorname{Beta}(\alpha,\beta) $$

with a mean near the expected background rate.

The prior predictive distribution should then be examined to ensure it produces reasonable numbers of events for the planned sample size.

Prior Choice for Count Outcomes

Suppose:

$$ Y_i\sim\operatorname{Poisson}(\lambda_i) $$

where \(\lambda_i>0\) is a rate. A natural prior might be specified on the log scale:

$$ \log(\lambda)\sim N(\mu,\sigma^2) $$

This guarantees positive values for \(\lambda\).

The prior can be interpreted by transforming it back to the count-rate scale.

Prior Choice for Bioequivalence and PK Parameters

Bayesian priors can also be useful for pharmacokinetic parameters such as clearance, volume of distribution, absorption rate, or interindividual variability.

For a positive parameter such as clearance:

$$ CL>0 $$

a log-scale model may be appropriate:

$$ \log(CL)\sim N(\mu_{CL},\sigma_{CL}^2) $$

The prior can then be calibrated using historical pharmacokinetic studies, mechanistic knowledge, or previous population-PK analyses.

This illustrates a broader principle: Prior specification should respect the biological scale and the natural geometry of the parameter.

Prior Choice in Dose-Response Models

Dose-response models often contain several parameters with different meanings and constraints.

For example, an Emax model can be written:

$$ E(D) = E_0+ \frac{ E_{\max}D }{ ED_{50}+D } $$

Here:

  • \(E_0\) is the baseline response.
  • \(E_{\max}\) is the maximum treatment effect.
  • \(ED_{50}\) is the dose producing half of the maximum effect.

Each parameter may require a different prior.

For example, \(ED_{50}\) must be positive and may be naturally modeled on a log scale.

The prior should reflect plausible dose-response behavior rather than simply assigning the same distribution to every parameter.

Prior Choice for Correlations

Correlation parameters must lie between \(-1\) and \(1\).

A simple normal prior on the correlation itself can therefore be problematic.

In multivariate Bayesian models, correlation structures are often modeled using specialized distributions such as the LKJ prior.

For example:

$$ R\sim\operatorname{LKJ}(\eta) $$

where \(R\) is a correlation matrix and \(\eta\) controls the degree of concentration around particular correlation structures.

The choice of \(\eta\) should be justified in the context of the expected correlation structure.

Prior Choice for Missing Data Models

Bayesian models can incorporate uncertainty about missing-data mechanisms.

For example, parameters governing a missing-not-at-random model may require priors that encode assumptions about departures from missing-at-random behavior.

These priors can be especially influential because the missing-data mechanism may be weakly identified by the observed data.

Sensitivity to missing-data priors can be more important than sensitivity to ordinary regression-coefficient priors. When parameters are weakly identified, prior assumptions can materially determine the posterior.

Prior Choice and Weak Identification

A parameter is weakly identified when the observed data contain limited information about its value.

Examples include:

  • Rare events
  • Small samples
  • Few clusters
  • Complex hierarchical models
  • Mixture models
  • Missing-data sensitivity parameters
  • Weakly informed subgroup effects

In these situations, the prior can have a large influence.

That is not necessarily a flaw. It simply means that the data alone do not determine the answer.

When the data are weakly informative, prior assumptions should become more—not less—transparent.

Prior Robustness and the Amount of Data

Suppose the same prior is used with two studies. Study A has:

$$ n=20 $$

while Study B has:

$$ n=2000 $$

The prior may have substantial influence in Study A and relatively little influence in Study B.

Therefore, the same prior can be appropriately weak in one setting and effectively informative in another.

Why "Noninformative" Is a Dangerous Goal

Analysts sometimes attempt to remove subjectivity by choosing an extremely broad prior.

But a broad prior can create unintended consequences.

For example, it may:

  • Generate implausible predictions
  • Produce numerical instability
  • Assign too much probability to extreme effects
  • Behave differently under reparameterization
  • Influence marginal likelihoods and Bayes factors

A more useful goal is often: Choose a prior that is weakly influential where appropriate but scientifically plausible.

Prior Choice for Estimation vs. Model Comparison

An important distinction is that prior sensitivity can be particularly strong for Bayes factors and other marginal-likelihood-based quantities.

The marginal likelihood is:

$$ p(D) = \int p(D\mid\theta)\pi(\theta)\,d\theta $$

A diffuse prior spreads probability over a very large parameter space. This can reduce the marginal likelihood even when the posterior estimate is not substantially affected.

Consequently, a prior that is relatively harmless for estimation may still have a major effect on model comparison.

Prior sensitivity should be evaluated with respect to the actual inferential target. A prior may have little effect on a posterior mean but substantial effect on a Bayes factor.

Prior Choice for Decision Thresholds

Suppose the clinical decision rule is:

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

Then prior selection should be evaluated with this decision rule in mind.

For example, if different reasonable priors lead to:

$$ 0.94,\quad0.96,\quad0.98 $$

the decision could change depending on the prior.

That would indicate meaningful prior sensitivity.

If all reasonable priors produce:

$$ P(\theta>\delta\mid D) \approx0.99 $$

then the decision is much more robust.

Prior Sensitivity Should Be Prespecified When Possible

In confirmatory clinical research, prior sensitivity analyses should ideally be defined before examining the final trial results.

The statistical analysis plan might specify:

  • Primary prior
  • Skeptical prior
  • Weakly informative prior
  • Alternative historical-data prior
  • Decision criterion

This reduces the risk that prior selection becomes an outcome-dependent analysis decision.

How to Justify a Prior in a Statistical Analysis Plan

A strong prior justification should answer at least five questions.

  1. What parameter does the prior describe?
  2. Where did the prior information come from?
  3. Why is that information relevant to the current population?
  4. How was the strength of the prior determined?
  5. Would reasonable alternative priors change the conclusion?

A weak justification might say:

"We selected a normal prior with mean 0 and variance 100 because it is noninformative."

A stronger justification would describe the scientifically plausible treatment effect range, explain why that range was chosen, show the prior predictive implications, and demonstrate robustness to alternative priors.

An Example of a Strong Prior Justification

Suppose \(\theta\) is a treatment difference in response probability. A statistical analysis plan might state conceptually:

Example: "The prior distribution for the treatment effect was centered at zero because no direction of treatment benefit was assumed before the trial. The prior standard deviation was selected such that approximately 95% of prior probability lies between -0.20 and 0.20, corresponding to treatment effects larger than 20 percentage points in absolute magnitude being considered unlikely based on previous trials and clinical experience. Prior predictive simulations were performed to verify that the implied response probabilities were clinically plausible. Sensitivity analyses using both a more skeptical and a less informative prior were prespecified."

This is much more informative than simply reporting the mathematical distribution.

What a Good Prior Justification Looks Like

Weak Justification Stronger Justification
"We used a vague prior." Explains plausible parameter range and prior predictive behavior.
"The prior is standard." Explains why the standard prior is appropriate for this parameter and scale.
"Historical data were used." Describes the source, comparability, and amount of information borrowed.
"The prior was noninformative." Explains the intended influence and demonstrates sensitivity.
"The prior was chosen by an expert." Documents the elicitation process and uncertainty.

Common Mistake: Choosing the Prior After Seeing the Data

Suppose the observed treatment effect is unexpectedly large. The analyst then selects a prior centered on a large positive effect.

This can make the resulting posterior appear more convincing.

But if the prior was chosen because of the observed current data, it is not a genuine pre-data prior in the usual sense.

The analysis should instead distinguish clearly between:

  • Pre-data scientific information
  • Current observed data
  • Any empirical or data-dependent prior estimation procedure

Common Mistake: Treating Historical Data as Literal Current Data

Suppose 100 historical patients are available. It is tempting to say: "The prior is equivalent to adding 100 patients to the current trial."

That is often too simplistic.

Historical observations may differ in quality, relevance, population, endpoint definition, and study conduct.

A Bayesian model can discount or partially pool historical information rather than treating it as perfectly exchangeable with the current trial.

Common Mistake: Using a Prior That Is Too Strong

A prior may be mathematically valid but overly concentrated.

Suppose the current data strongly favor:

$$ \theta=0.30 $$

but the prior is extremely concentrated around:

$$ \theta=0 $$

The posterior may remain near zero even with meaningful current evidence.

This is not automatically wrong. But it requires strong justification for the prior information.

Strong priors are legitimate when strong prior information exists. The issue is not whether a prior is influential. The issue is whether its influence is scientifically defensible.

Common Mistake: Using a Prior That Is Too Weak

The opposite problem can occur. An extremely broad prior may imply:

  • Implausibly large effects
  • Unrealistic future outcomes
  • Numerical instability
  • Poorly behaved posterior distributions

A prior should not be made broad merely to avoid having to justify it.

Common Mistake: Ignoring the Prior on the Outcome Scale

An analyst may choose a prior that looks reasonable for \(\beta\).

But after transformation:

$$ OR=e^\beta $$

the implied odds ratios may be extreme.

Always examine the induced prior on the clinically meaningful scale.

Common Mistake: Assuming a Uniform Prior Is Neutral

A uniform prior:

$$ \theta\sim\operatorname{Uniform}(a,b) $$
may seem objective.

But uniformity is parameterization-dependent.

A prior that is uniform on \(\theta\) will not generally be uniform on \(\log(\theta)\), \(1/\theta\), or another transformation.

Therefore, the choice should be justified in terms of the scientific parameterization.

Common Mistake: Calling Every Prior "Subjective"

Prior distributions can certainly encode expert judgment.

But they can also be constructed from:

  • Published datasets
  • Meta-analysis
  • Historical controls
  • Previous randomized trials
  • Mechanistic constraints
  • Physical or biological knowledge

The important issue is not whether the prior contains judgment. The important issue is whether the judgment is explicit, reproducible, and defensible.

Common Mistake: Assuming Bayesian Results Are Independent of the Prior

Bayesian inference is conditional on the prior.

The posterior is:

$$ \pi(\theta\mid D) \propto L(\theta;D)\pi(\theta) $$

Changing \(\pi(\theta)\) can change the posterior.

The question is not whether the prior has an effect. The question is whether that effect is appropriate and whether the conclusion is robust to reasonable alternatives.

Common Mistake: Performing Sensitivity Analysis With Implausible Priors

A sensitivity analysis should represent genuine uncertainty about prior specification.

Comparing the primary prior with a wildly unrealistic distribution does not provide useful evidence of robustness.

Alternative priors should generally be scientifically plausible.

R Example: Beta Prior

A Beta prior can be specified in R as follows:

alpha_prior <- 2
beta_prior  <- 8

prior_mean <-
  alpha_prior /
  (alpha_prior + beta_prior)

prior_mean
# 0.20

The prior standard deviation is:

prior_sd <-
  sqrt(
    alpha_prior * beta_prior /
    (
      (alpha_prior + beta_prior)^2 *
      (alpha_prior + beta_prior + 1)
    )
  )

prior_sd

Plot the Prior Distribution

The prior should be visualized whenever practical.

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

prior_density <-
  dbeta(
    p,
    shape1 = 2,
    shape2 = 8
  )

plot(
  p,
  prior_density,
  type = "l",
  xlab = "Response probability",
  ylab = "Prior density",
  main = "Beta(2, 8) Prior"
)

The plot makes it immediately clear that the prior centers on lower response probabilities and assigns relatively little probability to values near 1.

Calculate Prior Quantiles

Prior quantiles can communicate uncertainty more intuitively than the Beta parameters alone.

qbeta(
  c(0.025, 0.50, 0.975),
  shape1 = 2,
  shape2 = 8
)

These values provide the lower 2.5th percentile, median, and upper 97.5th percentile of the prior.

This is often useful when documenting prior assumptions.

Calculate Prior Probability of Clinical Benefit

Suppose the clinically important response-rate threshold is 25%. The prior probability that the response rate exceeds this threshold is:

1 - pbeta(
  0.25,
  shape1 = 2,
  shape2 = 8
)

This answers:

$$ P(p>0.25) $$

before observing the current trial data.

Bayesian Updating in R

Suppose 6 responses are observed among 20 patients.

x <- 6
n <- 20

alpha_post <-
  alpha_prior + x

beta_post <-
  beta_prior + n - x

alpha_post
beta_post

The posterior distribution is:

$$ p\mid D \sim \operatorname{Beta}(8,22) $$

because:

$$ 2+6=8 $$

and:

$$ 8+(20-6)=22 $$

The posterior mean is:

posterior_mean <-
  alpha_post /
  (alpha_post + beta_post)

posterior_mean
# 0.2667

Calculate a 95% Credible Interval

qbeta(
  c(0.025, 0.975),
  shape1 = alpha_post,
  shape2 = beta_post
)

The resulting values form an equal-tail 95% credible interval for the response probability.

The interpretation is conditional on the specified likelihood and prior: Given the model, prior, and observed data, 95% of the posterior probability lies within the interval.

Prior Sensitivity in R

Suppose we want to compare three priors with the same prior mean but different strength.

priors <- list(
  weak = c(1, 4),
  moderate = c(4, 16),
  strong = c(40, 160)
)

posterior_means <- sapply(
  priors,
  function(prior) {
    a <- prior[1] + 6
    b <- prior[2] + 14

    a / (a + b)
  }
)

posterior_means

All three priors have mean 20%, but their posterior influence differs because their effective sample sizes differ.

Prior Predictive Simulation in R

Prior predictive simulation is often the easiest way to understand what a prior implies.

set.seed(123)

nsim <- 10000
n <- 50

p_prior <-
  rbeta(
    nsim,
    shape1 = 2,
    shape2 = 8
  )

x_prior_predictive <-
  rbinom(
    nsim,
    size = n,
    prob = p_prior
  )

summary(x_prior_predictive)

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

This simulation asks: Before seeing the trial data, what number of responses does the model think could plausibly occur in 50 patients?

Prior Predictive Simulation for a Treatment Effect

Suppose:

$$ \theta\sim N(0,0.10^2) $$

We can simulate treatment effects from the prior:

set.seed(123)

theta_prior <-
  rnorm(
    10000,
    mean = 0,
    sd = 0.10
  )

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

The resulting quantiles provide a direct summary of the prior's implied plausible treatment-effect range.

A Practical Prior-Selection Workflow

1
Define the estimand and identify the exact parameter for which a prior is required.
2
Determine the parameter's natural scale and mathematical support.
3
Gather relevant external information, historical data, and scientific knowledge.
4
Define clinically plausible ranges and meaningful thresholds before choosing a distribution.
5
Select a probability distribution whose support and shape can represent the scientific information.
6
Calibrate the prior parameters using quantiles, moments, historical data, or expert elicitation.
7
Inspect the prior on the parameter, clinical-effect, and observable-data scales.
8
Perform prior predictive checks to determine whether the model generates plausible data.
9
Perform sensitivity analyses using scientifically reasonable alternative priors.
10
Prespecify and document the primary prior, rationale, diagnostics, and sensitivity analyses.

A Practical Checklist for Justifying a Prior

Question What to Document
What parameter? Exact estimand and model parameter
What scale? Original, log, logit, standardized, or other scale
What distribution? Mathematical family and support
What center? Prior mean, median, or other location measure
How much uncertainty? Standard deviation, quantiles, or equivalent concentration
Where did information come from? Historical studies, literature, experts, biological knowledge, etc.
How strong is the prior? Effective sample size or other meaningful measure
Does it generate plausible data? Prior predictive assessment
Is the conclusion robust? Prior sensitivity analysis
Was it prespecified? Protocol, SAP, or other analysis documentation

Prior Choice in Clinical Trial Design

Prior distributions can affect not only the final analysis but also the design of a Bayesian clinical trial.

Suppose a trial will declare success when:

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

The probability of satisfying this criterion depends on:

  • The sample size
  • The true treatment effect
  • The likelihood
  • The prior distribution
  • The decision threshold

Therefore, the trial's operating characteristics should be evaluated under multiple plausible true parameter values and, when appropriate, multiple plausible priors.

Prior Predictive Operating Characteristics

A Bayesian trial design can be evaluated by simulation. For each simulated trial:

  1. Draw a parameter value from the assumed data-generating distribution.
  2. Generate trial data conditional on that parameter.
  3. Apply the Bayesian analysis.
  4. Calculate the posterior probability or decision criterion.
  5. Record whether the trial meets the decision rule.

Repeating this process provides an estimate of the design's operating characteristics.

Do not choose a prior simply because it produces attractive design characteristics. The prior should be scientifically defensible first. Design performance should then be evaluated under realistic scenarios.

Prior Distribution vs. Sampling Distribution

The prior and sampling distributions answer different questions.

Distribution Question
Prior \(\pi(\theta)\) What parameter values are plausible before current data?
Likelihood \(p(D\mid\theta)\) How plausible are the observed data for a given parameter?
Posterior \(\pi(\theta\mid D)\) What parameter values are plausible after combining prior and data?
Prior predictive \(p(D)\) What datasets does the model predict before seeing current data?
Posterior predictive \(p(D_{\text{new}}\mid D)\) What future data are predicted after seeing current data?

Prior, Likelihood, and Posterior: A Mental Model

A useful conceptual sequence is:

$$ \boxed{ \text{Prior} \longrightarrow \text{Update with Data} \longrightarrow \text{Posterior} } $$

The likelihood determines how strongly different parameter values are supported by the current observations.

The prior determines how plausible those parameter values were before seeing the current observations.

The posterior combines both.

Simple memory aid: Prior = what was plausible before the current data. Likelihood = what the current data say about each parameter value. Posterior = what is plausible after combining both sources of information.

How Much Should the Prior Influence the Result?

There is no universal target such as: "The prior should contribute less than 5%."

The appropriate influence depends on the amount and quality of prior information.

If strong external evidence exists, a strong prior may be appropriate. If little external information exists, a weakly informative prior may be more appropriate.

The goal is not to minimize prior influence.

The goal is to represent prior uncertainty accurately.

Strong Prior Information Is Not a Statistical Error

Suppose decades of evidence establish that a biological parameter is within a narrow range. It would be inappropriate to deliberately use an extremely broad prior simply to appear objective.

Ignoring reliable prior information is itself a modeling decision.

The appropriate question is: How much reliable information exists, and how relevant is it to the current analysis?

When a Weak Prior Is Preferable

A weakly informative prior may be preferable when:

  • Little external information exists
  • The current trial is adequately powered
  • The purpose is primarily to regularize estimation
  • Strong prior information could introduce substantial bias
  • The parameter is being estimated primarily from current data

Even then, the prior should be checked for reasonable predictive implications.

When an Informative Prior Is Preferable

An informative prior may be preferable when:

  • High-quality historical data exist
  • The current sample size is small
  • The disease population is rare
  • Biological information is strong
  • Previous studies are highly relevant
  • The goal is explicitly to borrow external information

The stronger the external information, the more important it becomes to demonstrate comparability and transportability.

Prior Choice in Rare Disease Trials

Rare disease trials often have limited sample sizes.

A purely data-driven analysis may therefore produce highly uncertain estimates.

Reliable external information may be especially valuable.

Bayesian methods can incorporate such information through:

  • Historical controls
  • Hierarchical models
  • Informative priors
  • Robust mixture priors
  • Dynamic borrowing

However, the small sample size also means prior assumptions can be highly influential.

Rare disease settings create both the strongest motivation for informative priors and the greatest need for prior sensitivity analysis.

Prior Choice in Safety Analyses

Safety endpoints can also be analyzed using Bayesian priors.

For a rare adverse event probability \(p\), the prior should reflect the background event rate and uncertainty around that rate.

If the prior is too optimistic about safety, important adverse-event probabilities may be downweighted.

For this reason, safety priors may appropriately be conservative or skeptical when the cost of underestimating risk is high.

Prior Choice and Asymmetric Consequences

Not all decisions have symmetric consequences. Suppose falsely declaring a treatment effective is much more costly than failing to recognize a modest benefit.

A skeptical prior may be scientifically and operationally appropriate.

Conversely, if the goal is early signal detection for a rare disease with very limited alternatives, a different prior may be reasonable.

The prior should therefore be considered alongside the decision problem.

Prior Choice Is Not the Same as Decision Threshold

These concepts should be kept separate.

The prior describes uncertainty about the parameter before the current data. The decision threshold describes how much posterior evidence is required for a decision.

For example:

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

contains both a clinical threshold \(\delta\) and a posterior decision threshold of 0.95.

Changing either can change the decision, but they represent different assumptions.

Prior Choice and Clinical Meaning

A statistically convenient prior is not necessarily a clinically meaningful prior.

Whenever possible, prior assumptions should be translated into quantities that clinicians can understand.

For example:

  • Expected response rate
  • Expected hazard ratio
  • Probability of clinically important benefit
  • Probability of meaningful harm
  • Expected treatment difference
  • Expected variability

This makes prior elicitation and review much more transparent.

Common Mistakes

  1. Choosing a prior solely by convention. A familiar distribution is not automatically appropriate for the parameter.
  2. Reporting only the prior mean. Prior concentration or uncertainty must also be described.
  3. Calling an extremely broad prior "noninformative." Broad priors can have important consequences and may generate implausible predictions.
  4. Ignoring parameterization. A prior on \(\beta\) induces a different prior on \(e^\beta\).
  5. Failing to perform prior predictive checks. A prior can look reasonable on the parameter scale while producing unrealistic data.
  6. Using historical data without assessing comparability. Historical information must be relevant to the current population and estimand.
  7. Choosing the prior after seeing the current data. This can compromise the intended interpretation of the prior as pre-data information.
  8. Ignoring prior sensitivity. Small or weakly identified datasets can be highly sensitive to prior choice.
  9. Using implausible alternatives in sensitivity analyses. Sensitivity analyses should examine scientifically reasonable priors.
  10. Assuming a prior is equivalent to adding literal observations. Equivalent sample size is model-dependent and should not be overinterpreted.
  11. Ignoring the consequences of the prior for decision thresholds. A small change in the posterior can matter greatly when a decision threshold is near the posterior probability.
  12. Confusing prior uncertainty with posterior uncertainty. The prior describes pre-data uncertainty; the posterior incorporates the current evidence.

What Should Be Reported for a Bayesian Prior?

A reproducible Bayesian analysis should clearly report:

  • The parameter being assigned the prior
  • The parameterization
  • The prior distribution
  • The numerical hyperparameters
  • The source of prior information
  • The rationale for the prior center
  • The rationale for the prior concentration
  • The prior predictive assessment
  • The sensitivity analyses
  • Any historical-data borrowing
  • The method used to estimate the posterior

For complex hierarchical models, the complete prior structure should be reported, not merely the prior on the primary treatment effect.

What Should Be Included in a Statistical Analysis Plan?

For a prespecified Bayesian clinical trial analysis, the SAP should ideally describe:

  • Primary estimand
  • Likelihood
  • Prior distributions
  • Prior rationale
  • Historical-data sources
  • Borrowing methodology
  • Primary posterior summary
  • Credible interval definition
  • Decision threshold
  • Prior predictive checks
  • Prior sensitivity analyses
  • Model diagnostics
  • Missing-data assumptions
  • Simulation-based operating characteristics

Primary Prior vs. Sensitivity Priors

A clinical analysis may specify one prior as the primary analysis prior.

Additional priors can be used for sensitivity analyses.

Prior Purpose
Primary prior Best representation of prespecified scientific information
Skeptical prior Tests robustness against smaller effects
Weak prior Assesses dependence on external information
Alternative historical prior Assesses dependence on historical-data source
Robust prior Protects against prior-data conflict

How to Interpret Prior Sensitivity Results

Suppose three priors produce posterior probabilities:

$$ 0.96,\quad0.97,\quad0.98 $$

The decision appears robust.

Now suppose they produce:

$$ 0.71,\quad0.91,\quad0.99 $$

The conclusion is highly prior-sensitive.

That does not automatically invalidate the analysis. It means that the available current data do not completely resolve the question independently of prior assumptions.

Prior sensitivity is itself an important scientific result. It tells the reader how much the conclusion depends on assumptions about pre-data uncertainty.

The Relationship Between Prior and Sample Size

In conjugate models, prior strength can sometimes be described using an equivalent sample size. Suppose:

$$ p\sim\operatorname{Beta}(\alpha,\beta) $$

and the trial observes \(n\) patients. The prior-to-data information ratio can be approximated by:

$$ \frac{\alpha+\beta}{n} $$

when the Beta-binomial interpretation is appropriate.

For example, if:

$$ ESS_{\text{prior}}=20 $$

and:

$$ n=100 $$

the prior represents roughly one-fifth of the sample size under the conjugate interpretation.

If instead:

$$ n=20 $$

the prior and current data have comparable nominal information.

Why Effective Sample Size Is Only an Approximation

For simple conjugate models, equivalent sample size is intuitive.

For a complex hierarchical model, however, there may be no single meaningful number representing prior information.

For example, a prior can constrain:

  • A variance component
  • A correlation matrix
  • A regression coefficient
  • A mixture weight

in ways that cannot be represented as a simple number of observations.

Use effective sample size as a communication tool, not as a universal definition of prior information.

Prior Choice and Computational Stability

Priors can also influence computational behavior.

Extremely broad priors may create regions of parameter space with enormous likelihood contributions or weak identification.

This can lead to:

  • Divergent transitions in HMC
  • Poor mixing
  • Slow convergence
  • Numerical overflow
  • Weakly identified posterior distributions

A scientifically reasonable weakly informative prior can sometimes improve both inference and computation.

Computationally convenient does not automatically mean scientifically appropriate. The prior should first be justified statistically; computational benefits are an additional consideration.

Prior Choice and MCMC Diagnostics

When posterior inference is performed using MCMC, prior specification should be considered alongside computational diagnostics.

Important diagnostics can include:

  • \(\hat R\)
  • Effective sample size
  • Trace plots
  • Posterior geometry
  • Divergent transitions
  • Energy diagnostics

Poorly chosen priors can sometimes make the posterior difficult to explore.

However, good computational diagnostics do not prove that the prior is scientifically appropriate.

Prior Predictive Checking Before MCMC

Prior predictive simulation can often be performed before fitting the actual posterior model.

This is valuable because it separates two questions:

  1. Is the prior scientifically reasonable?
  2. Can the computational algorithm successfully fit the resulting model?

A model should ideally pass both checks.

A Practical Bayesian Prior Workflow for Clinical Research

1
Define the clinical estimand and identify the parameter to be modeled.
2
Determine the natural parameter scale and support.
3
Review relevant previous trials, external data, biological knowledge, and clinical expertise.
4
Elicit plausible values, quantiles, and clinically important thresholds.
5
Select and calibrate a prior distribution.
6
Examine the prior on the clinically interpretable scale.
7
Perform prior predictive simulations.
8
Evaluate prior-data compatibility after the current data are observed.
9
Perform prespecified sensitivity analyses.
10
Document the prior and its rationale in the protocol, SAP, or analysis report.

What Makes a Prior Defensible?

A defensible prior has several characteristics.

  • Scientifically plausible: it reflects realistic parameter values.
  • Transparent: its source and construction are documented.
  • Appropriately scaled: its concentration reflects actual uncertainty.
  • Relevant: external information applies to the current estimand and population.
  • Predictively reasonable: it produces plausible observable data.
  • Robust: reasonable alternatives do not produce unexplained contradictions.
  • Prespecified when appropriate: especially for confirmatory analyses.

What Makes a Prior Weak?

A prior is not weak merely because its variance is numerically large.

A prior is better described as weakly informative when it provides limited information relative to the current data while still imposing scientifically reasonable constraints.

This distinction is important.

Weakly informative is a relative concept. The same prior may be weak in a large dataset and influential in a small dataset.

What Makes a Prior Strong?

A prior is strong when it places substantial probability mass in a relatively narrow region of the parameter space.

This may be appropriate when:

  • Reliable external evidence exists
  • The parameter is biologically constrained
  • Historical evidence is highly relevant
  • The current study is small
  • The analysis intentionally borrows external information

The important issue is not avoiding strong priors. It is justifying them.

A Final Worked Conceptual Example

Suppose a new treatment has an unknown response probability \(p\). Historical studies suggest approximately 20% response. The investigator believes that 10% to 35% are plausible values, with values outside this range possible but increasingly unlikely.

A prior might be constructed around:

$$ E(p)\approx0.20 $$

with a concentration calibrated so that the desired prior quantiles match the historical uncertainty.

Before the current trial, the analyst performs prior predictive simulation. The simulations produce plausible numbers of responders.

The analyst then evaluates three priors:

  • Primary historical-information prior
  • Weakly informative prior
  • Skeptical prior

The resulting posterior probabilities of clinically meaningful benefit are similar.

This is a strong analysis because:

  • The prior has a scientific basis.
  • The uncertainty around the prior is explicit.
  • The prior generates plausible data.
  • The conclusion is robust to reasonable alternatives.

The Most Important Principle

The most important principle in prior selection is simple:

$$ \boxed{ \text{Choose the prior to represent genuine pre-data uncertainty} } $$

Do not choose a prior because it produces the result you want. Do not choose an extremely broad prior merely to avoid discussing assumptions. Do not borrow historical information simply because it is available. Do not report only the prior mean.

Instead:

  • Define the parameter.
  • Understand its scale.
  • Identify relevant prior information.
  • Quantify uncertainty.
  • Check the implied predictions.
  • Assess sensitivity.
  • Document the rationale.

Prior Selection: A Compact Decision Framework

Situation Potential Approach
Strong, high-quality external evidence Informative prior
Some reliable external information Moderately informative prior
Little external information but known plausible range Weakly informative prior
Need to discourage implausibly large effects Weakly informative or skeptical prior
Historical information may not be fully exchangeable Robust or hierarchical borrowing
Few observations and weak identification Careful prior specification + sensitivity analysis
Model comparison using Bayes factors Pay particular attention to prior sensitivity
Complex hierarchical model Specify priors for all key hierarchical parameters

Bottom Line

Prior distributions are not an afterthought in Bayesian statistics. They are a fundamental component of the statistical model. A good prior is chosen by first understanding the parameter, its scientific scale, and the available pre-data information. The analyst should then select an appropriate distribution, calibrate its center and uncertainty, examine its implications through prior predictive simulation, and evaluate the robustness of conclusions to reasonable alternative priors. In clinical research, historical data, expert knowledge, biological constraints, and weakly informative assumptions can all provide legitimate foundations for prior specification. The strongest Bayesian analyses make those assumptions explicit, clinically interpretable, prespecified when appropriate, and transparent enough that another investigator can understand exactly how the prior contributed to the final inference.

References

Gelman, A., Carlin, J.B., Stern, H.S., Dunson, D.B., Vehtari, A. & Rubin, D.B. (2013). Bayesian Data Analysis, 3rd ed. CRC Press.
Gelman, A., Simpson, D. & Betancourt, M. (2017). The prior can often only be understood in the context of the likelihood. Entropy, 19(10), 555.
Gelman, A. (2006). Prior distributions for variance parameters in hierarchical models. Bayesian Analysis, 1(3), 515–534.
McElreath, R. (2020). Statistical Rethinking: A Bayesian Course with Examples in R and Stan, 2nd ed. CRC Press.
Spiegelhalter, D.J., Abrams, K.R. & Myles, J.P. (2004). Bayesian Approaches to Clinical Trials and Health-Care Evaluation. Wiley.
Berry, D.A. (2006). Bayesian clinical trials. Nature Reviews Drug Discovery, 5, 27–36.
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.
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.
Kass, R.E. & Wasserman, L. (1996). The selection of prior distributions by formal rules. Journal of the American Statistical Association, 91(435), 1343–1370.
Gelman, A. (2006). Prior distributions for variance parameters in hierarchical models. Bayesian Analysis, 1(3), 515–534.