Minimax estimation of linear and quadratic functionals on sparsity classes

For the Gaussian sequence model, we obtain nonasymptotic minimax rates of estimation of the linear, quadratic and the l(2)-norm functionals on classes of sparse vectors and construct optimal estimators that attain these rates. The main object of interest is the class B-0(s) of s-sparse vectors theta = (theta(1) (,...,) theta(d)), for which we also provide completely adaptive estimators (independent of s and of the noise variance sigma) having logarithmically slower rates than the minimax ones. Furthermore, we obtain the minimax rates on the lq-balls Bq (r) = {theta is an element of R-d : parallel to theta parallel to(q) <= r} where 0 < q <= 2, and parallel to theta parallel to q = (Sigma(d)(i=1))1/q. This analysis shows that there are, in general, three zones in the rates of convergence that we call the sparse zone, the dense zone and the degenerate zone, while a fourth zone appears for estimation of the quadratic functional. We show that, as opposed to estimation of., the correct logarithmic terms in the optimal rates for the sparse zone scale as log(d/s(2)) and not as log(d/s). For the class B-0(s), the rates of estimation of the linear functional and of the l(2)-norm have a simple elbow at s = root d (boundary between the sparse and the dense zones) and exhibit similar performances, whereas the estimation of the quadratic functional Q(theta) reveals more complex effects: the minimax risk on B-0(s) is infinite and the sparseness assumption needs to be combined with a bound on the l(2)-norm. Finally, we apply our results on estimation of the l(2)-norm to the problem of testing against sparse alternatives. In particular, we obtain a nonasymptotic analog of the Ingster-Donoho-Jin theory revealing some effects that were not captured by the previous asymptotic analysis.