The authors would like to thank Projects “Agenda Petróleo” and “ECOS-Nord” for their financial support.
<ph f="cmbx">A Statistical view of Iterative Methods for Linear Inverse Problems</ph>

### Carenne Luden͂a

Escuela de Matematicas, Facultad de Ciencias, UCV, Av. Los Ilustres, Los Chaguaramos, Codigo Postal 1020-A, Caracas Venezuela. Telf.: (58)212-6051481.
Departamento de Matematicas, IVIC, Carretera Panamericana KM. 11, Aptdo. 21827, Codigo Postal 1020-A, Caracas Venezuela. Telf.: (58)212-5041412.
E-mail address : afermin@euler.ciens.ucv.ve E-mail address : cludena@ivic.ve
• Abstract. In this article we study the problem of recovering the unknown solution of a linear ill-posed problem, via iterative regularization methods. We review the problem of projection-regularization from a statistical point of view. A basic purpose of the paper is the consideration of adaptive model selection for determining regularization parameters. This article introduces a new regularized estimator which has the best possible adaptive properties for a wide range of linear functionals. We derive non asymptotic upper bounds for the mean square error of the estimator and give the optimal convergence rates.

1 Introduction

The area of mathematical inverse problems is quite broad and involves the qualitative and quantitative analysis of a wide variety of physical models. Moreover, a considerable number of problems arising in different scientific and technical fields belong to a class of ill-posed problems. For example, geophysicists scan the earth's subsurface by recording arrival times of waves reflected off different layers underneath the surface, and try to determine a meaningful solution and to understand which features in the solution are statistically significant.
From a statistical point of view, the problem can be seen as recovering an unobservable signal $\stackrel{~}{f}$  based on observations
 $\begin{array}{c}y\left({x}_{i}\right)=A\stackrel{~}{f}\left({x}_{i}\right)+{\varepsilon }_{i},\end{array}$ (1.1)
