Abstract:Gaussian approximations are routinely employed in Bayesian statistics to ease inference when the target posterior is intractable. Although these approximations are asymptotically justified by Bernstein-von Mises type results, in practice the expected Gaussian behavior may poorly represent the shape of the posterior, thus affecting approximation accuracy. Motivated by these considerations, we derive an improved class of closed-form approximations of posterior distributions which arise from a new treatment of a third-order version of the Laplace method yielding approximations in a tractable family of skew-symmetric distributions. Under general assumptions which account for misspecified models and non-i.i.d. settings, this family of approximations is shown to have a total variation distance from the target posterior whose rate of convergence improves by at least one order of magnitude the one established by the classical Bernstein-von Mises theorem. Specializing this result to the case of regular parametric models shows that the same improvement in approximation accuracy can be also derived for polynomially bounded posterior functionals. Unlike other higher-order approximations, our results prove that it is possible to derive closed-form and valid densities which are expected to provide, in practice, a more accurate, yet similarly-tractable, alternative to Gaussian approximations of the target posterior, while inheriting its limiting frequentist properties. We strengthen such arguments by developing a practical skew-modal approximation for both joint and marginal posteriors that achieves the same theoretical guarantees of its theoretical counterpart by replacing the unknown model parameters with the corresponding MAP estimate. Empirical studies confirm that our theoretical results closely match the remarkable performance observed in practice, even in finite, possibly small, sample regimes.

Variational Inference based on a Subclass of Closed Skew Normals

Structured Variational Approximations with Skew Normal Decomposable Graphical Models and Implicit Copulas

On mean and/or variance mixtures of normal distributions

Skewed Bernstein-von Mises theorem and skew-modal approximations

Informative Bayesian inference for the skew-normal distribution

On Bayesian predictive density estimation for skew-normal distributions

Comparison of Azzalini and Geometric Skew Normal Distributions Under Bayesian Paradigm

Statistical Inference for the Location and Scale Parameters of the Skew Normal Distribution

High-dimensional copula variational approximation through transformation

Variational Inference: A Review for Statisticians

Multivariate Geometric Skew-Normal Distribution

Mean mixtures of normal distributions: properties, inference and application

Statistical applications of the multivariate skew-normal distribution

Skew-symmetric approximations of posterior distributions

Statistical applications of the multivariate skew normal distribution

Clustering of non-Gaussian data by variational Bayes for normal inverse Gaussian mixture models

Sparse bayesian inference with regularized gaussian distributions

Variational Bayes Approximations for Clustering via Mixtures of Normal Inverse Gaussian Distributions

A flexible generalization of the skew normal distribution based on a weighted normal distribution

A Novel Skewed Generalized Normal Distribution: Properties, Statistical Inference, and Its Applications

Variational sparse inverse Cholesky approximation for latent Gaussian processes via double Kullback-Leibler minimization