Statistical Inverse Problems in Hilbert Scales

In this paper, we study the nonlinear operator equation

with the infinite dimensional separable real Hilbert spaces $H$ and $H^{\prime }$ with the inner products $⟨\cdot ,\cdot {⟩}_H$ and $⟨\cdot ,\cdot {⟩}_{H^{\prime }}$, respectively. Here, $H^{\prime }$ is the space of functions between a Polish space $X$ and a real separable Hilbert space $Y$. We observe the noisy values of the function $\hat{g}$ at the inputs $x_i$:

for $1\le i\le m$. Here, $m$ is the number of observations which is called the sample size. In contrary to the direct learning scheme where we estimate the function $\hat{g}$, here we aim to estimate the function $\hat{f}$ directly from the observations.

A common approach to stably approximate the solution of equation (1) is the Tikhonov regularization scheme. Sometimes, we have the information about the true solution, e.g., the true solution may be differentiable. To incorporate this information, we employ the Tikhonov regularization in Hilbert scales. This scheme consists of a functional which is a linear combination of a fidelity term measuring the fitness of the data and a penalty term in a stronger norm forcing the smoothness in the approximated solution. To define this scheme, we introduce a densely defined, unbounded, closed, linear, self-adjoint, strictly positive operator $L:D\left(L\right)\subset H\to H$ such that for some ${\ell }_L>0$,

Here, we observe that $L^{-1}:H\to H$ is a bounded operator from the strict positivity of the operator $L$.

The Tikhonov functional for the considered nonlinear inverse problem with sample $z={\left\{\right(x_i,y_i\left)\right\}}_{i=1}^m$ is given by

where $\overline{f}\in D\left(A\right)\cap D\left(L\right)$ is an initial guess. The regularization parameter $\lambda >0$ has to balance both terms appropriately. Then, the Tikhonov regularization scheme in Hilbert scales can be defined as

For the continuous and weakly sequentially closed operator $A$, there exists a global solution of the regularization scheme in (4). But it is not necessarily unique, since $A$ is nonlinear (see [14, Section 4.1.1]).

We consider the Hilbert scales $H_a$ generated by the operator $L$. Here the spaces $H_a:=D\left(L^a\right)$ are Hilbert spaces equipped with the inner product ${⟨f,g⟩}_{H_a}={⟨L^{a}f,L^{a}g⟩}_H,f,g\in H_a$. For the Hilbert scales, we have the well-known interpolation inequality

which holds for any $a<b<c$.

The regularization schemes in Hilbert scales have been well-studied and analysed under the different assumptions in the classical inverse problems [2, 5, 7, 14]. In learning theory, the general regularization in Hilbert scales is introduced for the linear inverse problems and established the rates of convergence [13]. Nicole et al. [10] studied the Stochastic Gradient Descent in Hilbert scales and the author provided different examples of Hilbert scales in learning. Further, the authors discussed the error estimates for the Stochastic Gradient Descent scheme for the direct learning problem. In the paper [11], the rates of convergence are established for nonlinear statistical inverse learning problems in the RKHS setting. The authors considered some assumptions on the nonlinearity of the operator $A$ such as Fréchet differentiability of the operator, Lipschitz continuity of the Fréchet derivative, and the link condition to transfer the smoothness in terms of the operator $L$ to the covariance operator.

Here, we consider the nonlinear inverse learning problem in Hilbert scales satisfying conditional stability estimates characterized by general concave index functions. We use the Tikhonov regularization schemes to obtain the stable approximate solution in the RKHS framework. Werner and Hofmann [16] illustrated the validity of the conditional stability estimates in different models and real-world situations. The authors showed that the derivative of $A$ is not always necessary for this condition.

For the regularization schemes in the RKHSs, generally, we consider the smoothness by the source condition in terms of the covariance operator which implies the rates of convergence. The covariance operator depends on the considered kernel and unknown probability measure. Therefore, the source condition cannot be verified practically. Moreover, the misspecified kernel affects the source condition and consequently, the rates of convergence for the regularization schemes. Here, we consider the smoothness of the true solution in terms of the known operator $L$ which can be checked in practice. We divide the smoothness into two cases: the regular case (i.e., $\hat{f}\in D\left(L\right)$) and the oversmoothing case (i.e., $\hat{f}\notin D\left(L\right)$). The oversmoothing case is very delicate. We consider that the regularized solution in Hilbert scales (4) belongs to $D\left(L\right)$. But the true solution does not belong to $D\left(L\right)$ in the oversmoothing case. The analysis is also tricky for nonlinear inverse problems since the Tikhonov regularization in Hilbert Scales does not have an explicit solution. The analysis starts with the step $E_{z,\lambda }\left(f_{z,\lambda }\right)\le E_{z,\lambda }\left(\hat{f}\right)$. But $E_{z,\lambda }\left(\hat{f}\right)$ is not well-defined in the oversmoothing case (since $\hat{f}\notin D\left(L\right)$). We will utilize the concept of distance functions to overcome this problem.

