Approximating Probability Densities by Iterated Laplace Approximations
View Researcher's Other CodesDisclaimer: The provided code links for this paper are external links. Science Nest has no responsibility for the accuracy, legality or content of these links. Also, by downloading this code(s), you agree to comply with the terms of use as set out by the author(s) of the code(s).
Authors | Bjorn Bornkamp |
Journal/Conference Name | Journal of Computational and Graphical Statistics |
Paper Category | Other |
Paper 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. |
Date of publication | 2011 |
Code Programming Language | R |
Comment |