Estimation of discrete distributions in relative entropy, and the deviations of the missing mass
Jaouad Mourtada
TL;DR
The paper provides sharp non-asymptotic high-probability guarantees for estimating discrete distributions under KL divergence. It first analyzes the Laplace (add-one) estimator, establishing a tight upper bound and proving its optimality among confidence-independent methods, while also deriving a matching lower bound. It then proves minimax-optimal high-probability guarantees for confidence-dependent smoothing and demonstrates that an adaptive, data-driven smoothing approach can beat fixed smoothing in sparse regimes. A major contribution is the introduction of effective sparsity concepts and adaptive estimators that leverage them to achieve near-optimal rates across varying support sizes, together with a sharp high-probability bound on missing mass. The results reveal a fundamental separation between asymptotic rates and uniform non-asymptotic guarantees, and they quantify the statistical-computational trade-offs in high-probability discrete distribution estimation.
Abstract
We study the problem of estimating a distribution over a finite alphabet from an i.i.d. sample, with accuracy measured in relative entropy (Kullback-Leibler divergence). While optimal bounds on the expected risk are known, high-probability guarantees remain less well-understood. First, we analyze the classical Laplace (add-one) estimator, obtaining matching upper and lower bounds on its performance and establishing its optimality among confidence-independent estimators. We then characterize the minimax-optimal high-probability risk and show that it is achieved by a simple confidence-dependent smoothing technique. Notably, the optimal non-asymptotic risk incurs an additional logarithmic factor compared to the ideal asymptotic rate. Next, motivated by regimes in which the alphabet size exceeds the sample size, we investigate methods that adapt to the sparsity of the underlying distribution. We introduce an estimator using data-dependent smoothing, for which we establish a high-probability risk bound depending on two effective sparsity parameters. As part of our analysis, we also derive a sharp high-probability upper bound on the missing mass.
