Approximations for the number of maxima and near-maxima in independent data
Fraser Daly
TL;DR
This work studies, for iid samples, the distributional error in approximating the number of observations attaining the sample maximum (discrete case) or within a threshold of a sample order statistic (continuous case). It develops explicit total variation bounds using Stein's method, introducing a novel logarithmic target distribution for the discrete setting and applying negative binomial approximations via mixed binomial representations for the continuous setting. The authors provide concrete bounds and illustrative examples (geometric, Gumbel, Uniform) and establish a Poisson approximation in the discrete case as a complementary tool. The results yield practical, quantitative error guarantees for simple, well-known distributions (logarithmic, NB, Poisson) as approximants to complex maxima-related statistics, with potential extensions to dependent data and broader applications in reliability and algorithmic contexts.
Abstract
In the setting where we have $n$ independent observations of a random variable $X$, we derive explicit error bounds in total variation distance when approximating the number of observations equal to the maximum of the sample (in the case where $X$ is discrete) or the number of observations within a given distance of an order statistic of the sample (in the case where $X$ is absolutely continuous). The logarithmic and Poisson distributions are used as approximations in the discrete case, with proofs which include the development of Stein's method for a logarithmic target distribution. In the absolutely continuous case our approximations are by the negative binomial distribution, and are established by considering negative binomial approximation for mixed binomials. The cases where $X$ is geometric, Gumbel and uniform are used as illustrative examples.
