27
2

Dimension-free Private Mean Estimation for Anisotropic Distributions

Abstract

We present differentially private algorithms for high-dimensional mean estimation. Previous private estimators on distributions over Rd\mathbb{R}^d suffer from a curse of dimensionality, as they require Ω(d1/2)\Omega(d^{1/2}) samples to achieve non-trivial error, even in cases where O(1)O(1) samples suffice without privacy. This rate is unavoidable when the distribution is isotropic, namely, when the covariance is a multiple of the identity matrix, or when accuracy is measured with respect to the affine-invariant Mahalanobis distance. Yet, real-world data is often highly anisotropic, with signals concentrated on a small number of principal components. We develop estimators that are appropriate for such signals\unicodex2013\unicode{x2013}our estimators are (ε,δ)(\varepsilon,\delta)-differentially private and have sample complexity that is dimension-independent for anisotropic subgaussian distributions. Given nn samples from a distribution with known covariance-proxy Σ\Sigma and unknown mean μ\mu, we present an estimator μ^\hat{\mu} that achieves error μ^μ2α\|\hat{\mu}-\mu\|_2\leq \alpha, as long as ntr(Σ)/α2+tr(Σ1/2)/(αε)n\gtrsim\mathrm{tr}(\Sigma)/\alpha^2+ \mathrm{tr}(\Sigma^{1/2})/(\alpha\varepsilon). In particular, when σ2=(σ12,,σd2)\pmb{\sigma}^2=(\sigma_1^2, \ldots, \sigma_d^2) are the singular values of Σ\Sigma, we have tr(Σ)=σ22\mathrm{tr}(\Sigma)=\|\pmb{\sigma}\|_2^2 and tr(Σ1/2)=σ1\mathrm{tr}(\Sigma^{1/2})=\|\pmb{\sigma}\|_1, and hence our bound avoids dimension-dependence when the signal is concentrated in a few principal components. We show that this is the optimal sample complexity for this task up to logarithmic factors. Moreover, for the case of unknown covariance, we present an algorithm whose sample complexity has improved dependence on the dimension, from d1/2d^{1/2} to d1/4d^{1/4}.

View on arXiv
Comments on this paper