Point Process-based Monte Carlo estimation
Abstract
This paper addresses the issue of estimating the expectation of a real-valued random variable of the form where is a deterministic function and can be a random finite- or infinite-dimensional vector. Using recent results on rare event simulation, we propose a unified framework for dealing with both probability and mean estimation for such random variables, \emph{i.e.} linking algorithms such as Tootsie Pop Algorithm (TPA) or Last Particle Algorithm with nested sampling. Especially, it extends nested sampling as follows: first the random variable does not need to be bounded any more: it gives the principle of an ideal estimator with an infinite number of terms that is unbiased and always better than a classical Monte Carlo estimator -- in particular it has a finite variance as soon as there exists such that . Moreover we address the issue of nested sampling termination and show that a random truncation of the sum can preserve unbiasedness while increasing the variance only by a factor up to 2 compared to the ideal case. We also build an unbiased estimator with fixed computational budget which supports a Central Limit Theorem and discuss parallel implementation of nested sampling, which can dramatically reduce its computational cost. Finally we extensively study the case where is heavy-tailed.
Cite
@article{arxiv.1412.6368,
title = {Point Process-based Monte Carlo estimation},
author = {Clément Walter},
journal= {arXiv preprint arXiv:1412.6368},
year = {2015}
}
Comments
13 pages + 4 pages of appendix, 7 figures