The main results of our paper can be summarized as follows:

• We discuss the rates of convergence for the Tikhonov regularization in Hilbert Scales under a conditional stability assumption for the inverse problem.

• We obtain the error estimates in the absence of the widely-considered source condition. We will use the concept of the distance functions for this.

• We establish the error bounds in both the regular case and the oversmoothing case for the appropriate benchmark smoothness.

The manuscript is organized as follows: In Section 2, we present the basic definitions, notation, and assumptions required in our analysis. In Section 3, we state and prove our main results. Here, we discuss the rates of convergence for Tikhonov regularization in Hilbert scales in the probabilistic sense. In Section 4, we present the explicit rates in terms of sample size by bounding the distance functions. In Appendix, we state the probabilistic estimate of perturbation inequalities.

Let the input space $X$ be a Polish space and the output space $(Y,{⟨\cdot ,\cdot ⟩}_Y)$ be a real separable Hilbert space. We consider the joint probability measure $\rho$ on the sample space $Z=X\times Y$. We denote the marginal distribution on $X$ by $\nu$ and the conditional distribution of $y$ given $x$ by $\rho \left(y|x\right)$. Therefore, the measure $\rho$ can be split as $\rho (x,y)=\rho \left(y|x\right)\nu \left(x\right)$.

For the probability measure $\rho$ on $X\times Y$, we assume that

For the considered model $y=\hat{g}\left(x\right)+\epsilon$ with centred noise $\epsilon$ we find ${\int }_{Y}yd\rho \left(y|x\right)=\hat{g}\left(x\right)$ provided that the conditional expectation w.r.t. $\rho$ of $y$ given $x$ exists (a.s.). This holds true under Assumption (6). This fact together with the operator equation (1) motivates us to consider the following assumption.

Here, $\hat{f}$ is the true solution of equation (1) which we aim at estimating. Here, we want to mention that the function $\hat{f}$ is also the minimizer of the expected risk considered in [11].

We consider a Bernstein-type assumption for the noise $\epsilon =y-A\left(\hat{f}\right)\left(x\right)$:

We want to utilize the properties of the Reproducing kernel Hilbert spaces (RKHSs) in our analysis. Therefore, we assume that $\mathrm{Ran}\left(A\right)$ is contained in a vector-valued Reproducing kernel Hilbert space (RKHSvv). The RKHSvv $H_K$ arises from the operator-valued positive semi-definite kernel $K:X\times X\to L\left(Y\right)$ [9]. Here, $L\left(Y\right)$ is the Banach space of bounded linear operators.

This assumption implies that $H^{\prime }\subset L^2(X,\nu ;Y)$. We denote the canonical injection map $H^{\prime }$ to $L^2(X,\nu ;Y)$ by $I_{\nu }$ and the corresponding covariance operator is $T_{\nu }:=I_{\nu }^{\ast }{I}_{\nu }$. From the above assumption, we see that the covariance operator is positive and trace class. The covariance operator is very important in our convergence analysis. We will need some regularity assumptions in terms of the covariance operator on the marginal probability measure $\nu$ to achieve the uniform convergence rates for the regularized solution (4).

The error estimates studied in our analysis are based on the smoothness of the true solution and the behaviour of the effective dimension. The error estimates and the optimal parameter choice depend on the effective dimension for the regularization methods in reproducing kernel Hilbert spaces [4, 3, 12]. To achieve the fast convergence rates, we introduce the concept of the effective dimension $N\left(\lambda \right)$ [17]:

The effective dimension is a continuous, decreasing function of $\lambda$. The effective dimension is finite, since the operator $T_{\nu }$ is a trace class, and we get

The different behaviours of the eigenvalues of the covariance operator lead to different decay rates of the effective dimension [8]. Under the different scenarios of the effective dimension, we will get the explicit convergence rates in the next section.

