Search papers, labs, and topics across Lattice.
This paper analyzes two existing Monte Carlo integration estimators that leverage determinantal point processes (DPPs) to reduce variance: a biased estimator using a fixed DPP and an unbiased estimator using a DPP tailored to the integrand. The authors generalize both estimators to continuous settings and provide practical sampling algorithms. Their analysis highlights the trade-off between bias and variance reduction when using DPPs for Monte Carlo integration.
DPP-based Monte Carlo integration can offer variance reduction, but choosing the right DPP鈥攆ixed vs. tailored to the integrand鈥攄etermines whether you get a biased but faster converging estimator or an unbiased but standard-rate estimator.
The standard Monte Carlo estimator $\widehat{I}_N^{\mathrm{MC}}$ of $\int fd蠅$ relies on independent samples from $蠅$ and has variance of order $1/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 $f$ and $蠅$. We examine two existing DPP-based estimators: one by Bardenet & Hardy (2020) with a rate of $\mathcal{O}(N^{-(1+1/d)})$ for smooth $f$, but relying on a fixed DPP. The other, by Ermakov & Zolotukhin (1960), is unbiased with rate of order $1/N$, like Monte Carlo, but its DPP is tailored to $f$. We revisit these estimators, generalize them to continuous settings, and provide sampling algorithms.
DeepMind