Bayesian Matching Estimator

The Bayesian Matching Estimator estimates average treatment effects in observational studies by combining classical nearest-neighbour or kernel matching with a Bayesian posterior over the treatment effect. It inherits matching's covariate-balancing logic while propagating uncertainty through a full posterior distribution rather than relying on asymptotic standard errors, yielding credible intervals that reflect both sampling variability and prior knowledge.