Quadratic functional estimation in inverse problems

Cristina Butucea (1) , Katia Meziani (2)

1 LPP - Laboratoire Paul Painlevé - UMR 8524
2 CEREMADE - CEntre de REcherches en MAthématiques de la DEcision

In this paper, we consider a Gaussian sequence of independent observations having a polynomially increasing variance. This model describes a large panel of inverse problems, such as the deconvolution of blurred images or the recovering of the fractional derivative of a signal. We estimate the sum of squares of the means of our observations. This quadratic functional has practical meanings, e.g. the energy of a signal, and it is often used for goodness-of-fit testing. We compute Pinsker estimators when the underlying signal has both a finite and infinite amount of smoothness. When the signal is sufficiently smoother than the difficulty of the inverse problem, we attain the parametric rate and the efficiency constant associated with it. Moreover, we give upper bounds of the second order term in the risk. Otherwise, when the parametric rate cannot be attained, we compute non parametric upper bounds of the risk.

Mots clés: second order risk; Gaussian sequence model; inverse problem; minimax upper bounds; parametric rate; Pinsker estimator; projection estimator; quadratic functional; second order risk.

Domaines: Statistiques [math.ST]; Théorie [stat.TH]

Cristina Butucea, Katia Meziani. Quadratic functional estimation in inverse problems. Statistical Methodology, 2010, 8 (1), pp.31-41. ⟨10.1016/j.stamet.2010.05.002⟩. ⟨hal-00634945⟩