68

Robust estimation of the mean with bounded relative standard deviation

Springer Proceedings in Mathematics & statistics (MS), 2018
Abstract

Many randomized approximation algorithms operate by giving a procedure for simulating a random variable XX which has mean μ\mu equal to the target answer, and a relative standard deviation bounded above by a known constant cc. Examples of this type of algorithm includes methods for approximating the number of satisfying assignments to 2-SAT or DNF, the volume of a convex body, and the partition function of a Gibbs distribution. Because the answer is usually exponentially large in the problem input size, it is typical to require an estimate μ^\hat \mu satisfy P(μ^/μ1>ϵ)δ\mathbb{P}(|\hat \mu/\mu - 1| > \epsilon) \leq \delta, where ϵ\epsilon and δ\delta are user specified nonnegative parameters. The current best algorithm uses 2c2ϵ2(1+ϵ)2ln(2/δ)\lceil 2c^2\epsilon^{-2}(1+\epsilon)^2 \ln(2/\delta) \rceil samples to achieve such an estimate. By modifying the algorithm in order to balance the tails, it is possible to improve this result to 2(c2ϵ2+1)/(1ϵ2)ln(2/δ)\lceil 2(c^2\epsilon^{-2} + 1)/(1-\epsilon^2)\ln(2/\delta) \rceil samples. Aside from the theoretical improvement, we also consider how to best implement this algorithm in practice. Numerical experiments show the behavior of the estimator on distributions where the relative standard deviation is unknown or infinite.

View on arXiv
Comments on this paper