Tutorials › Biostatistics › Posterior Probability vs. P-Value: A Comparison

Bayesian Inference & Hypothesis Testing

Posterior Probability vs. P-Value: A Comparison

A practical and mathematical guide to the difference between posterior probabilities and p-values, including Bayes' theorem, prior distributions, Bayesian updating, hypothesis probabilities, Bayes factors, clinical interpretation, worked examples, and implementation in R.

Intermediate 15 min read

What You'll Learn

  • What a p-value actually represents
  • What a posterior probability represents
  • Why a p-value is not the probability that the null hypothesis is true
  • How Bayes' theorem converts prior information into posterior probabilities
  • How the same clinical study can produce a p-value and a posterior probability
  • How to calculate both quantities in R and interpret them correctly

Introduction

Few concepts in clinical statistics are as frequently confused as the p-value and the posterior probability.

Both quantities are used to evaluate evidence from clinical data. Both can be expressed as probabilities. Both can be used in hypothesis-based decision making.

But they answer fundamentally different questions.

A p-value asks:

$$ P(\text{data as or more extreme}\mid H_0) $$

A posterior probability asks something much closer to the question clinicians often actually want to answer:

$$ P(H\mid\text{data}) $$

In other words, a p-value conditions on the hypothesis and asks about the data, whereas a posterior probability conditions on the observed data and asks about the hypothesis or parameter.

Key idea: A p-value is a probability about the data assuming a hypothesis is true. A posterior probability is a probability about a parameter or hypothesis after incorporating the observed data and a prior distribution.

The Fundamental Difference

Suppose a clinical trial compares an experimental treatment with a null hypothesis of no treatment effect.

The frequentist framework might calculate:

$$ p\text{-value} = P(\text{data as or more extreme}\mid H_0) $$

The Bayesian framework might instead calculate:

$$ P(H_0\mid\text{data}) $$

These expressions look similar, but their conditioning is reversed.

Quantity What It Conditions On What It Describes
P-value \(H_0\) Probability of data under the null
Posterior probability Observed data Probability of a hypothesis or parameter region after updating

The Most Important Equation in This Tutorial

The distinction can be summarized by two probabilities:

$$ P(\text{data}\mid H_0) \neq P(H_0\mid\text{data}) $$

The left side asks: How compatible are these data with the null hypothesis?

The right side asks: How plausible is the null hypothesis after seeing the data?

These are not the same probability.

Common error: A p-value of 0.03 does not mean there is a 3% probability that the null hypothesis is true.

What Exactly Is a P-Value?

Consider a null hypothesis:

$$ H_0:\theta=\theta_0 $$

Suppose the observed test statistic is \(T_{\text{obs}}\). For a one-sided test in which larger values are more extreme, the p-value is conceptually:

$$ p = P(T\ge T_{\text{obs}}\mid H_0) $$

For a two-sided test, the definition depends on the test statistic and the method used to measure extremeness, but the essential idea remains the same:

$$ p = P(\text{data at least as incompatible with }H_0 \mid H_0) $$

The p-value therefore assumes the null hypothesis is true for purposes of calculating the reference distribution.

What a P-Value Does Not Tell You

A p-value does not directly provide:

  • The probability that \(H_0\) is true
  • The probability that \(H_A\) is true
  • The probability that the treatment works
  • The probability that the observed effect will replicate
  • The probability that the result is clinically important
  • The probability that the study hypothesis is correct

For example, if:

$$ p=0.03 $$

the correct interpretation is that, under the null model and the assumptions of the test, data at least as extreme as those observed would have probability about 3%.

It is not correct to say:

Incorrect: "There is a 3% chance that the null hypothesis is true."

What Is a Posterior Probability?

Bayesian inference treats the unknown parameter as uncertain and represents that uncertainty probabilistically.

Let \(\theta\) be the parameter of interest. Before observing the data, the investigator specifies a prior distribution:

$$ \pi(\theta) $$

After observing data \(D\), Bayes' theorem gives the posterior distribution:

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

The posterior distribution describes uncertainty about the parameter after combining the prior information with the observed data.

A posterior probability can then be calculated for any region of interest. For example:

$$ P(\theta>\theta^\ast\mid D) $$

This can be interpreted directly as the posterior probability that the parameter exceeds the clinically relevant threshold \(\theta^\ast\), under the specified Bayesian model and prior.

Bayes' Theorem

For a discrete hypothesis \(H\), Bayes' theorem can be written:

$$ P(H\mid D) = \frac{P(D\mid H)P(H)} {P(D)} $$

where:

  • \(P(H)\) is the prior probability of the hypothesis.
  • \(P(D\mid H)\) is the likelihood of the observed data under the hypothesis.
  • \(P(D)\) is the marginal probability of the data.
  • \(P(H\mid D)\) is the posterior probability.

