Instance-Optimality in I/O-Efficient Sampling and Sequential Estimation

TL;DR

This paper develops instance-optimal algorithms for I/O-efficient sampling and estimation tasks, providing guarantees that adapt to the input's structure without assumptions, and extends these results to various statistical estimation problems.

Contribution

It introduces instance-optimal, order-oblivious algorithms for I/O-efficient sampling and estimation, with broad applications in statistical and distributional analysis.

Findings

Algorithms are instance-optimal among order-oblivious methods.

Results apply to estimating mean, histograms, quantiles, and distribution functions.

Provides non-parametric, instance-optimal bounds for several fundamental problems.

Abstract

Suppose we have a memory storing $0$s and $1$s and we want to estimate the frequency of $1$s by sampling. We want to do this I/O-efficiently, exploiting that each read gives a block of $B$ bits at unit cost; not just one bit. If the input consists of uniform blocks: either all 1s or all 0s, then sampling a whole block at a time does not reduce the number of samples needed for estimation. On the other hand, if bits are randomly permuted, then getting a block of $B$ bits is as good as getting $B$ independent bit samples. However, we do not want to make any such assumptions on the input. Instead, our goal is to have an algorithm with instance-dependent performance guarantees which stops sampling blocks as soon as we know that we have a probabilistically reliable estimate. We prove our algorithms to be instance-optimal among algorithms oblivious to the order of the blocks, which we argue is the strongest form of instance optimality we can hope for. We also present similar results for I/O-efficiently estimating mean with both additive and multiplicative error, estimating histograms, quantiles, as well as the empirical cumulative distribution function. We obtain our above results on I/O-efficient sampling by reducing to corresponding problems in the so-called sequential estimation. In this setting, one samples from an unknown distribution until one can provide an estimate with some desired error probability. We then provide non-parametric instance-optimal results for several fundamental problems: mean and quantile estimation, as well as learning mixture distributions with respect to $\ell_\infty$ and the so-called Kolmogorov-Smirnov distance.