The probability distribution of model parameters after updating prior beliefs with observed data, central to Bayesian inference.The posterior distribution is a fundamental concept in Bayesian statistics that represents the updated probability distribution of model parameters after incorporating observed data. It combines prior beliefs about the parameters with the likelihood of the observed data, yielding a complete probabilistic description of what is known about the parameters after seeing the evidence. This distribution is central to Bayesian inference, as it allows researchers to make probabilistic statements about parameters and quantify uncertainty in a principled way.
The concept of the posterior distribution dates back to Thomas Bayes' work in the 18th century, published posthumously in 1763, and was further formalized by Pierre-Simon Laplace in the early 19th century.
The posterior distribution is computed using Bayes' theorem: the posterior is proportional to the prior distribution multiplied by the likelihood function. In practice, for complex models, the posterior is often not analytically tractable, so computational methods such as Markov Chain Monte Carlo (MCMC) or variational inference are used to sample from or approximate it. The resulting distribution can then be summarized using statistics like the mean, median, or credible intervals, which provide a range of plausible parameter values given the data.
Posterior distributions are used in virtually every field that involves statistical inference, including astrophysics (estimating neutron star properties from gravitational wave signals), medical imaging (reconstructing images from noisy data), and machine learning (Bayesian neural networks).
In the context of gravitational wave astronomy, the posterior distribution for the masses of merging black holes can be so tightly constrained that it allows researchers to determine the redshift of the source with an uncertainty of only a few percent.