The Bayesian analysis therefore explicitly incorporates prior information.

The Direction of Conditioning

The easiest way to remember the distinction is to pay attention to what comes after the vertical bar.

Expression Meaning
\(P(D\mid H_0)\) Probability of data assuming \(H_0\)
\(P(H_0\mid D)\) Probability of \(H_0\) after observing data
\(P(\theta>0\mid D)\) Posterior probability that the treatment effect is positive
\(P(D\mid\theta=0)\) Probability of the data under zero treatment effect

This is the central conceptual distinction between the two approaches.

A Clinical Example

Suppose a Phase II study evaluates an experimental treatment using a binary response endpoint.

Assume the clinically uninteresting response rate is:

$$ p_0=0.30 $$

Suppose 12 responses are observed among 30 patients:

$$ X=12,\qquad n=30 $$

The observed response rate is:

$$ \hat p = \frac{12}{30} = 0.40 $$

Thus the observed response rate is 40%, compared with the null value of 30%.

The Frequentist Analysis

Consider the one-sided hypothesis test:

$$ H_0:p\le0.30 \qquad\text{vs.}\qquad H_A:p>0.30 $$

Using the exact binomial test, the one-sided p-value is the probability of observing at least 12 responses when the true response rate is 30%:

$$ p = P(X\ge12\mid p=0.30) $$

The exact value is approximately:

$$ \boxed{p\approx0.1593} $$

Therefore, using a conventional 5% significance threshold, the result would not be statistically significant.

Frequentist interpretation: If the true response rate were 30%, the probability of observing 12 or more responses among 30 patients would be approximately 15.9%.

The Bayesian Analysis

Now analyze the same data using a Bayesian model.

Assume a uniform prior for the response probability:

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

The binomial likelihood is:

$$ L(p\mid X) \propto p^X(1-p)^{n-X} $$

Because the beta distribution is conjugate to the binomial likelihood, the posterior distribution is:

$$ p\mid X \sim \operatorname{Beta}(1+X,1+n-X) $$

With \(X=12\) and \(n=30\):

$$ p\mid X \sim \operatorname{Beta}(13,19) $$

Posterior Probability of Clinical Benefit

Suppose the clinically relevant question is:

$$ P(p>0.30\mid X=12) $$

Under the Beta(13,19) posterior, this probability is approximately:

$$ \boxed{P(p>0.30\mid D)\approx0.8931} $$

Therefore:

$$ \boxed{P(p\le0.30\mid D)\approx0.1069} $$

The Bayesian statement is therefore quite different from the p-value.

Analysis Result Question Answered
Exact one-sided p-value 0.1593 How unusual are these data under \(p=0.30\)?
Posterior probability \(p>0.30\) 0.8931 How probable is \(p>0.30\) after observing the data?
Posterior probability \(p\le0.30\) 0.1069 How probable is the clinically uninteresting region?
These numbers should not be compared as if they were two versions of the same statistic. The p-value and posterior probability answer different questions and are based on different inferential frameworks.

Why Can They Give Different Answers?

The difference arises because Bayesian and frequentist inference use different probability statements.

The p-value evaluates the observed data relative to a null reference model.

The posterior probability updates uncertainty about the parameter using:

  • The prior distribution
  • The likelihood
  • The observed data

Consequently, two Bayesian analysts using different priors can obtain different posterior probabilities from the same data.

A conventional frequentist p-value does not depend on a prior distribution.

The Role of the Prior

The prior is one of the most important differences between Bayesian and frequentist inference.

Suppose an investigator believes before the trial that response rates near 30% are plausible. A weakly informative prior might be:

$$ p\sim\operatorname{Beta}(3,7) $$

This prior has mean:

$$ E(p)=\frac{3}{3+7}=0.30 $$

After observing 12 responses among 30 patients, the posterior becomes:

$$ p\mid D \sim \operatorname{Beta}(15,25) $$

The posterior probability that \(p>0.30\) is then different from the Beta(1,1) analysis because the prior contains information.

Prior Information Is Not "Cheating"

A common misconception is that Bayesian analysis manipulates the answer by choosing a prior.

That is not the appropriate way to think about prior distributions.

All statistical analyses make assumptions.

Bayesian analysis makes its assumptions about prior uncertainty explicit.

Frequentist analyses make assumptions through the sampling model, test statistic, reference distribution, stopping rule, and other aspects of the design.

Good Bayesian practice: The prior should be scientifically justified, documented before seeing the current data when possible, and subjected to sensitivity analyses when prior choice could materially affect the conclusion.

Posterior Probability of a Parameter Region

One major advantage of Bayesian inference is that clinically meaningful probability statements can be made directly.

For example:

$$ P(p>0.30\mid D)=0.893 $$

can be interpreted as:

