The Voronoigram: Minimax Estimation of Bounded Variation Functions From Scattered Data

TL;DR

This paper introduces the Voronoigram, a novel estimator for multivariate bounded variation functions using Voronoi diagrams and total variation regularization, achieving near-optimal rates in noisy, scattered data settings.

Contribution

The paper develops the Voronoigram estimator, analyzes its properties, and proves its minimax optimality for bounded variation functions in a scattered data context.

Findings

Voronoigram performs local averaging over Voronoi cells.

It is asymptotically equivalent to a continuum TV functional.

The estimator is minimax rate optimal up to log factors.

Abstract

We consider the problem of estimating a multivariate function $f_0$ of bounded variation (BV), from noisy observations $y_i = f_0(x_i) + z_i$ made at random design points $x_i \in \mathbb{R}^d$, $i=1,\ldots,n$. We study an estimator that forms the Voronoi diagram of the design points, and then solves an optimization problem that regularizes according to a certain discrete notion of total variation (TV): the sum of weighted absolute differences of parameters $\theta_i,\theta_j$ (which estimate the function values $f_0(x_i),f_0(x_j)$) at all neighboring cells $i,j$ in the Voronoi diagram. This is seen to be equivalent to a variational optimization problem that regularizes according to the usual continuum (measure-theoretic) notion of TV, once we restrict the domain to functions that are piecewise constant over the Voronoi diagram. The regression estimator under consideration hence performs (shrunken) local averaging over adaptively formed unions of Voronoi cells, and we refer to it as the Voronoigram, following the ideas in Koenker (2005), and drawing inspiration from Tukey's regressogram (Tukey, 1961). Our contributions in this paper span both the conceptual and theoretical frontiers: we discuss some of the unique properties of the Voronoigram in comparison to TV-regularized estimators that use other graph-based discretizations; we derive the asymptotic limit of the Voronoi TV functional; and we prove that the Voronoigram is minimax rate optimal (up to log factors) for estimating BV functions that are essentially bounded.