Bayesian vs Frequentist: Two Meanings of the Word Probability

6 minute read

Published:

TL;DR: Frequentists treat the parameter as a fixed unknown and the data as random, so probability statements attach to procedures. Bayesians treat the parameter as uncertain and the data as given, so probability statements attach to the parameter itself. That is why a credible interval can say "95% probability the parameter is in here" and a confidence interval cannot. With plenty of data and a mild prior the two intervals nearly coincide; with little data, or a real prior, they genuinely differ.

The disagreement

To a frequentist, \(\mathbb{P}(A)\) is the limiting relative frequency of \(A\) in repeated trials. A coin’s bias \(\theta\) is a fixed physical constant, not the outcome of a trial, so “\(\mathbb{P}(\theta > 0.5)\)” is not a meaningful expression — it is either true or false. Randomness lives entirely in the sampling, which is why every frequentist guarantee is a statement about what would happen across hypothetical repeated experiments.

To a Bayesian, \(\mathbb{P}(A)\) is a degree of belief obeying the probability axioms. Uncertainty about a fixed constant is still uncertainty, so \(\theta\) gets a distribution — before the data (the prior) and after (the posterior). The data, once observed, are not random at all; they are what you conditioned on.

Bayes’ theorem as inference

\[ \underbrace{p(\theta \mid x)}_{\text{posterior}} = \frac{\overbrace{p(x \mid \theta)}^{\text{likelihood}}\ \overbrace{p(\theta)}^{\text{prior}}}{\underbrace{p(x)}_{\text{evidence}}}, \qquad p(x) = \int p(x\mid\theta)\,p(\theta)\,d\theta \]

The denominator does not depend on \(\theta\), so for inference it is just the constant that makes the posterior integrate to one, and the working form is \(p(\theta\mid x) \propto p(x\mid\theta)\,p(\theta)\). The likelihood is the same object maximum likelihood maximises; the difference is that Bayes multiplies it by a prior and keeps the whole resulting function rather than its peak.

Beta–Binomial, worked

A conjugate prior is one that leaves the posterior in the same family, so the update is arithmetic on parameters rather than an integral. For binomial data the conjugate prior is the Beta.

Take \(\theta \sim \mathrm{Beta}(2,2)\) — symmetric, gently favouring values near \(0.5\), worth as much as two prior successes and two prior failures. Observe \(s = 7\) successes in \(n = 10\) trials. Since \(p(\theta) \propto \theta^{1}(1-\theta)^{1}\) and \(p(x\mid\theta)\propto\theta^{7}(1-\theta)^{3}\), the product is \(\theta^{8}(1-\theta)^{4}\):

\[ \theta \mid x \;\sim\; \mathrm{Beta}(\alpha + s,\ \beta + n - s) = \mathrm{Beta}(9,\ 5) \]

The posterior mean is \(9/14 = 0.643\), sitting between the prior mean \(0.5\) and the MLE \(0.7\). It is exactly a weighted average, with weights given by the counts:

\[ \frac{\alpha+\beta}{\alpha+\beta+n}\cdot 0.5 \;+\; \frac{n}{\alpha+\beta+n}\cdot 0.7 = \tfrac{4}{14}(0.5) + \tfrac{10}{14}(0.7) = 0.643 \]

The prior carries \(\alpha+\beta = 4\) pseudo-observations against the data’s 10, so it holds 29% of the weight. Double the data to 20 trials and the prior’s share falls to 17%; the prior is a fixed amount of evidence that real data eventually swamp.

Beta prior, likelihood and posterior for seven successes in ten trialsThree densities over theta from 0 to 1. The Beta(2,2) prior is a low broad hump peaking at 0.5 with height 1.5. The normalised likelihood peaks near 0.7. The Beta(9,5) posterior is the tallest and narrowest, peaking at 0.667 with height 3.1, and lies between the other two but much closer to the likelihood. A bar underneath marks the 95 percent equal-tailed credible interval from 0.386 to 0.861. 95% credible interval [0.386, 0.861] 00.20.40.60.81.0 prior Beta(2,2) likelihood posterior Beta(9,5) θ, the success probability — 7 successes in 10 trials
Notice that the posterior is taller and narrower than either input. Multiplying two densities concentrates: the posterior is more certain than the prior and than the data alone, because it uses both.

