Optimal estimation of high-dimensional Gaussian location mixtures

Abstract

This paper studies the optimal rate of estimation in a finite Gaussian location mixture model in high dimensions without separation conditions. We assume that the number of components k is bounded and that the centers lie in a ball of bounded radius, while allowing the dimension d to be as large as the sample size n. Extending the one-dimensional result of Heinrich and Kahn (Ann. Statist. 46 (2018) 2844–2870), we show that the minimax rate of estimating the mixing distribution in Wasserstein distance is $\mathrm{\Theta }({(\mathit{d}/\mathit{n})}^{1/4}+{\mathit{n}}^{-1/(4\mathit{k}-2)})$, achieved by an estimator computable in time $\mathit{O}(\mathit{n}{\mathit{d}}^2+{\mathit{n}}^{5/4})$. Furthermore, we show that the mixture density can be estimated at the optimal parametric rate $\mathrm{\Theta }(\sqrt{\mathit{d}/\mathit{n}})$ in Hellinger distance and provide a computationally efficient algorithm to achieve this rate in the special case of $\mathit{k}=2$.

Both the theoretical and methodological development rely on a careful application of the method of moments. Central to our results is the observation that the information geometry of finite Gaussian mixtures is characterized by the moment tensors of the mixing distribution, whose low-rank structure can be exploited to obtain a sharp local entropy bound.