Monte Carlo and Quasi-Monte Carlo Density Estimation via Conditioning
Abstract
Estimating the unknown density from which a given independent sample originates is more difficult than estimating the mean, in the sense that for the best popular non-parametric density estimators, the mean integrated square error converges more slowly than at the canonical rate of $\mathcal{O}(1/n)$. When the sample is generated from a simulation model and we have control over how this is done, we can do better. We examine an approach in which conditional Monte Carlo yields, under certain conditions, a random conditional density which is an unbiased estimator of the true density at any point. By averaging independent replications, we obtain a density estimator that converges at a faster rate than the usual ones. Moreover, combining this new type of estimator with randomized quasi-Monte Carlo to generate the samples typically brings a larger improvement on the error and convergence rate than for the usual estimators, because the new estimator is smoother as a function of the underlying uniform random numbers.
- Publication:
-
arXiv e-prints
- Pub Date:
- June 2019
- DOI:
- 10.48550/arXiv.1906.04607
- arXiv:
- arXiv:1906.04607
- Bibcode:
- 2019arXiv190604607L
- Keywords:
-
- Mathematics - Statistics Theory;
- 62G07;
- 65C05;
- 62G20
- E-Print:
- Main manuscript: 36 pages, 6 figures, 5 tables. Supplement: 11 pages, 5 figures, 5 tables. We are very thankful to the anonymous referees, whose comments were considered in this submission