cs.LGApr 21, 2026

On two ways to use determinantal point processes for Monte Carlo integration

Authors: Guillaume GautierRémi BardenetMichal Valko

Organizations: †Univ. Lille, CNRS, Centrale Lille, UMR 9189 – CRIStAL, 59651 Villeneuve d’Ascq, France · *Inria Lille-Nord Europe, 40 avenue Halley 59650 Villeneuve d’Ascq, France · ‡DeepMind Paris, 14 Rue de Londres, 75009 Paris, France · ∗Inria Lille-Nord Europe, 40 avenue Halley 59650 Villeneuve d’Ascq, France

Abstract

The standard Monte Carlo estimator I^NMC\widehat{I}_N^{\mathrm{MC}} of fdω\int fdω relies on independent samples from ωω and has variance of order 1/N1/N. Replacing the samples with a determinantal point process (DPP), a repulsive distribution, makes the estimator consistent, with variance rates that depend on how the DPP is adapted to ff and ωω. We examine two existing DPP-based estimators: one by Bardenet & Hardy (2020) with a rate of O(N(1+1/d))\mathcal{O}(N^{-(1+1/d)}) for smooth ff, but relying on a fixed DPP. The other, by Ermakov & Zolotukhin (1960), is unbiased with rate of order 1/N1/N, like Monte Carlo, but its DPP is tailored to ff. We revisit these estimators, generalize them to continuous settings, and provide sampling algorithms.

Explore similar work

CardsList
  1. Variance Reduction for Independent Metropolis

    Jun 25, 2024Siran Liu, Petros Dellaportas, Michalis K. TitsiasVariance Reduction