In order to establish the error estimate, we introduce the discrete operators for the samples. For the ordered set $(x{)}_i=x_i$, we define the Sampling Operator

We define the inner product space $Y^m$ with the inner product  for $(y{)}_i=y_i$ and $(y^{\prime }{)}_i=y_i^{\prime }$ for $1\le i\le m$. Then, we get the expression of its adjoint $S_x^{\ast }$ as

It can be easily checked that under Assumption 3, ${\parallel S_x\parallel }_{H^{\prime }\to Y^m}\le \kappa$.

We need to make some assumptions about the nonlinear structure of operator $A$. Following the work of Werner and Hofmann [16], we consider the following assumption on $A$, $D\left(A\right)$, and $\hat{f}$. To introduce this assumption, we define the closed balls $B_{\mu }^u\left(\hat{f}\right)=\{f\in H_u:{\Vert f-\hat{f}\Vert }_u\le \mu \}$ in $H_u(u\in R)$ with center $\hat{f}\in H_u$ and radius $\mu$ $(0<\mu \le 1)$ and their intersections with the domain of $A$, $D_{\mu }^u\left(\hat{f}\right):=B_{\mu }^u\left(\hat{f}\right)\cap D\left(A\right)$. For the simplicity, we will denote $B_{\mu }^0\left(\hat{f}\right)$ and $D_{\mu }^0\left(\hat{f}\right)$ and by $B_{\mu }\left(\hat{f}\right)$ and $D_{\mu }\left(\hat{f}\right)$.

Assumption 4 (iv) is called the conditional stability estimate which helps us to characterize the degree of ill-posedness of inverse problems. Here, we note that operator $A$ may not be differentiable, (see the examples in [16]).

The assertions about the convergence of Tikhonov-regularized solution $f_{z,\lambda }$ to the true solution $\hat{f}$ are formulated in this section. First of all, we introduce some standard quantities required to establish the error estimates. We denote

The probabilistic estimates of the above quantities are given in Appendix A. We will use the following standard assumption on the sample size $m$ and the regularization parameter $\lambda$ for our probabilistic estimates:

Now, we introduce the concept of the distance function (also known as ‘approximate source conditions’) which can be used in the absence of the source condition for $\hat{f}$ [1, 15]. It measures the violation of a benchmark smoothness of the true solution. It becomes very important in the ‘oversmoothing case’ $\hat{f}\notin D\left(L\right)$ for regularization in Hilbert Scales.

Here $q$ defines the benchmark smoothness. Let ${\hat{f}}^R$ be the minimizing element of the above problem. Here, we also denote the quantities $d_A\left(R\right)={\Vert I_{\nu }\left[A\right({\hat{f}}^R)-A(\hat{f}\left)\right]\Vert }_{L^2(X,\nu ;Y)}$ and $d^p\left(R\right)={\Vert {\hat{f}}^R-\hat{f}\Vert }_{H_{-p}}$.

Here, we note that when the true solution is of the form $\hat{f}=L^{-q}u$ with ${\parallel u\parallel }_H\le \overline{R}$, then the distance function $D\left(\overline{R}\right)=0$ and the minimizer ${\hat{f}}^{\overline{R}}=\hat{f}$.

The error analysis starts using the fact that $f_{z,\lambda }$ is the minimizer of the Tikhonov functional (4). We get the deterministic expressions (17), (18) for the quantities ${\Vert I_{\nu }\left[A\right(f_{z,\lambda })-A(\hat{f}\left)\right]\Vert }_{L^2(X,\nu ;Y)}$ and ${\Vert L(f_{z,\lambda }-\hat{f})\Vert }_H$ after some rearrangement, and using Cauchy-Schwarz inequality, Young’s inequality. After the simplification and using the probabilistic estimates from Proposition A.1 we get the error estimates in terms of the sample size $m$, the regularization parameter $\lambda$, and distance function by $R\left(\lambda \right)$. The distance function can be measured using the source condition for $\hat{f}$ (see Section 4). Consequently, we get the explicit dependency $\lambda \to R\left(\lambda \right)$. The estimates depend on the effective dimension which will be explicitly expressed in terms of $\lambda$ by using the different decay conditions on the effective dimension. Then, the bounds can be expressed explicitly in terms of $\lambda$ and $m$ for the given smoothness of the solution $\hat{f}$. In Section 4, the a-priori choice of regularization parameter will be obtained by balancing the terms in the error bounds.