Given the data and the specified Bayesian model and prior, the posterior probability that the response rate exceeds 30% is approximately 89.3%.

Similarly, an investigator might ask:

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

or:

$$ P(0.30

These questions are natural in a Bayesian framework.

Credible Intervals vs. Confidence Intervals

Posterior probabilities are closely related to Bayesian credible intervals.

A 95% posterior credible interval is a region \(C\) satisfying:

$$ P(\theta\in C\mid D)=0.95 $$

This is a probability statement about the parameter conditional on the data.

A frequentist 95% confidence interval has a different interpretation. Its repeated-sampling property is:

$$ P_\theta(\theta\in C(X))=0.95 $$

where the probability refers to the procedure over repeated samples.

Important distinction: A Bayesian 95% credible interval can be described probabilistically as containing the parameter with 95% posterior probability, given the model and prior. A frequentist 95% confidence interval does not have that same parameter-probability interpretation.

A Second Example: Stronger Evidence

Now suppose 18 responses are observed among 30 patients.

$$ X=18,\qquad n=30 $$

The observed response rate is:

$$ \hat p=\frac{18}{30}=0.60 $$

Continue to use:

$$ H_0:p\le0.30 $$

The exact one-sided binomial p-value is approximately:

$$ \boxed{p\approx0.000626} $$

Using the same uniform prior:

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

the posterior is:

$$ p\mid D \sim \operatorname{Beta}(19,13) $$

The posterior probability that the response rate exceeds 30% is approximately:

$$ \boxed{P(p>0.30\mid D)\approx0.99970} $$

Again, the two numbers are answering different questions.

Quantity Value
Observed response rate 60%
Exact one-sided p-value 0.000626
Posterior probability \(p>0.30\) 0.99970
Posterior probability \(p\le0.30\) 0.00030

Why a Small P-Value Does Not Mean "Small Probability of the Null"

The classic error becomes clearer mathematically. Suppose someone observes:

$$ p=0.001 $$

They might incorrectly conclude:

$$ P(H_0\mid D)=0.001 $$

But the p-value actually describes:

$$ P(D^\ast\ge D_{\text{obs}}\mid H_0) = 0.001 $$

The desired posterior probability would instead be:

$$ P(H_0\mid D) $$

To calculate that probability, one needs additional information about the prior probability of \(H_0\) and the competing hypotheses or parameter values.

Bayes' Theorem Shows Why the Difference Matters

For two competing hypotheses \(H_0\) and \(H_A\):

$$ P(H_0\mid D) = \frac{ P(D\mid H_0)P(H_0) }{ P(D\mid H_0)P(H_0) + P(D\mid H_A)P(H_A) } $$

Notice that the posterior probability depends on:

  • The probability of the data under \(H_0\)
  • The probability of the data under \(H_A\)
  • The prior probability of \(H_0\)
  • The prior probability of \(H_A\)

A p-value uses the first component but does not by itself provide the other ingredients needed to calculate a posterior probability.

Prior Odds and Posterior Odds

Bayesian updating can be expressed especially elegantly using odds.

The prior odds are:

$$ \frac{P(H_A)}{P(H_0)} $$

The posterior odds are:

$$ \frac{P(H_A\mid D)}{P(H_0\mid D)} $$

Bayes' theorem gives:

$$ \frac{P(H_A\mid D)}{P(H_0\mid D)} = \frac{P(D\mid H_A)}{P(D\mid H_0)} \times \frac{P(H_A)}{P(H_0)} $$

The first ratio is the Bayes factor:

$$ BF_{A0} = \frac{P(D\mid H_A)} {P(D\mid H_0)} $$

Therefore:

$$ \text{Posterior Odds} = \text{Bayes Factor} \times \text{Prior Odds} $$

Bayes Factors vs. P-Values

A Bayes factor and a p-value are also not interchangeable.

A Bayes factor compares how well two hypotheses predict the observed data:

$$ BF_{A0} = \frac{P(D\mid H_A)} {P(D\mid H_0)} $$

A p-value instead measures the extremeness of the data under a null reference distribution.

Quantity Primary Question
P-value How surprising are the data under \(H_0\)?
Bayes factor How much more does the data support one hypothesis than another?
Posterior probability What is the probability of a hypothesis after updating?

A Simple Discrete-Hypothesis Example

Suppose a trial considers two competing hypotheses:

$$ H_0:p=0.30 $$
and:

$$ H_A:p=0.50 $$

Suppose 18 responses are observed among 30 patients.

The likelihood under \(H_0\) is:

$$ P(X=18\mid H_0) = {30\choose18} (0.30)^{18}(0.70)^{12} $$

The likelihood under \(H_A\) is:

$$ P(X=18\mid H_A) = {30\choose18} (0.50)^{18}(0.50)^{12} $$

The combinatorial factor cancels when calculating the Bayes factor:

$$ BF_{A0} = \frac{ (0.50)^{18}(0.50)^{12} }{ (0.30)^{18}(0.70)^{12} } $$

The resulting Bayes factor quantifies how much the observed response count changes the relative evidence between these two specific hypotheses.

Important: A Bayes factor compares hypotheses. A posterior probability additionally requires prior probabilities for those hypotheses.

Posterior Probability Depends on the Prior

Suppose two analysts use the same data but different prior beliefs.

Analyst A assigns:

$$ P(H_A)=0.50 $$

while Analyst B assigns:

$$ P(H_A)=0.10 $$

Even if both analysts calculate exactly the same Bayes factor, their posterior probabilities will differ because their prior odds differ.

This is not a computational error. It is a mathematical consequence of Bayesian inference.

Why the Prior Matters Especially With Small Samples

When the sample size is small, the likelihood may contain limited information.

The prior can therefore contribute substantially to the posterior distribution.

As the amount of information in the data increases, the likelihood often dominates a reasonable prior.

Practical implication: Prior sensitivity analysis is particularly important when the sample size is small, the endpoint is noisy, or the prior is informative.

Posterior Probability of Treatment Benefit

In clinical research, one of the most useful Bayesian quantities is often:

$$ P(\theta>\theta_{\text{clinically relevant}}\mid D) $$

For example, if \(\theta\) is a treatment effect and zero represents no benefit:

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

Alternatively, if the clinically meaningful threshold is \(\delta\):

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

This allows the statistical analysis to focus directly on a clinically meaningful region.

Clinical Significance vs. Statistical Significance

Neither a p-value nor a posterior probability automatically establishes clinical importance.

Suppose a treatment effect is:

$$ \hat\theta=0.02 $$

with a very small p-value because the study is extremely large.

The result may be statistically convincing but clinically unimportant.

Bayesian analysis can explicitly address this distinction by calculating:

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

where \(\delta\) is a clinically meaningful threshold.

Key principle: Statistical evidence and clinical relevance are separate questions. A good analysis should make both explicit.

A Useful Mental Model

1
P-value: Start with the null hypothesis and ask how unusual the observed data are.
2
Bayesian analysis: Start with prior uncertainty about the parameter or hypotheses.
3
Update: Combine the prior with the likelihood from the observed data.
4
Posterior: Describe uncertainty about the parameter after seeing the data.

Posterior Probability Is Not Just a "Better P-Value"

It is tempting to think of Bayesian posterior probabilities as simply a more intuitive replacement for p-values.

That is not quite correct.

The two approaches are based on different inferential philosophies.

A p-value is a frequentist tail probability derived from a sampling distribution.

A posterior probability is derived from a posterior distribution that combines a prior distribution with the likelihood.

They therefore answer different questions and can legitimately lead to different numerical summaries of the same evidence.

Why P-Values Can Be Difficult to Interpret

The p-value requires conditioning on a null hypothesis that may be scientifically implausible or unrealistic.

For example:

$$ H_0:\theta=0 $$

A p-value asks how unusual the data would be if \(\theta\) were exactly zero.

But the scientific question may instead be:

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

or:

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

The Bayesian framework can answer these questions directly.

Why Posterior Probabilities Can Also Be Misused

Bayesian inference is not immune to statistical mistakes.

Common problems include:

  • Using an unjustified prior
  • Ignoring prior sensitivity
  • Choosing a prior after seeing the data
  • Interpreting a posterior probability without specifying the model
  • Ignoring model misspecification
  • Treating a posterior probability as objective and prior-free
  • Using an arbitrary clinical threshold

A posterior probability is conditional on the entire Bayesian model.

Important: A statement such as "there is a 95% probability that the treatment works" is incomplete unless the underlying parameter, model, prior, data, and meaning of "works" are defined.

Posterior Probability With a Beta-Binomial Model

For a binary endpoint, one of the simplest Bayesian models is:

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

with a beta prior:

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

The posterior is:

$$ p\mid X \sim \operatorname{Beta}(a+X,b+n-X) $$

This conjugate relationship makes the beta-binomial model especially useful for teaching and simple clinical examples.

Posterior Mean

The posterior mean is:

$$ E[p\mid X] = \frac{a+X}{a+b+n} $$

This combines the prior information and observed responses.

For the uniform prior \(Beta(1,1)\) and \(X=12,n=30\):

$$ E[p\mid D] = \frac{13}{32} = 0.40625 $$

Thus the posterior mean is approximately 40.6%.

The observed response rate was exactly 40%. The slight difference arises because the Bayesian posterior incorporates the prior.

Posterior Probability Above a Threshold

The posterior probability that \(p\) exceeds a clinically meaningful threshold \(p^\ast\) is:

$$ P(p>p^\ast\mid D) = 1-F_{\operatorname{Beta}} (p^\ast;a+X,b+n-X) $$

where \(F_{\operatorname{Beta}}\) is the beta cumulative distribution function.

For example:

$$ P(p>0.30\mid X=12,n=30) = 1-F_{\operatorname{Beta}}(0.30;13,19) $$

which gives approximately:

$$ \boxed{0.8931} $$

R Implementation: Exact P-Value

The exact one-sided binomial p-value can be calculated with binom.test():

x <- 12
n <- 30
p0 <- 0.30

binom.test(
  x = x,
  n = n,
  p = p0,
  alternative = "greater"
)

The resulting p-value is approximately:

p-value
# approximately 0.1593

R Implementation: Posterior Distribution

With a Beta(1,1) prior:

x <- 12
n <- 30

a <- 1
b <- 1

post_a <- a + x
post_b <- b + n - x

post_a
post_b

The posterior is:

# Posterior:
# Beta(13, 19)

R Implementation: Posterior Probability

The posterior probability that the response rate exceeds 30% is:

p_clinical <- 0.30

posterior_prob <-
  1 - pbeta(
    p_clinical,
    shape1 = post_a,
    shape2 = post_b
  )

posterior_prob

which gives approximately:

# approximately 0.8931

R Implementation: Posterior Probability of the Null Region

The posterior probability that the response rate is at or below 30% is:

posterior_null <-
  pbeta(
    p_clinical,
    shape1 = post_a,
    shape2 = post_b
  )

posterior_null

This gives approximately:

# approximately 0.1069

R Implementation: Posterior Mean

posterior_mean <-
  post_a / (post_a + post_b)

posterior_mean

# approximately 0.40625

R Implementation: Credible Interval

A 95% equal-tail posterior credible interval can be calculated using the beta quantile function:

qbeta(
  c(0.025, 0.975),
  shape1 = post_a,
  shape2 = post_b
)

This interval is a direct posterior probability statement about the response rate under the specified Bayesian model.

Frequentist Confidence Interval for the Same Data

The frequentist analysis can also produce a confidence interval:

binom.test(
  x = 12,
  n = 30,
  p = 0.30,
  alternative = "greater",
  conf.level = 0.95
)

The confidence interval and Bayesian credible interval may be numerically similar, particularly with weakly informative priors and moderate sample sizes.

However, their interpretations remain different.

Confidence Interval vs. Credible Interval

Interval Interpretation
95% Confidence Interval A procedure constructed to contain the true parameter in 95% of repeated samples
95% Credible Interval A posterior region containing 95% posterior probability

The numerical endpoints can sometimes be nearly identical even though the underlying inferential interpretation differs.

A Key Numerical Insight

In the worked example:

$$ p\text{-value}=0.1593 $$

while:

$$ P(p>0.30\mid D)=0.8931 $$

At first glance, these numbers may appear contradictory. They are not.

The p-value asks:

$$ P(X\ge12\mid p=0.30) $$

The posterior probability asks:

$$ P(p>0.30\mid X=12) $$

The conditioning is reversed.

The Inverse Probability Fallacy

The error of interpreting a p-value as the probability that the null hypothesis is true is sometimes described as an inverse probability fallacy.

The incorrect reasoning is:

$$ P(D\mid H_0) \Rightarrow P(H_0\mid D) $$

But this inference is not valid without applying Bayes' theorem and specifying the competing hypotheses and their prior probabilities.

Remember: Conditional probabilities generally cannot be reversed simply because the same two events appear in both expressions.

Why Prior Probability Can Change the Interpretation

Suppose a new treatment is being tested for a rare disease. If the prior probability that a genuinely effective treatment exists is very low, even a statistically unusual result under the null may not imply a high posterior probability of efficacy.

Conversely, if substantial prior evidence already supports efficacy, the same data may produce a considerably higher posterior probability.

This is a direct consequence of:

$$ \text{Posterior Odds} = \text{Bayes Factor} \times \text{Prior Odds} $$

Statistical Significance Is Not the Same as Evidence Probability

A conventional significance threshold such as:

$$ \alpha=0.05 $$
creates a decision rule for a frequentist testing procedure.

It does not mean:

$$ P(H_0\mid D)=0.05 $$

Likewise, a posterior probability of 0.95 is not automatically equivalent to a p-value of 0.05.

There are mathematical relationships between certain Bayesian and frequentist procedures under specific assumptions, but there is no general conversion:

$$ p=0.05 \not\Rightarrow P(H_0\mid D)=0.05 $$

Posterior Probability and Decision Making

Bayesian clinical decisions can be built around posterior probabilities.

For example, a development program might define a decision rule:

$$ P(p>0.30\mid D)>0.90 \quad\Rightarrow\quad \text{consider further development} $$

This is a Bayesian decision rule.

It directly expresses the probability that the parameter exceeds the clinically relevant threshold.

A frequentist trial might instead use:

$$ p<0.05 \quad\Rightarrow\quad \text{reject }H_0 $$

These are different decision frameworks.

Posterior Probability Can Incorporate Clinical Thresholds

Suppose a treatment effect of less than 5 percentage points is considered clinically unimportant. Then the relevant Bayesian question might be:

$$ P(\theta>0.05\mid D) $$

rather than simply:

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

This distinction is particularly important in large clinical trials, where tiny effects can produce very small p-values.

Multiple Testing and Posterior Probabilities

Bayesian posterior probabilities do not automatically solve every issue associated with multiple comparisons.

If many hypotheses, endpoints, subgroups, or treatments are analyzed, the Bayesian model should account for the multiplicity structure when appropriate.

Hierarchical models, shrinkage priors, multiplicity-adjusted decision criteria, and other approaches may be appropriate depending on the scientific question.

The important point is: switching from p-values to posterior probabilities does not eliminate the need for careful statistical design.

Posterior Probability in Adaptive Clinical Trials

Posterior probabilities are particularly useful in Bayesian adaptive designs.

For example, a trial might monitor:

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

and use prespecified thresholds for:

  • Stopping for futility
  • Stopping for efficacy
  • Increasing or decreasing allocation probabilities
  • Expanding a treatment arm
  • Dropping an ineffective treatment

The posterior distribution is updated as new data become available.

Important: A Bayesian adaptive trial still requires prespecified operating characteristics, simulation, and careful control of decision rules. Bayesian does not mean "anything can be changed after seeing the data."

Posterior Predictive Probability

Another Bayesian quantity commonly encountered in clinical development is the posterior predictive probability.

This is different from a posterior probability about the parameter.

For example, investigators might ask:

$$ P(\text{trial will ultimately succeed}\mid D) $$

This probability integrates uncertainty in the parameter with uncertainty in future observations.

It can therefore be useful for interim decision making.

Three Probabilities That Should Not Be Confused

Quantity Expression Question
P-value \(P(D^\ast\ge D\mid H_0)\) How unusual are data this extreme under the null?
Posterior probability \(P(H\mid D)\) How probable is the hypothesis after observing the data?
Posterior predictive probability \(P(D_{\text{future}}\mid D)\) What is the probability of future outcomes given current data?

These quantities can all be useful, but they answer different questions.

Common Misconceptions

  1. "A p-value is the probability that the null hypothesis is true."

    False. A p-value is calculated under a null hypothesis and describes the extremeness of the observed data under that hypothesis.

  2. "A p-value of 0.03 means there is a 97% chance the treatment works."

    False. That is not what the p-value represents.

  3. "A posterior probability is completely objective."

    Not necessarily. It depends on the specified Bayesian model and prior.

  4. "A 95% confidence interval means there is a 95% probability that the parameter is inside it."

    That is not the standard frequentist interpretation.

  5. "A 95% credible interval and a 95% confidence interval mean the same thing."

    They can have similar numerical endpoints but have different probability interpretations.

  6. "Bayesian and frequentist methods must give the same conclusion."

    Not necessarily. They can produce different numerical evidence summaries and different decisions, particularly with informative priors or small samples.

  7. "A Bayesian analysis does not need hypothesis testing."

    Bayesian analyses can use explicit posterior probability decision rules, Bayes factors, posterior predictive probabilities, or other criteria.

  8. "A small p-value proves the treatment is clinically important."

    Statistical evidence does not establish clinical relevance.

A Practical Comparison

Feature Frequentist P-Value Bayesian Posterior Probability
Uses a prior distribution? No Yes
Probability assigned to parameter? Not in the usual frequentist interpretation Yes
Probability of null hypothesis? Not given by the p-value Can be calculated under an appropriate Bayesian model
Depends on observed data? Yes Yes
Uses a sampling distribution? Yes Yes, through the likelihood/model
Directly supports \(P(\theta>\delta\mid D)\)? No Yes
Can incorporate prior evidence? Not through the p-value itself Yes
Can be used for adaptive decisions? Yes, with appropriate design Yes

When a P-Value Is Useful

P-values remain useful for many statistical purposes.

They can provide a standardized measure of how inconsistent observed data are with a specified null model.

They are particularly familiar in frequentist clinical trial designs, where type I error, power, confidence intervals, and prespecified testing procedures form an integrated framework.

The key is to interpret the p-value correctly.

When Posterior Probabilities Are Useful

Posterior probabilities are especially useful when the scientific question is naturally probabilistic.

Examples include:

  • Probability that a treatment effect exceeds a clinical threshold
  • Probability that a response rate exceeds a target
  • Probability that one treatment is better than another
  • Probability that a parameter lies within a clinically important range
  • Probability of future trial success
  • Bayesian adaptive trial decisions

A Bayesian Clinical Trial Statement

Instead of reporting only:

$$ p=0.03 $$

a Bayesian analysis might report:

$$ P(\theta>\delta\mid D)=0.96 $$

along with:

  • The prior distribution
  • The posterior distribution
  • A credible interval
  • The clinically relevant threshold \(\delta\)
  • A sensitivity analysis for alternative priors

This provides a direct probability statement about the parameter of interest.

How to Report a P-Value Correctly

A good frequentist report should connect the p-value to the hypothesis being tested.

For example:

Example: "Under the null hypothesis that the response probability is 30%, the exact one-sided binomial p-value was 0.159."

Avoid saying:

Incorrect: "There is a 15.9% probability that the null hypothesis is true."

How to Report a Posterior Probability Correctly

A Bayesian report should identify the prior and the parameter region.

Example: "Using a Beta(1,1) prior, the posterior probability that the response rate exceeded 30% was 89.3%."

This statement is interpretable because it identifies:

  • The parameter
  • The threshold
  • The posterior probability
  • The prior distribution

Sensitivity Analysis for the Prior

A strong Bayesian clinical analysis should often examine how much the posterior depends on reasonable alternative priors.

For example:

x <- 12
n <- 30
threshold <- 0.30

# Weak / uniform prior
post1 <- pbeta(
  threshold,
  1 + x,
  1 + n - x
)

# Moderately informative prior
post2 <- pbeta(
  threshold,
  3 + x,
  7 + n - x
)

# More concentrated prior centered near 30%
post3 <- pbeta(
  threshold,
  30 + x,
  70 + n - x
)

c(
  uniform = 1 - post1,
  moderate = 1 - post2,
  concentrated = 1 - post3
)

The exact posterior probabilities will differ because the priors encode different amounts of prior information.

If the substantive conclusion is unchanged across scientifically reasonable priors, the result is more robust to prior specification.

Posterior Probability and the Null Hypothesis

There is one additional subtlety.

In a continuous parameter model, the probability of an exact point such as:

$$ P(p=0.30\mid D) $$
is generally zero when the posterior is continuous.

Therefore, Bayesian analyses often formulate clinically meaningful regions:

$$ P(p\le0.30\mid D) $$
or:

$$ P(p>0.30\mid D) $$

Alternatively, a point null can be assigned positive prior probability using a mixture prior.

Point Nulls and Mixture Priors

Suppose an investigator explicitly wants to assign positive probability to:

$$ H_0:p=0.30 $$

A mixture prior can be constructed:

$$ \pi(p) = w\,\delta_{0.30}(p) + (1-w)g(p) $$

where \(\delta_{0.30}\) represents a point mass at 0.30 and \(g(p)\) represents a continuous alternative distribution.

This permits direct calculation of:

$$ P(H_0\mid D) $$

rather than merely calculating the posterior probability of a region such as \(p\le0.30\).

Why This Distinction Matters in Clinical Development

Clinical development decisions frequently involve questions such as:

  • Is the treatment sufficiently effective to continue?
  • What is the probability that the treatment provides clinically meaningful benefit?
  • What is the probability that the response rate exceeds a target?
  • How likely is the next study to succeed?
  • How much uncertainty remains about the treatment effect?

These questions are not always naturally answered by a p-value.

Bayesian posterior and predictive probabilities can provide direct answers, provided the Bayesian model and prior are appropriate.

The Same Data Can Support Multiple Summaries

For the example with 12 responses among 30 patients:

Summary Value
Observed response rate 40%
Exact one-sided p-value vs. 30% 15.9%
Posterior mean with Beta(1,1) prior 40.6%
Posterior probability \(p>30\%\) 89.3%
Posterior probability \(p\le30\%\) 10.7%

None of these quantities is "the correct number" independent of context. They summarize different aspects of the evidence.

A Useful Clinical Interpretation Framework

1
Estimate the effect: What was observed?
2
Assess uncertainty: How precise is the estimate?
3
Assess statistical evidence: How incompatible are the data with the null model?
4
Assess clinical relevance: Is the effect large enough to matter?
5
For Bayesian inference: What is the posterior probability of the clinically relevant region?

A Direct Side-by-Side Example

Question Frequentist Bayesian
Is the observed response rate unusual if \(p=0.30\)? Use a p-value Not the primary Bayesian question
What is the probability \(p>0.30\) after seeing the data? Not directly answered Use posterior probability
What parameter values remain plausible? Confidence interval / estimation Posterior distribution / credible interval
How much does data favor one hypothesis over another? Can use testing procedures Bayes factor
What is the probability of future success? Requires predictive calculations Posterior predictive probability

Common Reporting Errors

  1. Calling the p-value the probability of the null.
  2. Calling a 95% confidence interval a 95% probability interval.
  3. Reporting a posterior probability without specifying the prior.
  4. Ignoring prior sensitivity.
  5. Equating posterior probability with clinical importance.
  6. Assuming a p-value and posterior probability should numerically agree.
  7. Converting a p-value directly into a posterior probability.
  8. Using "probability the treatment works" without defining what "works" means.

Practical Workflow for Clinical Statisticians

1
Define the clinical parameter of interest.
2
Define what constitutes a clinically meaningful effect.
3
If using frequentist inference, specify the null hypothesis and testing procedure.
4
If using Bayesian inference, specify and justify the prior distribution.
5
Fit the appropriate likelihood/model.
6
Calculate the p-value, posterior distribution, or both as appropriate.
7
Do not interpret \(P(D\mid H_0)\) as \(P(H_0\mid D)\).
8
Evaluate clinical relevance separately from statistical evidence.
9
For Bayesian analyses, conduct prior sensitivity analyses when appropriate.
10
Report the inferential framework clearly so readers know what the probability means.

The Most Important Equations to Remember

The p-value is fundamentally:

$$ \boxed{ p\text{-value} = P(\text{data as or more extreme}\mid H_0) } $$

A posterior probability is:

$$ \boxed{ P(H\mid D) } $$

For a parameter region:

$$ \boxed{ P(\theta>\delta\mid D) } $$

Bayes' theorem connects prior and posterior probabilities:

$$ \boxed{ P(H\mid D) = \frac{P(D\mid H)P(H)} {P(D)} } $$

And the odds formulation is:

$$ \boxed{ \text{Posterior Odds} = \text{Bayes Factor} \times \text{Prior Odds} } $$

Final Comparison

Concept P-Value Posterior Probability
Basic form \(P(D\mid H_0)\) \(P(H\mid D)\)
Probability is about Data Parameter or hypothesis
Conditions on Null hypothesis Observed data
Uses prior information? No Yes
Can directly say \(P(\theta>\delta\mid D)\)? No Yes
Can directly say \(P(H_0\mid D)\)? No Yes, with an appropriate Bayesian model
Typical interpretation How unusual are these data under \(H_0\)? How probable is the parameter/hypothesis after observing the data?

Worked Example Summary

For 12 responses among 30 patients with a null response rate of 30%:

Component Result
Observed response rate 40%
Null response rate 30%
Exact one-sided p-value 0.1593
Prior Beta(1,1)
Posterior distribution Beta(13,19)
Posterior mean 40.625%
Posterior probability \(p>30\%\) 89.31%
Posterior probability \(p\le30\%\) 10.69%

The p-value and posterior probability are not contradictory. They simply answer different questions.

The Most Important Concept

The single most important concept in this tutorial is:

$$ \boxed{ P(D\mid H_0) \neq P(H_0\mid D) } $$

The first quantity is related to the logic of a p-value. The second is a posterior probability.

Confusing these two probabilities is one of the most common errors in statistical interpretation.

If a clinical statistician remembers only one sentence from this tutorial, it should be:

Bottom line: A p-value tells you how unusual the observed data would be under a specified null hypothesis. A posterior probability tells you how probable a parameter value or hypothesis is after combining the observed data with a prior distribution. The p-value is a probability about the data conditional on the hypothesis; the posterior probability is a probability about the hypothesis or parameter conditional on the data. They are therefore not interchangeable.

References

Bayes, T. (1763). An Essay towards solving a Problem in the Doctrine of Chances. Philosophical Transactions of the Royal Society of London, 53, 370–418.
Gelman, A., Carlin, J.B., Stern, H.S., Dunson, D.B., Vehtari, A. & Rubin, D.B. (2013). Bayesian Data Analysis, 3rd ed. CRC Press.
Berry, S.M., Carlin, B.P., Lee, J.J. & Muller, P. (2010). Bayesian Adaptive Methods for Clinical Trials. CRC Press.
Spiegelhalter, D.J., Abrams, K.R. & Myles, J.P. (2004). Bayesian Approaches to Clinical Trials and Health-Care Evaluation. Wiley.
Greenland, S. (2019). Valid P-values behave exactly as they should: Some misleading criticisms of P-values and their resolution. American Statistician, 73(sup1), 106–114.
Wasserstein, R.L., Schirm, A.L. & Lazar, N.A. (2019). Moving to a world beyond "p < 0.05". American Statistician, 73(sup1), 1–19.
Goodman, S.N. (1999). Toward evidence-based medical statistics. 1: The P value fallacy. Annals of Internal Medicine, 130(12), 995–1004.
Goodman, S.N. (1999). Toward evidence-based medical statistics. 2: The Bayes factor. Annals of Internal Medicine, 130(12), 1005–1013.