

Author: Bornkamp Björn
Publisher: Taylor & Francis Ltd
ISSN: 1061-8600
Source: Journal of Computational and Graphical Statistics, Vol.20, Iss.3, 2011-01, pp. : 656-669
Disclaimer: Any content in publications that violate the sovereignty, the constitution or regulations of the PRC is not accepted or approved by CNPIEC.
Abstract
The Laplace approximation is an old, but frequently used method to approximate integrals for Bayesian calculations. In this article we develop an extension of the Laplace approximation, by applying it iteratively to the residual, that is, the difference between the current approximation and the true function. The final approximation is thus a linear combination of multivariate normal densities, where the coefficients are chosen to achieve a good fit to the target distribution. We illustrate on real and artificial examples that the proposed procedure is a computationally efficient alternative to current approaches for approximation of multivariate probability densities. This article has supplementary material online.
Related content


96.51 Iterated dissections? Approximations to
The Mathematical Gazette, Vol. 96, Iss. 537, 2012-11 ,pp. :


Stability of precise Laplace's method under approximations;Applications
By Guionnet A.
ESAIM: Probability and Statistics, Vol. 3, Iss. issue, 2010-03 ,pp. :




Tail probability of random variable and Laplace transform
Applicable Analysis, Vol. 84, Iss. 5, 2005-05 ,pp. :


Uniformity of Double Saddlepoint Conditional Probability Approximations
By Kolassa J.E.
Journal of Multivariate Analysis, Vol. 64, Iss. 1, 1998-01 ,pp. :