Quantile Importance Sampling

Abstract

In Bayesian inference, the approximation of integrals of the form ψ = EF l(X) = χ l(x)dF(x) is a fundamental challenge. Such integrals are crucial for evidence estimation, which is important for various purposes, including model selection and numerical analysis. The existing strategies for evidence estimation are classified into four categories: deterministic approximation, density estimation, importance sampling, and vertical representation (Llorente et al., 2023). In this paper, we show that the Riemann sum estimator due to Yakowitz, Krimmel and Szidarovszky (1978) can be used in the context of nested sampling (Skilling, 2006) to achieve a O(n−4) rate of convergence, faster than the usual Ergodic Central Limit Theorem, under certain regularity conditions. We provide a brief overview of the literature on the Riemann sum estimators and the nested sampling algorithm and its connections to vertical likelihood Monte Carlo. We provide theoretical and numerical arguments to show how merging these two ideas may result in improved and more robust estimators for evidence estimation, especially in higher dimensional spaces. We also briefly discuss the idea of simulating the Lorenz curve that avoids the problem of intractable Λ functions, essential for the vertical representation and nested sampling.