Bernstein–von Mises theorem

From HandWiki
Revision as of 15:14, 6 February 2024 by ScienceGen (talk | contribs) (correction)
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Short description: Theorem in Bayesian inference

In Bayesian inference, the Bernstein–von Mises theorem provides the basis for using Bayesian credible sets for confidence statements in parametric models. It states that under some conditions, a posterior distribution converges in the limit of infinite data to a multivariate normal distribution centered at the maximum likelihood estimator with covariance matrix given by [math]\displaystyle{ n^{-1} I(\theta_0)^{-1} }[/math], where [math]\displaystyle{ \theta_0 }[/math] is the true population parameter and [math]\displaystyle{ I(\theta_0) }[/math] is the Fisher information matrix at the true population parameter value:[1]

[math]\displaystyle{ P(\theta|x_1,\dots x_n)= \mathcal{N}(\theta_0, n^{-1}I(\theta_0)^{-1}) \text{ for } n\to \infty. }[/math]

The Bernstein–von Mises theorem links Bayesian inference with frequentist inference. It assumes there is some true probabilistic process that generates the observations, as in frequentism, and then studies the quality of Bayesian methods of recovering that process, and making uncertainty statements about that process. In particular, it states that Bayesian credible sets of a certain credibility level [math]\displaystyle{ \alpha }[/math] will asymptotically be confidence sets of confidence level [math]\displaystyle{ \alpha }[/math], which allows for the interpretation of Bayesian credible sets.

Heuristic statement

In a model [math]\displaystyle{ (P_\theta: \theta \in \Theta) }[/math], under certain regularity conditions (finite-dimensional, well-specified, smooth, existence of tests), if the prior distribution [math]\displaystyle{ \Pi }[/math] on [math]\displaystyle{ \theta }[/math] has a density with respect to the Lebesgue measure which is smooth enough (near [math]\displaystyle{ \theta_0 }[/math] bounded away from zero), the total variation distance between the rescaled posterior distribution (by centring and rescaling to [math]\displaystyle{ \sqrt{n}(\theta - \theta_0) }[/math]) and a Gaussian distribution centred on any efficient estimator and with the inverse Fisher information as variance will converge in probability to zero.

Bernstein–von Mises and maximum likelihood estimation

In case the maximum likelihood estimator is an efficient estimator, we can plug this in, and we recover a common, more specific, version of the Bernstein–von Mises theorem.

Implications

The most important implication of the Bernstein–von Mises theorem is that the Bayesian inference is asymptotically correct from a frequentist point of view. This means that for large amounts of data, one can use the posterior distribution to make, from a frequentist point of view, valid statements about estimation and uncertainty.

History

The theorem is named after Richard von Mises and S. N. Bernstein, although the first proper proof was given by Joseph L. Doob in 1949 for random variables with finite probability space.[2] Later Lucien Le Cam, his PhD student Lorraine Schwartz, David A. Freedman and Persi Diaconis extended the proof under more general assumptions.[citation needed]

Limitations

In case of a misspecified model, the posterior distribution will also become asymptotically Gaussian with a correct mean, but not necessarily with the Fisher information as the variance. This implies that Bayesian credible sets of level [math]\displaystyle{ \alpha }[/math] cannot be interpreted as confidence sets of level [math]\displaystyle{ \alpha }[/math].[3]

In the case of nonparametric statistics, the Bernstein–von Mises theorem usually fails to hold with a notable exception of the Dirichlet process.

A remarkable result was found by Freedman in 1965: the Bernstein–von Mises theorem does not hold almost surely if the random variable has an infinite countable probability space; however, this depends on allowing a very broad range of possible priors. In practice, the priors used typically in research do have the desirable property even with an infinite countable probability space.

Different summary statistics such as the mode and mean may behave differently in the posterior distribution. In Freedman's examples, the posterior density and its mean can converge on the wrong result, but the posterior mode is consistent and will converge on the correct result.

Notes

  1. van der Vaart, A.W. (1998). "10.2 Bernstein–von Mises Theorem". Asymptotic Statistics. Cambridge University Press. ISBN 0-521-78450-6. 
  2. Doob, Joseph L. (1949). "Application of the theory of martingales". Colloq. Intern. Du C.N.R.S (Paris) 13: 23–27. 
  3. Kleijn, B.J.K.; van der Vaart, A.W. (2012). "The Bernstein-Von–Mises theorem under misspecification". Electronic Journal of Statistics 6: 354–381. doi:10.1214/12-EJS675. 

References

  • van der Vaart, A.W. (1998). "10.2 Bernstein–von Mises Theorem". Asymptotic Statistics. Cambridge University Press. ISBN 0-521-49603-9. 
  • Doob, Joseph L. (1949), Application of the theory of martingales. Colloq. Intern. du C.N.R.S (Paris), No. 13, pp. 23–27.
  • Freedman, David A. (1963). On the asymptotic behaviour of Bayes estimates in the discrete case I. The Annals of Mathematical Statistics, vol. 34, pp. 1386–1403.
  • Freedman, David A. (1965). On the asymptotic behaviour of Bayes estimates in the discrete case II. The Annals of Mathematical Statistics, vol. 36, pp. 454–456.
  • Le Cam, Lucien (1986). Asymptotic Methods in Statistical Decision Theory, Springer. ISBN:0-387-96307-3 (Pages 336 and 618–621).
  • Lorraine Schwartz (1965). On Bayes procedures. Z. Wahrscheinlichkeitstheorie, No. 4, pp. 10–26.