Credible versus confidence

The 95% equal-tailed credible interval for \(\mathrm{Beta}(9,5)\) runs from 0.386 to 0.861, and it means what people wish a confidence interval meant: given this prior and this data, there is a 95% probability that \(\theta\) lies in that range. You can also read off \(\mathbb{P}(\theta > 0.5 \mid x) = 0.867\) directly, a question a p-value cannot answer.

Method95% intervalWhat it claims
Credible, \(\mathrm{Beta}(9,5)\)[0.386, 0.861]95% posterior probability \(\theta\) is inside
Clopper–Pearson (exact)[0.348, 0.933]≥95% of such intervals cover \(\theta\)
Wilson score[0.397, 0.892]≈95% coverage
Wald, \(\hat p \pm 1.96\sqrt{\hat p(1-\hat p)/n}\)[0.416, 0.984]nominally 95%, actually much less
Interview trap — the two intervals answer different questions: a credible interval is a statement about \(\theta\) given the data and the prior. A confidence interval is a statement about the procedure's success rate over repeated samples; see confidence intervals. They are not interchangeable and they are not generally equal. Note also the Wald row above: it is the interval most people write down by reflex, and at \(n=10\) it reaches 0.984 — its coverage is far below the advertised 95% for \(p\) near 0 or 1, as Brown, Cai and DasGupta (2001) document in detail.

When it matters and when it does not

With a large sample and a prior that is not sharply informative, the Bernstein–von Mises theorem says the posterior converges to a Gaussian centred at the MLE with variance given by the inverse Fisher information — the same object the frequentist asymptotics produce. At \(n = 10\) the credible and Clopper–Pearson intervals above differ by a lot; at \(n = 1000\) they would agree to three digits. The choice is then a matter of interpretation, not of numbers.

The genuine divergences are: small samples, where the prior does real work; problems with sequential or optional stopping, where frequentist error rates depend on the stopping rule and the posterior does not; hierarchical models, where partial pooling is natural in one framework and awkward in the other; and any setting where you want a direct probability about a hypothesis rather than a tail probability under a null.

Key Insight — the prior is not the only assumption: the usual objection is that priors are subjective. But the likelihood is an assumption too, and it is a much stronger one — choosing "the data are i.i.d. Gaussian" constrains the analysis far more than choosing Beta(2,2) over Beta(1,1). Frequentist methods do not avoid subjectivity, they relocate it into the model, the test statistic, the stopping rule and the \(\alpha\) threshold. The honest difference is that a Bayesian analysis is obliged to write one of its assumptions down explicitly.

Recap

  • Frequentist: parameter fixed, data random, probability = long-run frequency. Bayesian: data given, parameter uncertain, probability = degree of belief.
  • Posterior ∝ likelihood × prior; conjugacy makes the update arithmetic — Beta(α,β) with \(s\) of \(n\) becomes Beta(α+s, β+n−s).
  • Beta(2,2) with 7/10 gives Beta(9,5): posterior mean 0.643, a weighted average of the prior mean 0.5 and the MLE 0.7 with weights 4/14 and 10/14.
  • Credible interval [0.386, 0.861] claims 95% posterior probability; a confidence interval claims 95% coverage of the procedure. Different claims.
  • Bernstein–von Mises: with enough data and a non-degenerate prior the two answers converge. Disagreement is a small-sample, stopping-rule or hierarchical phenomenon.

References

  1. Brown, L. D., Cai, T. T., & DasGupta, A. Interval Estimation for a Binomial Proportion. Statistical Science, 16(2), 101–133, 2001.
  2. Kass, R. E., & Wasserman, L. The Selection of Prior Distributions by Formal Rules. Journal of the American Statistical Association, 91(435), 1343–1370, 1996.
  3. van der Vaart, A. W. Asymptotic Statistics. Cambridge University Press, 1998.
  4. Efron, B., & Hastie, T. Computer Age Statistical Inference. Cambridge University Press, 2016.