Minimax Estimation of Functionals of Discrete Distributions

We propose a general methodology for the construction and analysis of essentially minimax estimators for a wide class of functionals of finite dimensional parameters, and elaborate on the case of discrete distributions, where the support size is unknown and may be comparable with or even much larger than the number of observations . We treat the respective regions where the functional is nonsmooth and smooth separately. In the nonsmooth regime, we apply an unbiased estimator for the best polynomial approximation of the functional whereas, in the smooth regime, we apply a bias-corrected version of the maximum likelihood estimator (MLE). We illustrate the merit of this approach by thoroughly analyzing the performance of the resulting schemes for estimating two important information measures: 1) the entropy and 2) . We obtain the minimax rates for estimating these functionals. In particular, we demonstrate that our estimator achieves the optimal sample complexity for entropy estimation. We also demonstrate that the sample complexity for estimating , is , which can be achieved by our estimator but not the MLE. For - /inline-formula>, we show the minimax rate for estimating is for infinite support size, while the maximum rate for the MLE is . For all the above cases, the behavior of the minimax rate-optimal estimators with samples is essentially that of the MLE (plug-in rule) with samples, which we term “effective sample size enlargement.” We highlight the practical advantages of our schemes for the estimation of entropy and mutual information. We compare our performance with various existing approaches, and demonstrate that our approach reduces running time and boosts the accuracy. Moreover, we show that the minimax rate-optimal mutual information estimator yielded by our framework leads to significant performance boosts over the Chow–Liu algorithm in learning graphical models. The wide use of information measure estimation suggests that the insights and estimators obtained in this paper could be broadly applicable.