38
3

High-Order Statistical Functional Expansion and Its Application To Some Nonsmooth Problems

Abstract

Let \bxj=\btheta+\bepj,j=1,...,n\bx_j = \btheta +\bep_j, j=1,...,n, be observations of an unknown parameter \btheta\btheta in a Euclidean or separable Hilbert space \scrH\scrH, where \bepj\bep_j are noises as random elements in \scrH\scrH from a general distribution. We study the estimation of f(\btheta)f(\btheta) for a given functional f:\scrH\RRf:\scrH\rightarrow \RR based on \bxj\bx_j's. The key element of our approach is a new method which we call High-Order Degenerate Statistical Expansion. It leverages the use of classical multivariate Taylor expansion and degenerate UU-statistic and yields an elegant explicit formula. In the univariate case of \scrH=R\scrH=\R, the formula expresses the error of the proposed estimator as a sum of order kk degenerate UU-products of the noises with coefficient f(k)(\btheta)/k!f^{(k)}(\btheta)/k! and an explicit remainder term in the form of the Riemann-Liouville integral as in the Taylor expansion around the true \btheta\btheta. For general \scrH\scrH, the formula expresses the estimation error in terms of the inner product of f(k)(\btheta)/k!f^{(k)}(\btheta)/k! and the average of the tensor products of kk noises with distinct indices and a parallel extension of the remainder term from the univariate case. This makes the proposed method a natural statistical version of the classical Taylor expansion. The proposed estimator can be viewed as a jackknife estimator of an ideal degenerate expansion of f()f(\cdot) around the true \btheta\btheta with the degenerate UU-product of the noises, and can be approximated by bootstrap. Thus, the jackknife, bootstrap and Taylor expansion approaches all converge to the proposed estimator. We develop risk bounds for the proposed estimator and a central limit theorem under a second moment condition (even in expansions of higher than the second order). We apply this new method to generalize several existing results with smooth and nonsmooth ff to universal \bepj\bep_j's with only minimum moment constraints.

View on arXiv
Comments on this paper