where $A:F\to Y$  is some known compact linear operator defined over a separable Hilbert space $F$  , with values in a separable Hilbert space $Y$  and ${x}_{i},i=1,\dots ,n$  is a fixed observation scheme. We assume that the observations $y\left({x}_{i}\right)\in IR$  and that the observation noise ${\varepsilon }_{i}$  are i.i.d. realizations of a certain random variable $\varepsilon$  . Throughout the paper, we shall denote $\mathbf{y}=\left(y\left({x}_{i}\right){\right)}_{i=1}^{n}$  . In this article we study the problem of estimating $\stackrel{~}{f}$  using fixed iterative methods.
The best possible accuracy, regardless of any discretization and noise corruption is determined by some a priori smoothness assumption on the exact solution $\stackrel{~}{f}$  . Here, smoothness is given in terms of some index function $\eta$  on the spectrum de ${A}^{*}A$  by ${A}_{\eta ,\rho }=\left\{f\in F,f=\eta \left({A}^{*}A\right)\omega ,\parallel \omega \parallel \le \rho \right\}$  where ${A}_{\eta ,\rho }$  is called a source condition. For classical Hilbert scales, the smoothness is measured in terms of powers $\eta \left(t\right):={t}^{\mu }$  with $0\le \mu \le {\mu }_{0}$  , ${\mu }_{0}>0$  .
In a deterministic framework, the statistical model  1.1 is formulated as the problem of finding the best-approximate solution of $Af=y,$  in the situation where only perturbed data ${y}^{\delta }$  are available with $\parallel y-{y}^{\delta }\parallel \le \delta .$  Here, $\delta$  is called the noise level. It is important to remark that whereas in this case consistency of the estimators depends on the approximation parameter $\delta$  , in ( 1.1 ) it depends on the number of observations $n$  .
In general, the best ${L}^{2}$  approximation ${A}^{†}y$  , where ${A}^{†}$  is the Moore-Penrose (generalized) inverse of A, does not depend continuously on the left-hand side $y$  . We define the Moore-Penrose inverse in an operator-theoretic way by restricting the domain and range of $A$  in such a way that the resulting restricted operator is invertible; its inverse will then be extended to its maximal domain $\mathcal{D}\left({A}^{†}\right)=\mathcal{ℛ}\left(A\right)+\mathcal{ℛ}\left(A{\right)}^{\perp }$  , with $\mathcal{ℛ}\left(A\right)$  the range of the operator $A$  and $\mathcal{ℛ}\left(A{\right)}^{\perp }$  the orthogonal complement of the range of $A$  .
The inverse problems that we study in this article are called ill-posed problems because the operator $A$  is compact and consequently equation  1.1 can not be inverted directly since ${A}^{-1}$  is not a bounded operator. Ill-posed problems are usually treated by applying some linear regularization procedure, often based on a singular value decomposition; see Tikhonov and Arsenin in [23. An interesting early survey of the statistical perspective on ill-posed problems is studied in great detail by O'Sullivan in [21.
In practice however  1.1 is hardly ever considered. Instead, we project the problem onto a smaller dimensional space ${Y}_{m}$  of $Y$  . This yields a sequence of closed subspaces ${Y}_{m}$  indexed by $m\in {M}_{n}$  , a collection of index sets. Clearly, an important problem is thus how to choose subspace ${Y}_{m}$  based on the data. This can be done by selection of a cutoff point or by threshold methods. Choosing the right subspace will be called model selection.
Sometimes this projection provides enough regularization to produce a good approximate solution, but often additional regularization is needed. Regularization methods replace an ill-posed problem by a family of well-posed problems, their solution, called regularized solutions, are used as approximations to the desired solution of the inverse problem. These methods always involve some parameter measuring the closeness of the regularized and the original (unregularized) inverse problem, rules (and algorithms) for the choice of these regularization parameters as well as convergence properties of the regularized solutions are central points in the theory of these methods, since they allow to find the right balance between stability and accuracy. The general principles of regularization for ill-posed problems are known. In particular, such principles have been established by A.N. Tikhonov. The literature on various regularization methods based on these general principles is extensive ( Engl, Hanke, and Neubauer [9, Gilyazov and Gol'dman [11).
In statistic, regularization, is associated to penalty based methods or thresholding methods or more generality to “smoothing” techniques. In applications, regularization offers a unifying perspective for many diverse ill-posed inverse problems, a wide range of problems concerned with recovering information from indirect and usually noisy measurements, arising in geophysics, tomography and econometrics. One of the most important, but still insufficiently developed, topics in the theory of ill posed problems is connected with iteration regularization [11; i.e, with the utilization of iteration methods of any form for the stable approximate solution of ill-posed problems. Iterative regularization methods tend to be more attractive in terms of numerical cost and implementation, but a number of open questions remain in their theoretical analysis.
In this article we propose an iterative regularized estimator for linear ill-posed problems.
Necessary conditions for convergence are established. These conditions connect the choice of the regularization parameter (i.e., the iteration index) with the projection dimension.
Moreover, we prove that the iterative regularized estimator is optimal in the sense that the estimator achieves the best rate of convergence among all the regularized estimators. A recent work in this direction is developed by Loubes and Luden͂a in [16, which discusses the problem of estimating inverse nonlinear ill-posed problems with different types of complexity penalties leading either to a model selection estimator or to a regularized estimator.
The choice of the regularization sequence is here crucial, and a lot of work associated with the selection of a good regularization parameter can be found in the literature [10, [13. When using iterative methods the problem is finding a good stopping criteria for terminating the iteration procedure. In this article we will use tools developed in the context of model selection via penalization, [2,[1, based on the use of concentration inequalities.
Our article is organized as follows. Section 2 presents basic assumptions and a statement of the discretized inverse problem. In section 3 we discuss regularization methods and we prove consistency of the estimator when the regularization parameter is known. In section 4 we present our main result, prove optimality of an adaptive regularized estimator and give its rate of convergence. Finally, in the last section we introduce regularization by iterative methods for the solution of inverse problems and provide some examples to explain the properties of iterative regularization methods.

2 Preliminaries

2.1 Formulation of the problem and basic assumptions

We assume that the inverse problem is given by $y=A\stackrel{~}{f}+\varepsilon .$  where $\varepsilon$  is a centered random variable satisfying the moment condition $IE\left(|\varepsilon {|}^{p}/{\sigma }^{p}\right)\le p!/2$  and $IE\left({\varepsilon }^{2}\right)={\sigma }^{2}$  .
We also need some notations concerning the fixed design settings, ${x}_{i},i=1,\dots ,n$  . Define the empirical measure:
${P}_{n}=\frac{1}{n}{\sum }_{i=1}^{n}{\delta }_{{x}_{i}}.$  and the associated empirical norm $\parallel y{\parallel }_{n}^{2}=\parallel y{\parallel }_{{P}_{n}}^{2}=\frac{1}{n}{\sum }_{i=1}^{n}\left(y\left({x}_{i}\right){\right)}^{2}$  as well as the empirical scalar product $}_{n}=\frac{1}{n}{\sum }_{i=1}^{n}{\varepsilon }_{i}y\left({t}_{i}\right).$  We assume $A:F\to Y$  , $F,Y$  separable Hilbert spaces. Let $〈,〉$  stands for the inner product defined over $F$  . We assume that the range of the operator $A$  , $\mathcal{ℛ}\left(A\right),$  is closed in the sense that $\stackrel{~}{f}\in \mathcal{N}\left(A{\right)}^{\perp }$  , where $\mathcal{N}\left(A{\right)}^{\perp }$  is the orthogonal complement of the null space of the operator $A$  .
With this notation let $J\left(f\right)=\parallel y-Af{\parallel }_{n}^{2}$  the quadratic risk function. We will denote by $\stackrel{^}{f}$  the function that minimizes the risk (which may not be unique), defined as
 $\begin{array}{c}\stackrel{^}{f}=arg{min}_{f\in F}J\left(f\right),\end{array}$ (2.1)
where the minimum is taken over all functions from $F$  to $Y$  .
The solution of the problem ${min}_{f\in F}J\left(f\right)$  exists if and only if $f$  is a solution of the normal equation
 $\begin{array}{c}{A}^{*}Af={A}^{*}y\end{array}$ (2.2)
where ${A}^{*}:Y\to F$  is the adjoint operator of $A$  (introduced via the requirement that for all $f\in F$  and $y\in Y$  , $〈Ax,y{〉}_{n}=\frac{1}{n}〈x,{A}^{*}y〉$  holds). It is important to remark that the operator ${A}^{*}$  actually depends on the observation sequence ${x}_{i},i=1,\dots n$  . If $Y$  is generated by $\left\{{\phi }_{j}{\right\}}_{j=1}^{m}$  and is such that this basis is orthonormal with respect to the ${L}^{2}\left({P}_{n}\right)$  norm over $Y$  , and $A$  is the identity then ${A}^{*}=\frac{1}{n}\left({\phi }_{j}\left({x}_{i}\right){\right)}_{i,j}$  .
It is necessary to mention that the convergence rates can thus be given only over subsets of $F$  , i.e., under a-priori assumptions on the exact solution $\stackrel{~}{f}$  . We will formulate such a-priori assumptions, encountered typically in the inverse problem literature, in terms of the exact solutions by considering subsets of $F$  given by some source condition of the form ${A}_{\mu ,\rho }=\left\{f\in F,f=\left({A}^{*}A{\right)}^{\mu }\omega ,\parallel \omega \parallel \le \rho \right\}$  where $0\le \mu \le {\mu }_{0}$  , ${\mu }_{0}>0$  and use the notation
 $\begin{array}{c}{A}_{\mu }={\bigcup }_{\rho >0}{A}_{\mu ,\rho }=\mathcal{ℛ}\left(\left({A}^{*}A{\right)}^{\mu }\right)\end{array}$ (2.3)
These sets are usually called source sets, $f\in {A}_{\mu ,\rho }$  is said to have a source representation.
The requirement that $f$  be in ${A}_{\mu ,\rho }$  can be considered as an smoothness condition.

2.2 Projection methods

For numerical calculations, we have to approximate the space $F$  by a finite-dimensional subspace. Estimating over all $F$  is in general not possible. One approach in this direction is regularization by projection, where the regularization is achieved by a finite-dimensional approximation through projection.
Let ${M}_{n}$  be a collection of index sets ( $m\in {M}_{n},m=\left\{{j}_{1}\dots ,{j}_{{d}_{m}}\right\}$  ). We give a sequence ${Y}_{1}\subset {Y}_{2}\dots \subset {Y}_{m}\dots \subset Y$  whose of union is dense in $Y$  . We assume $dim\left({Y}_{m}\right)={d}_{m}.$  Let ${\Pi }_{{Y}_{m}}^{n}$  be the orthogonal projector in the empirical norm over the subspace ${Y}_{m}$  and let ${A}_{m}={\Pi }_{{Y}_{m}}^{n}A$  . Define ${F}_{m}={A}^{*}{Y}_{m}$  , with ${A}_{m}^{*}:{Y}_{m}\to F$  the adjoint operator of ${A}_{m}$  , and ${\Pi }_{{F}_{m}}$  to be the orthogonal projector onto the subspace ${F}_{m}$  . Then, by construction ${\Pi }_{{F}_{m}}=\left({\Pi }_{{Y}_{m}}^{n}A{\right)}^{+}{\Pi }_{{Y}_{m}}^{n}A.$  Thus, we shall assume that data are give through an orthogonal design, corresponding to an orthogonal projection ${\Pi }_{{Y}_{m}}^{n}$  as
 $\begin{array}{c}{\Pi }_{{Y}_{m}}^{n}y={\Pi }_{{Y}_{m}}^{n}A\stackrel{~}{f}+{\Pi }_{{Y}_{m}}^{n}\varepsilon .\end{array}$ (2.4)
With this notation we have that the best-approximate ${L}^{2}$  solution has the expression ${\Pi }_{{F}_{m}}\stackrel{~}{f}={A}_{m}^{†}{y}_{m}.$  for ${y}_{m}={\Pi }_{{Y}_{m}}^{n}y$  in the domain of ${A}_{m}^{†}$  . In the following we shall denote ${\stackrel{~}{f}}_{m}={\Pi }_{{F}_{m}}\stackrel{~}{f}$  .
Our goal is to find the solution of the equation  1.1 in the finite-dimensional subspace ${F}_{m}$  of $F$  . We have that for projection without regularization the choice of ${F}_{m}$  and of ${A}_{m}$  has many advantages. For noisy data and severely ill-posed problems the dimension of the subspace has to be rather low to keep total error estimate small, since for these problems the smallest singular value of ${A}_{m}$  decreases rapidly as $m$  increases. To be able to use larger dimensions we have to combine the projection method with additional regularization methods,such as iterative methods [9,[11.

2.3 Singular value decomposition

As often ${A}_{m}$  is not of full rank, the singular value decomposition (SVD) of the operator ${A}_{m}$  is then a useful tool. Let $\left({\sigma }_{j};{\phi }_{j},{\phi }_{j}{\right)}_{j\in m}$  be a singular system for a linear operator ${A}_{m}$  , that is, ${A}_{m}{\phi }_{j}={\sigma }_{j}{\phi }_{j}$  and ${A}_{m}^{*}{\phi }_{j}={\sigma }_{j}{\phi }_{j}$  ; where $\left\{{\sigma }_{j}^{2}{\right\}}_{j\in m}$  are the nonzero eigenvalues of the selfadjoint operator ${A}_{m}^{*}{A}_{m}$  (and also of ${A}_{m}{A}_{m}^{*}$  ), considered in decreasing order. Furthermore, $\left\{{\phi }_{j}{\right\}}_{j\in m}$  and $\left\{{\phi }_{j}{\right\}}_{j\in m}$  are a corresponding complete orthonormal system of eigenvectors of ${A}_{m}^{*}{A}_{m}$  and ${A}_{m}{A}_{m}^{*}$  , respectively. For general linear operators with an SVD decomposition, we can write
 $\begin{array}{c}{A}_{m}f={\sum }_{j\in m}{\sigma }_{j}〈f,{\phi }_{j}〉{\phi }_{j}\end{array}$ (2.5)
 $\begin{array}{c}{A}_{m}^{*}{y}_{m}={\sum }_{j\in m}{\sigma }_{j}〈{y}_{m},{\phi }_{j}〉{\phi }_{j}.\end{array}$ (2.6)
For ${y}_{m}$  in the domain of ${A}_{m}^{†}$  , $\mathcal{D}\left({A}_{m}^{†}\right)$  , the best-approximate ${L}^{2}$  solution hast the expression ${A}_{m}^{†}{y}_{m}={\sum }_{j\in m}\frac{〈{y}_{m},{\phi }_{j}〉}{{\sigma }_{j}}{\phi }_{j}={\sum }_{j\in m}\frac{〈{A}_{m}^{*}{y}_{m},{\phi }_{j}〉}{{\sigma }_{j}^{2}}{\phi }_{j}.$  Note that for large $j$  , the term $1/{\sigma }_{j}$  grows to infinity. Thus, the high frequency errors are strongly amplified. We will asume that ${\sigma }_{j}=\mathcal{O}\left({j}^{-p}\right)$  for some $p>1/2$  , which is clearly related to the ill-posedness of the operator $A$  and the approximation properties of ${Y}_{m}$  . For the construction and analysis of regularization methods, we will require some general notation for functions of the operators ${A}_{m}^{*}{A}_{m}$  and ${A}_{m}{A}_{m}^{*}$  .
Let ${E}_{\lambda }$  be the spectral decomposition of ${A}_{m}^{*}{A}_{m}$  given by ${E}_{\lambda }\left(\cdot \right)={\sum }_{{\sigma }^{2}<{\lambda }_{j},j\in m}〈\cdot ,{\phi }_{j}〉{\phi }_{j}$  and ${H}_{\lambda }$  the spectral decomposition of ${A}_{m}{A}_{m}^{*}$  . Then ${E}_{\lambda }$  is an orthogonal projector and projects onto $span\left\{{\phi }_{j}|j\in m,{\sigma }^{2}<\lambda \right\}.$  Since $\left({\sigma }_{j}^{2};{\phi }_{j}\right)$  is an eigensystem for the selfadjoint compact operator ${A}_{m}^{*}{A}_{m}$  , ${A}_{m}^{*}{A}_{m}f={\sum }_{j\in m}{\sigma }_{j}^{2}〈f,{\phi }_{j}〉{\phi }_{j}$  holds, which will be written (using the definition of the integral below) as ${A}_{m}^{*}{A}_{m}f=\int \lambda d{E}_{\lambda }f$  for $f\in \mathcal{D}\left({A}_{m}\right)$  . Here the limits of integration could be 0 and $\parallel {A}_{m}{\parallel }^{2}+\epsilon$  for any $\epsilon >0$  . We sometimes omit the limits of integration.
This, motivates the definition
 $\begin{array}{c}G\left({A}_{m}^{*}{A}_{m}\right):=\int G\left(\lambda \right)d{E}_{\lambda }:={\sum }_{{\sigma }^{2}={\lambda }_{j},j\in m}G\left({\sigma }_{j}^{2}\right)〈\cdot ,{\phi }_{j}〉{\phi }_{j}\end{array}$ (2.7)
of a (piecewise) continuous function $G$  of a selfadjoint linear operator on ${F}_{m}$  . If ${A}^{*}A$  is continuously invertible, then $\left({A}^{*}A{\right)}^{-1}=\int \frac{1}{\lambda }d{E}_{\lambda }$  .
In this case the best-approximate ${L}^{2}$  solution, for ${y}_{m}$  in the domain of ${A}_{m}^{†}$  , can be characterized by the equation
 $\begin{array}{c}{\stackrel{~}{f}}_{m}={A}_{m}^{†}{y}_{m}=\int \frac{1}{\lambda }d{E}_{\lambda }{A}_{m}^{*}{y}_{m}.\end{array}$ (2.8)
If $G\left({A}_{m}^{*}{A}_{m}\right)$  is defined via  2.7 , then for $f\in \mathcal{D}\left({G}_{1}\left({A}_{m}^{*}{A}_{m}\right)\right)$  and $g\in \mathcal{D}\left({G}_{2}\left({A}_{m}^{*}{A}_{m}\right)\right)$
 $\begin{array}{c}〈{G}_{1}\left({A}_{m}^{*}{A}_{m}\right)f,{G}_{2}\left({A}_{m}^{*}{A}_{m}\right)g〉=\int {G}_{1}\left(\lambda \right){G}_{2}\left(\lambda \right)d〈{E}_{\lambda }f,g〉\end{array}$ (2.9)
and
 $\begin{array}{c}\parallel G\left({A}_{m}^{*}{A}_{m}\right)f{\parallel }^{2}=\int G\left(\lambda \right)d\parallel {E}_{\lambda }f{\parallel }^{2}.\end{array}$ (2.10)
The source set, ${A}_{\mu }$   2.3 , can be characterized via the singular values as follows:
$\left({A}_{m}^{*}{A}_{m}{\right)}^{\mu }\omega =\int {\lambda }^{\mu }d{E}_{\lambda }\omega ={\sum }_{j\in m}{\sigma }_{j}^{2\mu }〈\omega ,{\phi }_{j}〉{\phi }_{j}$

3 Regularization methods

After the general considerations of the last section, we now explain the construction of a regularization method for the important special case of selfadjoint linear operators. The basic idea for deriving a regularization method is to replace the amplification factors $1/{\lambda }_{j}$  by a filtered version $Q\left({\lambda }_{j},\alpha \right)$  , where the filter function is a piecewise continuous, nonnegative and nonincreasing function of $\lambda$  on the segment $\left[0,\parallel {A}_{m}{\parallel }^{2}\right]$  for a regularization parameter $\alpha >0$  .
The assumptions over the regularizing coefficients ${Q}_{\alpha }\left(\lambda \right)$  are technical and are given in order to control fluctuations over set $\left[0,\parallel {A}_{m}{\parallel }^{2}\right]$  .
The filter family $\left\{Q\left({\lambda }_{j},\alpha \right){\right\}}_{j\in m}$  approximates the function ${\lambda }^{-1}$  for $\alpha \to \infty$  . Intuitively, a regularization on ${A}_{m}^{†}$  should then be the replacement of the ill conditioned operator ${A}_{m}^{†}$  by a family $\left\{R\left({\lambda }_{j},\alpha \right){\right\}}_{j\in m}:{Y}_{m}\to {F}_{m}$  of continuous operators. Throughout all the article, we shall denote $\left\{Q\left({\lambda }_{j},\alpha \right){\right\}}_{j\in m}$  and $\left\{R\left({\lambda }_{j},\alpha \right){\right\}}_{j\in m}$  by ${Q}_{\alpha }$  and ${R}_{\alpha }$  , respectively. Obviously, for all $\alpha >0,$  ${R}_{\alpha }$  is bounded.
As the approximation of ${\stackrel{~}{f}}_{m}$  , we then take ${f}_{m,\alpha }={Q}_{\alpha }\left({A}_{m}^{*}{A}_{m}\right){A}_{m}^{*}{y}_{m}={R}_{\alpha }{y}_{m},$  where ${R}_{\alpha }:=\int {Q}_{\alpha }\left(\lambda \right)d{E}_{\lambda }{A}_{m}^{*}$  .
Remark 3.1. Note that with the above notation
 $\begin{array}{c}{f}_{m,\alpha }={R}_{\alpha }{y}_{m}={R}_{\alpha }{\Pi }_{{Y}_{m}}^{n}A\stackrel{~}{f}+{R}_{\alpha }{\Pi }_{{Y}_{m}}^{n}\varepsilon .\end{array}$ (3.1)
Also that we can write ${R}_{\alpha }={Q}_{\alpha }\left({A}_{m}^{*}{A}_{m}\right){A}_{m}^{*}.$
The next theorem gives conditions under which the first term in  3.1 converges to ${\stackrel{~}{f}}_{m}={\Pi }_{{F}_{m}}\stackrel{~}{f}$  . The proof follows that of [9, but we include it for the sake of completeness.
Theorem 3.2. Let, for all $\alpha >0$  , ${Q}_{\alpha }:\left[0,\parallel {A}_{m}{\parallel }^{2}\right]\to IR$  be a piecewise continuous and nonincreasing function of $\lambda$  on the segment $\left[0,\parallel {A}_{m}{\parallel }^{2}\right]$  . Assume also that there is a $C>0$  such that $|\lambda {Q}_{\alpha }\left(\lambda \right)|\le C,$  and ${lim}_{\alpha \to \infty }{Q}_{\alpha }\left(\lambda \right)=\frac{1}{\lambda }$  for all $\lambda \in \left(0,\parallel {A}_{m}{\parallel }^{2}\right).$  Then, for all $y\in D\left({A}_{m}^{†}\right)$  , ${lim}_{\alpha \to \infty }{Q}_{\alpha }\left({A}_{m}^{*}{A}_{m}\right){A}_{m}^{*}{\Pi }_{{Y}_{m}}^{n}A\stackrel{~}{f}={\stackrel{~}{f}}_{m}$  holds with ${\stackrel{~}{f}}_{m}={A}_{m}^{†}{y}_{m}.$
Remark 3.3. In order to assume convergence as $\alpha \to \infty$  , it is necessary to choose ${Q}_{\alpha }$  such that it approximates $1/\lambda$  for all $\lambda \in \left(0,\parallel {A}_{m}{\parallel }^{2}\right]$  . Also, note that the condition $|\lambda {Q}_{\alpha }\left(\lambda \right)|\le C$  implies that $\parallel {A}_{m}{R}_{\alpha }\parallel =\parallel {A}_{m}{A}_{m}^{*}{Q}_{\alpha }\left({A}_{m}^{*}{A}_{m}\right)\parallel \le C$  , i.e, $\parallel {A}_{m}{R}_{\alpha }\parallel$  is uniformly bounded.
• Proof. As in [9, if ${\stackrel{~}{f}}_{m}$  is defined by  2.8 , then by  2.2 the residual norm has the representation $\parallel {\stackrel{~}{f}}_{m}-{Q}_{\alpha }\left({A}_{m}^{*}{A}_{m}\right){A}_{m}^{*}{A}_{m}\stackrel{~}{f}{\parallel }^{2}=\parallel \left(I-{Q}_{\alpha }\left({A}_{m}^{*}{A}_{m}\right){A}_{m}^{*}{A}_{m}\right){\stackrel{~}{f}}_{m}{\parallel }^{2}$  From the formula  2.10 , it follows that $\parallel {\stackrel{~}{f}}_{m}-{Q}_{\alpha }\left({A}_{m}^{*}{A}_{m}\right){A}_{m}^{*}{A}_{m}\stackrel{~}{f}{\parallel }^{2}={\int }_{0}^{\parallel {A}_{m}{\parallel }^{2}+}\left(1-\lambda {Q}_{\alpha }\left(\lambda \right){\right)}^{2}d\parallel {E}_{\lambda }{\stackrel{~}{f}}_{m}{\parallel }^{2}.$  Since $\left(1-\lambda {Q}_{\alpha }\left(\lambda \right){\right)}^{2}$  is bounded by the constant $\left(1+C{\right)}^{2}$  , which is integrable with respect to the measure $d\parallel {E}_{\lambda }{\stackrel{~}{f}}_{m}{\parallel }^{2},$  then by the Dominated Convergence Theorem,  $\begin{array}{c}{lim}_{\alpha \to \infty }{\int }_{0}^{\parallel {A}_{m}{\parallel }^{2}+}\left(1-\lambda {Q}_{\alpha }\left(\lambda \right){\right)}^{2}d\parallel {E}_{\lambda }{\stackrel{~}{f}}_{m}{\parallel }^{2}={\int }_{0}^{\parallel {A}_{m}{\parallel }^{2}+}{lim}_{\alpha \to \infty }\left(1-\lambda {Q}_{\alpha }\left(\lambda \right){\right)}^{2}d\parallel {E}_{\lambda }{\stackrel{~}{f}}_{m}{\parallel }^{2}.\end{array}$ (3.2)
Since for $\lambda >0$  , ${lim}_{\alpha \to \infty }\left(1-\lambda {Q}_{\alpha }\left(\lambda \right)\right)=0$  then the integral on the right-hand side of  3.2 equals to $0$  . On the other hand, if $\lambda =0$  , ${lim}_{\alpha \to \infty }\left(1-\lambda {Q}_{\alpha }\left(\lambda \right)\right)=1$  then the equation  3.2 has the form  $\begin{array}{c}{lim}_{\alpha \to \infty }{\int }_{0}^{\parallel {A}_{m}{\parallel }^{2}+}\left(1-\lambda {Q}_{\alpha }\left(\lambda \right){\right)}^{2}d\parallel {E}_{\lambda }{\stackrel{~}{f}}_{m}{\parallel }^{2}={lim}_{\lambda \to {0}^{+}}\parallel {E}_{\lambda }{\stackrel{~}{f}}_{m}{\parallel }^{2}-\parallel {E}_{0}{\stackrel{~}{f}}_{m}{\parallel }^{2}\end{array}$ (3.3)
which is equal the jump of $\parallel {E}_{\lambda }{\stackrel{~}{f}}_{m}{\parallel }^{2}$  at $\lambda =0$  . Since ${\stackrel{~}{f}}_{m}\in \mathcal{N}\left({A}_{m}{\right)}^{\perp }$  , the term on the right-hand side of  3.3 equals to 0. Thus, ${R}_{\alpha }{A}_{m}\stackrel{~}{f}$  converges to ${\stackrel{~}{f}}_{m}$  as $\alpha \to \infty$  for ${y}_{m}\in D\left({A}_{m}^{†}\right),$  which ends the proof.
Let $Tr\left(B\right)$  the trace of the selfadjoint operator ${B}^{t}B$  for any square matriz $B$  , which is defined by $Tr\left(B\right)=\frac{1}{n}{\sum }_{j\in m}{b}_{j}$  for ${b}_{j}$  eigenvalues of ${B}^{t}B$  .
We then have the following result,
Theorem 3.4. Let ${Q}_{\alpha }$  be as in theorem  3.2 . Let $\mu ,\rho >0$  and let ${\omega }_{\mu }:\left(0,{\alpha }_{0}\right)\to R$  be such that for all $\alpha \in \left(0,{\alpha }_{0}\right)$  and $\lambda \in \left[0,{\sigma }_{1}^{2}\right]$  ${sup}_{0\le \lambda \le {\sigma }_{1}^{2}}{\lambda }^{\mu }|1-\lambda {Q}_{\alpha }\left(\lambda \right)|\le {\omega }_{\mu }\left(\alpha \right)$  holds. Then for $\stackrel{~}{f}\in {A}_{\mu ,\rho },$  the following inequality holds true
 $\begin{array}{c}IE\parallel {\stackrel{~}{f}}_{m}-{f}_{\alpha ,m}{\parallel }^{2}\le 2{\omega }_{\mu }\left(\alpha {\right)}^{2}{\rho }^{2}+2{\sigma }^{2}Tr\left({Q}_{\alpha }^{2}\left({A}_{m}^{*}{A}_{m}\right){A}_{m}{A}_{m}^{*}\right).\end{array}$ (3.4)
• Proof. The proof of this inequality is based on the definition of the estimator ${f}_{m,\alpha }$  and on the assumptions over this function. We have that the ${L}^{2}-$  norm of the difference between the regularized function and the true data function can be bounded by  $\begin{array}{c}IE\parallel {\stackrel{~}{f}}_{m}-{f}_{m,\alpha }{\parallel }^{2}\le 2IE\parallel {\stackrel{~}{f}}_{m}-{R}_{\alpha }{A}_{m}\stackrel{~}{f}{\parallel }^{2}+2IE\parallel {R}_{\alpha }{A}_{m}\stackrel{~}{f}-{f}_{m,\alpha }{\parallel }^{2}\end{array}$ (3.5)
where ${f}_{m,\alpha }={R}_{\alpha }{y}_{m}$  .
This is the typical bias-variance decomposition. The first term on the right-hand side is an approximation error, which corresponds to the bias, whereas the second term, variance, is a stability bound on the regularizing operator ${R}_{\alpha }$  . Note that by the Theorem  3.2 , the first term in  3.5 goes to 0 if ${y}_{m}\in \mathcal{D}\left({T}_{m}^{†}\right)$  .
Let $\omega \in {F}_{m}$  with $\parallel \omega \parallel \le \rho$  . Since $\stackrel{~}{f}\in {F}_{m}$  then ${\Pi }_{{F}_{m}}\stackrel{~}{f}=\left({A}_{m}^{*}{A}_{m}{\right)}^{\mu }\omega$  . On the other hand, ${\lambda }^{\mu }{sup}_{\lambda }|1-\lambda {Q}_{\alpha }\left(\lambda \right)|\le {\omega }_{\mu }\left(\alpha \right)$  , then the first term in this equation can be bounded by
 $\begin{array}{cc}IE\parallel {\stackrel{~}{f}}_{m}-{R}_{\alpha }{A}_{m}\stackrel{~}{f}{\parallel }^{2}& =IE\parallel {\stackrel{~}{f}}_{m}-{Q}_{\alpha }\left({A}_{m}^{*}{A}_{m}\right){A}_{m}^{*}{A}_{m}{\stackrel{~}{f}}_{m}{\parallel }^{2}\end{array}$
 $\begin{array}{cc}& =IE\parallel \left(I-{Q}_{\alpha }\left({A}_{m}^{*}{A}_{m}\right){A}_{m}^{*}{A}_{m}\right){\stackrel{~}{f}}_{m}{\parallel }^{2}\end{array}$
 $\begin{array}{cc}& \le {\omega }_{\mu }^{2}{\rho }^{2}.\end{array}$
 $\begin{array}{}\end{array}$
In order to control the term corresponding to the variance we used that the data perturbation is white noise. Thus,
 $\begin{array}{cc}IE\parallel {R}_{\alpha }{A}_{m}\stackrel{~}{f}-{f}_{m,\alpha }{\parallel }^{2}& =IE〈\varepsilon ,\left({Q}_{\alpha }\left({A}_{m}^{*}{A}_{m}\right){A}_{m}^{*}{\right)}^{*}{Q}_{\alpha }\left({A}_{m}^{*}{A}_{m}\right){A}_{m}^{*}\varepsilon 〉\end{array}$
 $\begin{array}{cc}& =IE〈\varepsilon ,{Q}_{\alpha }^{2}\left({A}_{m}^{*}{A}_{m}\right){A}_{m}{A}_{m}^{*}\varepsilon 〉\end{array}$
 $\begin{array}{cc}& ={\sigma }^{2}Tr\left({Q}_{\alpha }^{2}\left({A}_{m}^{*}{A}_{m}\right){A}_{m}{A}_{m}^{*}\right)\end{array}$
 $\begin{array}{}\end{array}$
which yields the desired result.
The next result will be useful when studying iterative methods.
Theorem 3.5. Let ${Q}_{\alpha }$  be as in theorem  3.2 . Assume also that ${Q}_{\alpha }$  is continuously differentiable and that the function $|1-\lambda {Q}_{\alpha }\left(\lambda \right){|}^{\prime }|\lambda {Q}_{\alpha }\left(\lambda \right)-1{|}^{-1}$  doest not decrease. Then the estimates are valid ${sup}_{0\le \lambda \le {\sigma }_{1}^{2}}|{Q}_{\alpha }\left(\lambda \right)|={Q}_{\alpha }\left(0\right),$  and ${sup}_{0\le \lambda \le {\sigma }_{1}^{2}}{\lambda }^{\mu }|1-\lambda {Q}_{\alpha }\left(\lambda \right)|<{\mu }^{\mu }\left(\mu +1{\right)}^{-1}{\omega }_{\mu }\left(\alpha \right)$  where ${\omega }_{\mu }\left(\alpha \right)={Q}_{\alpha }\left(0{\right)}^{-\mu }.$
• Proof. The proof can be carried out by standard techniques. A proof of this result can be found in [11.

4 Rates of convergence for the regularized estimator

In any regularization method, the regularization parameter $\alpha$  plays a crucial role. For choosing the parameter, there are general methods of parameter selection. For example, the Discrepancy Principe [20, Cross-Validation [7and the L-curve [10. They differ in the amount of a priori information required as well as in the decision criteria. The appropriate choice of regularization parameter is a difficult problem. We would like too choose $\alpha$  , based on the data in such a way that optimal rates are maintained. This choice should not depend on a priori regularity assumptions.
Our goal is to introduce adaptive methods in the context of statistical inverse problems.
In this section we introduce our adaptive estimator, for a fixed $m={m}_{0}$  . We choose ${m}_{0}$  such that $\parallel \stackrel{~}{f}-{\Pi }_{{F}_{{m}_{0}}}\stackrel{~}{f}{\parallel }^{2}$  satisfies the optimal rates with high probability since we know $\parallel \stackrel{~}{f}-{\Pi }_{{F}_{{m}_{0}}}\stackrel{~}{f}{\parallel }^{2}<\parallel I-{\Pi }_{{Y}_{{m}_{0}}}{\parallel }^{4\mu }=\mathcal{O}\left({d}_{{m}_{0}}^{-4\mu p}\right)$  for a certain $p$  and $0<\mu \le 1/2$  . It is satisfied if the dimension of the set is such that
 $\begin{array}{c}{d}_{{m}_{0}}\ge {n}^{\frac{1}{2p+1}}.\end{array}$ (4.1)
This leads to the rate $\parallel \stackrel{~}{f}-{f}_{m,\alpha }{\parallel }^{2}=\mathcal{O}\left({n}^{-\frac{4\mu p}{4\mu p+2p+1}}\right).$  Analogous results are obtained in the case of Hilbert scales ([12,[16).
Adaptive model selection is a technique which penalizes the regularization parameter, in such a way that we choose ${\stackrel{^}{f}}_{{m}_{0},{\alpha }_{\stackrel{^}{k}}}$  by minimizing $arg{min}_{k\in \mathcal{K},f\in {F}_{m}}\left(\parallel {R}_{{\alpha }_{k}}\left({y}_{m}-{A}_{m}f\right){\parallel }^{2}+pen\left({\alpha }_{k}\right)\right)$  where $\stackrel{^}{k}=arg{min}_{k\in \mathcal{K}}\left(\parallel {R}_{{\alpha }_{k}}\left({y}_{m}-{A}_{m}f\right){\parallel }^{2}+pen\left({\alpha }_{k}\right)\right)$  and $pen\left({\alpha }_{k}\right)=r{\sigma }^{2}\left(1+{L}_{k}\right)\left[Tr\left({R}_{{\alpha }_{k}}^{t}{R}_{{\alpha }_{k}}\right)+{\rho }^{2}\left({R}_{{\alpha }_{k}}\right)\right],$  with $r>2$  and ${L}_{k}$  is a sequence which is incorporated in order to control the complexity of the set $\mathcal{K}=\left\{1,2,\dots ,{k}_{n}\right\}$  , of all possible index up to ${k}_{n}$  . Here ${\rho }^{2}\left(B\right)=\rho \left({B}^{t}B\right)$  is the spectral radius of the selfadjoint operator ${B}^{t}B$  for any square matriz $B$  , which is defined by ${\rho }^{2}\left(B\right)=\frac{1}{n}{max}_{j\in m}{b}_{j}$  for ${b}_{j}$  eigenvalues of ${B}^{t}B$  .
Thus, $\stackrel{^}{k}$  is selected by minimizing
 $\begin{array}{c}arg{min}_{k\in \mathcal{K}}\left(\parallel {R}_{{\alpha }_{k}}\left({y}_{m}-{A}_{m}f\right){\parallel }^{2}+\frac{r{\sigma }^{2}\left(1+{L}_{k}\right)}{n}\left[{\sum }_{j\in m}{Q}^{2}\left({\lambda }_{j}\right){\lambda }_{j}+{max}_{j\in m}{Q}^{2}\left({\lambda }_{j}\right){\lambda }_{j}\right]\right).\end{array}$ (4.2)
The strategy as proposed in this article automatically provides the optimal order of accuracy.
The regularized estimator has a rate of convergence less or equal than the best rate achieved by the best estimator for a selected model. We have the following result,
Theorem 4.1. For any $f\in {F}_{m}$  and any ${\alpha }_{k}$  the following inequality holds true for $d$  a positive constant that depends on r (as in Lemma  4.4 ),
 $\begin{array}{c}IE\parallel {\stackrel{~}{f}}_{m}-{\stackrel{^}{f}}_{{\alpha }_{\stackrel{^}{k}}}{\parallel }^{2}\le \frac{1}{\left(1-\nu \right)}{inf}_{k\in \mathcal{K}}\left[C\left(1+\nu \right)\parallel {\stackrel{~}{f}}_{m}-{f}_{{\alpha }_{k}}{\parallel }^{2}+2pen\left({\alpha }_{k}\right)\right]+\frac{{C}_{1}\left(d\right)}{n}\end{array}$ (4.3)
where ${C}_{1}\left(d\right)=4{\sigma }^{2}{\sum }_{k}\frac{n{\rho }^{2}\left({R}_{{\alpha }_{k}}\right)}{d}\left[\sqrt{dr{L}_{k}\left[\frac{Tr\left({R}_{{\alpha }_{k}}^{t}{R}_{{\alpha }_{k}}\right)}{{\rho }^{2}\left({R}_{{\alpha }_{k}}\right)}+1\right]}+1\right]{e}^{-\sqrt{dr{L}_{k}\left[\frac{Tr\left({R}_{{\alpha }_{k}}^{t}{R}_{{\alpha }_{k}}\right)}{{\rho }^{2}\left({R}_{{\alpha }_{k}}\right)}+1\right]}}.$
Remark 4.2. An important issue is that equation  4.3 is non asymptotic. The goodness of fit of the estimator is defined by trace, $Tr\left({R}_{\alpha }^{t}{R}_{\alpha }\right)$  , and spectral radius, ${\rho }^{2}\left({R}_{\alpha }\right)$  . Also, the estimator is optimal in the sense that the adaptive estimator achieves the best rate of convergence among all the regularized estimators.
Remark 4.3. Remark that under our assumptions, namely that the basis is orthonormal for the fixed design, both $n{\rho }^{2}\left({R}_{k}\right)$  and $Tr\left({R}_{k}^{t}{R}_{k}\right)/{\rho }^{2}\left({R}_{k}\right)$  do not depend on n.
• Proof. For any ${f}_{{\alpha }_{k}}$  and any $k\in IN$   $\begin{array}{ccc}& \parallel {R}_{{\alpha }_{\stackrel{^}{k}}}\left({y}_{m}-{A}_{m}{\stackrel{^}{f}}_{{\alpha }_{\stackrel{^}{k}}}\right){\parallel }^{2}+pen\left({\alpha }_{\stackrel{^}{k}}\right)\le \parallel {R}_{{\alpha }_{k}}\left({y}_{m}-{A}_{m}{f}_{{\alpha }_{k}}\right){\parallel }^{2}+pen\left({\alpha }_{k}\right)& \end{array}$
and  $\begin{array}{ccc}& \parallel {R}_{{\alpha }_{k}}\left({y}_{m}-{A}_{m}{f}_{{\alpha }_{k}}\right){\parallel }^{2}=\parallel {R}_{{\alpha }_{k}}{A}_{m}\left(\stackrel{~}{f}-{f}_{{\alpha }_{k}}\right){\parallel }^{2}+2〈{R}_{{\alpha }_{k}}{A}_{m}\left(\stackrel{~}{f}-{f}_{{\alpha }_{k}}\right),{R}_{{\alpha }_{k}}{\Pi }_{{Y}_{m}}^{n}\varepsilon 〉+\parallel {R}_{{\alpha }_{k}}{\Pi }_{{Y}_{m}}^{n}\varepsilon {\parallel }^{2}& \end{array}$
Thus, following standard arguments we have  $\begin{array}{ccc}& & \parallel {R}_{{\alpha }_{\stackrel{^}{k}}}{A}_{m}\left(\stackrel{~}{f}-{\stackrel{^}{f}}_{{\alpha }_{\stackrel{^}{k}}}\right){\parallel }^{2}\end{array}$
 $\begin{array}{ccc}& \le & \parallel {R}_{{\alpha }_{k}}{A}_{m}\left(\stackrel{~}{f}-{f}_{{\alpha }_{k}}\right){\parallel }^{2}-2<{R}_{{\alpha }_{\stackrel{^}{k}}}{A}_{m}\left(\stackrel{~}{f}-{\stackrel{^}{f}}_{{\alpha }_{\stackrel{^}{k}}}\right),{R}_{{\alpha }_{\stackrel{^}{k}}}{\Pi }_{{Y}_{m}}^{n}\varepsilon >\end{array}$
 $\begin{array}{ccc}& & +2<{R}_{{\alpha }_{k}}{A}_{m}\left(\stackrel{~}{f}-{f}_{{\alpha }_{k}}\right),{R}_{{\alpha }_{k}}{\Pi }_{{Y}_{m}}^{n}\varepsilon >-\parallel {R}_{{\alpha }_{\stackrel{^}{k}}}{\Pi }_{{Y}_{m}}^{n}\varepsilon {\parallel }^{2}+\parallel {R}_{{\alpha }_{k}}{\Pi }_{{Y}_{m}}^{n}\varepsilon {\parallel }^{2}+pen\left({\alpha }_{k}\right)+pen\left({\alpha }_{\stackrel{^}{k}}\right).\end{array}$
Let $0<\nu <1$  . Since the algebraic inequality $2ab\le \nu {a}^{2}+\frac{1}{\nu }{b}^{2}$  holds for all $a,b\in IR$  , we find that  $\begin{array}{ccc}& & \left(1-\nu \right)\parallel {R}_{{\alpha }_{\stackrel{^}{k}}}{A}_{m}\left(\stackrel{~}{f}-{\stackrel{^}{f}}_{{\alpha }_{\stackrel{^}{k}}}\right){\parallel }^{2}\end{array}$
 $\begin{array}{ccc}& \le & \left(1+\nu \right)\parallel {R}_{{\alpha }_{k}}{A}_{m}\left(\stackrel{~}{f}-{f}_{{\alpha }_{k}}\right){\parallel }^{2}+2pen\left({\alpha }_{k}\right)+2{sup}_{{\alpha }_{k}}\left\{\frac{1}{\nu }\parallel {R}_{{\alpha }_{k}}{\Pi }_{{Y}_{m}}^{n}\varepsilon {\parallel }^{2}-pen\left({\alpha }_{k}\right)\right\},\end{array}$
holds for any $k$  and ${f}_{{\alpha }_{k}}\in {F}_{m}$  .
On the other hand, using that is $1\le \parallel {R}_{{\alpha }_{k}}A\parallel \le C$  , we have that for any ${f}_{{\alpha }_{k}}\in {F}_{{m}_{0}}$  and any $k\in IN$  ,  $\begin{array}{ccc}& & \left(1-\nu \right)\parallel {\stackrel{~}{f}}_{m}-{\stackrel{^}{f}}_{{\alpha }_{\stackrel{^}{k}}}{\parallel }^{2}\le C\left(1+\nu \right)\parallel {\stackrel{~}{f}}_{m}-{f}_{{\alpha }_{k}}{\parallel }^{2}\end{array}$
 $\begin{array}{ccc}& +& 2pen\left({\alpha }_{k}\right)+2{C}_{1}{sup}_{{\alpha }_{k}}\left\{\parallel {R}_{{\alpha }_{k}}{\Pi }_{{Y}_{m}}^{n}\varepsilon {\parallel }^{2}-pen\left({\alpha }_{k}\right)\right\}.\end{array}$
The proof then follows directly from the following technical lemma ([3,[16) which characterizes the supremum of an empirical process by the regularization family.
Lemma 4.4. Let $\eta \left(A\right)=\sqrt{{\varepsilon }^{t}{A}^{t}A\varepsilon }=\parallel A\varepsilon \parallel .$  Then, there exists a positive constant $d$  that depends on $r/2$  such that the following inequality holds  $\begin{array}{ccc}& & P\left({\eta }^{2}\left(A\right)\ge {\sigma }^{2}\left[Tr\left({A}^{t}A\right)+\rho \left({A}^{t}A\right)\right]r/2\left(1+L\right)+{\sigma }^{2}u\right)\end{array}$ (4.4)
 $\begin{array}{ccc}& & \le exp\left\{-\sqrt{d\left(1/\rho \left({A}^{t}A\right)u+r/2L\left[Tr\left({A}^{t}A\right)/\rho \left({A}^{t}A\right)+1\right]\right)}\right\}.\end{array}$
With the above notation, $\eta \left({R}_{{\alpha }_{k}}\right)=\parallel {R}_{k}{\varepsilon }_{m}\parallel$  where ${\varepsilon }_{m}={\Pi }_{{Y}_{m}}^{n}\varepsilon$  .
Now, with this lemma we have  $\begin{array}{ccc}& & P\left({sup}_{{\alpha }_{k}}\parallel {R}_{{\alpha }_{k}}{\varepsilon }_{m}{\parallel }^{2}-pen\left({\alpha }_{k}\right)>{\sigma }^{2}x\right)\end{array}$
 $\begin{array}{ccc}& \le & {\sum }_{k}P\left[{\eta }^{2}\left({R}_{{\alpha }_{k}}\right)\ge r{\sigma }^{2}\left(1+{L}_{k}\right)\left[Tr\left({R}_{{\alpha }_{k}}^{t}{R}_{{\alpha }_{k}}\right)+{\rho }^{2}\left({R}_{{\alpha }_{k}}\right)\right]+{\sigma }^{2}x\right]\end{array}$
 $\begin{array}{ccc}& \le & {\sum }_{k}exp\left\{-\sqrt{d\left(1/{\rho }^{2}\left({R}_{{\alpha }_{k}}\right)x+r{L}_{k}\left[Tr\left({R}_{{\alpha }_{k}}^{t}{R}_{{\alpha }_{k}}\right)/{\rho }^{2}\left({R}_{{\alpha }_{k}}\right)+1\right]\right)}\right\}\end{array}$
Since for $X$  positive $IEX={\int }_{0}^{\infty }P\left(X>u\right)du$  , we then have that  $\begin{array}{ccc}& & IE\left[{sup}_{{\alpha }_{k}}\parallel {R}_{{\alpha }_{k}}{\varepsilon }_{m}{\parallel }^{2}-pen\left({\alpha }_{k}\right)\right]={\int }_{0}^{\infty }P\left[{sup}_{{\alpha }_{k}}\parallel {R}_{{\alpha }_{k}}{\varepsilon }_{m}{\parallel }^{2}-pen\left({\alpha }_{k}\right)\ge x\right]dx\end{array}$
 $\begin{array}{ccc}& =& {\sigma }^{2}{\int }_{0}^{\infty }P\left[{sup}_{{\alpha }_{k}}\parallel {R}_{{\alpha }_{k}}{\varepsilon }_{m}{\parallel }^{2}-pen\left({\alpha }_{k}\right)\ge {\sigma }^{2}u\right]du\end{array}$
 $\begin{array}{ccc}& \le & {\sigma }^{2}{\sum }_{k}{\int }_{0}^{\infty }exp\left\{-\sqrt{{k}_{1}u+{k}_{2}}\right\}du.\end{array}$
where ${k}_{1}=d/{\rho }^{2}\left({R}_{{\alpha }_{k}}\right)$  and ${k}_{2}=dr{L}_{k}\left[Tr\left({R}_{{\alpha }_{k}}^{T}{R}_{{\alpha }_{k}}\right)/{\rho }^{2}\left({R}_{{\alpha }_{k}}\right)+1\right]$  .
Let $w={k}_{1}u+{k}_{2},$  then
 $\begin{array}{cc}IE\left[{sup}_{{\alpha }_{k}}\parallel {R}_{{\alpha }_{k}}{\varepsilon }_{m}{\parallel }^{2}-pen\left({\alpha }_{k}\right)\right]& \le {\sigma }^{2}{\sum }_{k}{\int }_{{k}_{2}}^{\infty }\frac{1}{{k}_{1}}exp\left\{-\sqrt{w}\right\}dw\end{array}$
 $\begin{array}{cc}& ={\sigma }^{2}{\sum }_{k}\frac{2}{{k}_{1}}\left[-\sqrt{{k}_{2}}+1\right]exp\left\{-\sqrt{{k}_{2}}\right\}\end{array}$
 $\begin{array}{}\end{array}$
Finally, we have the desired result.
 $\begin{array}{cc}IE\parallel {\stackrel{~}{f}}_{m}-{\stackrel{^}{f}}_{{\alpha }_{\stackrel{^}{k}}}{\parallel }^{2}\le \frac{1}{\left(1-\nu \right)}{inf}_{k\in \mathcal{K}}\left[C\left(1+\nu \right)\parallel {\stackrel{~}{f}}_{m}-{f}_{{\alpha }_{k}}{\parallel }^{2}+2pen\left({\alpha }_{k}\right)\right]+\frac{{C}_{1}\left(d\right)}{n}& \end{array}$
 $\begin{array}{}\end{array}$
where ${C}_{1}\left(d\right)=4{\sigma }^{2}{\sum }_{k}\frac{n{\rho }^{2}\left({R}_{{\alpha }_{k}}\right)}{d}\left[\sqrt{dr{L}_{k}\left[\frac{Tr\left({R}_{{\alpha }_{k}}^{t}{R}_{{\alpha }_{k}}\right)}{{\rho }^{2}\left({R}_{{\alpha }_{k}}\right)}+1\right]}+1\right]{e}^{-\sqrt{dr{L}_{k}\left[\frac{Tr\left({R}_{{\alpha }_{k}}^{t}{R}_{{\alpha }_{k}}\right)}{{\rho }^{2}\left({R}_{{\alpha }_{k}}\right)}+1\right]}},$

5 Regularization by iterative methods

Iterative regularization methods, are very competitive methods for linear inverse problems.
In iterative regularization, one picks an initial guess ${f}_{0}$  for the unknown $\stackrel{~}{f}$  , and then one iteratively constructs updated approximations via a regularization scheme. The regularization parameter associated with iterative regularization is thus the $“$  stopping point”of the iterative sequence, and an important part of the mathematical theory is the development of stopping criteria for terminating the iteration. In other words, the iteration index plays the role of the regularization parameter $\alpha$  , and the stopping criteria plays of the parameter selection method.

5.1 Descent Methods for Linear Inverse Problems

As an example of iterative regularization, we consider descent methods. Descendent methods have become quite popular in the last years for the solution of linear inverse problems and for nonlinear inverse problems [11. In this subsection we consider two examples.
As an approximation of ${\stackrel{~}{f}}_{m}$  we will choose ${f}_{m,\alpha }$  such that
 $\begin{array}{c}{f}_{m,\alpha }=\left[I-{A}_{m}^{*}{A}_{m}{Q}_{\alpha }\left({A}_{m}^{*}{A}_{m}\right)\right]{f}_{0}+{Q}_{\alpha }\left({A}_{m}^{*}{A}_{m}\right){A}_{m}^{*}{\Pi }_{{Y}_{m}}^{n}y\end{array}$ (5.1)
where ${f}_{0}\in {F}_{m}$  is an initial approach and this ${f}_{0}\in \mathcal{N}\left({A}_{m}{\right)}^{\perp }$  [11.
Most iterative methods for approximating $\stackrel{~}{f}$  are based on a transformation of the normal equation into equivalent fixed point equations like $f=f+{A}_{m}^{*}\left({A}_{m}f-y\right)$  If $\parallel {A}_{m}{\parallel }^{2}<2$  then the corresponding fixed point operator $I-{A}_{m}^{*}{A}_{m}$  is nonexpansive and one may apply the method of successive approximations. It must be emphasized that $I-{A}_{m}^{*}{A}_{m}$  is no contradiction if our inverse problem is ill-posed, since the spectrum of ${A}_{m}^{*}{A}_{m}$  clusters at the origin.

5.2 Landweber iteration

In this subsection we presented the well-known Landweber iteration, which arises from converting the necessary conditions for minimizing  2.1 into a fixed point iteration. Much development in the last few years has taken place in advancing the theory of Landweber iteration for linear and nonlinear inverse problems.
Using the terminology of the last sections, we introduce the function
 $\begin{array}{c}{Q}_{k}\left(\lambda \right)={\sum }_{j=0}^{k-1}\left(1-\lambda {\right)}^{j}={\lambda }^{-1}\left(1-\left(1-\lambda {\right)}^{k}\right)\end{array}$ (5.2)
We call ${Q}_{k}$  the iteration polynomial of degree $k-1$  . Associated with it is the polynomial ${r}_{k}\left(\lambda \right)=1-\lambda {Q}_{k}\left(\lambda \right)=\left(1-\lambda {\right)}^{k}$  of degree $k$  , which is called the residual polynomial since it determines the residual $y-{A}_{m}{f}_{m,k}$  .
Thus, inserting the equation  5.2 in  5.1 we obtain recursively,
 $\begin{array}{c}{f}_{m,k+1}={f}_{m,k}-{A}_{m}^{*}\left({A}_{m}{f}_{m,k}-{y}_{m}\right),k=0,1,\dots \end{array}$ (5.3)
starting from an initial guess ${f}_{0}$  . This is a steepest descent method called the linear version of Landweber's iteration. Each step of the iterative process  5.3 is carried out along the direction opposite to the direction of the gradient of the quadratic functional $J\left(f\right)$  in  2.1 . It is known that there is the greatest decrease of the functional along this direction.
If $\parallel {A}_{m}\parallel \le 1$  , we considerer $\lambda \in \left(0,1\right]$  such that in this interval $\lambda {Q}_{k}\left(\lambda \right)$  is uniformly bounded and since ${Q}_{k}\left(\lambda \right)$  converge to $1/\lambda$  as $k\to \infty$  then according to Theorem  3.2 the sequences ${f}_{m,k}$  converge to ${\stackrel{~}{f}}_{m}$  when $y\in \mathcal{D}\left({A}_{m}^{†}\right)$  . If $\parallel {A}_{m}\parallel$  is not bounded by one, then we introduce a relaxation parameter $0<\tau <\parallel {A}_{m}{\parallel }^{-2}$  in front of ${A}_{m}^{*}$  in  5.3 , i.e, we would iterate
 $\begin{array}{c}{f}_{m,k+1}={f}_{m,k}-\tau {A}_{m}^{*}\left({A}_{m}{f}_{m,k}-y\right),k=0,1,\dots \end{array}$ (5.4)
If $\tau \equiv {\tau }_{k}$  , one can obtain various variants of the method of steepest descent depending on a choice of the sequence ${\tau }_{k}$  . The Landweber iteration  5.4 is usually called a method of simple iteration.
In the following we derive a simple estimate for the error propagation in the Landweber iteration. We then have the following result,
Corollary 5.1. Let $\tau =1/\left(2\parallel {A}_{m}{\parallel }^{2}\right)<1/{\lambda }_{1}$  . If $y\in \mathcal{ℛ}\left({A}_{m}\right)$  , then the Landweber iteration is an order optimal regularization method, i.e, $\parallel {\stackrel{~}{f}}_{m}-{f}_{k\left(m\right)}{\parallel }^{2}\le 2{c}_{1}{k}^{-2\mu }+2{c}_{2}\frac{{\sigma }^{2}}{n}\left(\tau k{\right)}^{\left(2p+1\right)/2p},$  where ${c}_{1}={\rho }^{2}\left(\frac{\mu }{\tau e}{\right)}^{\mu }$  and ${c}_{2}=\frac{1}{2p+1}\left(\frac{2p+1}{2p-1}{\right)}^{\left(2p+1\right)/4p}.$
• Proof. To apply Theorem  3.4 we have to study the terms of the bias, $IE\parallel {\stackrel{~}{f}}_{m}-{R}_{\alpha }{A}_{m}\stackrel{~}{f}{\parallel }^{2}$  , and variance $IE\parallel {R}_{\alpha }{A}_{m}\stackrel{~}{f}-{f}_{k\left(m\right)}\parallel$  . By  5.2 we have ${\stackrel{~}{f}}_{m}-{R}_{\alpha }{A}_{m}\stackrel{~}{f}=\left(I-{A}_{m}^{*}{A}_{m}{Q}_{k}\left({A}_{m}^{*}{A}_{m}\right)\right){\stackrel{~}{f}}_{m}=\left(I-{A}_{m}^{*}{A}_{m}{\right)}^{k}{\stackrel{~}{f}}_{m}$  We have to study the residual polynomial ${r}_{k}\left(\lambda \right)=\left(1-\lambda {\right)}^{k}$  of the Landweber iteration.
For $0\le \lambda \le \parallel {A}_{m}{\parallel }^{2}$  the function ${\lambda }^{\mu }|1-\lambda {Q}_{k}\left(\lambda \right)|$  assumes its maximum for $\lambda ={\tau }^{-1}\mu \left(\mu +k{\right)}^{-1}$  .
Thus, we have
 $\begin{array}{cc}{\lambda }^{\mu }|1-\lambda {Q}_{k}\left(\lambda \right)|& \le {max}_{0\le \lambda \le \parallel {A}_{m}{\parallel }^{2}}{\lambda }^{\mu }|1-\lambda {Q}_{k}\left(\lambda \right)|\end{array}$
 $\begin{array}{cc}& <\frac{{\mu }^{\mu }}{{\tau }^{\mu }\left(\mu +k{\right)}^{\mu }}\frac{{k}^{k}}{\left(\mu +k{\right)}^{k}}\end{array}$
 $\begin{array}{cc}& <{\left(\frac{\mu }{\tau e}\right)}^{\mu }{k}^{-\mu }\end{array}$
 $\begin{array}{}\end{array}$
This leads to numbers ${\omega }_{\mu }\left(k\right)$  as introduced in Theorem  3.4  ${\omega }_{\mu }\left(k\right)={\left(\frac{\mu }{\tau e}\right)}^{\mu }{k}^{-\mu }$  Thus, the term corresponding to the bias is bounded by $\parallel {\stackrel{~}{f}}_{m}-{R}_{\alpha }{A}_{m}\stackrel{~}{f}{\parallel }^{2}\le {\rho }^{2}\left(\frac{\mu }{\tau e}{\right)}^{2\mu }{k}^{-2\mu }.$  Next, we establish bounds for the variance term. By assumption, the singular values satisfy ${\lambda }_{j}\approx {j}^{-2p}.$  Note that for small values of ${\lambda }_{j}$  we have
${Q}_{k}^{2}\left(\lambda \right)\le \left(\tau k{\right)}^{2}$  $\forall$  $j>{m}^{\prime }$
and for big values of ${\lambda }_{j}$  ( ${\lambda }_{j}\approx {\lambda }_{1}$  )
${Q}_{k}^{2}\left({\lambda }_{j}\right)\le {\lambda }_{j}^{-2}$  $\forall$  $j<{m}^{\prime }$  .
Consequently,
 $\begin{array}{cc}nTr\left({Q}_{k}^{2}\left({A}_{m}{A}_{m}^{*}\right){A}_{m}{A}_{m}^{*}\right)={\sum }_{j=1}^{m}{Q}_{k}^{2}\left({\lambda }_{j}\right){\lambda }_{j}& \le {\sum }_{j=1}^{{m}^{\prime }}{\lambda }_{j}^{-1}+{\sum }_{j>{m}^{\prime }}\left(\tau k{\right)}^{2}{\lambda }_{j}\end{array}$
 $\begin{array}{cc}& \le {\int }_{0}^{{m}^{\prime }}{s}^{2p}ds+\left(\tau k{\right)}^{2}{\int }_{0}^{{m}^{\prime }}{s}^{-2p}ds\end{array}$
 $\begin{array}{}\end{array}$
This suggest searching ${m}^{\prime }\approx c\left(\tau k{\right)}^{1/2p}$  for $p>1/2$  , where $c=\left(\frac{2p+1}{2p+1}{\right)}^{1/4p}$  . Hence we have,
 $\begin{array}{cc}IE\parallel {R}_{\alpha }{A}_{m}\stackrel{~}{f}-{f}_{k\left(m\right)}{\parallel }^{2}& =Tr\left({Q}_{k}^{2}\left({A}_{m}{A}_{m}^{*}\right){A}_{m}{A}_{m}^{*}\right)\end{array}$
 $\begin{array}{cc}& \le \frac{{c}^{2p+1}}{2p+1}\frac{\left(\tau k{\right)}^{\left(2p+1\right)/2p}}{n}.\end{array}$
 $\begin{array}{}\end{array}$
Finally, this implies $IE\parallel {\stackrel{~}{f}}_{m}-{f}_{k\left(m\right)}{\parallel }^{2}\le 2{c}_{1}{k}^{-2\mu }+2{c}_{2}\frac{{\sigma }^{2}}{n}\left(\tau k{\right)}^{\left(2p+1\right)/2p},$  where ${c}_{1}={\rho }^{2}\left(\frac{\mu }{\tau e}{\right)}^{\mu }$  and ${c}_{2}=\frac{1}{2p+1}\left(\frac{2p+1}{2p-1}{\right)}^{\left(2p+1\right)/4p}$  .
Remark 5.2. Note that under the above inequality is satisfied if the dimension of the set is such that ${d}_{{m}_{0}}\approx {n}^{\frac{1}{4\mu p+2p+1}}$  . Here, the optimal choice of regularization sequence, depending on $p$  and $\mu$  . The optimal rates are of order $IE\parallel \stackrel{~}{f}-{f}_{k\left(m\right)}{\parallel }^{2}=\mathcal{O}\left({n}^{-\frac{4\mu p}{4\mu p+2p+1}}\right)$  . Analogous results are obtained in the ill-posed problem literature, see for example [5, where typically in a Hilbert scale setting optimal rates are of order $\mathcal{O}\left({n}^{-\frac{2s}{2s+2p+1}}\right)$  , with $s=2\mu p$  .
We are ready to state our main result for the Landweber iteration, which bounds the mean squared error of the select estimate ${\stackrel{^}{f}}_{\stackrel{^}{k}}$  basically by the smallest mean squared error among the estimates ${f}_{k}$  plus a remainder term of order $1/n$  . The result follows from Theorem  4.1 .
Corollary 5.3. Let $\tau =1/\left(2\parallel {A}_{m}{\parallel }^{2}\right)<1/{\lambda }_{1}$  . Next assume $\stackrel{^}{k}$  as in  4.2 and ${d}_{{m}_{0}}$  as in  4.1 . If $y\in \mathcal{ℛ}\left({A}_{m}\right)$  then for any $f\in {F}_{m}$  and any $k$  , the following inequality holds true
 $\begin{array}{cc}IE\parallel {\stackrel{~}{f}}_{m}-{\stackrel{^}{f}}_{\stackrel{^}{k}}{\parallel }^{2}& \le \frac{1}{\left(1-\nu \right)}{inf}_{k\in \mathcal{K}}\left[C\left(1+\nu \right)\parallel {\stackrel{~}{f}}_{m}-{f}_{k}{\parallel }^{2}+\frac{2r{\sigma }^{2}\left(1+{L}_{k}\right)\left(c\left(\tau k{\right)}^{\frac{2p+1}{2p}}+\tau k\right)}{n}\right]\end{array}$
 $\begin{array}{cc}& +\frac{4{\sigma }^{2}}{n}{\sum }_{k}\frac{\tau k}{d}\left[\sqrt{dr{L}_{k}\left[c\left(\tau k{\right)}^{1/2p}+1\right]}+1\right]{e}^{-\sqrt{dr{L}_{k}\left[c\left(\tau k{\right)}^{1/2p}+1\right]}},\end{array}$
 $\begin{array}{}\end{array}$
for some $C>0$  and $c=\frac{1}{2p+1}\left(\frac{2p+1}{2p-1}{\right)}^{\left(2p+1\right)/4p}$  .
• Proof. For fixed ${\lambda }_{j}$  and $m$  we have that the terms of the trace and spectral radius are bounded by the follows expression  $\begin{array}{c}Tr\left({R}_{k}^{t}{R}_{k}\right)={\sum }_{j=1}^{m}{Q}_{k}^{2}\left({\lambda }_{j}\right){\lambda }_{j}\le {\sum }_{j=1}^{{m}^{\prime }}{\lambda }_{j}^{-1}+{\sum }_{j>{m}^{\prime }}\left(\tau k{\right)}^{2}{\lambda }_{j}\end{array}$ (5.5)
and  $\begin{array}{c}{\rho }^{2}\left({R}_{k}^{t}{R}_{k}\right)={max}_{j\in m}{Q}_{k}^{2}\left({\lambda }_{j}\right){\lambda }_{j}\le {max}_{j\le {m}^{\prime }}{\lambda }_{j}^{-1}+{max}_{j>{m}^{\prime }}\left(\tau k{\right)}^{2}{\lambda }_{j}.\end{array}$ (5.6)
Balancing both terms in  5.5 and  5.6 gives the optimal choice of the trace and the spectral radius, respectively. Thus, we have $Tr\left({R}_{k}^{t}{R}_{k}\right)\approx \frac{1}{2p+1}{\left(\frac{2p+1}{2p-1}\right)}^{\frac{\left(2p+1\right)}{4p}}\frac{\left(\tau k{\right)}^{\frac{2p+1}{2p}}}{n}$  and ${\rho }^{2}\left({R}_{k}^{t}{R}_{k}\right)\approx \frac{\tau k}{n}$  Note that the penalization term is roughly proportional to $\frac{1}{n}\left[\frac{1}{2p+1}{\left(\frac{2p+1}{2p-1}\right)}^{\frac{\left(2p+1\right)}{4p}}\left(\tau k{\right)}^{\frac{2p+1}{2p}}+\tau k\right]$  On the other hand $\frac{Tr\left({R}_{\alpha }^{t}{R}_{\alpha }\right)}{{\rho }^{2}\left({R}_{\alpha }\right)}=\frac{1}{2p+1}{\left(\frac{2p+1}{2p-1}\right)}^{\left(2p+1\right)/4p}\left(\tau k{\right)}^{1/2p}.$  The result then follows directly from Theorem  4.1 .

5.3 Nonlinear multistep iterative process

Many approximate methods widely used in practice are nonlinear. We cite a important example of nonlinear approximate method. We considerer a nonlinear multistep iterative process, which have error residual $1-\lambda {Q}_{k}\left(\lambda \right){=}^{k}{\prod }_{i=1}\left(1-{\tau }_{ik}^{-1}\lambda \right)$  with ${\tau }_{ik}={\tau }_{ik}\left({f}_{0},A,y\right)>0,0<{\tau }_{1k}\le {\tau }_{2k}\dots \le {\tau }_{kk}\le {\lambda }_{1}$  . Then for $\lambda >0$  , ${Q}_{k}\left(\lambda \right)$  have the following representation ${Q}_{k}\left(\lambda \right)={\lambda }^{-1}\left[1{-}^{k}{\prod }_{i=1}\left(1-{\tau }_{ik}^{-1}\lambda \right)\right]$  The following corollary is established.
Corollary 5.4. Let ${\tau }_{ik}={\tau }_{ik}\left({f}_{0},A,y\right)>0,$  with $0<{\tau }_{1k}\le {\tau }_{2k}\dots \le {\tau }_{kk}\le {\lambda }_{1}$  . If $y\in \mathcal{ℛ}\left({A}_{m}\right)$  , then the nonlinear multistep iterative process is an order optimal regularization method, i.e, $IE\parallel {\stackrel{~}{f}}_{m}-{f}_{k\left(m\right)}{\parallel }^{2}\le 2{c}_{1}\left({\sum }_{i=1}^{k}{\tau }_{{i}_{k}}^{-1}{\right)}^{-2\mu }+2{c}_{2}\frac{{\sigma }^{2}}{n}\left({\sum }_{i=1}^{k}{\tau }_{{i}_{k}}^{-1}{\right)}^{\left(2p+1\right)/2p},$  where ${c}_{1}={\rho }^{2}{\mu }^{\mu }\left(\mu +1{\right)}^{-1}$  and ${c}_{2}=\frac{1}{2p+1}\left(\frac{2p+1}{2p-1}{\right)}^{\left(2p+1\right)/4p}.$
• Proof. As before, we investigate the behavior of the bias and the variance.
In the relation ${\lambda }^{\mu }\left(1-\lambda {Q}_{k}\left(\lambda \right)\right)={{\lambda }^{\mu }}^{k}{\prod }_{i=1}\left(1-{\tau }_{ik}^{-1}\lambda \right)$  the least upper can not be reached at the points $\lambda =0$  and $\lambda ={\tau }_{1k}$  , since the estimated function is not equal to zero identically.
On the other hand  $\begin{array}{c}\left[1-\lambda {Q}_{k}\left(\lambda \right){\right]}^{\prime }\left[\lambda {Q}_{k}\left(\lambda \right)-1{\right]}^{-1}={\sum }_{i=1}^{k}\frac{1}{{\tau }_{ik}-\lambda }\end{array}$ (5.7)
Since the function in the right-hand of  5.7 does not decrease as a function of $\lambda$  on the half-interval $\left[0,{\tau }_{1k}\right)$  then, the estimates of the Theorem  3.5 are valid.
Thus, for $0\le \lambda \le {\tau }_{1k}$  , we have
 $\begin{array}{cc}{sup}_{0\le \lambda \le {\tau }_{1k}}|{Q}_{k}\left(\lambda \right)|& ={Q}_{k}\left(0\right)={\sum }_{i=1}^{k}{\tau }_{{i}_{k}}^{-1}\end{array}$
 $\begin{array}{}\end{array}$
and
 $\begin{array}{cc}{sup}_{0\le \lambda \le {\tau }_{1k}}{\lambda }^{\mu }|1-\lambda {Q}_{k}\left(\lambda \right)|& <{\mu }^{\mu }\left(\mu +1{\right)}^{-1}\left({\sum }_{i=1}^{k}{\tau }_{{i}_{k}}^{-1}{\right)}^{-\mu }\end{array}$
 $\begin{array}{}\end{array}$
Note that ${\omega }_{\mu }\left(k\right)=\left({\sum }_{i=1}^{k}{\tau }_{{i}_{k}}^{-1}{\right)}^{-\mu }.$  Thus, the bias is bounded by $\parallel {\stackrel{~}{f}}_{m}-{R}_{\alpha }{A}_{m}\stackrel{~}{f}{\parallel }^{2}\le {c}_{1}\left({\sum }_{i=1}^{k}{\tau }_{{i}_{k}}^{-1}{\right)}^{-2\mu },$  where ${c}_{1}={\rho }^{2}{\mu }^{2\mu }\left(\mu +1{\right)}^{-2}$  .
On the order hand, it is not difficult to see that
 $\begin{array}{cc}Tr\left({Q}_{k}^{2}\left({A}_{m}{A}_{m}^{*}\right){A}_{m}{A}_{m}^{*}\right)& \le \frac{1}{n}\left[{\sum }_{1\le j\le {m}^{\prime }}{j}^{2p}+\left({\sum }_{i=1}^{k}{\tau }_{ik}^{-1}{\right)}^{2}{\sum }_{j>{m}^{\prime }}{j}^{-2p}\right]\le \frac{1}{n}\left[\frac{{{m}^{\prime }}^{2p+1}}{2p+1}+\left({\sum }_{i=1}^{k}{\tau }_{ik}^{-1}{\right)}^{2}\frac{{{m}^{\prime }}^{-2p+1}}{2p-1}\right]\end{array}$
 $\begin{array}{}\end{array}$
This suggest searching ${m}^{\prime }\approx c\left({\sum }_{i=1}^{k}{\tau }_{{i}_{k}}^{-1}{\right)}^{1/2p}$  with $c=\left(\frac{2p+1}{2p+1}{\right)}^{1/4p}$  .
Thus, we have that the term variance is bounded by $IE\parallel {R}_{\alpha }{A}_{m}\stackrel{~}{f}-{f}_{k\left(m\right)}\parallel ={\sigma }^{2}Tr\left({Q}_{k}^{2}\left({A}_{m}^{*}{A}_{m}\right){A}_{m}^{*}{A}_{m}\right)\le {c}_{2}\frac{{\sigma }^{2}}{n}\left({\sum }_{i=1}^{k}{\tau }_{{i}_{k}}^{-1}{\right)}^{\left(2p+1\right)/2p}$  where ${c}_{2}=\frac{1}{2p+1}\left(\frac{2p+1}{2p-1}{\right)}^{\left(2p+1\right)/4p}$  Finally we have $IE\parallel {\stackrel{~}{f}}_{m}-{f}_{k\left(m\right)}{\parallel }^{2}\le 2{c}_{1}\left({\sum }_{i=1}^{k}{\tau }_{{i}_{k}}^{-1}{\right)}^{-2\mu }+2{c}_{2}\frac{{\sigma }^{2}}{n}\left({\sum }_{i=1}^{k}{\tau }_{{i}_{k}}^{-1}{\right)}^{\left(2p+1\right)/2p}$  Balancing the bias and variance terms gives the optimal choice $IE\parallel {\stackrel{~}{f}}_{m}-{f}_{k\left(m\right)}{\parallel }^{2}=\mathcal{O}\left({n}^{-\frac{4\mu p}{4\mu p+2p+1}}\right).$
We have the following result.
Corollary 5.5. Let ${\tau }_{ik}$  be as in corollary  5.4 . Next assume $\stackrel{^}{k}$  as in  4.2 and ${d}_{{m}_{0}}$  as in  4.1 . If $y\in \mathcal{ℛ}\left({A}_{m}\right)$  then for any $f\in {F}_{m}$  and any $k$  , the following inequality holds true
 $\begin{array}{cc}IE\parallel {\stackrel{~}{f}}_{m}-{\stackrel{^}{f}}_{\stackrel{^}{k}}{\parallel }^{2}& \le \frac{1}{\left(1-\nu \right)}{inf}_{k\in \mathcal{K}}\left[C\left(1+\nu \right)\parallel {\stackrel{~}{f}}_{m}-{f}_{k}{\parallel }^{2}+\frac{2r{\sigma }^{2}\left(1+{L}_{k}\right)\left(c\left({\sum }_{i=1}^{k}{\tau }_{{i}_{k}}^{-1}{\right)}^{\frac{2p+1}{2p}}+{\sum }_{i=1}^{k}{\tau }_{{i}_{k}}^{-1}\right)}{n}\right]\end{array}$
 $\begin{array}{cc}& +\frac{4{\sigma }^{2}}{n}{\sum }_{k}\frac{{\sum }_{i=1}^{k}{\tau }_{{i}_{k}}^{-1}}{d}\left[\sqrt{dr{L}_{k}\left[c\left({\sum }_{i=1}^{k}{\tau }_{{i}_{k}}^{-1}{\right)}^{1/2p}+1\right]}+1\right]{e}^{-\sqrt{dr{L}_{k}\left[c\left({\sum }_{i=1}^{k}{\tau }_{{i}_{k}}^{-1}{\right)}^{1/2p}+1\right]}},\end{array}$
 $\begin{array}{}\end{array}$
for some $C>0$  and $c=\frac{1}{2p+1}\left(\frac{2p+1}{2p-1}{\right)}^{\left(2p+1\right)/4p}$  .
• Proof. First observe that $Tr\left({R}_{k}^{t}{R}_{k}\right)\approx \frac{1}{2p+1}{\left(\frac{2p+1}{2p-1}\right)}^{\frac{\left(2p+1\right)}{4p}}\frac{\left({\sum }_{i=1}^{k}{\tau }_{{i}_{k}}^{-1}{\right)}^{\frac{2p+1}{2p}}}{n}$  and ${\rho }^{2}\left({R}_{k}^{t}{R}_{k}\right)\approx \frac{{\sum }_{i=1}^{k}{\tau }_{{i}_{k}}^{-1}}{n}$  Consequently $\frac{Tr\left({R}_{k}^{t}{R}_{k}\right)}{\rho \left({R}_{k}\right)}\approx \frac{1}{2p+1}{\left(\frac{2p+1}{2p-1}\right)}^{\frac{\left(2p+1\right)}{4p}}\left({\sum }_{i=1}^{k}{\tau }_{{i}_{k}}^{-1}{\right)}^{\frac{1}{2p}}$  Note that both $n{\rho }^{2}\left({R}_{k}\right)$  and $Tr\left({R}_{k}^{t}{R}_{k}\right)/{\rho }^{2}\left({R}_{k}\right)$  do not depend on n. The proof then follows directly from theorem  4.1 .
References

1. Barron A. et al., (1999). “Risk Bounds for Model Selection via Penalization ”. Probab. Theory and Related Fields. 113, pp. 467-493.
2. Birg L., and Massart P., (1978). “Minimum Contrast Estimators on Sieves: Exponential Bounds and Rates of Convergence ”. Bernoulli, Vol4.N3, pp 329-395.
3. Bousquet O., (2002). “ Concentration Inequalities for Sub-Additive Functions Using of Entropy Method.”
4. Burger, M., (2001). “A level set method for inverse problems ”. Inverse Problems 17, pp. 13271356.
5. Cavalier L., Golubev G., Picard D., and Tsybakov A., (2002). “Oracle inequalities for inverse problem ”. Ann. Statist, Vol30. N3, pp 843-874.
6. Deuflhard P., Engl H. and Scherzer O., (1998) “A convergence analysis of iterative methods for the solution of nonlinear ill-posed problems under affinely invariant conditions ”. Inverse Problems Vol14, pp. 1081-1106.
7. Dey A.K. et al, (1996) “Cross-Validation for Parameter Selection in inverse estimation problems ”. Scand. J. Statist., Vol23. N4, pp. 609-620.
8. Engl H., (1993) “Regularization methods for the stable solution of inverse problems ”. Surveys on Mathematics for Industry 3, pp. 71-143.
9. Engl H., Hanke M. and Neubauer A., (1996) “Reguralization of Inverse Problems ”. Kluwer Academic Publishers.
10. Engl H. and Grever W., (1994) “Using the L-curve for Determinig Optimal Regularization parameters ”. Numer. MAth. Vol. 69, pp. 25-31.
11. Gilyazov S.F., and Gol'dman N.L., (2000). “Reguralization of Ill-Posed Problems by Iteration Methods”. Kluwer Academic Publishers.
12. Cohen A., Hoffmann M., and Reiss M. “Adaptive Wavelet Galerkin Methods for Linear Inverse Problems”. SIAM J. Numer. Anal. Vol. 42, No. 4, pp. 1479-1501.
13. Kilmer M.E, and O'leary D.P., (2001). “Choosing Reguralization Parameters in iterative Methods for Ill-Posed Problems”. SIAM J. MATRIX ANAL. APPL. Vol. 22, N4, pp. 1204-1221.
14. Lamm P.K., (1999).“Some Recent Developments and Open Problems in Solution Methods for Mathematical Inverse Problems”. Department of Mathematics, Michihan state University, USA.
15. Ledoux, M., and Talagrand, (1996). “Deviation Inequalities for Product Measures ”. ESAIM: Probabilities and Statistics 1, pp. 63-87.
16. Loubes J.-M. and Lude $\stackrel{~}{n}$  a C., (2004). “Penalized Estimators for Nonlinear Inverse Problems. ”.
17. Loubes J.-M., and Van De Geer S., (2002). “Adaptive estimation in regression, using soft thresholding type penalties ”. Statistica Neerlandica, 56, pp 453-478.
18. Lude $\stackrel{~}{n}$  a C., and Rios,(2003).“Penalized Model Selection for Ill-posed Linear Problems. ”.
19. Mathé P., and Pereverezev S.V., (2003) “Discretization Strategy for Linear Ill-posed problems in variable Hilbert Scales ”. Inverse Problems, Vol.19, N6, pp. 1263-1279.
20. Morozov V.A., (1966). “On the Solution of Functional Equations by the Method of Regularization ”. Soviet Math. Dokl.,7, pp. 414-417.
21. F. O'sullivan., (1986).“A Statistical Perspective on Ill-Posed Inverse Problems ”. Statistical Science. Vol. 1, N4, pp. 502-527.
22. Pereverzev S. and Schock E., (2003). “On the adaptive selection of the parameter in regularization of ill-posed problems ”.
23. Tikhonov, A.,and Arsenin, V., (1977).“Solutions of Ill-Posed Problems ”. Wiley, New York.

Escuela de Matematicas, Facultad de Ciencias, UCV, Av. Los Ilustres, Los Chaguaramos, Codigo Postal 1020-A, Caracas Venezuela. Telf.: (58)212-6051481.
Departamento de Matematicas, IVIC, Carretera Panamericana KM. 11, Aptdo. 21827, Codigo Postal 1020-A, Caracas Venezuela. Telf.: (58)212-5041412.
E-mail address : afermin@euler.ciens.ucv.ve E-mail address : cludena@ivic.ve