MSC2020 subject classifications: 62H30, 41A25, 41A35, 68Q32.
Introduction
Learning theory studies how to infer the structure of an unknown function from sample data. In high-dimensional spaces, learning algorithms often face the curse of dimensionality, as the sample size n n n required to achieve a given accuracy can grow rapidly with the ambient dimension d d d . This motivates the manifold assumption underlying many manifold learning methods, namely that the data lie on or near a low-dimensional manifold embedded in the ambient space R d \mathbb{R}^{d} R d . A central question is then whether learning rates can be governed by the intrinsic dimension ϱ \varrho ϱ rather than the ambient dimension d d d .
Kernel methods, including support vector machines (SVMs), are a classical family of machine learning algorithms with a solid mathematical foundation in reproducing kernel Hilbert space (RKHS) theory [2 , 15 ] . The choice of kernel determines the hypothesis space and its approximation properties. In this paper, we focus on the Gaussian kernel, which measures similarity through Euclidean distances and provides a nonparametric framework for nonlinear prediction. Classical results on Euclidean domains establish fast learning rates for Gaussian SVMs. For hinge-loss classification, the work [17 ] combines Tsybakov’s noise condition with a geometric noise condition to control estimation and approximation errors. Classification with varying Gaussian kernels and general convex losses is studied in [27 ] , which combines Sobolev approximation, RKHS covering estimates, and comparison inequalities for excess misclassification risk. For least-squares and conditional quantile regression, the work [8 ] derives oracle inequalities and essentially minimax optimal learning rates under suitable smoothness and distributional assumptions. The work [14 ] considers correntropy loss and uses varying Gaussian kernels to achieve the minimax rate of convergence.
Despite the dependence on the ambient metric, existing studies show that Gaussian kernel methods can overcome the curse of ambient dimensionality under suitable manifold and regularity assumptions, with learning rates governed by the intrinsic dimension [10 , 13 , 21 , 28 , 29 , 30 ] .
The early work in [29 , 30 ] develops the connection between Gaussian approximation and learning on manifolds and establishes approximation estimates and learning bounds governed by the manifold dimension. However, the approximation results in [29 , 30 ] identify a saturation issue. For the direct Gaussian integral approximation, increasing the regularity α \alpha α of the target beyond α = 2 \alpha= 2 α = 2 does not in general improve the approximation order. On a manifold, the error includes contributions from the target function, the volume density, and the difference between Euclidean and geodesic distances. Estimating these contributions separately without correcting their combined effect leaves a second-order error even for smoother target functions.
Subsequent work in [28 ] studies Gaussian process regression on manifolds and derives intrinsic-dimensional RKHS entropy estimates, but the accompanying approximation result remains restricted to smoothness at most two, leaving the second-order saturation issue unresolved.
The work [10 ] introduces a more general box-counting dimension condition. By expressing the low-dimensional structure through covering numbers of the input domain, their work obtains learning rates for Gaussian SVMs on supports that need not be smooth manifolds. Related ideas are developed for pairwise learning in our earlier work [13 ] , where oracle inequalities and a box-counting dimension assumption yield fast rates for Gaussian kernel estimators. In their squared-loss approximation analyses, both works use signed combinations of Gaussian kernels with different bandwidths to cancel lower-order approximation terms. Under their smoothness assumptions on extensions to the ambient Euclidean space, these constructions achieve arbitrarily high approximation orders and thus overcome second-order saturation of the approximation error. However, the use of ambient-space extensions yields squared RKHS norm bounds of order γ − d \gamma^{-d} γ − d for regression, where γ \gamma γ denotes the Gaussian bandwidth. As we explain in Remark 4.2, this dependence on the ambient dimension would also enter the learning rates obtained by combining these bounds with our oracle inequality for logistic loss.
Our principal inspiration is the recent work [21 ] on adaptive Bayesian regression with low intrinsic dimensionality. The results in [21 ] , Lemma 4.1 and Proposition 4.2 establish high-order Gaussian expansions and remove the smoothness restriction by recursively correcting the input function. They provide approximation of intrinsically defined Hölder functions of arbitrary positive order, with error measured in the supremum norm. By constructing the approximant directly on the manifold without extending the target function to the ambient Euclidean space, they also obtain a squared RKHS norm bound of order γ − ϱ \gamma^{-\varrho} γ − ϱ whose growth as the bandwidth decreases is governed by the intrinsic dimension ϱ \varrho ϱ . Their key technique is a high-order expansion in geodesic normal coordinates that includes the target function, the squared Euclidean distance, and the volume density. Uniform bounds on the remainder and on the Hölder norms of the coefficient functions then justify recursive corrections to the input that cancel the lower-order bias terms.
Building on the high-order Gaussian expansion in [21 ] , we extend the approximation framework from intrinsic Hölder classes to the broader Besov spaces B 2 , ∞ α ( X ) B_{2,\infty}^{\alpha}(X) B 2 , ∞ α ( X ) for every α > 0 \alpha> 0 α > 0 , with the approximation error measured in L 2 ( X ) L_{2}(X) L 2 ( X ) . Our construction shares their principle of canceling lower-order terms in the Gaussian expansion, but implements this cancellation through a finite signed combination of Gaussian integral operators at different bandwidths. The resulting operators act directly on f f f , without evaluating its derivatives or explicitly computing geometric correction coefficients, and are bounded on L 2 ( X ) L_2(X) L 2 ( X ) uniformly in the bandwidth. Combining this boundedness with a Sobolev operator expansion, we obtain the Besov approximation bound by real interpolation. The approximants also satisfy a squared RKHS norm bound of the form C γ − ϱ ∥ f ∥ L 2 ( X ) 2 C\gamma^{-\varrho}\|f\|_{L_2(X)}^2 C γ − ϱ ∥ f ∥ L 2 ( X ) 2 , whose bandwidth exponent is governed by the intrinsic dimension and whose dependence on the target is through its L 2 ( X ) L_2(X) L 2 ( X ) norm alone. At the same smoothness order, this extension admits a larger class of targets, including potentially unbounded functions at low smoothness, and supplies the approximation estimates needed for logistic risk when uniform approximation is unavailable. Related approximation and regression guarantees for deep neural networks over Besov classes are also established in [19 , 20 , 23 ] .
In this paper, we study the statistical learning problem of binary classification with logistic loss
ϕ ( t ) = log ( 1 + e − t ) . \phi(t) = \log(1 + e^{-t}). ϕ ( t ) = log ( 1 + e − t ) .
The Bayes function
f ∗ = log η 1 − η f^* = \log\frac{\eta}{1-\eta} f ∗ = log 1 − η η
can be unbounded when the conditional class probability η \eta η approaches zero or one. Standard oracle inequalities based on fixed-level clipping, such as [15 Theorem 7.22] , are not directly applicable to logistic loss. Some existing analyses, such as [9 ] , impose a boundedness assumption on the Bayes function f ∗ f^* f ∗ , which is equivalent to requiring the conditional class probability to stay uniformly away from zero and one. A fundamental difficulty in the analysis of logistic loss is therefore to obtain effective approximation and concentration bounds without imposing such a boundedness assumption on the Bayes function f ∗ f^* f ∗ .
Logistic classification with Gaussian kernels is studied directly in [26 ] . It combines Sobolev approximation with a truncation level that can increase with the sample size, and its general analysis already addresses potentially unbounded targets. Its stochastic analysis does not exploit variance bounds that shrink with the excess risk. Combined with its approximation and truncation estimates, this yields a slow rate for the misclassification error. Previous works [1 , 9 ] consider the variance estimate, but with an exponential dependence on the supremum bound for logistic loss. Hence their bound can deteriorate rapidly when regularization decreases.
The concentration argument in this paper follows the localization strategy developed in [16 ] . That analysis uses the regularized excess risk to control both the norm of the candidate functions and the variance of their loss differences, and then combines local complexity estimates with concentration. We adapt this framework using a logistic Bernstein-type inequality whose constant grows only linearly with the supremum norm of the predictor. This variance control permits an oracle analysis without a uniform bound on the Bayes function f ∗ f^* f ∗ and connects the statistical estimate to our manifold approximation bounds. We further use this linear variance dependence to obtain faster rates by clipping predictions at a level that grows logarithmically with the sample size.
The unboundedness issue also motivates the work [32 ] on classification with deep neural networks and logistic loss. Their result [32 Theorem 2.1] introduces an oracle-type inequality with a flexible auxiliary reference loss. A carefully chosen bounded reference preserves the Bayes risk in expectation while controlling the second moments needed for concentration. They also provide an efficient ReLU approximation of the logarithm near zero, with only polylogarithmic dependence on the inverse distance to the singularity in the network complexity. See [32 Theorem 2.4] . These tools yield learning rates under regularity assumptions on the conditional class probability, without requiring a bounded Bayes function f ∗ f^* f ∗ .
In this paper, we study norm-regularized Gaussian kernel estimators on a compact Riemannian manifold X X X without boundary and of intrinsic dimension ϱ \varrho ϱ . Our main contributions are as follows.
We establish an oracle inequality for logistic loss without assuming a uniform bound on the Bayes function f ∗ f^* f ∗ . The proof uses a Bernstein-type bound whose constant grows only linearly with the supremum norm of the predictor. Theorem 3.2 applies to separable normed function spaces, with constants uniform in the relevant entropy exponent. Corollary 3.3 further allows covering estimates that are available only above a positive radius b b b , retaining the additional error b 2 / λ b^2/\lambda b 2 / λ .
We prove approximation bounds for Besov functions on the manifold by constructing a signed multiscale Gaussian approximation S γ , s f S_{\gamma,s}f S γ , s f that cancels the lower-order bias terms in the Sobolev expansion. For every α > 0 \alpha>0 α > 0 , Theorems 4.1 and 4.3 give
∥ S γ , s f ∥ H γ ( X ) ≤ C γ − ϱ / 2 ∥ f ∥ L 2 ( X ) , \|S_{\gamma,s}f\|_{\mathcal{H}_\gamma(X)} \le C\gamma^{-\varrho/2}\|f\|_{L_2(X)}, ∥ S γ , s f ∥ H γ ( X ) ≤ C γ − ϱ /2 ∥ f ∥ L 2 ( X ) ,
∥ S γ , s f − f ∥ L 2 ( X ) ≤ C γ α ∥ f ∥ B 2 , ∞ α ( X ) . \|S_{\gamma,s}f-f\|_{L_2(X)} \le C\gamma^\alpha\|f\|_{B^\alpha_{2,\infty}(X)}. ∥ S γ , s f − f ∥ L 2 ( X ) ≤ C γ α ∥ f ∥ B 2 , ∞ α ( X ) .
Thus the construction removes second-order saturation while retaining a squared RKHS norm bound of order γ − ϱ \gamma^{-\varrho} γ − ϱ . This approximation result is of independent interest and warrants special attention beyond its application to logistic classification.
We combine the oracle and approximation bounds with the manifold covering estimate in [28 Lemma 5.4] to derive, for a density g ∈ L p ( X ) g \in L_p(\mathcal{X}) g ∈ L p ( X ) with 2 ≤ p ≤ ∞ 2 \le p \le\infty 2 ≤ p ≤ ∞ , the logistic excess-risk rate
( ( log n ) ℓ + 1 n ) 4 α ( 1 − 1 / p ) 4 α ( 1 − 1 / p ) + 3 ℓ \left(\frac{(\log n)^{\ell+1}}{n}\right)^{\frac{4\alpha(1-1/p)}{4\alpha(1-1/p)+3\ell}} ( n ( log n ) ℓ + 1 ) 4 α ( 1 − 1/ p ) + 3 ℓ 4 α ( 1 − 1/ p )
Its polynomial and logarithmic exponents depend on the intrinsic dimension ℓ \ell ℓ , the smoothness α \alpha α of the Bayes function f ∗ f^* f ∗ , and the density integrability p p p , but not on the ambient dimension d d d . See Theorem 4.7 .
( ( log n ) ℓ + 2 n ) 2 α ( 1 − 1 / p ) 2 α ( 1 − 1 / p ) + ℓ \left(\frac{(\log n)^{\ell+2}}{n}\right)^{\frac{2\alpha(1-1/p)}{2\alpha(1-1/p)+\ell}} ( n ( log n ) ℓ + 2 ) 2 α ( 1 − 1/ p ) + ℓ 2 α ( 1 − 1/ p )
For bounded densities, the logistic excess-risk bound for the clipped predictor has the same polynomial exponent 2 α / ( 2 α + ℓ ) 2\alpha/(2\alpha+\ell) 2 α / ( 2 α + ℓ ) as the minimax optimal rate for squared-loss regression discussed in Remark 4.2, and the excess misclassification risk also matches the minimax optimal rate n − α / ( 2 α + ℓ ) n^{-\alpha/(2\alpha+\ell)} n − α / ( 2 α + ℓ ) for binary classification without an additional noise condition.
The rest of this paper is organized as follows. § 2 introduces the problem setting and related notations. § 3 presents the Bernstein-type bound and the oracle inequalities for logistic loss. § 4 establishes Gaussian approximation bounds for Besov functions on manifolds and derives learning rates for both the original and clipped predictions. § 5 summarizes the results and discusses further questions. All proofs are collected in Appendix § 6 .
Preliminaries
This section reviews some basic notions of the logistic classification problem, function spaces on manifolds, Gaussian RKHSs and complexity measures of function classes.
The Logistic Classification Problem
We consider the binary classification problem where we regard a separable metric space X \mathcal{X} X as the input space and Y = { − 1 , 1 } \mathcal{Y}=\{-1,1\} Y = { − 1 , 1 } as the label space. Let P P P be a probability distribution on X × Y \mathcal{X}\times\mathcal{Y} X × Y . We denote by P X P_{\mathcal{X}} P X the marginal distribution of P P P on X \mathcal{X} X , and we fix a measurable version of the conditional probability η ( x ) : = P ( Y = 1 ∣ X = x ) \eta(x):=P(Y=1\mid X=x) η ( x ) := P ( Y = 1 ∣ X = x ) . The logistic loss function is defined by
ϕ ( t ) : = log ( 1 + e − t ) . \phi(t):=\log(1+e^{-t}). ϕ ( t ) := log ( 1 + e − t ) .
We define the logistic risk of a measurable function f : X → R f:\mathcal{X}\to\mathbb{R} f : X → R by
R ( f ) : = E P ϕ ( Y f ( X ) ) . \mathcal{R}(f):=\mathbb{E}_{P}\phi(Yf(X)). R ( f ) := E P ϕ ( Y f ( X )) .
With the convention 0 log 0 = 0 0\log0=0 0 log 0 = 0 , the Bayes logistic risk is
R ∗ : = inf f : X → R f m e a s u r a b l e R ( f ) = ∫ X ( − η log η − ( 1 − η ) log ( 1 − η ) ) d P X . \mathcal{R}^*:=\inf_{\substack{f:\mathcal{X}\to\mathbb{R}\\ f\ \mathrm{measurable}}}\mathcal{R}(f)=\int_{\mathcal{X}}\left(-\eta\log\eta-(1-\eta)\log(1-\eta)\right)dP_{\mathcal{X}}. R ∗ := f : X → R f measurable inf R ( f ) = ∫ X ( − η log η − ( 1 − η ) log ( 1 − η ) ) d P X .
When 0 < η ( x ) < 1 0<\eta(x)<1 0 < η ( x ) < 1 , the pointwise minimizer is
f ∗ ( x ) : = log η ( x ) 1 − η ( x ) . f^*(x):=\log\frac{\eta(x)}{1-\eta(x)}. f ∗ ( x ) := log 1 − η ( x ) η ( x ) .
We do not regard f ∗ f^* f ∗ as a real-valued function at points where η \eta η is 0 or 1. Instead, define the measurable Bayes reference loss
ψ ( x , y ) : = { ϕ ( y f ∗ ( x ) ) , 0 < η ( x ) < 1 , 0 , η ( x ) ∈ { 0 , 1 } . \psi(x,y):=
\begin{cases}
\phi(yf^*(x)), & 0<\eta(x)<1,\\
0, & \eta(x)\in\{0,1\}.
\end{cases} ψ ( x , y ) := { ϕ ( y f ∗ ( x )) , 0 , 0 < η ( x ) < 1 , η ( x ) ∈ { 0 , 1 } .
Then E P ψ ( X , Y ) = R ∗ \mathbb{E}_{P}\psi(X,Y)=R^* E P ψ ( X , Y ) = R ∗ , including when P X { η ∈ { 0 , 1 } } > 0 P_X\{\eta\in\{0,1\}\}>0 P X { η ∈ { 0 , 1 }} > 0 .
To learn a classifier from n n n i.i.d. samples D = ( ( X 1 , Y 1 ) , … , ( X n , Y n ) ) D=((X_1,Y_1),\ldots,(X_n,Y_n)) D = (( X 1 , Y 1 ) , … , ( X n , Y n )) , we minimize the empirical risk on a hypothesis class F \mathcal{F} F :
f ^ D ∈ arg min f ∈ F 1 n ∑ i = 1 n ϕ ( Y i f ( X i ) ) . \widehat{f}_D\in\underset{f\in\mathcal{F}}{\arg\min}\ \frac{1}{n}\sum_{i=1}^{n}\phi(Y_i f(X_i)). f D ∈ f ∈ F arg min n 1 i = 1 ∑ n ϕ ( Y i f ( X i )) .
or, usually, a regularized version:
f ^ D , λ ∈ arg min f ∈ F { 1 n ∑ i = 1 n ϕ ( Y i f ( X i ) ) + λ ∥ f ∥ F 2 } . (1) \widehat{f}_{D,\lambda}\in\underset{f\in\mathcal{F}}{\arg\min}\left\{\frac{1}{n}\sum_{i=1}^{n}\phi(Y_i f(X_i))+\lambda\lVert f\rVert_{\mathcal{F}}^2\right\}. \tag*{(1)} f D , λ ∈ f ∈ F arg min { n 1 i = 1 ∑ n ϕ ( Y i f ( X i )) + λ ∥ f ∥ F 2 } . (1)
where λ > 0 \lambda>0 λ > 0 is a regularization parameter and ∥ ⋅ ∥ F \lVert\cdot\rVert_{\mathcal{F}} ∥ ⋅ ∥ F is a norm on F \mathcal{F} F . We aim to establish an oracle inequality to bound the excess risk R ( f ^ D , λ ) − R ∗ R(\widehat{f}_{D,\lambda})-R^* R ( f D , λ ) − R ∗ in terms of the sample size n n n , the regularization parameter λ \lambda λ , and the complexity of the hypothesis class F \mathcal{F} F . In particular, by choosing the regularization parameter λ n \lambda_n λ n and the hypothesis class F n \mathcal{F}_n F n appropriately depending on n n n , we can obtain an asymptotic convergence rate for the excess risk as n → ∞ n\to\infty n → ∞ , commonly called a learning rate.
The logistic loss serves as a convex surrogate for the 0–1 classification loss to facilitate optimization. For a measurable decision function f : X → R f:\mathcal{X}\to\mathbb{R} f : X → R , the misclassification risk and the Bayes risk are defined by
R 01 ( f ) : = E P 1 { sgn ( f ( X ) ) ≠ Y } , R_{01}(f):=\mathbb{E}_{P}\mathbf{1}_{\{\operatorname{sgn}(f(X))\ne Y\}}, R 01 ( f ) := E P 1 { sgn ( f ( X )) = Y } ,
R 01 ∗ : = ∫ X min { η ( x ) , 1 − η ( x ) } d P X ( x ) . R^*_{01}:=\int_{\mathcal{X}}\min\{\eta(x),1-\eta(x)\}\,dP_X(x). R 01 ∗ := ∫ X min { η ( x ) , 1 − η ( x )} d P X ( x ) .
Here we agree that sgn ( 0 ) = 1 \operatorname{sgn}(0)=1 sgn ( 0 ) = 1 . R 01 ( f ) R_{01}(f) R 01 ( f ) directly reflects the performance of classifier f f f . We therefore use comparison or calibration inequalities to translate bounds on the excess logistic risk into bounds on the excess misclassification risk. Such inequalities are widely studied in statistical learning theory. See [4 , 7 , 15 , 31 ] . For the logistic loss ϕ \phi ϕ , by [4 ] , the calibration inequality takes the form
R 01 ( f ) − R 01 ∗ ≤ C R ( f ) − R ∗ . (2) R_{01}(f)-R^*_{01}\le C\sqrt{R(f)-R^*}. \tag*{(2)} R 01 ( f ) − R 01 ∗ ≤ C R ( f ) − R ∗ . (2)
This calibration inequality can be improved when the distribution P P P satisfies the Tsybakov noise condition [24 ] . See a detailed discussion in Remark 4.1.
Function Spaces on Manifolds
For a data domain X ⊂ R d \mathcal{X}\subset\mathbb{R}^d X ⊂ R d , it is well-known that typical learning rates are subject to the curse of dimensionality, and the sample size n n n required to achieve a given accuracy can grow rapidly with the dimension d d d . For example, the minimax nonparametric squared-loss regression rate for f ∗ ∈ C α ( [ 0 , 1 ] d ) f^*\in C^\alpha([0,1]^d) f ∗ ∈ C α ([ 0 , 1 ] d ) is n − 2 α / ( 2 α + d ) n^{-2\alpha/(2\alpha+d)} n − 2 α / ( 2 α + d ) , cf. [18 ] . To address the curse of dimensionality, a common assumption is that the marginal distribution P X P_X P X is supported on a low-dimensional manifold X ⊂ R d \mathcal{X}\subset\mathbb{R}^d X ⊂ R d with intrinsic dimension ϱ ≪ d \varrho\ll d ϱ ≪ d . This motivates the study of learning and approximation on manifolds.
We now specify the manifold assumption and introduce the function spaces used to measure regularity on it. Suppose that X ⊂ R d \mathcal{X}\subset\mathbb{R}^d X ⊂ R d is a connected compact C ∞ C^\infty C ∞ embedded submanifold without boundary, of intrinsic dimension 1 ≤ ϱ < d 1\le\varrho<d 1 ≤ ϱ < d , and is equipped with the metric induced by the Euclidean embedding. We write d X d_{\mathcal{X}} d X , d V dV d V , E x \mathcal{E}_x E x , and inj ( X ) \operatorname{inj}(\mathcal{X}) inj ( X ) for its geodesic distance, volume measure, exponential map at x ∈ X x\in\mathcal{X} x ∈ X , and injectivity radius. Compactness implies inj ( X ) > 0 \operatorname{inj}(\mathcal{X})>0 inj ( X ) > 0 and provides a finite collection of relatively compact coordinate neighborhoods carrying smooth orthonormal frames. The Sobolev and Besov spaces provide a regularity framework for stating our approximation bounds and learning rates. For an integer m ≥ 1 m\ge1 m ≥ 1 , let H m ( X ) H^m(\mathcal{X}) H m ( X ) be the Sobolev space equipped with the equivalent norm
∥ f ∥ H m ( X ) 2 = ∑ ℓ = 0 m ∫ X ∣ ∇ ℓ f ( x ) ∣ 2 d V ( x ) . \lVert f\rVert_{H^m(\mathcal{X})}^2=\sum_{\ell=0}^{m}\int_{\mathcal{X}}\lvert\nabla^\ell f(x)\rvert^2\,dV(x). ∥ f ∥ H m ( X ) 2 = ℓ = 0 ∑ m ∫ X ∣ ∇ ℓ f ( x ) ∣ 2 d V ( x ) .
Here, ∇ \nabla ∇ denotes the Levi-Civita covariant derivative, and all pointwise tensor norms are induced by the Riemannian metric on X X X . For f ∈ L 2 ( X ) + H m ( X ) f \in L_2(X) + H^m(X) f ∈ L 2 ( X ) + H m ( X ) , define the Peetre K K K -functional by
K ( t , f ; L 2 ( X ) , H m ( X ) ) : = inf f = f 0 + f 1 f 0 ∈ L 2 ( X ) , f 1 ∈ H m ( X ) { ∥ f 0 ∥ L 2 ( X ) + t ∥ f 1 ∥ H m ( X ) } . K(t,f;L_2(X),H^m(X)):=\inf_{\substack{f=f_0+f_1\\ f_0\in L_2(X),\ f_1\in H^m(X)}}\{\lVert f_0\rVert_{L_2(X)}+t\lVert f_1\rVert_{H^m(X)}\}. K ( t , f ; L 2 ( X ) , H m ( X )) := f = f 0 + f 1 f 0 ∈ L 2 ( X ) , f 1 ∈ H m ( X ) inf {∥ f 0 ∥ L 2 ( X ) + t ∥ f 1 ∥ H m ( X ) } .
Given α > 0 \alpha>0 α > 0 , choose any integer m > α m>\alpha m > α and define
B 2 , ∞ α ( X ) : = ( L 2 ( X ) , H m ( X ) ) α / m , ∞ , B^\alpha_{2,\infty}(X):=(L_2(X),H^m(X))_{\alpha/m,\infty}, B 2 , ∞ α ( X ) := ( L 2 ( X ) , H m ( X ) ) α / m , ∞ ,
∥ f ∥ B 2 , ∞ α ≍ sup t > 0 t − α / m K ( t , f ; L 2 ( X ) , H m ( X ) ) . (3) \lVert f\rVert_{B^\alpha_{2,\infty}}\asymp\sup_{t>0}t^{-\alpha/m}K(t,f;L_2(X),H^m(X)). \tag*{(3)} ∥ f ∥ B 2 , ∞ α ≍ t > 0 sup t − α / m K ( t , f ; L 2 ( X ) , H m ( X )) . (3)
The real-interpolation identity in (3) and the independence, up to equivalence of norms, of the choice of m > α m>\alpha m > α are standard. See [5 Chapter 3] and [22 Chapter 7] .
We also recall the chartwise definition of the classical Hölder spaces. Fix a finite smooth atlas { ( U ν , ϕ ν ) } ν = 1 J \{(U_\nu,\phi_\nu)\}_{\nu=1}^J {( U ν , ϕ ν ) } ν = 1 J and a smooth subordinate partition of unity { χ ν } ν = 1 J \{\chi_\nu\}_{\nu=1}^J { χ ν } ν = 1 J such that supp χ ν ⋐ U ν \operatorname{supp}\chi_\nu\Subset U_\nu supp χ ν ⋐ U ν . For k ∈ N 0 k\in\mathbb{N}_0 k ∈ N 0 and 0 < β ≤ 1 0<\beta\leq1 0 < β ≤ 1 , set
∥ f ∥ C k , β ( X ) : = ∑ ν = 1 J ∥ ( χ ν f ) ∘ ϕ ν − 1 ∥ C k , β ( R ϱ ) . \lVert f\rVert_{C^{k,\beta}(X)}:=\sum_{\nu=1}^{J}\lVert(\chi_\nu f)\circ\phi_\nu^{-1}\rVert_{C^{k,\beta}(\mathbb{R}^{\varrho})}. ∥ f ∥ C k , β ( X ) := ν = 1 ∑ J ∥( χ ν f ) ∘ ϕ ν − 1 ∥ C k , β ( R ϱ ) .
where each localized function is extended by zero from ϕ ν ( U ν ) \phi_\nu(U_\nu) ϕ ν ( U ν ) to R ϱ \mathbb{R}^{\varrho} R ϱ and
∥ u ∥ C k , β ( R ϱ ) : = ∑ ∣ a ∣ ≤ k ∥ D a u ∥ ∞ + ∑ ∣ a ∣ = k sup z ≠ z ′ ∣ D a u ( z ) − D a u ( z ′ ) ∣ ∣ z − z ′ ∣ β . \lVert u\rVert_{C^{k,\beta}(\mathbb{R}^{\varrho})}:=\sum_{\lvert a\rvert\leq k}\lVert D^a u\rVert_\infty+\sum_{\lvert a\rvert=k}\sup_{z\neq z'}\frac{\lvert D^a u(z)-D^a u(z')\rvert}{\lvert z-z'\rvert^\beta}. ∥ u ∥ C k , β ( R ϱ ) := ∣ a ∣ ≤ k ∑ ∥ D a u ∥ ∞ + ∣ a ∣ = k ∑ z = z ′ sup ∣ z − z ′ ∣ β ∣ D a u ( z ) − D a u ( z ′ )∣ .
Thus C k , β ( X ) C^{k,\beta}(X) C k , β ( X ) consists of the functions for which this norm is finite. Different finite atlases and subordinate partitions give equivalent norms.
There is also a chartwise difference characterization of the Besov space B ∞ , ∞ α B^\alpha_{\infty,\infty} B ∞ , ∞ α at every order α > 0 \alpha>0 α > 0 . Let
k : = ⌈ α ⌉ − 1 , β : = α − k ∈ ( 0 , 1 ] , k:=\lceil\alpha\rceil-1,\qquad\beta:=\alpha-k\in(0,1], k := ⌈ α ⌉ − 1 , β := α − k ∈ ( 0 , 1 ] ,
and, for u ∈ C k ( R ϱ ) u\in C^k(\mathbb{R}^{\varrho}) u ∈ C k ( R ϱ ) , set
∥ u ∥ B ∞ , ∞ α ( R ϱ ) : = ∑ ∣ a ∣ ≤ k ∥ D a u ∥ ∞ + ∑ ∣ a ∣ = k sup h ≠ 0 ∣ h ∣ − β ∥ Δ h 2 D a u ∥ ∞ , \lVert u\rVert_{B^\alpha_{\infty,\infty}(\mathbb{R}^{\varrho})}:=\sum_{\lvert a\rvert\leq k}\lVert D^a u\rVert_\infty+\sum_{\lvert a\rvert=k}\sup_{h\neq0}\lvert h\rvert^{-\beta}\lVert\Delta_h^2D^a u\rVert_\infty, ∥ u ∥ B ∞ , ∞ α ( R ϱ ) := ∣ a ∣ ≤ k ∑ ∥ D a u ∥ ∞ + ∣ a ∣ = k ∑ h = 0 sup ∣ h ∣ − β ∥ Δ h 2 D a u ∥ ∞ ,
Δ h 2 u ( z ) : = u ( z + 2 h ) − 2 u ( z + h ) + u ( z ) . \Delta_h^2u(z):=u(z+2h)-2u(z+h)+u(z). Δ h 2 u ( z ) := u ( z + 2 h ) − 2 u ( z + h ) + u ( z ) .
Then, using the same atlas and partition of unity, define
∥ f ∥ B ∞ , ∞ α ( X ) : = ∑ ν = 1 J ∥ ( χ ν f ) ∘ ϕ ν − 1 ∥ B ∞ , ∞ α ( R ϱ ) . \lVert f\rVert_{B^\alpha_{\infty,\infty}(X)}:=\sum_{\nu=1}^{J}\lVert(\chi_\nu f)\circ\phi_\nu^{-1}\rVert_{B^\alpha_{\infty,\infty}(\mathbb{R}^{\varrho})}. ∥ f ∥ B ∞ , ∞ α ( X ) := ν = 1 ∑ J ∥( χ ν f ) ∘ ϕ ν − 1 ∥ B ∞ , ∞ α ( R ϱ ) .
with zero extension. This defines B ∞ , ∞ α ( X ) B^\alpha_{\infty,\infty}(X) B ∞ , ∞ α ( X ) . Up to equivalence of norms, it is independent of the atlas and partition. See [22 Section 1.2.2, Definition 7.2.2, Theorem 7.2.3, and Theorem 7.3.1(ii)] . If α = k + β \alpha=k+\beta α = k + β with 0 < β < 1 0<\beta<1 0 < β < 1 , this is the classical C k , β ( X ) C^{k,\beta}(X) C k , β ( X ) space with an equivalent norm. If α = k + 1 \alpha=k+1 α = k + 1 , however, B ∞ , ∞ k + 1 ( X ) B^{k+1}_{\infty,\infty}(X) B ∞ , ∞ k + 1 ( X ) is the Zygmund class which contains C k , 1 ( X ) C^{k,1}(X) C k , 1 ( X ) continuously and is generally larger.
Using B ∞ , ∞ α ( X ) B^\alpha_{\infty,\infty}(X) B ∞ , ∞ α ( X ) as a bridge, we can prove that the Besov class B 2 , ∞ α ( X ) B^\alpha_{2,\infty}(X) B 2 , ∞ α ( X ) is indeed very rich and contains the Hölder class C k , β ( X ) C^{k,\beta}(X) C k , β ( X ) continuously.
Proposition 2.1 (Hölder–Besov embedding). Let k ∈ N 0 k\in\mathbb{N}_0 k ∈ N 0 , β ∈ ( 0 , 1 ] \beta\in(0,1] β ∈ ( 0 , 1 ] , and α = k + β \alpha=k+\beta α = k + β . Then
C k , β ( X ) ↪ B 2 , ∞ α ( X ) . C^{k,\beta}(X)\hookrightarrow B^\alpha_{2,\infty}(X). C k , β ( X ) ↪ B 2 , ∞ α ( X ) . More precisely, there is a constant C = C ( X , k , β ) C=C(X,k,\beta) C = C ( X , k , β ) such that
∥ f ∥ B 2 , ∞ α ( X ) ≤ C ∥ f ∥ C k , β ( X ) for all f ∈ C k , β ( X ) . \lVert f\rVert_{B^\alpha_{2,\infty}(X)}\leq C\lVert f\rVert_{C^{k,\beta}(X)}\quad\text{for all }f\in C^{k,\beta}(X). ∥ f ∥ B 2 , ∞ α ( X ) ≤ C ∥ f ∥ C k , β ( X ) for all f ∈ C k , β ( X ) .
Gaussian RKHSs
For a bandwidth γ > 0 \gamma>0 γ > 0 , let H γ ( X ) \mathcal{H}_\gamma(X) H γ ( X ) be the RKHS of the Gaussian kernel
k γ ( x , y ) : = exp ( − ∥ x − y ∥ R d 2 γ 2 ) , x , y ∈ X . k_\gamma(x,y):=\exp\left(-\frac{\|x-y\|_{\mathbb{R}^d}^2}{\gamma^2}\right),\qquad x,y\in X. k γ ( x , y ) := exp ( − γ 2 ∥ x − y ∥ R d 2 ) , x , y ∈ X .
H γ ( X ) \mathcal{H}_\gamma(X) H γ ( X ) is separable and consists of continuous functions. Since k γ ( x , x ) = 1 k_\gamma(x,x)=1 k γ ( x , x ) = 1 , the reproducing property gives
∥ f ∥ ∞ ≤ ∥ f ∥ H γ ( X ) , f ∈ H γ ( X ) . \|f\|_\infty\le\|f\|_{\mathcal{H}_\gamma(X)},\qquad f\in\mathcal{H}_\gamma(X). ∥ f ∥ ∞ ≤ ∥ f ∥ H γ ( X ) , f ∈ H γ ( X ) .
For a detailed analysis of H γ ( X ) \mathcal{H}_\gamma(X) H γ ( X ) , see [15 ] , Section 4.4. Its separability and the norm bound above verify the function-space assumptions used in the oracle inequality developed in § 3 below. For fixed n n n , λ > 0 \lambda>0 λ > 0 , and γ > 0 \gamma>0 γ > 0 , the Gaussian regularized empirical objective is coercive and weakly lower semicontinuous, so it admits a minimizer. Strong convexity makes this minimizer unique, and continuity of the kernel implies continuous dependence on the sample. The resulting estimator f ^ D , λ , γ \hat f_{D,\lambda,\gamma} f ^ D , λ , γ is therefore measurable with respect to D D D .
Complexity Measures of Function Classes
To describe the complexity of the function class, we introduce dyadic entropy numbers, covering numbers and empirical Rademacher complexity. For a bounded linear operator T : E → F T:E\to F T : E → F , its i i i th dyadic entropy number is
e i ( T ) : = inf { δ > 0 : T ( B E ) is covered by 2 i − 1 balls of radius δ in F } . e_i(T):=\inf\{\delta>0:T(B_E)\text{ is covered by }2^{i-1}\text{ balls of radius }\delta\text{ in }F\}. e i ( T ) := inf { δ > 0 : T ( B E ) is covered by 2 i − 1 balls of radius δ in F } .
Here B E B_E B E is the unit ball in E E E . For an arbitrary subset m a t h c a l Q ⊂ F mathcal{Q}\subset F ma t h c a l Q ⊂ F , we use the analogous notation
e i ( Q , F ) : = inf { δ > 0 : Q is covered by 2 i − 1 balls of radius δ in F } . e_i(\mathcal{Q},F):=\inf\{\delta>0:\mathcal{Q}\text{ is covered by }2^{i-1}\text{ balls of radius }\delta\text{ in }F\}. e i ( Q , F ) := inf { δ > 0 : Q is covered by 2 i − 1 balls of radius δ in F } .
For an arbitrary subset Q ⊂ F \mathcal{Q}\subset F Q ⊂ F , define covering numbers
N ( δ , Q , ∥ ⋅ ∥ F ) : = min { n ≥ 1 : ∃ q 1 , … , q n ∈ F such that Q ⊂ ⋃ i = 1 n B ( q i , δ ) } . N(\delta,\mathcal{Q},\|\cdot\|_F):=\min\left\{n\ge1:\exists q_1,\ldots,q_n\in F\text{ such that }\mathcal{Q}\subset\bigcup_{i=1}^{n}B(q_i,\delta)\right\}. N ( δ , Q , ∥ ⋅ ∥ F ) := min { n ≥ 1 : ∃ q 1 , … , q n ∈ F such that Q ⊂ i = 1 ⋃ n B ( q i , δ ) } .
For a class Q \mathcal{Q} Q of measurable functions on X × Y X\times Y X × Y , define its empirical Rademacher complexity on the sample D D D by
Rad D ( Q ) : = E σ sup q ∈ Q ∣ 1 n ∑ i = 1 n σ i q ( X i , Y i ) ∣ . \operatorname{Rad}_D(\mathcal{Q}):=\mathbb{E}_{\sigma}\sup_{q\in\mathcal{Q}}\left|\frac{1}{n}\sum_{i=1}^{n}\sigma_iq(X_i,Y_i)\right|. Rad D ( Q ) := E σ q ∈ Q sup n 1 i = 1 ∑ n σ i q ( X i , Y i ) .
where σ 1 , … , σ n \sigma_1,\ldots,\sigma_n σ 1 , … , σ n are independent Rademacher variables, each taking the values − 1 -1 − 1 and 1 1 1 with equal probability, and are independent of D D D .
Our oracle inequality imposes a polynomial decay expected empirical dyadic entropy condition on the hypothesis class F \mathcal{F} F . In particular, for each fixed bandwidth γ > 0 \gamma> 0 γ > 0 , the dyadic entropy numbers of the Gaussian RKHS H γ ( X ) \mathcal{H}_{\gamma}(X) H γ ( X ) exhibit super-polynomial decay. This flexibility allows us to optimize the entropy exponent to obtain the best learning rate for the Gaussian kernel estimator.
An Oracle Inequality for Logistic Loss with Potentially Unbounded Bayes Functions
To derive an oracle inequality for regularized empirical risk minimization, typical arguments combine a Bernstein-type variance condition with a capacity assumption for the hypothesis class and concentration inequalities such as Talagrand’s inequality. The variance condition controls the second moment of the excess loss in terms of the excess risk. For logistic loss, obtaining such a condition with a useful constant is nontrivial. The results in [9 Lemma 8] and [1 Proposition 8.1] yield variance bounds with constants that grow exponentially with the uniform bound ∥ f ∥ ∞ \|f\|_{\infty} ∥ f ∥ ∞ . This dependence can lead to poor statistical guarantees when ∥ f ∥ ∞ \|f\|_{\infty} ∥ f ∥ ∞ increases. This is particularly relevant for logistic loss, since the Bayes function f ∗ ( x ) → − ∞ f^{*}(x) \to-\infty f ∗ ( x ) → − ∞ as η ( x ) → 0 \eta(x) \to0 η ( x ) → 0 and f ∗ ( x ) → + ∞ f^{*}(x) \to+\infty f ∗ ( x ) → + ∞ as η ( x ) → 1 \eta(x) \to1 η ( x ) → 1 .
The following lemma gives a sharp Bernstein-type bound with only linear dependence on ∥ f ∥ ∞ \|f\|_{\infty} ∥ f ∥ ∞ , without requiring the conditional probability η \eta η to be uniformly bounded away from 0 and 1.
Lemma 3.1 (Bernstein-type bound for logistic loss). For every bounded measurable f : X → R f : X \to\mathbb{R} f : X → R ,
E P ( ϕ ( Y f ( X ) ) − ψ ( X , Y ) ) 2 ≤ ( ∥ f ∥ ∞ + 2 ) ( R ( f ) − R ∗ ) . \mathbb{E}_{P}\left(\phi(Yf(X))-\psi(X,Y)\right)^{2} \leq\left(\|f\|_{\infty}+2\right)\left(\mathcal{R}(f)-\mathcal{R}^{*}\right). E P ( ϕ ( Y f ( X )) − ψ ( X , Y ) ) 2 ≤ ( ∥ f ∥ ∞ + 2 ) ( R ( f ) − R ∗ ) .
We now state an oracle inequality for regularized empirical risk minimization over a separable normed function space. Using the Bernstein-type bound in Lemma 3.1, we do not require the conditional probability η \eta η to be bounded away from 0 and 1. For D = ( ( X 1 , Y 1 ) , … , ( X n , Y n ) ) D=((X_{1},Y_{1}),\ldots,(X_{n},Y_{n})) D = (( X 1 , Y 1 ) , … , ( X n , Y n )) , let L 2 ( D ) L_{2}(D) L 2 ( D ) be equipped with the empirical norm
∥ g ∥ L 2 ( D ) 2 : = 1 n ∑ i = 1 n ∣ g ( X i , Y i ) ∣ 2 , g : X × Y → R . \|g\|_{L_{2}(D)}^{2} := \frac{1}{n}\sum_{i=1}^{n}|g(X_{i},Y_{i})|^{2},\qquad g:X\times Y\to\mathbb{R}. ∥ g ∥ L 2 ( D ) 2 := n 1 i = 1 ∑ n ∣ g ( X i , Y i ) ∣ 2 , g : X × Y → R .
Similarly, for D X = ( X 1 , … , X n ) D_{X}=(X_{1},\ldots,X_{n}) D X = ( X 1 , … , X n ) , let L 2 ( D X ) L_{2}(D_{X}) L 2 ( D X ) be equipped with the empirical norm
∥ f ∥ L 2 ( D X ) 2 : = 1 n ∑ i = 1 n ∣ f ( X i ) ∣ 2 , f : X → R . \|f\|_{L_{2}(D_{X})}^{2} := \frac{1}{n}\sum_{i=1}^{n}|f(X_{i})|^{2},\qquad f:X\to\mathbb{R}. ∥ f ∥ L 2 ( D X ) 2 := n 1 i = 1 ∑ n ∣ f ( X i ) ∣ 2 , f : X → R .
Theorem 3.2 (Oracle inequality for logistic loss). Let F \mathcal{F} F be a separable normed linear space of bounded measurable functions f : X → R f:X\to\mathbb{R} f : X → R such that
∥ f ∥ ∞ ≤ ∥ f ∥ F , f ∈ F . \|f\|_{\infty}\leq\|f\|_{\mathcal{F}},\qquad f\in\mathcal{F}. ∥ f ∥ ∞ ≤ ∥ f ∥ F , f ∈ F . For every n ≥ 1 n\geq1 n ≥ 1 and 0 < λ ≤ 1 0<\lambda\leq1 0 < λ ≤ 1 , assume that the regularized empirical risk minimization problem (1) admits a minimizer for every sample D D D , and fix a choice f ^ D , λ ∈ F \widehat{f}_{D,\lambda}\in\mathcal{F} f D , λ ∈ F that is measurable as a function of D D D . Define the approximation error by
a F ( λ ) : = inf f ∈ F { λ ∥ f ∥ F 2 + R ( f ) − R ∗ } . (4) a_{\mathcal{F}}(\lambda):=\inf_{f\in\mathcal{F}}\left\{\lambda\|f\|_{\mathcal{F}}^{2}+\mathcal{R}(f)-\mathcal{R}^{*}\right\}. \tag*{(4)} a F ( λ ) := f ∈ F inf { λ ∥ f ∥ F 2 + R ( f ) − R ∗ } . (4) Furthermore, suppose that, for some p ∈ ( 0 , 1 / 2 ] p\in(0,1/2] p ∈ ( 0 , 1/2 ] and A ≥ 1 A\geq1 A ≥ 1 , for every n ≥ 1 n\geq1 n ≥ 1 ,
E D X ∼ P X n e i ( id : F → L 2 ( D X ) ) ≤ A 1 / ( 2 p ) i − 1 / ( 2 p ) , i ≥ 1. (5) \mathbb{E}_{D_{X}\sim P_{X}^{n}}e_{i}\left(\operatorname{id}:\mathcal{F}\to L_{2}(D_{X})\right)\leq A^{1/(2p)}i^{-1/(2p)},\qquad i\geq1. \tag*{(5)} E D X ∼ P X n e i ( id : F → L 2 ( D X ) ) ≤ A 1/ ( 2 p ) i − 1/ ( 2 p ) , i ≥ 1. (5) Then there is a numerical constant C ≥ 1 C\geq1 C ≥ 1 , uniform over p ∈ ( 0 , 1 / 2 ] p\in(0,1/2] p ∈ ( 0 , 1/2 ] , such that, for every n ≥ 1 n\geq1 n ≥ 1 , 0 < λ ≤ 1 0<\lambda\leq1 0 < λ ≤ 1 , and t ≥ 1 t\geq1 t ≥ 1 , with probability at least 1 − e − t 1-e^{-t} 1 − e − t ,
R ( f ^ D , λ ) − R ∗ ≤ C [ a F ( λ ) + λ + λ − 1 ( A n ) 2 / ( 1 + p ) + λ − 1 t 2 n 2 ] . (6) \mathcal{R}(\widehat{f}_{D,\lambda})-\mathcal{R}^{*}\leq C\left[a_{\mathcal{F}}(\lambda)+\lambda+\lambda^{-1}\left(\frac{A}{n}\right)^{2/(1+p)}+\lambda^{-1}\frac{t^{2}}{n^{2}}\right]. \tag*{(6)} R ( f D , λ ) − R ∗ ≤ C [ a F ( λ ) + λ + λ − 1 ( n A ) 2/ ( 1 + p ) + λ − 1 n 2 t 2 ] . (6)
The entropy assumption in Theorem 3.2 requires control of every entropy number. In our manifold application, the available covering estimate holds only above a positive radius, as stated in Proposition 4.6 . The following corollary considers this restriction by requiring the covering estimate only for b ≤ δ ≤ 1 b \le\delta\le1 b ≤ δ ≤ 1 .
Corollary 3.3 (Oracle inequality under a truncated entropy condition). Suppose that the assumptions of Theorem 3.2 hold with the entropy condition (5) replaced by the following covering number condition: for some p ∈ ( 0 , 1 / 2 ] p \in(0,1/2] p ∈ ( 0 , 1/2 ] , A ≥ 1 A \ge1 A ≥ 1 , and b ∈ ( 0 , 1 / 4 ] b \in(0,1/4] b ∈ ( 0 , 1/4 ] ,
log N ( δ , B F , ∥ ⋅ ∥ ∞ ) ≤ A δ − 2 p , b ≤ δ ≤ 1. (7) \log N(\delta, B_F, \lVert\cdot\rVert_{\infty}) \le A\delta^{-2p}, \qquad b \le\delta\le1. \tag*{(7)} log N ( δ , B F , ∥ ⋅ ∥ ∞ ) ≤ A δ − 2 p , b ≤ δ ≤ 1. (7) Then there is a numerical constant C ≥ 1 C \ge1 C ≥ 1 , independent of p p p , A A A , b b b , such that, for every n ≥ 1 n \ge1 n ≥ 1 , 0 < λ ≤ 1 0 < \lambda\le1 0 < λ ≤ 1 , and t ≥ 1 t \ge1 t ≥ 1 , with probability at least 1 − e − t 1-e^{-t} 1 − e − t ,
R ( f ^ D , λ ) − R ∗ ≤ C [ a F ( λ ) + λ + λ − 1 ( A n ) 2 / ( 1 + p ) + b 2 λ + λ − 1 t 2 n 2 ] . (8) \mathcal{R}(\widehat{f}_{D,\lambda})-\mathcal{R}^{*} \le C\left[a_F(\lambda)+\lambda+\lambda^{-1}\left(\frac{A}{n}\right)^{2/(1+p)}+\frac{b^2}{\lambda}+\lambda^{-1}\frac{t^2}{n^2}\right]. \tag*{(8)} R ( f D , λ ) − R ∗ ≤ C [ a F ( λ ) + λ + λ − 1 ( n A ) 2/ ( 1 + p ) + λ b 2 + λ − 1 n 2 t 2 ] . (8)
Remark 3.1 (Regularization dependence). The factor λ − 1 \lambda^{-1} λ − 1 in (6) comes from the dependence of both the localized envelope and the variance bound on the predictor norm. In the proof of Theorem 3.2 , functions in a localized class at level ϵ \epsilon ϵ have supremum norm bounded by C ( ϵ / λ ) 1 / 2 C(\epsilon/\lambda)^{1/2} C ( ϵ / λ ) 1/2 . Lemma 3.1 then gives the second-moment bound C ϵ 3 / 2 λ − 1 / 2 C\epsilon^{3/2}\lambda^{-1/2} C ϵ 3/2 λ − 1/2 . Together, these estimates yield the entropy contribution C λ − 1 ( A / n ) 2 / ( 1 + p ) C\lambda^{-1}(A/n)^{2/(1+p)} C λ − 1 ( A / n ) 2/ ( 1 + p ) through the local Rademacher bound.
In the squared-loss settings of [10 , 13 ] , fixed-level clipping gives a bounded loss class and a Bernstein constant independent of λ \lambda λ . Hence, for its localized class, the variance bound is of order ϵ \epsilon ϵ and gives the entropy contribution C λ − p A / n C\lambda^{-p}A/n C λ − p A / n .
This difference in the dependence on λ \lambda λ explains why approximation by ambient-space extension can yield an intrinsic-dimensional polynomial learning rate in bounded squared-loss regression, whereas the same norm bounds do not suffice in our logistic analysis. We will discuss this distinction further in Remark 4.2.
Gaussian Approximation and Learning Rates on Manifolds
We now combine the oracle inequality with Gaussian approximation bounds to derive learning rates on the manifold X X X . Throughout this section and the proofs of its results, we assume that X X X satisfies the manifold assumptions specified in § 2.2 . We construct Gaussian RKHS approximants that fully exploit the Besov smoothness α \alpha α while retaining RKHS norm bounds governed by the intrinsic dimension ϱ \varrho ϱ . The construction uses a signed combination of Gaussian integral operators at different bandwidths. Together with covering bounds for the Gaussian RKHS, we establish fast learning rates for logistic classification, and faster rates for clipped predictions.
Gaussian Approximation on Besov Spaces
Define the Gaussian integral operator
T γ f ( x ) : = ( 1 γ π ) ϱ ∫ X exp ( − ∥ x − y ∥ R d 2 γ 2 ) f ( y ) d V ( y ) . T_{\gamma}f(x):=\left(\frac{1}{\gamma\sqrt{\pi}}\right)^{\varrho}\int_X \exp\left(-\frac{\lVert x-y\rVert_{\mathbb{R}^{d}}^2}{\gamma^2}\right)f(y)\,dV(y). T γ f ( x ) := ( γ π 1 ) ϱ ∫ X exp ( − γ 2 ∥ x − y ∥ R d 2 ) f ( y ) d V ( y ) .
Results in [21 Lemma 4.1 and Proposition 4.2] use a high-order correction scheme to establish a one-bandwidth expansion in the C k , β C^{k,\beta} C k , β setting. Their construction recursively corrects the input function. The signed multiscale construction below instead cancels the common expansion coefficients across the bandwidths j γ j\gamma j γ .
Given α > 0 \alpha> 0 α > 0 , set
s : = ⌈ α ⌉ + 1 , c j : = ( − 1 ) j + 1 ( s j ) , 1 ≤ j ≤ s . s:=\lceil\alpha\rceil+1,\qquad c_j:=(-1)^{j+1}\binom{s}{j},\qquad1\le j\le s. s := ⌈ α ⌉ + 1 , c j := ( − 1 ) j + 1 ( j s ) , 1 ≤ j ≤ s .
and define the multiscale approximation operator
S γ , s f : = ∑ j = 1 s c j T j γ f . (9) S_{\gamma,s}f:=\sum_{j=1}^{s}c_jT_{j\gamma}f. \tag*{(9)} S γ , s f := j = 1 ∑ s c j T j γ f . (9)
The coefficients satisfy the finite-difference identities
∑ j = 1 s c j = 1 , ∑ j = 1 s c j j ℓ = 0 , 1 ≤ ℓ < s . (10) \sum_{j=1}^{s} c_j = 1,\qquad\sum_{j=1}^{s} c_j j^\ell= 0,\qquad1 \le\ell< s. \tag*{(10)} j = 1 ∑ s c j = 1 , j = 1 ∑ s c j j ℓ = 0 , 1 ≤ ℓ < s . (10)
The following theorem bounds the RKHS norm of S γ , s f S_{\gamma,s}f S γ , s f .
Theorem 4.1 (RKHS norm bound). For every f ∈ L 2 ( X ) f \in L_2(X) f ∈ L 2 ( X ) and γ > 0 \gamma> 0 γ > 0 ,
∥ S γ , s f ∥ H γ ( X ) ≤ C α , X , d γ − ϱ / 2 ∥ f ∥ L 2 ( X ) . \left\|S_{\gamma,s}f\right\|_{H_\gamma(X)} \le C_{\alpha,X,d}\gamma^{-\varrho/2}\left\|f\right\|_{L_2(X)}. ∥ S γ , s f ∥ H γ ( X ) ≤ C α , X , d γ − ϱ /2 ∥ f ∥ L 2 ( X ) .
We next establish an operator expansion in which all operator coefficients are independent of the bandwidth. This result is a Sobolev counterpart of [21 Lemma 4.1] . They establish a uniform expansion for functions in C k , β ( X ) C^{k,\beta}(X) C k , β ( X ) with the remainder controlled in the supremum norm. Here, we establish an L 2 ( X ) L_2(X) L 2 ( X ) remainder bound of order γ 2 N + 2 \gamma^{2N+2} γ 2 N + 2 controlled by ∥ f ∥ H 2 N + 2 ( X ) \left\|f\right\|_{H^{2N+2}(X)} ∥ f ∥ H 2 N + 2 ( X ) , and identify the differential operator coefficients. This Sobolev formulation enables the interpolation argument used below to obtain approximation bounds on Besov spaces.
Lemma 4.2 (Uniform Sobolev expansion). Let Diff m ( X ) \operatorname{Diff}^{m}(X) Diff m ( X ) denote the space of linear differential operators on X X X of order at most m m m with smooth coefficients. Then, there exist differential operators A r A_r A r , r ≥ 0 r \ge0 r ≥ 0 , depending only on X X X , such that A 0 = I A_0 = I A 0 = I , A r ∈ Diff 2 r ( X ) A_r \in\operatorname{Diff}^{2r}(X) A r ∈ Diff 2 r ( X ) , and, for every N ∈ N 0 N \in\mathbb{N}_0 N ∈ N 0 , there exists a constant C N , X , d > 0 C_{N,X,d} > 0 C N , X , d > 0 satisfying
∥ T γ f − ∑ r = 0 N γ 2 r A r f ∥ L 2 ( X ) ≤ C N , X , d γ 2 N + 2 ∥ f ∥ H 2 N + 2 ( X ) (11) \left\|T_\gamma f-\sum_{r=0}^{N}\gamma^{2r}A_r f\right\|_{L_2(X)}
\le C_{N,X,d}\gamma^{2N+2}\left\|f\right\|_{H^{2N+2}(X)} \tag*{(11)} T γ f − r = 0 ∑ N γ 2 r A r f L 2 ( X ) ≤ C N , X , d γ 2 N + 2 ∥ f ∥ H 2 N + 2 ( X ) (11) for all 0 < γ ≤ 1 0 < \gamma\le1 0 < γ ≤ 1 and f ∈ H 2 N + 2 ( X ) f \in H^{2N+2}(X) f ∈ H 2 N + 2 ( X ) . Moreover,
A r = 1 4 r r ! Δ X r + B r , B r ∈ Diff 2 r − 2 ( X ) , r ≥ 1 , A_r=\frac{1}{4^r r!}\Delta_X^r+B_r,\qquad B_r\in\operatorname{Diff}^{2r-2}(X),\qquad r\ge1, A r = 4 r r ! 1 Δ X r + B r , B r ∈ Diff 2 r − 2 ( X ) , r ≥ 1 , where Δ X = div ∇ \Delta_X=\operatorname{div}\nabla Δ X = div ∇ .
The identities (10) cancel the lower-order terms in this expansion when the operators are combined in S γ , s S_{\gamma,s} S γ , s . The resulting Sobolev approximation bound and uniform L 2 ( X ) L_2(X) L 2 ( X ) boundedness yield the following Besov approximation estimate by real interpolation.
Theorem 4.3 (Besov approximation). Let α > 0 \alpha> 0 α > 0 and let S γ , s S_{\gamma,s} S γ , s be defined by (9) , with s = ⌊ α ⌋ + 1 s=\lfloor\alpha\rfloor+1 s = ⌊ α ⌋ + 1 . Then there is a constant C α > 0 C_\alpha>0 C α > 0 , independent of γ \gamma γ and f f f , such that, for every f ∈ B 2 , ∞ α ( X ) f\in B^\alpha_{2,\infty}(X) f ∈ B 2 , ∞ α ( X ) ,
∥ S γ , s f − f ∥ L 2 ( X ) ≤ C α γ α ∥ f ∥ B 2 , ∞ α ( X ) , 0 < γ ≤ 1. (12) \left\|S_{\gamma,s}f-f\right\|_{L_2(X)}
\le C_\alpha\gamma^\alpha\left\|f\right\|_{B^\alpha_{2,\infty}(X)},\qquad0<\gamma\le1. \tag*{(12)} ∥ S γ , s f − f ∥ L 2 ( X ) ≤ C α γ α ∥ f ∥ B 2 , ∞ α ( X ) , 0 < γ ≤ 1. (12)
Recall the approximation error defined in (4) . For the bandwidth-dependent Gaussian RKHS H γ ( X ) H_\gamma(X) H γ ( X ) , we write
a γ ( λ ) : = inf f ∈ H γ ( X ) { λ ∥ f ∥ H γ ( X ) 2 + R ( f ) − R ∗ } . (13) a_\gamma(\lambda):=\inf_{f\in H_\gamma(X)}
\left\{\lambda\left\|f\right\|_{H_\gamma(X)}^2+\mathcal{R}(f)-\mathcal{R}^*\right\}. \tag*{(13)} a γ ( λ ) := f ∈ H γ ( X ) inf { λ ∥ f ∥ H γ ( X ) 2 + R ( f ) − R ∗ } . (13)
Theorem 4.4 (Approximation error). Suppose that 0 < η < 1 0<\eta<1 0 < η < 1 holds P X P_X P X -almost surely and the Bayes decision function
f ∗ = log η 1 − η f^*=\log\frac{\eta}{1-\eta} f ∗ = log 1 − η η has a measurable extension in B 2 , ∞ α ( X ) B^\alpha_{2,\infty}(X) B 2 , ∞ α ( X ) for some α > 0 \alpha>0 α > 0 . The extension is again denoted by f ∗ f^* f ∗ . Suppose also that P X P_X P X has a density g ∈ L p ( X ) g\in L_p(X) g ∈ L p ( X ) with respect to d V dV d V for a fixed 2 ≤ p ≤ ∞ 2\le p\le\infty 2 ≤ p ≤ ∞ . Then, for every λ > 0 \lambda>0 λ > 0 and 0 < γ ≤ 1 0<\gamma\le1 0 < γ ≤ 1 ,
a γ ( λ ) ≤ C α , p , X , d γ 2 α ( 1 − 1 / p ) ∥ g ∥ L p ( X ) ∥ f ∗ ∥ B 2 , ∞ α ( X ) 2 ( 1 − 1 / p ) + C α , X , d λ γ − ϱ ∥ f ∗ ∥ L 2 ( X ) 2 . (14) a_\gamma(\lambda)\le C_{\alpha,p,X,d}\gamma^{2\alpha(1-1/p)}
\lVert g\rVert_{L_p(X)}\lVert f^*\rVert_{B^\alpha_{2,\infty}(X)}^{2(1-1/p)}
+C_{\alpha,X,d}\lambda\gamma^{-\varrho}\lVert f^*\rVert_{L_2(X)}^2. \tag*{(14)} a γ ( λ ) ≤ C α , p , X , d γ 2 α ( 1 − 1/ p ) ∥ g ∥ L p ( X ) ∥ f ∗ ∥ B 2 , ∞ α ( X ) 2 ( 1 − 1/ p ) + C α , X , d λ γ − ϱ ∥ f ∗ ∥ L 2 ( X ) 2 . (14)
The Besov assumption in Theorem 4.4 can be verified through regularity and endpoint conditions on η \eta η . The small value bound in [6 Definition 3.1] controls the probability that a conditional class probability is close to zero. For binary classification, it controls both η \eta η near zero and η \eta η near one. This condition is combined with low-order Hölder regularity in [25 Assumption 3.1] for kernel logistic regression. Unlike the usual Tsybakov noise condition, which concerns η \eta η near the decision level 1 / 2 1/2 1/2 , it controls the regions where f ∗ f^* f ∗ can become large. Motivated by these works, we give two sufficient conditions for the Besov regularity of f ∗ f^* f ∗ , covering low orders under Hölder assumptions and arbitrary positive orders under stronger smoothness and endpoint conditions. We impose the small value condition with respect to d V dV d V , since our approximation space is defined relative to manifold volume.
Proposition 4.5 (Besov regularity of the Bayes function). Let α > 0 \alpha> 0 α > 0 , and suppose that the conditional probability admits a version η : X → [ 0 , 1 ] \eta: X \to[0,1] η : X → [ 0 , 1 ] . Each of the following conditions ensures that 0 < η < 1 0 < \eta< 1 0 < η < 1 holds d V dV d V -almost everywhere and that the Bayes function f ∗ : = log ( η / ( 1 − η ) ) f^* := \log(\eta/(1-\eta)) f ∗ := log ( η / ( 1 − η )) , extended arbitrarily on the remaining null set, belongs to B 2 , ∞ α ( X ) B^\alpha_{2,\infty}(X) B 2 , ∞ α ( X ) .
(i) Let 0 < α < 1 0 < \alpha< 1 0 < α < 1 and α ≤ β ≤ 1 \alpha\le\beta\le1 α ≤ β ≤ 1 . If η ∈ C 0 , β ( X ) \eta\in C^{0,\beta}(X) η ∈ C 0 , β ( X ) and
∫ X [ η ( x ) ( 1 − η ( x ) ) ] − 2 α / β d V ( x ) < ∞ , (15) \int_X [\eta(x)(1-\eta(x))]^{-2\alpha/\beta}\,dV(x) < \infty, \tag*{(15)} ∫ X [ η ( x ) ( 1 − η ( x )) ] − 2 α / β d V ( x ) < ∞ , (15) then f ∗ ∈ B 2 , ∞ α ( X ) f^* \in B^\alpha_{2,\infty}(X) f ∗ ∈ B 2 , ∞ α ( X ) and
∥ f ∗ ∥ B 2 , ∞ α ( X ) ≤ C X , α , β ( 1 + ∥ η ∥ C 0 , β ( X ) α / β ) ( ∫ X [ η ( 1 − η ) ] − 2 α / β d V ) 1 / 2 . (16) \|f^*\|_{B^\alpha_{2,\infty}(X)} \le C_{X,\alpha,\beta}(1+\|\eta\|_{C^{0,\beta}(X)}^{\alpha/\beta})\left(\int_X [\eta(1-\eta)]^{-2\alpha/\beta}\,dV\right)^{1/2}. \tag*{(16)} ∥ f ∗ ∥ B 2 , ∞ α ( X ) ≤ C X , α , β ( 1 + ∥ η ∥ C 0 , β ( X ) α / β ) ( ∫ X [ η ( 1 − η ) ] − 2 α / β d V ) 1/2 . (16) (ii) Let m ≥ 2 m \ge2 m ≥ 2 be an integer with α ≤ m \alpha\le m α ≤ m . If η ∈ C m ( X ) \eta\in C^m(X) η ∈ C m ( X ) and
∫ X [ η ( x ) ( 1 − η ( x ) ) ] − m d V ( x ) < ∞ , (17) \int_X [\eta(x)(1-\eta(x))]^{-m}\,dV(x) < \infty, \tag*{(17)} ∫ X [ η ( x ) ( 1 − η ( x )) ] − m d V ( x ) < ∞ , (17) then f ∗ ∈ H m ( X ) f^* \in H^m(X) f ∗ ∈ H m ( X ) and
∥ f ∗ ∥ H m ( X ) ≤ C ( ∫ X [ η ( 1 − η ) ] − m d V ) 1 / 2 , (18) \|f^*\|_{H^m(X)} \le C\left(\int_X [\eta(1-\eta)]^{-m}\,dV\right)^{1/2}, \tag*{(18)} ∥ f ∗ ∥ H m ( X ) ≤ C ( ∫ X [ η ( 1 − η ) ] − m d V ) 1/2 , (18) where C C C depends only on X X X , m m m , and ∥ η ∥ C m ( X ) \|\eta\|_{C^m(X)} ∥ η ∥ C m ( X ) .
In particular, suppose that, for some C 0 > 0 C_0 > 0 C 0 > 0 and κ > 0 \kappa> 0 κ > 0 ,
V { x ∈ X : min { η ( x ) , 1 − η ( x ) } ≤ u } ≤ C 0 u κ , 0 < u ≤ 1 2 . (19) V\{x \in X : \min\{\eta(x),1-\eta(x)\} \le u\} \le C_0u^\kappa,\qquad0 < u \le\frac{1}{2}. \tag*{(19)} V { x ∈ X : min { η ( x ) , 1 − η ( x )} ≤ u } ≤ C 0 u κ , 0 < u ≤ 2 1 . (19) Then (15) follows when κ > 2 α / β \kappa> 2\alpha/\beta κ > 2 α / β , and (17) follows when κ > m \kappa> m κ > m .
It is worth noting that the small value condition in [6 Definition 3.1] is imposed under P X P_X P X . This does not by itself imply (19) . If one additionally assumes that the density g ≥ c > 0 g \ge c > 0 g ≥ c > 0 almost everywhere, the corresponding bound under P X P_X P X implies (19) .
Learning Rates
We now combine the Gaussian approximation bounds with the oracle inequalities to derive learning rates on the manifold. The remaining ingredient is a covering estimate for the Gaussian RKHS whose dependence on the bandwidth is governed by the intrinsic dimension.
Proposition 4.6 (Gaussian covering numbers above a fixed scale). Fix an integer m ≥ 2 m \ge2 m ≥ 2 . There are constants C m ≥ 1 C_m \ge1 C m ≥ 1 and γ m ∈ ( 0 , 1 ) \gamma_m \in(0,1) γ m ∈ ( 0 , 1 ) , depending only on m m m and X X X , such that γ m m − 1 < 1 / 4 \gamma_m^{m-1} < 1/4 γ m m − 1 < 1/4 and, for every 0 < γ ≤ γ m 0 < \gamma\le\gamma_m 0 < γ ≤ γ m ,
log N ( δ , B H γ ( X ) , ∥ ⋅ ∥ ∞ ) ≤ C m γ − ϱ ( log 2 δ ) ϱ + 1 , γ m − 1 < δ ≤ 1. (20) \log N(\delta,B_{H_\gamma}(X),\|\cdot\|_\infty) \le C_m\gamma^{-\varrho}\left(\log\frac{2}{\delta}\right)^{\varrho+1},\qquad\gamma^{m-1}<\delta\le1. \tag*{(20)} log N ( δ , B H γ ( X ) , ∥ ⋅ ∥ ∞ ) ≤ C m γ − ϱ ( log δ 2 ) ϱ + 1 , γ m − 1 < δ ≤ 1. (20) This implies, uniformly for the entropy exponent q ∈ ( 0 , 1 / 2 ] q \in(0,1/2] q ∈ ( 0 , 1/2 ] ,
log N ( δ , B H γ ( X ) , ∥ ⋅ ∥ ∞ ) ≤ C m q − ( ϱ + 1 ) γ − ϱ δ − 2 q , γ m − 1 < δ ≤ 1. (21) \log N(\delta,B_{H_\gamma}(X),\|\cdot\|_\infty) \le C_mq^{-(\varrho+1)}\gamma^{-\varrho}\delta^{-2q},\qquad\gamma^{m-1}<\delta\le1. \tag*{(21)} log N ( δ , B H γ ( X ) , ∥ ⋅ ∥ ∞ ) ≤ C m q − ( ϱ + 1 ) γ − ϱ δ − 2 q , γ m − 1 < δ ≤ 1. (21)
The proof of Proposition 4.6 is based on [28 Lemma 5.4] at the fixed smoothness order m m m . The manifold improvement in [10 Theorem A.10] also uses [28 Lemma 5.4] . Here we retain the strict lower limit δ > γ m − 1 \delta> \gamma^{m-1} δ > γ m − 1 in that covering estimate and apply Corollary 3.3 at a covering scale above this limit. The resulting additional term is of order λ − 1 γ 2 ( m − 1 ) \lambda^{-1}\gamma^{2(m-1)} λ − 1 γ 2 ( m − 1 ) . By choosing a sufficiently large fixed order m m m , we can absorb this term and obtain the following learning rate.
Theorem 4.7 (Learning rate). Under the assumptions of Theorem 4.4 , write f ^ n , λ , γ \widehat{f}_{n,\lambda,\gamma} f n , λ , γ for the regularized empirical minimizer over H γ ( X ) \mathcal{H}_{\gamma}(X) H γ ( X ) . Set β : = 2 α ( 1 − 1 / p ) \beta:= 2\alpha(1-1/p) β := 2 α ( 1 − 1/ p ) , and define
γ n : = ( ( log n ) ϱ + 1 n ) 2 2 β + 3 ϱ , λ n : = ( ( log n ) ϱ + 1 n ) 2 ( β + ϱ ) 2 β + 3 ϱ . \gamma_n := \left(\frac{(\log n)^{\varrho+1}}{n}\right)^{\frac{2}{2\beta+3\varrho}}, \qquad\lambda_n := \left(\frac{(\log n)^{\varrho+1}}{n}\right)^{\frac{2(\beta+\varrho)}{2\beta+3\varrho}}. γ n := ( n ( log n ) ϱ + 1 ) 2 β + 3 ϱ 2 , λ n := ( n ( log n ) ϱ + 1 ) 2 β + 3 ϱ 2 ( β + ϱ ) . For every fixed t ≥ 1 t \ge1 t ≥ 1 and all sufficiently large n n n , with probability at least 1 − e − t 1-e^{-t} 1 − e − t ,
R ( f ^ n , λ n , γ n ) − R ∗ ≤ C ( ( log n ) ϱ + 1 n ) 4 α ( 1 − 1 / p ) 4 α ( 1 − 1 / p ) + 3 ϱ . (22) \mathcal{R}(\widehat{f}_{n,\lambda_n,\gamma_n})-\mathcal{R}^{*}\le C\left(\frac{(\log n)^{\varrho+1}}{n}\right)^{\frac{4\alpha(1-1/p)}{4\alpha(1-1/p)+3\varrho}}. \tag*{(22)} R ( f n , λ n , γ n ) − R ∗ ≤ C ( n ( log n ) ϱ + 1 ) 4 α ( 1 − 1/ p ) + 3 ϱ 4 α ( 1 − 1/ p ) . (22) Here C C C may depend on t , α , p , X , ∥ g ∥ L p ( X ) t,\alpha,p,X,\lVert g\rVert_{L_p(X)} t , α , p , X , ∥ g ∥ L p ( X ) , and ∥ f ∗ ∥ B 2 , ∞ α ( X ) \lVert f^{*}\rVert_{B^{\alpha}_{2,\infty}(X)} ∥ f ∗ ∥ B 2 , ∞ α ( X ) , but is independent of n n n .
Remark 4.1 (From logistic risk to misclassification risk). Recall the calibration inequality (2) . This calibration inequality can be improved when the distribution P P P satisfies the Tsybakov noise condition [24 ] . Specifically, assume that there exist 0 ≤ q < ∞ 0\le q<\infty 0 ≤ q < ∞ and C q > 0 C_q>0 C q > 0 such that
P X { x ∈ X : ∣ 2 η ( x ) − 1 ∣ ≤ u } ≤ C q u q , u > 0. P_X\{x\in X:|2\eta(x)-1|\le u\}\le C_q u^q,\qquad u>0. P X { x ∈ X : ∣2 η ( x ) − 1∣ ≤ u } ≤ C q u q , u > 0. For q = ∞ q=\infty q = ∞ , we instead assume that there exists u 0 > 0 u_0>0 u 0 > 0 such that P X { x ∈ X : ∣ 2 η ( x ) − 1 ∣ ≤ u 0 } = 0 P_X\{x\in X:|2\eta(x)-1|\le u_0\}=0 P X { x ∈ X : ∣2 η ( x ) − 1∣ ≤ u 0 } = 0 . Then, by [26 (2.4)] , the excess misclassification risk satisfies
R 01 ( f ) − R 01 ∗ ≤ C ( R ( f ) − R ∗ ) ( q + 1 ) / ( q + 2 ) (23) \mathcal{R}_{01}(f)-\mathcal{R}_{01}^{*}\le C\left(\mathcal{R}(f)-\mathcal{R}^{*}\right)^{(q+1)/(q+2)} \tag*{(23)} R 01 ( f ) − R 01 ∗ ≤ C ( R ( f ) − R ∗ ) ( q + 1 ) / ( q + 2 ) (23) for every measurable f : X → R f:X\to\mathbb{R} f : X → R . For finite q q q , the constant C C C depends only on q q q and C q C_q C q . For q = ∞ q=\infty q = ∞ , the exponent ( q + 1 ) / ( q + 2 ) (q+1)/(q+2) ( q + 1 ) / ( q + 2 ) is interpreted as 1 1 1 , and C C C depends on u 0 u_0 u 0 .
Under the assumptions and with the parameter choices of Theorem 4.7 , if the above noise condition holds, then, for every fixed t ≥ 1 t\ge1 t ≥ 1 and all sufficiently large n n n , the calibration inequality (23) gives, with probability at least 1 − e − t 1-e^{-t} 1 − e − t ,
R 01 ( f ^ n , λ n , γ n ) − R 01 ∗ ≤ C ( ( log n ) ϱ + 1 n ) 4 α ( 1 − 1 / p ) ( q + 1 ) ( 4 α ( 1 − 1 / p ) + 3 ϱ ) ( q + 2 ) . \mathcal{R}_{01}(\widehat{f}_{n,\lambda_n,\gamma_n})-\mathcal{R}_{01}^{*}\le C\left(\frac{(\log n)^{\varrho+1}}{n}\right)^{\frac{4\alpha(1-1/p)(q+1)}{(4\alpha(1-1/p)+3\varrho)(q+2)}}. R 01 ( f n , λ n , γ n ) − R 01 ∗ ≤ C ( n ( log n ) ϱ + 1 ) ( 4 α ( 1 − 1/ p ) + 3 ϱ ) ( q + 2 ) 4 α ( 1 − 1/ p ) ( q + 1 ) .
Remark 4.2 (Intrinsic norm control and comparison with squared-loss regression). Throughout this comparison, we suppress logarithmic factors, write β : = 2 α ( 1 − 1 / p ) \beta:=2\alpha(1-1/p) β := 2 α ( 1 − 1/ p ) , and use q q q for the entropy exponent. As we have mentioned in Remark 3.1, in the bounded squared-loss settings of [10 , 13 ] , fixed-level clipping gives an envelope and a Bernstein constant independent of λ \lambda λ . At localized risk level ε \varepsilon ε , the second moment is of order ε \varepsilon ε , and the leading Rademacher term has the form
( ε λ ) q / 2 ε ( 1 − q ) / 2 ( A n ) 1 / 2 = ε 1 / 2 λ − q / 2 ( A n ) 1 / 2 . \left(\frac{\varepsilon}{\lambda}\right)^{q/2}\varepsilon^{(1-q)/2}\left(\frac{A}{n}\right)^{1/2}
=\varepsilon^{1/2}\lambda^{-q/2}\left(\frac{A}{n}\right)^{1/2}. ( λ ε ) q /2 ε ( 1 − q ) /2 ( n A ) 1/2 = ε 1/2 λ − q /2 ( n A ) 1/2 . For logistic loss, the second-moment bound instead gives
( ε λ ) q / 2 ( ε 3 / 2 λ − 1 / 2 ) ( 1 − q ) / 2 ( A n ) 1 / 2 = ε ( 3 − q ) / 4 λ − ( 1 + q ) / 4 ( A n ) 1 / 2 . \left(\frac{\varepsilon}{\lambda}\right)^{q/2}\left(\varepsilon^{3/2}\lambda^{-1/2}\right)^{(1-q)/2}\left(\frac{A}{n}\right)^{1/2}
=\varepsilon^{(3-q)/4}\lambda^{-(1+q)/4}\left(\frac{A}{n}\right)^{1/2}. ( λ ε ) q /2 ( ε 3/2 λ − 1/2 ) ( 1 − q ) /2 ( n A ) 1/2 = ε ( 3 − q ) /4 λ − ( 1 + q ) /4 ( n A ) 1/2 . Requiring each term to be bounded by a sufficiently small ε \varepsilon ε gives ε ≥ C λ − q A / n \varepsilon\ge C\lambda^{-q}A/n ε ≥ C λ − q A / n for squared loss, and ε ≥ C λ − 1 ( A / n ) 2 / ( 1 + q ) \varepsilon\ge C\lambda^{-1}(A/n)^{2/(1+q)} ε ≥ C λ − 1 ( A / n ) 2/ ( 1 + q ) for logistic loss.
This distinction matters when the squared RKHS norm of an ambient-space Gaussian approximant is bounded by C γ − d C\gamma^{-d} C γ − d , where γ \gamma γ is the Gaussian bandwidth. For any fixed κ > 0 \kappa>0 κ > 0 and sufficiently large n n n , consider
λ n = n − κ , q n = 1 log n , \lambda_n=n^{-\kappa},\qquad q_n=\frac{1}{\log n}, λ n = n − κ , q n = log n 1 , λ n − q n = e κ , λ n − 1 = n κ . \lambda_n^{-q_n}=e^{\kappa},\qquad\lambda_n^{-1}=n^{\kappa}. λ n − q n = e κ , λ n − 1 = n κ . In the squared-loss bounds of [10 , 13 ] , the remaining dependence on q n q_n q n contributes only logarithmic factors. One can therefore choose κ \kappa κ large enough to control the norm penalty λ γ − d \lambda\gamma^{-d} λ γ − d of an ambient-space approximant without introducing d d d into the polynomial exponent of the stochastic term. With intrinsic-dimensional covering
bounds, balancing the squared-loss risk bias γ 2 α \gamma^{2\alpha} γ 2 α against the stochastic term γ − ϱ / n \gamma^{-\varrho}/n γ − ϱ / n gives the rate n − 2 α / ( 2 α + ϱ ) n^{-2\alpha/(2\alpha+\varrho)} n − 2 α / ( 2 α + ϱ ) , whose exponent matches the minimax rate in [18 ] .
However, for logistic loss, Theorem 4.4 gives a risk bias of order γ β \gamma^\beta γ β . The Gaussian covering estimate makes the stochastic contribution of polynomial order λ − 1 γ − 2 ϱ / n 2 \lambda^{-1}\gamma^{-2\varrho}/n^2 λ − 1 γ − 2 ϱ / n 2 . With the same risk bias, an ambient-space squared RKHS norm bound of order γ − d \gamma^{-d} γ − d and the intrinsic bound of Theorem 4.1 lead to the respective balances
inf λ > 0 { λ γ − d + λ − 1 γ − 2 ϱ n 2 } = 2 γ − ( d + 2 ϱ ) / 2 n , inf λ > 0 { λ γ − ϱ + λ − 1 γ − 2 ϱ n 2 } = 2 γ − 3 ϱ / 2 n . \inf_{\lambda>0}\left\{\lambda\gamma^{-d}+\lambda^{-1}\frac{\gamma^{-2\varrho}}{n^2}\right\}=\frac{2\gamma^{-(d+2\varrho)/2}}{n},
\qquad
\inf_{\lambda>0}\left\{\lambda\gamma^{-\varrho}+\lambda^{-1}\frac{\gamma^{-2\varrho}}{n^2}\right\}=\frac{2\gamma^{-3\varrho/2}}{n}. λ > 0 inf { λ γ − d + λ − 1 n 2 γ − 2 ϱ } = n 2 γ − ( d + 2 ϱ ) /2 , λ > 0 inf { λ γ − ϱ + λ − 1 n 2 γ − 2 ϱ } = n 2 γ − 3 ϱ /2 .
Balancing each term with γ β \gamma^\beta γ β gives the ambient rate n − 2 β / ( 2 β + 2 ϱ + d ) n^{-2\beta/(2\beta+2\varrho+d)} n − 2 β / ( 2 β + 2 ϱ + d ) and the intrinsic rate n − 2 β / ( 2 β + 3 ϱ ) n^{-2\beta/(2\beta+3\varrho)} n − 2 β / ( 2 β + 3 ϱ ) in Theorem 4.7 . Thus the intrinsic norm bound removes the ambient dimension from the polynomial exponent. This comparison explains the role of the intrinsic RKHS norm bound within the present oracle analysis.
The resulting rate remains slower than the squared-loss regression rate n − 2 α / ( 2 α + ϱ ) n^{-2\alpha/(2\alpha+\varrho)} n − 2 α / ( 2 α + ϱ ) in [10 ] , whose dependence on smoothness and dimension matches the minimax rate in [18 ] . Two differences in the analysis explain this gap. The first is, as discussed above, the norm-dependent variance bound in Lemma 3.1 . The second is the density-dependent risk bias in Theorem 4.4 , which is of order γ 2 α ( 1 − 1 / p ) \gamma^{2\alpha(1-1/p)} γ 2 α ( 1 − 1/ p ) for g ∈ L p ( X ) g\in L_p(X) g ∈ L p ( X ) .
In view of Remark 4.2, we now seek to close the gap between the logistic learning rate and the squared-loss regression rate by clipping the predictions at a level that increases with the sample size. The idea of clipping is to control the supremum norm and the Bernstein constant through the clipping level M M M . As we will show in the proof of Theorem 4.8 , clipping only increases the logistic risk by at most e − M e^{-M} e − M . Therefore, a clipping level that grows logarithmically with n n n can make the clipping bias negligible. The following theorem gives a faster rate for the clipped predictions under the same Besov and density assumptions, with a new choice of bandwidth and regularization parameter. For bounded densities, it recovers the optimal polynomial exponent of the regression rate n − 2 α / ( 2 α + ϱ ) n^{-2\alpha/(2\alpha+\varrho)} n − 2 α / ( 2 α + ϱ ) . Since clipping preserves the classification rule, it also improves the misclassification bound for the original estimator.
Theorem 4.8 (Learning rates for clipped predictions). Under the assumptions of Theorem 4.4 , define the clipping operator
π M ( u ) = max { − M , min { u , M } } , M > 0. \pi_M(u)=\max\{-M,\min\{u,M\}\},\qquad M>0. π M ( u ) = max { − M , min { u , M }} , M > 0. Set β : = 2 α ( 1 − 1 / p ) \beta:=2\alpha(1-1/p) β := 2 α ( 1 − 1/ p ) , and define
γ n : = ( ( log n ) ϱ + 2 n ) 1 / ( β + ϱ ) , λ n : = ( log n ) ϱ + 2 n , M n : = 3 log n . \gamma_n:=\left(\frac{(\log n)^{\varrho+2}}{n}\right)^{1/(\beta+\varrho)},
\qquad
\lambda_n:=\frac{(\log n)^{\varrho+2}}{n},
\qquad
M_n:=3\log n. γ n := ( n ( log n ) ϱ + 2 ) 1/ ( β + ϱ ) , λ n := n ( log n ) ϱ + 2 , M n := 3 log n . For every fixed t ≥ 1 t\geq1 t ≥ 1 and all sufficiently large n n n , with probability at least 1 − e − t 1-e^{-t} 1 − e − t ,
R ( π M n f ^ n , λ n , γ n ) − R ∗ ≤ C ( ( log n ) ϱ + 2 n ) β / ( β + ϱ ) , (24) \mathcal{R}(\pi_{M_n}\widehat f_{n,\lambda_n,\gamma_n})-\mathcal{R}^*
\leq C\left(\frac{(\log n)^{\varrho+2}}{n}\right)^{\beta/(\beta+\varrho)}, \tag*{(24)} R ( π M n f n , λ n , γ n ) − R ∗ ≤ C ( n ( log n ) ϱ + 2 ) β / ( β + ϱ ) , (24) R 01 ( f ^ n , λ n , γ n ) − R 01 ∗ ≤ C ( ( log n ) ϱ + 2 n ) β / ( 2 ( β + ϱ ) ) . \mathcal{R}_{01}(\widehat f_{n,\lambda_n,\gamma_n})-\mathcal{R}_{01}^*
\leq C\left(\frac{(\log n)^{\varrho+2}}{n}\right)^{\beta/(2(\beta+\varrho))}. R 01 ( f n , λ n , γ n ) − R 01 ∗ ≤ C ( n ( log n ) ϱ + 2 ) β / ( 2 ( β + ϱ )) . Here C C C may depend on t , α , p , X , ∥ g ∥ L p ( X ) t,\alpha,p,X,\|g\|_{L_p(X)} t , α , p , X , ∥ g ∥ L p ( X ) , and ∥ f ∗ ∥ B 2 , ∞ α ( X ) \|f^*\|_{B^\alpha_{2,\infty}(X)} ∥ f ∗ ∥ B 2 , ∞ α ( X ) , but is independent of n n n .
For p = ∞ p=\infty p = ∞ , the misclassification bound in (24) has polynomial order n − α / ( 2 α + ϱ ) n^{-\alpha/(2\alpha+\varrho)} n − α / ( 2 α + ϱ ) . In the classical setting of binary classification with an α \alpha α -Hölder conditional class probability on a d d d -dimensional Euclidean domain, under standard density assumptions and without an additional margin condition, the minimax rate of expected excess misclassification risk is n − α / ( 2 α + d ) n^{-\alpha/(2\alpha+d)} n − α / ( 2 α + d ) . See [3 (1.5)] . Thus, up to logarithmic factors, our bound also matches the minimax optimal rate for binary classification, with the ambient dimension d d d replaced by the intrinsic dimension ϱ \varrho ϱ .
Conclusions
We have studied regularized Gaussian logistic classification on a compact smooth embedded manifold. The analysis combines a Bernstein-type bound with linear dependence on the prediction bound, an oracle inequality based on regularized-risk localization, and Gaussian approximation of intrinsic Besov functions. oracle analysis does not require class probabilities to be separated from zero and one. The signed multiscale Gaussian approximation construction cancels the geometric and differential contributions to the smoothing bias and achieves arbitrary positive approximation order while retaining the intrinsic-dimensional RKHS norm bound. A version of the oracle inequality based on covering estimates above a positive radius makes these approximation results compatible with the available manifold entropy bound. Under Besov regularity of the Bayes function f ∗ f^{*} f ∗ and a density g ∈ L p ( X ) g \in L_{p}(X) g ∈ L p ( X ) with 2 ≤ p ≤ ∞ 2 \le p \le\infty 2 ≤ p ≤ ∞ , we prove the logistic excess-risk rate ( ( log n ) ϱ + 1 / n ) 4 α ( 1 − 1 / p ) / ( 4 α ( 1 − 1 / p ) + 3 ϱ ) ((\log n)^{\varrho+1}/n)^{4\alpha(1-1/p)/(4\alpha(1-1/p)+3\varrho)} (( log n ) ϱ + 1 / n ) 4 α ( 1 − 1/ p ) / ( 4 α ( 1 − 1/ p ) + 3 ϱ ) in Theorem 4.7 . Its exponents depend on the intrinsic dimension ϱ \varrho ϱ , regularity α \alpha α , and density integrability p p p , and are independent of the ambient dimension d d d . Theorem 4.8 further improves the learning rate by clipping predictions. Under the same assumptions, the clipped logistic excess risk has rate ( ( log n ) ϱ + 2 / n ) 2 α ( 1 − 1 / p ) / ( 2 α ( 1 − 1 / p ) + ϱ ) ((\log n)^{\varrho+2}/n)^{2\alpha(1-1/p)/(2\alpha(1-1/p)+\varrho)} (( log n ) ϱ + 2 / n ) 2 α ( 1 − 1/ p ) / ( 2 α ( 1 − 1/ p ) + ϱ ) . For bounded densities, this bound recovers the polynomial exponent 2 α / ( 2 α + ϱ ) 2\alpha/(2\alpha+\varrho) 2 α / ( 2 α + ϱ ) of the minimax optimal squared-loss regression rate, without assuming a uniformly bounded Bayes function. The corresponding excess misclassification risk also matches, up to logarithmic factors, the classical minimax optimal rate n − α / ( 2 α + ϱ ) n^{-\alpha/(2\alpha+\varrho)} n − α / ( 2 α + ϱ ) for binary classification without an additional margin condition.
Several questions remain open. First, the bandwidth and regularization choices depend on the smoothness, intrinsic dimension, and density integrability. A data-dependent selection rule that retains the same rate would connect the present analysis more closely with the adaptive procedures in [10 , 21 ] . Second, for 2 ≤ p < ∞ 2 \le p < \infty 2 ≤ p < ∞ , it remains open whether the dependence of the clipped learning rates on density integrability is optimal. Matching upper and lower bounds would determine whether this dependence is intrinsic or reflects a limitation of the present analysis. Third, it would be useful to extend the analysis to other surrogate loss functions, and manifolds with boundary or lower regularity.
Proofs
This appendix collects the proofs of the paper. Throughout the proofs, for a P P P -integrable measurable function g : X × Y → R g : X \times Y \to\mathbb{R} g : X × Y → R , we write
E P g : = ∫ X × Y g ( x , y ) d P ( x , y ) , \mathbb{E}_{P}g := \int_{X \times Y} g(x,y)\,\mathrm{d}P(x,y), E P g := ∫ X × Y g ( x , y ) d P ( x , y ) ,
E D g : = 1 n ∑ i = 1 n g ( X i , Y i ) . \mathbb{E}_{D}g := \frac{1}{n}\sum_{i=1}^{n}g(X_i,Y_i). E D g := n 1 i = 1 ∑ n g ( X i , Y i ) .
Thus E D \mathbb{E}_{D} E D denotes the average over the observed samples. Expectations over the random samples are written as E D ∼ P n \mathbb{E}_{D\sim P^{n}} E D ∼ P n .
Proof of Proposition 2.1
If 0 < β < 1 0 < \beta< 1 0 < β < 1 , the Euclidean Hölder–Zygmund theorem [22 Theorem 1.2] , applied in a finite atlas after multiplication by a subordinate partition of unity whose supports are compactly contained in the corresponding charts, gives C k , β ( X ) = B ∞ , ∞ β ( X ) C^{k,\beta}(X)=B^{\beta}_{\infty,\infty}(X) C k , β ( X ) = B ∞ , ∞ β ( X ) with equivalent norms. If β = 1 \beta=1 β = 1 , the same argument gives the continuous inclusion C k , 1 ( X ) ↪ B ∞ , ∞ k ( X ) C^{k,1}(X)\hookrightarrow B^{k}_{\infty,\infty}(X) C k , 1 ( X ) ↪ B ∞ , ∞ k ( X ) .
It remains to pass from L ∞ L_{\infty} L ∞ to L 2 L_{2} L 2 . Choose 0 < r ≤ inj ( X ) 0<r\le\operatorname{inj}(X) 0 < r ≤ inj ( X ) , ε > 0 \varepsilon>0 ε > 0 sufficiently small, and take the local-means order N ∈ N N\in\mathbb{N} N ∈ N with N > max { α , ϱ + 6 } N>\max\{\alpha,\varrho+6\} N > max { α , ϱ + 6 } . Then [22 Theorem 7.3.2] applies, with the same local means k 0 ( ε , f ) k_{0}(\varepsilon,f) k 0 ( ε , f ) and k ( t , f ) k(t,f) k ( t , f ) , both for ( p , q ) = ( ∞ , ∞ ) (p,q)=(\infty,\infty) ( p , q ) = ( ∞ , ∞ ) and for ( p , q ) = ( 2 , ∞ ) (p,q)=(2,\infty) ( p , q ) = ( 2 , ∞ ) . Compactness gives
∥ u ∥ L 2 ( X ) ≤ Vol ( X ) 1 / 2 ∥ u ∥ L ∞ ( X ) . \|u\|_{L_{2}(X)}\le\operatorname{Vol}(X)^{1/2}\|u\|_{L_{\infty}(X)}. ∥ u ∥ L 2 ( X ) ≤ Vol ( X ) 1/2 ∥ u ∥ L ∞ ( X ) .
Applying this inequality to every local mean in the two equivalent quasi-norms yields
∥ f ∥ B 2 , ∞ α ( X ) ≤ C α , X ∥ f ∥ B ∞ , ∞ α ( X ) ≤ C ∥ f ∥ C k , β ( X ) , \|f\|_{B^{\alpha}_{2,\infty}(X)}\le C_{\alpha,X}\|f\|_{B^{\alpha}_{\infty,\infty}(X)}\le C\|f\|_{C^{k,\beta}(X)}, ∥ f ∥ B 2 , ∞ α ( X ) ≤ C α , X ∥ f ∥ B ∞ , ∞ α ( X ) ≤ C ∥ f ∥ C k , β ( X ) ,
which completes the proof.
Proof of Lemma 3.1 . At a point where η ∈ { 0 , 1 } \eta\in\{0,1\} η ∈ { 0 , 1 } , since 0 ≤ ϕ ( t ) ≤ ∣ t ∣ + log 2 ≤ ∣ t ∣ + 2 0\le\phi(t)\le|t|+\log2\le|t|+2 0 ≤ ϕ ( t ) ≤ ∣ t ∣ + log 2 ≤ ∣ t ∣ + 2 , one has ϕ ( t ) 2 ≤ ( ∣ t ∣ + 2 ) ϕ ( t ) \phi(t)^2\le(|t|+2)\phi(t) ϕ ( t ) 2 ≤ ( ∣ t ∣ + 2 ) ϕ ( t ) , which proves the inequality at either endpoint. Fix now η ∈ ( 0 , 1 ) \eta\in(0,1) η ∈ ( 0 , 1 ) and t ∈ R t\in\mathbb{R} t ∈ R , and set
A : = ϕ ( t ) + log η , B : = ϕ ( − t ) + log ( 1 − η ) . A:=\phi(t)+\log\eta,\qquad B:=\phi(-t)+\log(1-\eta). A := ϕ ( t ) + log η , B := ϕ ( − t ) + log ( 1 − η ) .
It is enough to prove
η A 2 + ( 1 − η ) B 2 ≤ ( ∣ t ∣ + 2 ) ( η A + ( 1 − η ) B ) . (25) \eta A^2 + (1-\eta)B^2 \le(|t|+2)(\eta A+(1-\eta)B). \tag*{(25)} η A 2 + ( 1 − η ) B 2 ≤ ( ∣ t ∣ + 2 ) ( η A + ( 1 − η ) B ) . (25)
Let
F ( η , t ) : = ( ∣ t ∣ + 2 ) ( η A + ( 1 − η ) B ) − ( η A 2 + ( 1 − η ) B 2 ) . F(\eta,t):=(|t|+2)(\eta A+(1-\eta)B)-(\eta A^2+(1-\eta)B^2). F ( η , t ) := ( ∣ t ∣ + 2 ) ( η A + ( 1 − η ) B ) − ( η A 2 + ( 1 − η ) B 2 ) .
Direct differentiation gives
∂ η F ( η , t ) = ( ∣ t ∣ − A − B ) ( A − B ) . \partial_{\eta}F(\eta,t)=(|t|-A-B)(A-B). ∂ η F ( η , t ) = ( ∣ t ∣ − A − B ) ( A − B ) .
Note that
∣ t ∣ − ϕ ( t ) − ϕ ( − t ) = − 2 log ( 1 + e − ∣ t ∣ ) ≥ − 2 log 2 ≥ log η + log ( 1 − η ) . |t|-\phi(t)-\phi(-t)=-2\log(1+e^{-|t|})\ge-2\log2\ge\log\eta+\log(1-\eta). ∣ t ∣ − ϕ ( t ) − ϕ ( − t ) = − 2 log ( 1 + e − ∣ t ∣ ) ≥ − 2 log 2 ≥ log η + log ( 1 − η ) .
Hence ∣ t ∣ − A − B ≥ 0 |t|-A-B\ge0 ∣ t ∣ − A − B ≥ 0 . Moreover, A − B A-B A − B is strictly increasing in η \eta η and vanishes exactly at η ∗ = e t / ( 1 + e t ) \eta^*=e^t/(1+e^t) η ∗ = e t / ( 1 + e t ) . At this point A = B = 0 A=B=0 A = B = 0 , so F ( η ∗ , t ) = 0 F(\eta^*,t)=0 F ( η ∗ , t ) = 0 . Thus F ( η , t ) ≥ 0 F(\eta,t)\ge0 F ( η , t ) ≥ 0 , proving (A.1). Applying this inequality with t = f ( x ) t=f(x) t = f ( x ) and integrating over P X P_X P X completes the proof. □
Proof of Theorem 3.2 . By the definition of the infimum in (4) and comparison with the zero function, choose f λ ∈ F f_\lambda\in\mathcal{F} f λ ∈ F such that
a : = λ ∥ f λ ∥ F 2 + R ( f λ ) − R ∗ ≤ min { a F ( λ ) + λ , log 2 − R ∗ } . a:=\lambda\lVert f_\lambda\rVert_{\mathcal F}^2+\mathcal R(f_\lambda)-\mathcal R^*
\le\min\{a_{\mathcal F}(\lambda)+\lambda,\log2-\mathcal R^*\}. a := λ ∥ f λ ∥ F 2 + R ( f λ ) − R ∗ ≤ min { a F ( λ ) + λ , log 2 − R ∗ } .
For f ∈ F f\in\mathcal F f ∈ F , set
g f ( z ) : = λ ∥ f ∥ F 2 + ϕ ( y f ( x ) ) − λ ∥ f λ ∥ F 2 − ϕ ( y f λ ( x ) ) , z = ( x , y ) . g_f(z):=\lambda\lVert f\rVert_{\mathcal F}^2+\phi(yf(x))-\lambda\lVert f_\lambda\rVert_{\mathcal F}^2-\phi(yf_\lambda(x)),\qquad z=(x,y). g f ( z ) := λ ∥ f ∥ F 2 + ϕ ( y f ( x )) − λ ∥ f λ ∥ F 2 − ϕ ( y f λ ( x )) , z = ( x , y ) .
We see that
E P g f = λ ∥ f ∥ F 2 + R ( f ) − λ ∥ f λ ∥ F 2 − R ( f λ ) = λ ∥ f ∥ F 2 + R ( f ) − R ∗ − a . \mathbb E_Pg_f=\lambda\lVert f\rVert_{\mathcal F}^2+\mathcal R(f)-\lambda\lVert f_\lambda\rVert_{\mathcal F}^2-\mathcal R(f_\lambda)
=\lambda\lVert f\rVert_{\mathcal F}^2+\mathcal R(f)-\mathcal R^*-a. E P g f = λ ∥ f ∥ F 2 + R ( f ) − λ ∥ f λ ∥ F 2 − R ( f λ ) = λ ∥ f ∥ F 2 + R ( f ) − R ∗ − a .
Comparison of the empirical objective with its value at zero gives λ ∥ f ^ D , λ ∥ F 2 ≤ log 2 \lambda\lVert\widehat f_{D,\lambda}\rVert_{\mathcal F}^2\le\log2 λ ∥ f D , λ ∥ F 2 ≤ log 2 . The choice of f λ f_\lambda f λ and the nonnegativity of the excess risk give λ ∥ f λ ∥ F 2 ≤ a ≤ log 2 \lambda\lVert f_\lambda\rVert_{\mathcal F}^2\le a\le\log2 λ ∥ f λ ∥ F 2 ≤ a ≤ log 2 . Thus both f λ f_\lambda f λ and f ^ D , λ \widehat f_{D,\lambda} f D , λ belong to λ − 1 / 2 B F \lambda^{-1/2}B_{\mathcal F} λ − 1/2 B F .
We only consider f ∈ λ − 1 / 2 B F f\in\lambda^{-1/2}B_{\mathcal F} f ∈ λ − 1/2 B F , and define the slice
F λ , ε : = { f ∈ λ − 1 / 2 B F : E P g f = ε } . \mathcal F_{\lambda,\varepsilon}:=\{f\in\lambda^{-1/2}B_{\mathcal F}:\mathbb E_Pg_f=\varepsilon\}. F λ , ε := { f ∈ λ − 1/2 B F : E P g f = ε } .
The map f ↦ E P g f f\mapsto\mathbb E_Pg_f f ↦ E P g f is continuous with respect to the norm ∥ ⋅ ∥ F \lVert\cdot\rVert_{\mathcal F} ∥ ⋅ ∥ F because the logistic loss is 1-Lipschitz and ∥ f ∥ ∞ ≤ ∥ f ∥ F \lVert f\rVert_\infty\le\lVert f\rVert_{\mathcal F} ∥ f ∥ ∞ ≤ ∥ f ∥ F . Hence the slice F λ , ε \mathcal F_{\lambda,\varepsilon} F λ , ε is a closed separable subset of F \mathcal F F .
For each fixed sample set D D D , the map f ↦ E D g f f\mapsto\mathbb E_Dg_f f ↦ E D g f is convex. Indeed, f ↦ ∥ f ∥ F 2 f\mapsto\lVert f\rVert_{\mathcal F}^2 f ↦ ∥ f ∥ F 2 is convex, and each f ↦ ϕ ( Y i f ( X i ) ) f\mapsto\phi(Y_if(X_i)) f ↦ ϕ ( Y i f ( X i )) is convex as the composition of the convex loss ϕ \phi ϕ with a linear functional.
Let ε > 0 \varepsilon>0 ε > 0 and suppose that f ∈ λ − 1 / 2 B F f\in\lambda^{-1/2}B_{\mathcal F} f ∈ λ − 1/2 B F satisfies E D g f ≤ 0 \mathbb E_Dg_f\le0 E D g f ≤ 0 and E P g f ≥ ε \mathbb E_Pg_f\ge\varepsilon E P g f ≥ ε . For 0 ≤ θ ≤ 1 0\le\theta\le1 0 ≤ θ ≤ 1 , set f θ = ( 1 − θ ) f λ + θ f f_\theta=(1-\theta)f_\lambda+\theta f f θ = ( 1 − θ ) f λ + θ f . Since both endpoints belong to the convex ball λ − 1 / 2 B F \lambda^{-1/2}B_{\mathcal F} λ − 1/2 B F , so does every f θ f_\theta f θ . The convexity of f ↦ E D g f f\mapsto\mathbb E_Dg_f f ↦ E D g f and g f λ = 0 g_{f_\lambda}=0 g f λ = 0 give
E D g f θ ≤ ( 1 − θ ) E D g f λ + θ E D g f = θ E D g f ≤ 0. \mathbb E_Dg_{f_\theta}\le(1-\theta)\mathbb E_Dg_{f_\lambda}+\theta\mathbb E_Dg_f=\theta\mathbb E_Dg_f\le0. E D g f θ ≤ ( 1 − θ ) E D g f λ + θ E D g f = θ E D g f ≤ 0.
Moreover, the continuous map θ ↦ E P g f θ \theta\mapsto\mathbb E_Pg_{f_\theta} θ ↦ E P g f θ takes the value 0 at θ = 0 \theta=0 θ = 0 and the value E P g f ≥ ε \mathbb E_Pg_f\ge\varepsilon E P g f ≥ ε at θ = 1 \theta=1 θ = 1 . The intermediate value theorem therefore gives a θ ∈ [ 0 , 1 ] \theta\in[0,1] θ ∈ [ 0 , 1 ] for which f θ ∈ F λ , ε f_\theta\in\mathcal F_{\lambda,\varepsilon} f θ ∈ F λ , ε and E D g f θ ≤ 0 \mathbb E_Dg_{f_\theta}\le0 E D g f θ ≤ 0 . To apply this reduction to the estimator, note that E D g f ^ D , λ ≤ 0 \mathbb E_Dg_{\widehat f_{D,\lambda}}\le0 E D g f D , λ ≤ 0 and f ^ D , λ ∈ λ − 1 / 2 B F \widehat f_{D,\lambda}\in\lambda^{-1/2}B_{\mathcal F} f D , λ ∈ λ − 1/2 B F . The preceding argument shows that whenever E P g f ^ D , λ ≥ ε \mathbb E_Pg_{\widehat f_{D,\lambda}}\ge\varepsilon E P g f D , λ ≥ ε , there is an f ∈ F λ , ε f\in\mathcal F_{\lambda,\varepsilon} f ∈ F λ , ε with E D g f ≤ 0 \mathbb E_Dg_f\le0 E D g f ≤ 0 . Consequently,
P D ∼ P n { E P g f ^ D , λ ≥ ε } ≤ P D ∼ P n { E D g f ≤ 0 for some f ∈ F λ , ε } . \mathbb P_{D\sim P^n}\{\mathbb E_Pg_{\widehat f_{D,\lambda}}\ge\varepsilon\}
\le\mathbb P_{D\sim P^n}\{\mathbb E_Dg_f\le0\text{ for some }f\in\mathcal F_{\lambda,\varepsilon}\}. P D ∼ P n { E P g f D , λ ≥ ε } ≤ P D ∼ P n { E D g f ≤ 0 for some f ∈ F λ , ε } .
Moreover, the definitions of a a a and g f g_f g f give
R ( f ^ D , λ ) − R ∗ = a + E P g f ^ D , λ − λ ∥ f ^ D , λ ∥ F 2 ≤ a + E P g f ^ D , λ . \mathcal R(\widehat f_{D,\lambda})-\mathcal R^*
=a+\mathbb E_Pg_{\widehat f_{D,\lambda}}-\lambda\lVert\widehat f_{D,\lambda}\rVert_{\mathcal F}^2
\le a+\mathbb E_Pg_{\widehat f_{D,\lambda}}. R ( f D , λ ) − R ∗ = a + E P g f D , λ − λ ∥ f D , λ ∥ F 2 ≤ a + E P g f D , λ .
Thus, bounding the probability on the right by e − t e^{-t} e − t ensures that the excess risk is at most a + ε a+\varepsilon a + ε with probability at least 1 − e − t 1-e^{-t} 1 − e − t . It therefore suffices to control the event on the slice F λ , ε \mathcal F_{\lambda,\varepsilon} F λ , ε .
Assume ε ≥ a + λ \varepsilon\ge a+\lambda ε ≥ a + λ . For every f ∈ F λ , ε f\in\mathcal F_{\lambda,\varepsilon} f ∈ F λ , ε ,
∥ f ∥ F 2 ≤ a + ε λ ≤ 2 ε λ , ∥ f λ ∥ F 2 ≤ a λ ≤ ε λ . (26) \lVert f\rVert_{\mathcal F}^2\le\frac{a+\varepsilon}{\lambda}\le\frac{2\varepsilon}{\lambda},\qquad
\lVert f_\lambda\rVert_{\mathcal F}^2\le\frac{a}{\lambda}\le\frac{\varepsilon}{\lambda}. \tag*{(26)} ∥ f ∥ F 2 ≤ λ a + ε ≤ λ 2 ε , ∥ f λ ∥ F 2 ≤ λ a ≤ λ ε . (26)
Define the loss difference
q f ( x , y ) : = ϕ ( y f ( x ) ) − ϕ ( y f λ ( x ) ) . q_f(x,y):=\phi(yf(x))-\phi(yf_\lambda(x)). q f ( x , y ) := ϕ ( y f ( x )) − ϕ ( y f λ ( x )) .
and the auxiliary class
Q λ , ε : = { q f : f ∈ F λ , ε } . \mathcal{Q}_{\lambda,\varepsilon}:=\{q_f:f\in\mathcal{F}_{\lambda,\varepsilon}\}. Q λ , ε := { q f : f ∈ F λ , ε } .
Since f , f λ ∈ λ − 1 / 2 B F f,f_\lambda\in\lambda^{-1/2}B_{\mathcal{F}} f , f λ ∈ λ − 1/2 B F and 0 < λ ≤ 1 0<\lambda\le1 0 < λ ≤ 1 ,
ε = E P g f ≤ λ ∥ f ∥ F 2 + R ( f ) ≤ 1 + ∥ f ∥ ∞ + log 2 ≤ 1 + λ − 1 / 2 + log 2 ≤ 3 / λ . \varepsilon=\mathbb{E}_{P}g_f\le\lambda\|f\|_{\mathcal{F}}^2+\mathcal{R}(f)\le1+\|f\|_\infty+\log2\le1+\lambda^{-1/2}+\log2\le3/\lambda. ε = E P g f ≤ λ ∥ f ∥ F 2 + R ( f ) ≤ 1 + ∥ f ∥ ∞ + log 2 ≤ 1 + λ − 1/2 + log 2 ≤ 3/ λ .
We first establish boundedness in the supremum norm. The 1-Lipschitz property of ϕ \phi ϕ and (A.2) give
∥ q f ∥ ∞ ≤ ∥ f − f λ ∥ ∞ ≤ ∥ f ∥ F + ∥ f λ ∥ F ≤ C ε 1 / 2 λ − 1 / 2 . (27) \|q_f\|_\infty\le\|f-f_\lambda\|_\infty\le\|f\|_{\mathcal{F}}+\|f_\lambda\|_{\mathcal{F}}\le C\varepsilon^{1/2}\lambda^{-1/2}. \tag*{(27)} ∥ q f ∥ ∞ ≤ ∥ f − f λ ∥ ∞ ≤ ∥ f ∥ F + ∥ f λ ∥ F ≤ C ε 1/2 λ − 1/2 . (27)
Next, we use the Bernstein condition in Lemma 3.1 to bound the second moment. It gives
E P ( q f 2 ) ≤ 2 E P [ ( ϕ ( Y f ) − ψ ) 2 ] + 2 E P [ ( ϕ ( Y f λ ) − ψ ) 2 ] ≤ 2 ( ∥ f ∥ ∞ + 2 ) ( R ( f ) − R ∗ ) + 2 ( ∥ f λ ∥ ∞ + 2 ) ( R ( f λ ) − R ∗ ) ≤ C ε 3 / 2 λ − 1 / 2 . (28) \begin{aligned}
\mathbb{E}_{P}(q_f^2)&\le2\mathbb{E}_{P}\left[(\phi(Yf)-\psi)^2\right]+2\mathbb{E}_{P}\left[(\phi(Yf_\lambda)-\psi)^2\right]\\
&\le2(\|f\|_\infty+2)(\mathcal{R}(f)-\mathcal{R}^*)+2(\|f_\lambda\|_\infty+2)(\mathcal{R}(f_\lambda)-\mathcal{R}^*)\\
&\le C\varepsilon^{3/2}\lambda^{-1/2}. \tag*{(28)}
\end{aligned} E P ( q f 2 ) ≤ 2 E P [ ( ϕ ( Y f ) − ψ ) 2 ] + 2 E P [ ( ϕ ( Y f λ ) − ψ ) 2 ] ≤ 2 ( ∥ f ∥ ∞ + 2 ) ( R ( f ) − R ∗ ) + 2 ( ∥ f λ ∥ ∞ + 2 ) ( R ( f λ ) − R ∗ ) ≤ C ε 3/2 λ − 1/2 . (28)
The last inequality uses (A.2) and the bounds R ( f ) − R ∗ ≤ a + ε ≤ 2 ε \mathcal{R}(f)-\mathcal{R}^*\le a+\varepsilon\le2\varepsilon R ( f ) − R ∗ ≤ a + ε ≤ 2 ε and R ( f λ ) − R ∗ ≤ a ≤ ε \mathcal{R}(f_\lambda)-\mathcal{R}^*\le a\le\varepsilon R ( f λ ) − R ∗ ≤ a ≤ ε .
The regularization terms are constants and cancel in ( E P − E D ) g f (\mathbb{E}_{P}-\mathbb{E}_{D})g_f ( E P − E D ) g f , so that ( E P − E D ) g f = ( E P − E D ) q f (\mathbb{E}_{P}-\mathbb{E}_{D})g_f=(\mathbb{E}_{P}-\mathbb{E}_{D})q_f ( E P − E D ) g f = ( E P − E D ) q f . Hence we only need to estimate the complexity of the auxiliary class Q λ , ε \mathcal{Q}_{\lambda,\varepsilon} Q λ , ε . Fix a sample set D D D and define Φ D : L 2 ( D X ) → L 2 ( D ) \Phi_D:L_2(D_X)\to L_2(D) Φ D : L 2 ( D X ) → L 2 ( D ) on the sample points by
( Φ D u ) ( X j , Y j ) : = ϕ ( Y j u ( X j ) ) − ϕ ( Y j f λ ( X j ) ) , 1 ≤ j ≤ n . (\Phi_Du)(X_j,Y_j):=\phi(Y_ju(X_j))-\phi(Y_jf_\lambda(X_j)),\qquad1\le j\le n. ( Φ D u ) ( X j , Y j ) := ϕ ( Y j u ( X j )) − ϕ ( Y j f λ ( X j )) , 1 ≤ j ≤ n .
The 1-Lipschitz property of ϕ \phi ϕ implies that, for any u , v ∈ L 2 ( D X ) u,v\in L_2(D_X) u , v ∈ L 2 ( D X ) ,
∥ Φ D u − Φ D v ∥ L 2 ( D ) 2 = 1 n ∑ j = 1 n ∣ ϕ ( Y j u ( X j ) ) − ϕ ( Y j v ( X j ) ) ∣ 2 ≤ 1 n ∑ j = 1 n ∣ u ( X j ) − v ( X j ) ∣ 2 = ∥ u − v ∥ L 2 ( D X ) 2 . \begin{aligned}
\|\Phi_Du-\Phi_Dv\|_{L_2(D)}^2
&=\frac{1}{n}\sum_{j=1}^{n}\left|\phi(Y_ju(X_j))-\phi(Y_jv(X_j))\right|^2\\
&\le\frac{1}{n}\sum_{j=1}^{n}|u(X_j)-v(X_j)|^2\\
&=\|u-v\|_{L_2(D_X)}^2.
\end{aligned} ∥ Φ D u − Φ D v ∥ L 2 ( D ) 2 = n 1 j = 1 ∑ n ∣ ϕ ( Y j u ( X j )) − ϕ ( Y j v ( X j )) ∣ 2 ≤ n 1 j = 1 ∑ n ∣ u ( X j ) − v ( X j ) ∣ 2 = ∥ u − v ∥ L 2 ( D X ) 2 .
In particular, Φ D \Phi_D Φ D does not increase distances. Hence, by (A.2) and (5) ,
e i ( Q λ , ε , L 2 ( D ) ) ≤ e i ( F λ , ε , L 2 ( D X ) ) ≤ ( 2 ε λ ) 1 / 2 e i ( id : F → L 2 ( D X ) ) . e_i(\mathcal{Q}_{\lambda,\varepsilon},L_2(D))\le e_i(\mathcal{F}_{\lambda,\varepsilon},L_2(D_X))\le\left(\frac{2\varepsilon}{\lambda}\right)^{1/2}e_i(\operatorname{id}:\mathcal{F}\to L_2(D_X)). e i ( Q λ , ε , L 2 ( D )) ≤ e i ( F λ , ε , L 2 ( D X )) ≤ ( λ 2 ε ) 1/2 e i ( id : F → L 2 ( D X )) .
Consequently, (5) shows that
E D ∼ P n e i ( Q λ , ε , L 2 ( D ) ) ≤ C ( ε λ ) 1 / 2 A 1 / ( 2 p ) i − 1 / ( 2 p ) . (29) \mathbb{E}_{D\sim P^n}e_i(\mathcal{Q}_{\lambda,\varepsilon},L_2(D))\le C\left(\frac{\varepsilon}{\lambda}\right)^{1/2}A^{1/(2p)}i^{-1/(2p)}. \tag*{(29)} E D ∼ P n e i ( Q λ , ε , L 2 ( D )) ≤ C ( λ ε ) 1/2 A 1/ ( 2 p ) i − 1/ ( 2 p ) . (29)
Since the map ( z , f ) ↦ q f ( z ) (z,f)\mapsto q_f(z) ( z , f ) ↦ q f ( z ) is measurable in z z z and continuous in f f f , the class Q λ , ε \mathcal{Q}_{\lambda,\varepsilon} Q λ , ε is a separable Carathéodory class. Applying [15 Theorem 7.16] with (A.3), (A.4) and (A.5) gives the local Rademacher bound
E D ∼ P n Rad D ( Q λ , ε ) ≤ C [ ε ( 3 − p ) / 4 λ − ( 1 + p ) / 4 ( A n ) 1 / 2 + ε 1 / 2 λ − 1 / 2 ( A n ) 1 / ( 1 + p ) ] . (30) \mathbb{E}_{D\sim P^n}\operatorname{Rad}_D(\mathcal{Q}_{\lambda,\varepsilon})\le C\left[\varepsilon^{(3-p)/4}\lambda^{-(1+p)/4}\left(\frac{A}{n}\right)^{1/2}+\varepsilon^{1/2}\lambda^{-1/2}\left(\frac{A}{n}\right)^{1/(1+p)}\right]. \tag*{(30)} E D ∼ P n Rad D ( Q λ , ε ) ≤ C [ ε ( 3 − p ) /4 λ − ( 1 + p ) /4 ( n A ) 1/2 + ε 1/2 λ − 1/2 ( n A ) 1/ ( 1 + p ) ] . (30)
The constant C C C is uniform for p ∈ ( 0 , 1 / 2 ] p\in(0,1/2] p ∈ ( 0 , 1/2 ] .
Let h f : = E P q f − q f h_f:=\mathbb{E}_{P}q_f-q_f h f := E P q f − q f . Since ( E P − E D ) g f = ( E P − E D ) q f (\mathbb{E}_{P}-\mathbb{E}_{D})g_f=(\mathbb{E}_{P}-\mathbb{E}_{D})q_f ( E P − E D ) g f = ( E P − E D ) q f , one has
P D ∼ P n { E D g f ≤ 0 for some f ∈ F λ , ε } ≤ P D ∼ P n { sup f ∈ F λ , ε E D h f ≥ ε } . \mathbb{P}_{D\sim P^n}\{ \mathbb{E}_{D}g_f\le0\text{ for some }f\in\mathcal{F}_{\lambda,\varepsilon}\}
\le\mathbb{P}_{D\sim P^n}\left\{\sup_{f\in\mathcal{F}_{\lambda,\varepsilon}}\mathbb{E}_{D}h_f\ge\varepsilon\right\}. P D ∼ P n { E D g f ≤ 0 for some f ∈ F λ , ε } ≤ P D ∼ P n { f ∈ F λ , ε sup E D h f ≥ ε } .
By symmetrization [15 Proposition 7.10] ,
E D ∼ P n sup f ∈ F λ , ε E D h f ≤ 2 E D ∼ P n Rad D ( Q λ , ε ) . \mathbb{E}_{D\sim P^n}\sup_{f\in\mathcal{F}_{\lambda,\varepsilon}}\mathbb{E}_{D}h_f
\le2\mathbb{E}_{D\sim P^n}\operatorname{Rad}_D(\mathcal{Q}_{\lambda,\varepsilon}). E D ∼ P n f ∈ F λ , ε sup E D h f ≤ 2 E D ∼ P n Rad D ( Q λ , ε ) .
For every f ∈ F λ , ε f \in\mathcal{F}_{\lambda,\varepsilon} f ∈ F λ , ε , we have E P h f = 0 \mathbb{E}_{P}h_f=0 E P h f = 0 , ∥ h f ∥ ∞ ≤ 2 ∥ q f ∥ ∞ \|h_f\|_{\infty} \le2\|q_f\|_{\infty} ∥ h f ∥ ∞ ≤ 2∥ q f ∥ ∞ , and E P ( h f 2 ) ≤ E P ( q f 2 ) \mathbb{E}_{P}(h_f^2) \le\mathbb{E}_{P}(q_f^2) E P ( h f 2 ) ≤ E P ( q f 2 ) . Thus (A.3) and (A.4) provide the supremum and variance bounds required by the simplified Talagrand’s inequality [15 Theorem 7.5] . Together with the preceding symmetrization bound, it gives, for every t ≥ 1 t \ge1 t ≥ 1 , with probability at least 1 − e − t 1-e^{-t} 1 − e − t ,
sup f ∈ F λ , ε E D h f < C [ E D ∼ P n Rad D ( Q λ , ε ) + ε 3 / 4 λ − 1 / 4 ( t n ) 1 / 2 + ε 1 / 2 λ − 1 / 2 t n ] . \sup_{f\in\mathcal{F}_{\lambda,\varepsilon}} \mathbb{E}_{D}h_f<C\left[\mathbb{E}_{D\sim P^n}\operatorname{Rad}_{D}(\mathcal{Q}_{\lambda,\varepsilon})+\varepsilon^{3/4}\lambda^{-1/4}\left(\frac{t}{n}\right)^{1/2}+\varepsilon^{1/2}\lambda^{-1/2}\frac{t}{n}\right]. f ∈ F λ , ε sup E D h f < C [ E D ∼ P n Rad D ( Q λ , ε ) + ε 3/4 λ − 1/4 ( n t ) 1/2 + ε 1/2 λ − 1/2 n t ] .
Substituting (A.6), we obtain
P D ∼ P n { sup f ∈ F λ , ε E D h f ≥ ε } ≤ e − t . \mathbb{P}_{D\sim P^n}\left\{\sup_{f\in\mathcal{F}_{\lambda,\varepsilon}}\mathbb{E}_{D}h_f\ge\varepsilon\right\}\le e^{-t}. P D ∼ P n { f ∈ F λ , ε sup E D h f ≥ ε } ≤ e − t .
provided that ε ≥ a + λ \varepsilon\ge a+\lambda ε ≥ a + λ , as required for the localized estimates, and
ε ≥ C [ ε ( 3 − p ) / 4 λ − ( 1 + p ) / 4 ( A n ) 1 / 2 + ε 1 / 2 λ − 1 / 2 ( A n ) 1 / ( 1 + p ) + ε 3 / 4 λ − 1 / 4 ( t n ) 1 / 2 + ε 1 / 2 λ − 1 / 2 t n ] . \varepsilon\ge C\left[\varepsilon^{(3-p)/4}\lambda^{-(1+p)/4}\left(\frac{A}{n}\right)^{1/2}+\varepsilon^{1/2}\lambda^{-1/2}\left(\frac{A}{n}\right)^{1/(1+p)}+\varepsilon^{3/4}\lambda^{-1/4}\left(\frac{t}{n}\right)^{1/2}+\varepsilon^{1/2}\lambda^{-1/2}\frac{t}{n}\right]. ε ≥ C [ ε ( 3 − p ) /4 λ − ( 1 + p ) /4 ( n A ) 1/2 + ε 1/2 λ − 1/2 ( n A ) 1/ ( 1 + p ) + ε 3/4 λ − 1/4 ( n t ) 1/2 + ε 1/2 λ − 1/2 n t ] .
Choose
ε : = C 0 [ a + λ + λ − 1 ( A n ) 2 / ( 1 + p ) + λ − 1 t 2 n 2 ] . \varepsilon:=C_0\left[a+\lambda+\lambda^{-1}\left(\frac{A}{n}\right)^{2/(1+p)}+\lambda^{-1}\frac{t^2}{n^2}\right]. ε := C 0 [ a + λ + λ − 1 ( n A ) 2/ ( 1 + p ) + λ − 1 n 2 t 2 ] .
with C 0 C_0 C 0 a sufficiently large numerical constant, so that both conditions hold. If F λ , ε \mathcal{F}_{\lambda,\varepsilon} F λ , ε is empty, the earlier event inclusion gives E P g f ^ D , λ < ε \mathbb{E}_{P}g_{\widehat{f}_{D,\lambda}}<\varepsilon E P g f D , λ < ε for every D D D . Otherwise, the probability bound above and the same event inclusion show that E P g f ^ D , λ < ε \mathbb{E}_{P}g_{\widehat{f}_{D,\lambda}}<\varepsilon E P g f D , λ < ε with probability at least 1 − e − t 1-e^{-t} 1 − e − t .
Finally, note that a ≤ a F ( λ ) + λ a\le a_{\mathcal{F}}(\lambda)+\lambda a ≤ a F ( λ ) + λ . We therefore have, when E P g f ^ D , λ < ε \mathbb{E}_{P}g_{\widehat{f}_{D,\lambda}}<\varepsilon E P g f D , λ < ε ,
R ( f ^ D , λ ) − R ∗ = a + E P g f ^ D , λ − λ ∥ f ^ D , λ ∥ F 2 ≤ a + ε ≤ C [ a F ( λ ) + λ + λ − 1 ( A n ) 2 / ( 1 + p ) + λ − 1 t 2 n 2 ] . \begin{aligned}
\mathcal{R}(\widehat{f}_{D,\lambda})-\mathcal{R}^{*}=a+\mathbb{E}_{P}g_{\widehat{f}_{D,\lambda}}-\lambda\|\widehat{f}_{D,\lambda}\|_{\mathcal{F}}^{2} \\
&\le a+\varepsilon\\
&\le C\left[a_{\mathcal{F}}(\lambda)+\lambda+\lambda^{-1}\left(\frac{A}{n}\right)^{2/(1+p)}+\lambda^{-1}\frac{t^2}{n^2}\right].
\end{aligned} R ( f D , λ ) − R ∗ = a + E P g f D , λ − λ ∥ f D , λ ∥ F 2 ≤ a + ε ≤ C [ a F ( λ ) + λ + λ − 1 ( n A ) 2/ ( 1 + p ) + λ − 1 n 2 t 2 ] .
This completes the proof of (6) . □
Proof of Corollary 3.3. Use the function f λ f_\lambda f λ , the number a a a , and the classes F λ , ε \mathcal{F}_{\lambda,\varepsilon} F λ , ε and Q λ , ε \mathcal{Q}_{\lambda,\varepsilon} Q λ , ε from the proof of Theorem 3.2 . The localization argument and the bounds (A.2), (A.3), and (A.4) do not use the entropy assumption. Fix ε ≥ a + λ \varepsilon\ge a+\lambda ε ≥ a + λ and first suppose that the slice F λ , ε \mathcal{F}_{\lambda,\varepsilon} F λ , ε is nonempty. By (A.2), the slice is contained in ( 2 ε / λ ) 1 / 2 B F (2\varepsilon/\lambda)^{1/2}B_{\mathcal{F}} ( 2 ε / λ ) 1/2 B F . The map f ↦ q f f\mapsto q_f f ↦ q f is 1-Lipschitz in the supremum norm. Thus (7) gives
log N ( δ ( 2 ε λ ) 1 / 2 , Q λ , ε , ∥ ⋅ ∥ ∞ ) ≤ A δ − 2 p , b ≤ δ ≤ 1. \log N\left(\delta\left(\frac{2\varepsilon}{\lambda}\right)^{1/2},\mathcal{Q}_{\lambda,\varepsilon},\|\cdot\|_{\infty}\right)\le A\delta^{-2p},\qquad b\le\delta\le1. log N ( δ ( λ 2 ε ) 1/2 , Q λ , ε , ∥ ⋅ ∥ ∞ ) ≤ A δ − 2 p , b ≤ δ ≤ 1.
Choose such a cover at radius b ( 2 ε / λ ) 1 / 2 b(2\varepsilon/\lambda)^{1/2} b ( 2 ε / λ ) 1/2 . Since the centers of the covering balls need not belong to Q λ , ε \mathcal{Q}_{\lambda,\varepsilon} Q λ , ε , we select one member from each ball that intersects it. These representatives form a finite set N ⊂ Q λ , ε \mathcal{N}\subset\mathcal{Q}_{\lambda,\varepsilon} N ⊂ Q λ , ε such that
sup q ∈ Q λ , ε min v ∈ N ∥ q − v ∥ ∞ ≤ 2 b ( 2 ε λ ) 1 / 2 , log ∣ N ∣ ≤ A b − 2 p . (31) \sup_{q\in\mathcal{Q}_{\lambda,\varepsilon}}\min_{v\in\mathcal{N}}\|q-v\|_{\infty}\le2b\left(\frac{2\varepsilon}{\lambda}\right)^{1/2},\qquad\log|\mathcal{N}|\le Ab^{-2p}. \tag*{(31)} q ∈ Q λ , ε sup v ∈ N min ∥ q − v ∥ ∞ ≤ 2 b ( λ 2 ε ) 1/2 , log ∣ N ∣ ≤ A b − 2 p . (31)
Set N 0 : = N ∪ { 0 } \mathcal{N}_0:=\mathcal{N}\cup\{0\} N 0 := N ∪ { 0 } and consider a δ ( 2 ε / λ ) 1 / 2 \delta(2\varepsilon/\lambda)^{1/2} δ ( 2 ε / λ ) 1/2 -cover of it. For b ≤ δ ≤ 1 b\le\delta\le1 b ≤ δ ≤ 1 , we can cover N \mathcal{N} N using the corresponding δ ( 2 ε / λ ) 1 / 2 \delta(2\varepsilon/\lambda)^{1/2} δ ( 2 ε / λ ) 1/2 -cover of Q λ , ε \mathcal{Q}_{\lambda,\varepsilon} Q λ , ε and add one ball for zero. For 0 < δ < b 0<\delta<b 0 < δ < b , use one ball for each member of N 0 \mathcal{N}_0 N 0 . Both cases give
log N ( δ ( 2 ε λ ) 1 / 2 , N 0 , ∥ ⋅ ∥ ∞ ) ≤ 2 A δ − 2 p , 0 < δ ≤ 1. \log N\left(\delta\left(\frac{2\varepsilon}{\lambda}\right)^{1/2},\mathcal{N}_0,\|\cdot\|_{\infty}\right)\le2A\delta^{-2p},\qquad0<\delta\le1. log N ( δ ( λ 2 ε ) 1/2 , N 0 , ∥ ⋅ ∥ ∞ ) ≤ 2 A δ − 2 p , 0 < δ ≤ 1.
By (A.3), every member of N 0 \mathcal{N}_0 N 0 has supremum norm at most 2 ( 2 ε / λ ) 1 / 2 2(2\varepsilon/\lambda)^{1/2} 2 ( 2 ε / λ ) 1/2 . For 1 < δ ≤ 2 1 < \delta\le2 1 < δ ≤ 2 , the ( 2 ε / λ ) 1 / 2 (2\varepsilon/\lambda)^{1/2} ( 2 ε / λ ) 1/2 -cover gives the same estimate with 2 A 2A 2 A replaced by 4 A 4A 4 A . Inverting this covering bound and using ∥ ⋅ ∥ L 2 ( D ) ≤ ∥ ⋅ ∥ ∞ \|\cdot\|_{L_2(D)} \le\|\cdot\|_\infty ∥ ⋅ ∥ L 2 ( D ) ≤ ∥ ⋅ ∥ ∞ yields, for a numerical constant C 0 ≥ 1 C_0 \ge1 C 0 ≥ 1 ,
e i ( N 0 , L 2 ( D ) ) ≤ ( 2 ε λ ) 1 / 2 ( C 0 A ) 1 / ( 2 p ) i − 1 / ( 2 p ) , i ≥ 1. e_i(\mathcal{N}_0,L_2(D)) \le\left(\frac{2\varepsilon}{\lambda}\right)^{1/2}(C_0A)^{1/(2p)}i^{-1/(2p)}, \qquad i \ge1. e i ( N 0 , L 2 ( D )) ≤ ( λ 2 ε ) 1/2 ( C 0 A ) 1/ ( 2 p ) i − 1/ ( 2 p ) , i ≥ 1.
Apply [15 Theorem 7.16] to the finite class N 0 \mathcal{N}_0 N 0 , using (A.3) and (A.4). As in (A.6), this gives
E D ∼ P n Rad D ( N 0 ) ≤ C [ ε ( 3 − p ) / 4 λ − ( 1 + p ) / 4 ( A n ) 1 / 2 + ε 1 / 2 λ − 1 / 2 ( A n ) 1 / ( 1 + p ) ] . \mathbb{E}_{D\sim P^n}\operatorname{Rad}_D(\mathcal{N}_0) \le C\left[\varepsilon^{(3-p)/4}\lambda^{-(1+p)/4}\left(\frac{A}{n}\right)^{1/2}+\varepsilon^{1/2}\lambda^{-1/2}\left(\frac{A}{n}\right)^{1/(1+p)}\right]. E D ∼ P n Rad D ( N 0 ) ≤ C [ ε ( 3 − p ) /4 λ − ( 1 + p ) /4 ( n A ) 1/2 + ε 1/2 λ − 1/2 ( n A ) 1/ ( 1 + p ) ] .
with C C C uniform for p ∈ ( 0 , 1 / 2 ] p \in(0,1/2] p ∈ ( 0 , 1/2 ] . For every sample set D D D , (A.7) also implies
Rad D ( Q λ , ε ) ≤ Rad D ( N 0 ) + 2 b ( 2 ε λ ) 1 / 2 . \operatorname{Rad}_D(\mathcal{Q}_{\lambda,\varepsilon}) \le\operatorname{Rad}_D(\mathcal{N}_0)+2b\left(\frac{2\varepsilon}{\lambda}\right)^{1/2}. Rad D ( Q λ , ε ) ≤ Rad D ( N 0 ) + 2 b ( λ 2 ε ) 1/2 .
Consequently, the estimate (A.6) remains valid with an additional term C b ( ε / λ ) 1 / 2 Cb(\varepsilon/\lambda)^{1/2} C b ( ε / λ ) 1/2 on the right-hand side. Then, applying symmetrization and the simplified Talagrand’s inequality and following the argument in the proof of Theorem 3.2 yields (8) , with the additional term contributing C b 2 / λ Cb^2/\lambda C b 2 / λ to the bound.
Proof of Theorem 4.1 . The compact embedded manifold X X X satisfies the uniform extrinsic volume-growth bound [11 ]
sup y ∈ X V ( X ∩ B R d ( y , r ) ) ≤ C X r ϱ , 0 < r ≤ 1. \sup_{y\in X}V(X\cap B_{\mathbb{R}^d}(y,r))\le C_Xr^\varrho,\qquad0<r\le1. y ∈ X sup V ( X ∩ B R d ( y , r )) ≤ C X r ϱ , 0 < r ≤ 1.
This gives the uniform integral bound
sup y ∈ X ∫ X exp ( − ∥ x − y ∥ R d 2 γ 2 ) d V ( x ) ≤ C X γ ϱ . \sup_{y\in X}\int_X\exp\left(-\frac{\|x-y\|_{\mathbb{R}^d}^2}{\gamma^2}\right)dV(x)\le C_X\gamma^\varrho. y ∈ X sup ∫ X exp ( − γ 2 ∥ x − y ∥ R d 2 ) d V ( x ) ≤ C X γ ϱ .
By the reproducing property,
∥ T γ f ∥ H γ 2 ( X ) 2 = ( γ π ) − 2 ϱ ∫ X ∫ X f ( x ) f ( y ) ⟨ k γ ( ⋅ , x ) , k γ ( ⋅ , y ) ⟩ H γ ( X ) d V ( x ) d V ( y ) = ( γ π ) − 2 ϱ ∫ X ∫ X exp ( − ∥ x − y ∥ R d 2 γ 2 ) f ( x ) f ( y ) d V ( x ) d V ( y ) ≤ 1 2 ( γ π ) 2 ϱ ∫ X ∫ X ( ∣ f ( x ) ∣ 2 + ∣ f ( y ) ∣ 2 ) exp ( − ∥ x − y ∥ R d 2 γ 2 ) d V ( x ) d V ( y ) ≤ C X γ − ϱ ∥ f ∥ L 2 ( X ) 2 . \begin{aligned}
\|T_\gamma f\|_{\mathcal{H}_{\gamma^2}(X)}^2=(\gamma\sqrt{\pi})^{-2\varrho}\int_X\int_X f(x)f(y)\langle k_\gamma(\cdot,x),k_\gamma(\cdot,y)\rangle_{\mathcal{H}_\gamma(X)}dV(x)dV(y) \\
&=(\gamma\sqrt{\pi})^{-2\varrho}\int_X\int_X\exp\left(-\frac{\|x-y\|_{\mathbb{R}^d}^2}{\gamma^2}\right)f(x)f(y)dV(x)dV(y) \\
&\le\frac{1}{2(\gamma\sqrt{\pi})^{2\varrho}}\int_X\int_X\left(|f(x)|^2+|f(y)|^2\right)\exp\left(-\frac{\|x-y\|_{\mathbb{R}^d}^2}{\gamma^2}\right)dV(x)dV(y) \\
&\le C_X\gamma^{-\varrho}\|f\|_{L_2(X)}^2.
\end{aligned} ∥ T γ f ∥ H γ 2 ( X ) 2 = ( γ π ) − 2 ϱ ∫ X ∫ X f ( x ) f ( y ) ⟨ k γ ( ⋅ , x ) , k γ ( ⋅ , y ) ⟩ H γ ( X ) d V ( x ) d V ( y ) = ( γ π ) − 2 ϱ ∫ X ∫ X exp ( − γ 2 ∥ x − y ∥ R d 2 ) f ( x ) f ( y ) d V ( x ) d V ( y ) ≤ 2 ( γ π ) 2 ϱ 1 ∫ X ∫ X ( ∣ f ( x ) ∣ 2 + ∣ f ( y ) ∣ 2 ) exp ( − γ 2 ∥ x − y ∥ R d 2 ) d V ( x ) d V ( y ) ≤ C X γ − ϱ ∥ f ∥ L 2 ( X ) 2 .
The Gaussian RKHS inclusion criterion [15 Proposition 4.46] gives
H j γ ( X ) ⊂ H γ ( X ) , ∥ u ∥ H γ ( X ) ≤ j d / 2 ∥ u ∥ H j γ ( X ) , j ≥ 1. \mathcal{H}_{j\gamma}(X)\subset\mathcal{H}_{\gamma}(X),\qquad\|u\|_{\mathcal{H}_{\gamma}(X)}\le j^{d/2}\|u\|_{\mathcal{H}_{j\gamma}(X)},\qquad j\ge1. H j γ ( X ) ⊂ H γ ( X ) , ∥ u ∥ H γ ( X ) ≤ j d /2 ∥ u ∥ H j γ ( X ) , j ≥ 1.
Applying these two estimates to every term in (9) proves the result.
Proof of Lemma 4.2. We first prove the estimate for f ∈ C ∞ ( X ) f\in C^\infty(X) f ∈ C ∞ ( X ) and then extend it by density. Write L = 2 N + 2 L=2N+2 L = 2 N + 2 . Choose 0 < δ < inj ( X ) / 4 0<\delta<\operatorname{inj}(X)/4 0 < δ < inj ( X ) /4 and χ 0 ∈ C ∞ ( [ 0 , ∞ ) ) \chi_0\in C^\infty([0,\infty)) χ 0 ∈ C ∞ ([ 0 , ∞ )) with 0 ≤ χ 0 ≤ 1 0\le\chi_0\le1 0 ≤ χ 0 ≤ 1 such that χ 0 = 1 \chi_0=1 χ 0 = 1 on [ 0 , δ 2 ] [0,\delta^2] [ 0 , δ 2 ] and χ 0 = 0 \chi_0=0 χ 0 = 0 on [ 4 δ 2 , ∞ ) [4\delta^2,\infty) [ 4 δ 2 , ∞ ) . Define
χ ( x , y ) : = χ 0 ( d X ( x , y ) 2 ) , \chi(x,y):=\chi_0(d_X(x,y)^2), χ ( x , y ) := χ 0 ( d X ( x , y ) 2 ) ,
T γ = T γ l o c + T γ f a r . T_\gamma=T_\gamma^{\mathrm{loc}}+T_\gamma^{\mathrm{far}}. T γ = T γ loc + T γ far .
T γ l o c f ( x ) : = ( γ π ) − ϱ ∫ X exp ( − ∥ x − y ∥ R d 2 γ 2 ) χ ( x , y ) f ( y ) d V ( y ) , T_\gamma^{\mathrm{loc}}f(x):=(\gamma\sqrt{\pi})^{-\varrho}\int_X\exp\left(-\frac{\|x-y\|_{\mathbb{R}^d}^2}{\gamma^2}\right)\chi(x,y)f(y)dV(y), T γ loc f ( x ) := ( γ π ) − ϱ ∫ X exp ( − γ 2 ∥ x − y ∥ R d 2 ) χ ( x , y ) f ( y ) d V ( y ) ,
T γ f a r f ( x ) : = ( γ π ) − ϱ ∫ X exp ( − ∥ x − y ∥ R d 2 γ 2 ) ( 1 − χ ( x , y ) ) f ( y ) d V ( y ) . T_\gamma^{\mathrm{far}}f(x):=(\gamma\sqrt{\pi})^{-\varrho}\int_X\exp\left(-\frac{\|x-y\|_{\mathbb{R}^d}^2}{\gamma^2}\right)(1-\chi(x,y))f(y)dV(y). T γ far f ( x ) := ( γ π ) − ϱ ∫ X exp ( − γ 2 ∥ x − y ∥ R d 2 ) ( 1 − χ ( x , y )) f ( y ) d V ( y ) .
Since X X X is compact and its smooth embedding in R d \mathbb{R}^{d} R d induces the Riemannian metric, there is a c δ > 0 c_{\delta}>0 c δ > 0 such that
d X ( x , y ) ≥ δ ⟹ ∥ x − y ∥ R d ≥ c δ . d_{X}(x,y)\ge\delta\Longrightarrow\lVert x-y\rVert_{\mathbb{R}^{d}}\ge c_{\delta}. d X ( x , y ) ≥ δ ⟹ ∥ x − y ∥ R d ≥ c δ .
Let K γ f a r ( x , y ) K_{\gamma}^{\mathrm{far}}(x,y) K γ far ( x , y ) denote the integral kernel of T γ f a r T_{\gamma}^{\mathrm{far}} T γ far . Since 1 − χ ( x , y ) 1-\chi(x,y) 1 − χ ( x , y ) vanishes for d X ( x , y ) ≤ δ d_{X}(x,y)\le\delta d X ( x , y ) ≤ δ and 0 ≤ χ ≤ 1 0\le\chi\le1 0 ≤ χ ≤ 1 , we have
sup x ∈ X ∫ X ∣ K γ f a r ( x , y ) ∣ d V ( y ) ≤ ( γ π ) − ϱ V ( X ) exp ( − c δ 2 / γ 2 ) . \sup_{x\in X}\int_{X}\left|K_{\gamma}^{\mathrm{far}}(x,y)\right|dV(y)\le(\gamma\sqrt{\pi})^{-\varrho}V(X)\exp\left(-c_{\delta}^{2}/\gamma^{2}\right). x ∈ X sup ∫ X K γ far ( x , y ) d V ( y ) ≤ ( γ π ) − ϱ V ( X ) exp ( − c δ 2 / γ 2 ) .
By symmetry, the same bound holds with x x x and y y y interchanged. The Schur test therefore gives
∥ T γ f a r ∥ L 2 → L 2 ≤ C X , δ γ − ϱ exp ( − c δ 2 / γ 2 ) . (32) \left\|T_{\gamma}^{\mathrm{far}}\right\|_{L^{2}\to L^{2}}\le C_{X,\delta}\gamma^{-\varrho}\exp\left(-c_{\delta}^{2}/\gamma^{2}\right). \tag*{(32)} T γ far L 2 → L 2 ≤ C X , δ γ − ϱ exp ( − c δ 2 / γ 2 ) . (32)
For the local term, let E x : T x X → X \mathcal{E}_{x}:T_{x}X\to X E x : T x X → X denote the exponential map at x x x . Since χ ( x , y ) \chi(x,y) χ ( x , y ) vanishes for d X ( x , y ) ≥ 2 δ d_{X}(x,y)\ge2\delta d X ( x , y ) ≥ 2 δ and δ < inj ( X ) / 4 \delta<\operatorname{inj}(X)/4 δ < inj ( X ) /4 , we can put y = E x ( u ) y=\mathcal{E}_{x}(u) y = E x ( u ) , and write
q x ( u ) : = ∥ x − E x ( u ) ∥ R d 2 , d V ( E x ( u ) ) = J x ( u ) d u x . q_{x}(u):=\left\|x-\mathcal{E}_{x}(u)\right\|_{\mathbb{R}^{d}}^{2},\qquad dV(\mathcal{E}_{x}(u))=J_{x}(u)\,du_{x}. q x ( u ) := ∥ x − E x ( u ) ∥ R d 2 , d V ( E x ( u )) = J x ( u ) d u x .
By [21 Lemma D.12] , uniformly in x x x and for ∥ u ∥ < 2 δ \lVert u\rVert<2\delta ∥ u ∥ < 2 δ ,
q x ( u ) = ∥ u ∥ 2 + O ( ∥ u ∥ 4 ) , J x ( u ) = 1 + O ( ∥ u ∥ 2 ) , c ∥ u ∥ 2 ≤ q x ( u ) ≤ C ∥ u ∥ 2 . (33) \begin{aligned}
q_{x}(u)&=\lVert u\rVert^{2}+O(\lVert u\rVert^{4}),\\
J_{x}(u)&=1+O(\lVert u\rVert^{2}),\\
c\lVert u\rVert^{2}&\le q_{x}(u)\le C\lVert u\rVert^{2}. \tag*{(33)}
\end{aligned} q x ( u ) J x ( u ) c ∥ u ∥ 2 = ∥ u ∥ 2 + O (∥ u ∥ 4 ) , = 1 + O (∥ u ∥ 2 ) , ≤ q x ( u ) ≤ C ∥ u ∥ 2 . (33)
Rescale the tangent coordinates by the bandwidth via u = γ z u=\gamma z u = γ z . We then define
Ψ x ( t , z ) : = { q x ( t z ) / t 2 , t ≠ 0 , ∥ t z ∥ < 3 δ , ∥ z ∥ 2 , t = 0 , \Psi_{x}(t,z):=
\begin{cases}
q_{x}(tz)/t^{2}, & t\ne0,\ \lVert tz\rVert<3\delta,\\
\lVert z\rVert^{2}, & t=0,
\end{cases} Ψ x ( t , z ) := { q x ( t z ) / t 2 , ∥ z ∥ 2 , t = 0 , ∥ t z ∥ < 3 δ , t = 0 ,
and extend the product
exp ( − Ψ x ( t , z ) ) J x ( t z ) χ 0 ( t 2 ∥ z ∥ 2 ) \exp\left(-\Psi_{x}(t,z)\right)J_{x}(tz)\chi_{0}\left(t^{2}\lVert z\rVert^{2}\right) exp ( − Ψ x ( t , z ) ) J x ( t z ) χ 0 ( t 2 ∥ z ∥ 2 )
by zero when ∥ t z ∥ ≥ 3 δ \lVert tz\rVert\ge3\delta ∥ t z ∥ ≥ 3 δ . Since q x q_{x} q x and its first u u u -derivative vanish at the origin, applying Taylor’s formula to the scalar function θ ↦ q x ( θ t z ) \theta\mapsto q_{x}(\theta tz) θ ↦ q x ( θ t z ) and dividing by t 2 t^{2} t 2 for t ≠ 0 t\ne0 t = 0 gives
Ψ x ( t , z ) = ∫ 0 1 ( 1 − θ ) D 2 q x ( θ t z ) [ z , z ] d θ . (34) \Psi_{x}(t,z)=\int_{0}^{1}(1-\theta)D^{2}q_{x}(\theta tz)[z,z]\,d\theta. \tag*{(34)} Ψ x ( t , z ) = ∫ 0 1 ( 1 − θ ) D 2 q x ( θ t z ) [ z , z ] d θ . (34)
For fixed x , z x,z x , z , the right-hand side is smooth in t t t near zero and takes the value ∥ z ∥ 2 \lVert z\rVert^{2} ∥ z ∥ 2 there, so this representation also holds at zero. The cutoff restricts the nonzero terms to ∥ t z ∥ < 2 δ \lVert tz\rVert<2\delta ∥ t z ∥ < 2 δ . On this region, Ψ x ( t , z ) ≥ c ∥ z ∥ 2 \Psi_{x}(t,z)\ge c\lVert z\rVert^{2} Ψ x ( t , z ) ≥ c ∥ z ∥ 2 , and its integral representation bounds each t t t -derivative by a polynomial in ∥ z ∥ \lVert z\rVert ∥ z ∥ . The same holds for the derivatives of J x ( t z ) J_{x}(tz) J x ( t z ) and the cutoff. The chain rule and Leibniz’s rule therefore give, for ∣ t ∣ ≤ 1 \lvert t\rvert\le1 ∣ t ∣ ≤ 1 and 0 ≤ ℓ ≤ L 0\le\ell\le L 0 ≤ ℓ ≤ L ,
∣ ∂ t ℓ [ exp ( − Ψ x ( t , z ) ) J x ( t z ) χ 0 ( t 2 ∥ z ∥ 2 ) ] ∣ ≤ C ℓ ( 1 + ∥ z ∥ 2 ) M ℓ / 2 exp ( − c ∥ z ∥ 2 ) . (35) \left|\partial_{t}^{\ell}\left[\exp\left(-\Psi_{x}(t,z)\right)J_{x}(tz)\chi_{0}\left(t^{2}\lVert z\rVert^{2}\right)\right]\right|\le C_{\ell}(1+\lVert z\rVert^{2})^{M_{\ell}/2}\exp(-c\lVert z\rVert^{2}). \tag*{(35)} ∂ t ℓ [ exp ( − Ψ x ( t , z ) ) J x ( t z ) χ 0 ( t 2 ∥ z ∥ 2 ) ] ≤ C ℓ ( 1 + ∥ z ∥ 2 ) M ℓ /2 exp ( − c ∥ z ∥ 2 ) . (35)
Here, C ℓ > 0 C_{\ell}>0 C ℓ > 0 , M ℓ ∈ N 0 M_{\ell}\in\mathbb{N}_{0} M ℓ ∈ N 0 , and c > 0 c>0 c > 0 are independent of x , t , z x,t,z x , t , z . Thus each derivative is integrable over T x X T_{x}X T x X . Let I t f ( x ) \mathcal{I}_{t}f(x) I t f ( x ) denote the rescaled integral
I t f ( x ) = π − ϱ / 2 ∫ T x X exp ( − Ψ x ( t , z ) ) J x ( t z ) χ 0 ( t 2 ∥ z ∥ 2 ) f ( E x ( t z ) ) d z x . \mathcal{I}_{t}f(x)=\pi^{-\varrho/2}\int_{T_{x}X}\exp\left(-\Psi_{x}(t,z)\right)J_{x}(tz)\chi_{0}\left(t^{2}\lVert z\rVert^{2}\right)f(\mathcal{E}_{x}(tz))\,dz_{x}. I t f ( x ) = π − ϱ /2 ∫ T x X exp ( − Ψ x ( t , z ) ) J x ( t z ) χ 0 ( t 2 ∥ z ∥ 2 ) f ( E x ( t z )) d z x .
Then I t = T t l o c \mathcal{I}_{t}=T_{t}^{\mathrm{loc}} I t = T t loc for t ≠ 0 t\ne0 t = 0 and I 0 = I \mathcal{I}_{0}=I I 0 = I .
We now apply Taylor’s formula at the origin to the integrand as a function of t t t , retaining terms up to order L − 1 L-1 L − 1 . For fixed x x x and z z z , write the full integrand as the scalar function
F x , z ( s ) = exp ( − Ψ x ( s , z ) ) J x ( s z ) χ 0 ( s 2 ∥ z ∥ 2 ) f ( E x ( s z ) ) . F_{x,z}(s)=\exp\left(-\Psi_{x}(s,z)\right)J_{x}(sz)\chi_{0}\left(s^{2}\lVert z\rVert^{2}\right)f(\mathcal{E}_{x}(sz)). F x , z ( s ) = exp ( − Ψ x ( s , z ) ) J x ( sz ) χ 0 ( s 2 ∥ z ∥ 2 ) f ( E x ( sz )) .
Let F x , z ( j ) F_{x,z}^{(j)} F x , z ( j ) denote differentiation with respect to s s s . For the factor f ( E x ( t z ) ) f(\mathcal{E}_{x}(tz)) f ( E x ( t z )) , differentiation along the radial geodesic gives, for 1 ≤ ℓ ≤ L 1\le\ell\le L 1 ≤ ℓ ≤ L and ∥ t z ∥ < 2 δ \lVert tz\rVert<2\delta ∥ t z ∥ < 2 δ ,
d ℓ d t ℓ f ( E x ( t z ) ) = ( ∇ ℓ f ) E x ( t z ) ( P x , E x ( t z ) z , … , P x , E x ( t z ) z ⏟ ℓ a r g u m e n t s ) . \frac{d^{\ell}}{dt^{\ell}}f(\mathcal{E}_{x}(tz))=(\nabla^{\ell}f)_{\mathcal{E}_{x}(tz)}\left(\underbrace{P_{x,\mathcal{E}_{x}(tz)}z,\ldots,P_{x,\mathcal{E}_{x}(tz)}z}_{\ell\ \mathrm{arguments}}\right). d t ℓ d ℓ f ( E x ( t z )) = ( ∇ ℓ f ) E x ( t z ) ℓ arguments P x , E x ( t z ) z , … , P x , E x ( t z ) z .
Since parallel transport preserves lengths, this yields
∣ d ℓ d t ℓ f ( E x ( t z ) ) ∣ ≤ ∥ z ∥ ℓ ∣ ∇ ℓ f ∣ ( E x ( t z ) ) . (36) \left|\frac{d^\ell}{dt^\ell}f(\mathcal{E}_x(tz))\right| \leq\|z\|^\ell|\nabla^\ell f|(\mathcal{E}_x(tz)). \tag*{(36)} d t ℓ d ℓ f ( E x ( t z )) ≤ ∥ z ∥ ℓ ∣ ∇ ℓ f ∣ ( E x ( t z )) . (36)
By Leibniz’s rule, each term in F x , z ( L ) ( s ) F_{x,z}^{(L)}(s) F x , z ( L ) ( s ) has ℓ \ell ℓ derivatives on f ( E x ( s z ) ) f(\mathcal{E}_x(sz)) f ( E x ( sz )) and L − ℓ L-\ell L − ℓ derivatives on the other factor. Combining (A.12) with (A.11), and using ∥ z ∥ ℓ ≤ ( 1 + ∥ z ∥ 2 ) ℓ / 2 \|z\|^\ell\leq(1+\|z\|^2)^{\ell/2} ∥ z ∥ ℓ ≤ ( 1 + ∥ z ∥ 2 ) ℓ /2 , gives
∣ F x , z ( L ) ( s ) ∣ ≤ C ( 1 + ∥ z ∥ 2 ) M / 2 exp ( − c ∥ z ∥ 2 ) 1 { ∥ s z ∥ ≤ 2 δ } ∑ ℓ = 0 L ∣ ∇ ℓ f ∣ ( E x ( s z ) ) , (37) |F_{x,z}^{(L)}(s)| \leq C(1+\|z\|^2)^{M/2}\exp(-c\|z\|^2)\mathbf{1}_{\{\|sz\|\leq2\delta\}}\sum_{\ell=0}^{L}|\nabla^\ell f|(\mathcal{E}_x(sz)), \tag*{(37)} ∣ F x , z ( L ) ( s ) ∣ ≤ C ( 1 + ∥ z ∥ 2 ) M /2 exp ( − c ∥ z ∥ 2 ) 1 { ∥ sz ∥ ≤ 2 δ } ℓ = 0 ∑ L ∣ ∇ ℓ f ∣ ( E x ( sz )) , (37)
M : = max 0 ≤ ℓ ≤ L ( M L − ℓ + ℓ ) . M:=\max_{0\leq\ell\leq L}(M_{L-\ell}+\ell). M := 0 ≤ ℓ ≤ L max ( M L − ℓ + ℓ ) .
Applying Taylor’s formula to F x , z F_{x,z} F x , z at zero and then integrating over z z z gives the following expansion, with remainder operator R L , t R_{L,t} R L , t :
I t f ( x ) = π − ϱ / 2 ∑ j = 0 L − 1 t j j ! ∫ T x X F x , z ( j ) ( 0 ) d z x + R L , t f ( x ) , \mathcal{I}_t f(x)=\pi^{-\varrho/2}\sum_{j=0}^{L-1}\frac{t^j}{j!}\int_{T_xX}F_{x,z}^{(j)}(0)\,dz_x+R_{L,t}f(x), I t f ( x ) = π − ϱ /2 j = 0 ∑ L − 1 j ! t j ∫ T x X F x , z ( j ) ( 0 ) d z x + R L , t f ( x ) ,
R L , t f ( x ) : = π − ϱ / 2 t L ( L − 1 ) ! ∫ 0 1 ( 1 − θ ) L − 1 ∫ T x X F x , z ( L ) ( θ t ) d z x d θ . R_{L,t}f(x):=\frac{\pi^{-\varrho/2}t^L}{(L-1)!}\int_0^1(1-\theta)^{L-1}\int_{T_xX}F_{x,z}^{(L)}(\theta t)\,dz_x\,d\theta. R L , t f ( x ) := ( L − 1 )! π − ϱ /2 t L ∫ 0 1 ( 1 − θ ) L − 1 ∫ T x X F x , z ( L ) ( θ t ) d z x d θ .
Applying (A.13), we obtain
∣ R L , t f ( x ) ∣ ≤ C ∣ t ∣ L ∑ ℓ = 0 L ∫ 0 1 ( 1 − θ ) L − 1 ∫ T x X ( 1 + ∥ z ∥ 2 ) M / 2 exp ( − c ∥ z ∥ 2 ) × 1 { ∥ θ t z ∥ ≤ 2 δ } ∣ ∇ ℓ f ∣ ( E x ( θ t z ) ) d z x d θ . (38) |R_{L,t}f(x)|\leq C|t|^L\sum_{\ell=0}^{L}\int_0^1(1-\theta)^{L-1}\int_{T_xX}(1+\|z\|^2)^{M/2}\exp(-c\|z\|^2)
\times\mathbf{1}_{\{\|\theta tz\|\leq2\delta\}}|\nabla^\ell f|(\mathcal{E}_x(\theta tz))\,dz_x\,d\theta. \tag*{(38)} ∣ R L , t f ( x ) ∣ ≤ C ∣ t ∣ L ℓ = 0 ∑ L ∫ 0 1 ( 1 − θ ) L − 1 ∫ T x X ( 1 + ∥ z ∥ 2 ) M /2 exp ( − c ∥ z ∥ 2 ) × 1 { ∥ θ t z ∥ ≤ 2 δ } ∣ ∇ ℓ f ∣ ( E x ( θ t z )) d z x d θ . (38)
Choose finitely many pairs of coordinate neighborhoods U a ⋐ V a U_a\Subset V_a U a ⋐ V a such that a partition of unity is supported in the U a U_a U a , each V a V_a V a carries a smooth orthonormal frame, and φ a ( V a ) \varphi_a(V_a) φ a ( V a ) is convex (for example, an open Euclidean ball). For each a a a , let e 1 , … , e ϱ e_1,\ldots,e_\varrho e 1 , … , e ϱ denote the frame on V a V_a V a . For fixed θ ∈ [ 0 , 1 ] \theta\in[0,1] θ ∈ [ 0 , 1 ] , ∣ t ∣ ≤ 1 |t|\leq1 ∣ t ∣ ≤ 1 , and z ∈ R ϱ z\in\mathbb{R}^{\varrho} z ∈ R ϱ with ∥ θ t z ∥ ≤ 2 δ \|\theta tz\|\leq2\delta ∥ θ t z ∥ ≤ 2 δ , define
Θ θ , t , z ( x ) = E x ( θ t ∑ i = 1 ϱ z i e i ( x ) ) , x ∈ V a . \Theta_{\theta,t,z}(x)=\mathcal{E}_x\left(\theta t\sum_{i=1}^{\varrho}z_i e_i(x)\right),\qquad x\in V_a. Θ θ , t , z ( x ) = E x ( θ t i = 1 ∑ ϱ z i e i ( x ) ) , x ∈ V a .
After decreasing δ \delta δ , we have Θ θ , t , z ( U a ) ⊂ V a \Theta_{\theta,t,z}(U_a)\subset V_a Θ θ , t , z ( U a ) ⊂ V a , and the coordinate representation φ a ∘ Θ θ , t , z ∘ φ a − 1 \varphi_a\circ\Theta_{\theta,t,z}\circ\varphi_a^{-1} φ a ∘ Θ θ , t , z ∘ φ a − 1 is uniformly C 1 C^1 C 1 -close to the identity on a convex neighborhood of φ a ( U a ‾ ) \varphi_a(\overline{U_a}) φ a ( U a ) . Consequently, Θ θ , t , z \Theta_{\theta,t,z} Θ θ , t , z is injective on U a U_a U a and satisfies
0 < c ≤ ∣ det D x Θ θ , t , z ∣ ≤ C , x ∈ U a . 0<c\leq|\det D_x\Theta_{\theta,t,z}|\leq C,\qquad x\in U_a. 0 < c ≤ ∣ det D x Θ θ , t , z ∣ ≤ C , x ∈ U a .
Here ∣ det D x Θ θ , t , z ∣ |\det D_x\Theta_{\theta,t,z}| ∣ det D x Θ θ , t , z ∣ denotes the Jacobian with respect to Riemannian volume on the domain and image. The constants are independent of a , x , θ , t , z a,x,\theta,t,z a , x , θ , t , z . For every 0 ≤ ℓ ≤ L 0\leq\ell\leq L 0 ≤ ℓ ≤ L , the lower Jacobian bound and change of variables under the injective map Θ θ , t , z \Theta_{\theta,t,z} Θ θ , t , z give
∫ U a ∣ ∇ ℓ f ∣ 2 ( Θ θ , t , z ( x ) ) d V ( x ) ≤ c − 1 ∫ U a ∣ ∇ ℓ f ∣ 2 ( Θ θ , t , z ( x ) ) ∣ det D x Θ θ , t , z ∣ d V ( x ) = c − 1 ∫ Θ θ , t , z ( U a ) ∣ ∇ ℓ f ∣ 2 ( y ) d V ( y ) ≤ c − 1 ∫ V a ∣ ∇ ℓ f ∣ 2 ( y ) d V ( y ) . \begin{aligned}
\int_{U_a}|\nabla^\ell f|^2(\Theta_{\theta,t,z}(x))\,dV(x)
\leq c^{-1}\int_{U_a}|\nabla^\ell f|^2(\Theta_{\theta,t,z}(x))|\det D_x\Theta_{\theta,t,z}|\,dV(x) \\
&=c^{-1}\int_{\Theta_{\theta,t,z}(U_a)}|\nabla^\ell f|^2(y)\,dV(y)
\leq c^{-1}\int_{V_a}|\nabla^\ell f|^2(y)\,dV(y).
\end{aligned} ∫ U a ∣ ∇ ℓ f ∣ 2 ( Θ θ , t , z ( x )) d V ( x ) ≤ c − 1 ∫ U a ∣ ∇ ℓ f ∣ 2 ( Θ θ , t , z ( x )) ∣ det D x Θ θ , t , z ∣ d V ( x ) = c − 1 ∫ Θ θ , t , z ( U a ) ∣ ∇ ℓ f ∣ 2 ( y ) d V ( y ) ≤ c − 1 ∫ V a ∣ ∇ ℓ f ∣ 2 ( y ) d V ( y ) .
Taking square roots yields
∥ ∇ ℓ f ∘ Θ θ , t , z ∥ L 2 ( U a ) ≤ C ∥ ∇ ℓ f ∥ L 2 ( V a ) . (39) \|\nabla^\ell f\circ\Theta_{\theta,t,z}\|_{L^2(U_a)}\leq C\|\nabla^\ell f\|_{L^2(V_a)}. \tag*{(39)} ∥ ∇ ℓ f ∘ Θ θ , t , z ∥ L 2 ( U a ) ≤ C ∥ ∇ ℓ f ∥ L 2 ( V a ) . (39)
On each U a U_a U a , the orthonormal frame identifies integration over T x X T_xX T x X with Lebesgue integration over R ϱ \mathbb{R}^{\varrho} R ϱ . In these coordinates, the Gaussian weight and the condition ∥ θ t z ∥ ≤ 2 δ \|\theta tz\|\leq2\delta ∥ θ t z ∥ ≤ 2 δ are independent of x x x . Applying Minkowski’s integral inequality to the pointwise remainder bound moves the L 2 ( U a ) L_2(U_a) L 2 ( U a ) norm inside the integrals in θ \theta θ and z z z . The estimates (A.14) and (A.15) then give
∥ R L , t f ∥ L 2 ( U a ) ≤ C ∣ t ∣ L ∑ ℓ = 0 L ∫ 0 1 ( 1 − θ ) L − 1 ∫ { z ∈ R ϱ : ∥ θ t z ∥ ≤ 2 δ } ( 1 + ∥ z ∥ 2 ) M / 2 exp ( − c ∥ z ∥ 2 ) × ∥ ∇ ℓ f ∘ Θ θ , t , z ∥ L 2 ( U a ) d z d θ ≤ C ∣ t ∣ L ∑ ℓ = 0 L ∥ ∇ ℓ f ∥ L 2 ( V a ) ∫ 0 1 ( 1 − θ ) L − 1 ∫ R ϱ ( 1 + ∥ z ∥ 2 ) M / 2 exp ( − c ∥ z ∥ 2 ) d z d θ . \begin{aligned}
\lVert R_{L,t}f\rVert_{L_2(U_a)}
&\le C|t|^L\sum_{\ell=0}^{L}\int_0^1(1-\theta)^{L-1}\int_{\{z\in\mathbb{R}^{\varrho}:\lVert\theta tz\rVert\le2\delta\}}(1+\lVert z\rVert^2)^{M/2}\exp(-c\lVert z\rVert^2)\\
&\qquad\times\lVert\nabla^\ell f\circ\Theta_{\theta,t,z}\rVert_{L_2(U_a)}\,dz\,d\theta\\
&\le C|t|^L\sum_{\ell=0}^{L}\lVert\nabla^\ell f\rVert_{L_2(V_a)}\int_0^1(1-\theta)^{L-1}\int_{\mathbb{R}^{\varrho}}(1+\lVert z\rVert^2)^{M/2}\exp(-c\lVert z\rVert^2)\,dz\,d\theta.
\end{aligned} ∥ R L , t f ∥ L 2 ( U a ) ≤ C ∣ t ∣ L ℓ = 0 ∑ L ∫ 0 1 ( 1 − θ ) L − 1 ∫ { z ∈ R ϱ : ∥ θ t z ∥ ≤ 2 δ } ( 1 + ∥ z ∥ 2 ) M /2 exp ( − c ∥ z ∥ 2 ) × ∥ ∇ ℓ f ∘ Θ θ , t , z ∥ L 2 ( U a ) d z d θ ≤ C ∣ t ∣ L ℓ = 0 ∑ L ∥ ∇ ℓ f ∥ L 2 ( V a ) ∫ 0 1 ( 1 − θ ) L − 1 ∫ R ϱ ( 1 + ∥ z ∥ 2 ) M /2 exp ( − c ∥ z ∥ 2 ) d z d θ .
Both scalar integrals are finite and independent of t t t , so they can be absorbed into C C C . Since the finitely many sets U a U_a U a cover X X X , summing these local bounds yields
∥ R L , t f ∥ L 2 ( X ) ≤ ∑ a ∥ R L , t f ∥ L 2 ( U a ) ≤ C ∣ t ∣ L ∑ ℓ = 0 L ∥ ∇ ℓ f ∥ L 2 ( X ) ≤ C ∣ t ∣ L ∥ f ∥ H L ( X ) . (40) \lVert R_{L,t}f\rVert_{L_2(X)}\le\sum_a\lVert R_{L,t}f\rVert_{L_2(U_a)}\le C|t|^L\sum_{\ell=0}^{L}\lVert\nabla^\ell f\rVert_{L_2(X)}\le C|t|^L\lVert f\rVert_{H^L(X)}. \tag*{(40)} ∥ R L , t f ∥ L 2 ( X ) ≤ a ∑ ∥ R L , t f ∥ L 2 ( U a ) ≤ C ∣ t ∣ L ℓ = 0 ∑ L ∥ ∇ ℓ f ∥ L 2 ( X ) ≤ C ∣ t ∣ L ∥ f ∥ H L ( X ) . (40)
Note that T t = T ∣ t ∣ l o c \mathcal{T}_t=T^{\mathrm{loc}}_{|t|} T t = T ∣ t ∣ loc is even in t t t . Thus its odd Taylor coefficients vanish. Define
A r f : = 1 ( 2 r ) ! ∂ t 2 r T t f ∣ t = 0 . A_r f:=\frac{1}{(2r)!}\left.\partial_t^{2r}\mathcal{T}_t f\right|_{t=0}. A r f := ( 2 r )! 1 ∂ t 2 r T t f t = 0 .
This definition is independent of the truncation order N N N . It is also independent of the auxiliary cutoff as χ 0 \chi_0 χ 0 is identically one near zero.
To compute the term where all 2 r 2r 2 r derivatives fall on f ( E x ( t z ) ) f(\mathcal{E}_x(tz)) f ( E x ( t z )) , fix x x x and use geodesic normal coordinates u u u . The radial derivative formula identifies ( ∇ 2 r f ) x ( z , … , z ) (\nabla^{2r}f)_x(z,\ldots,z) ( ∇ 2 r f ) x ( z , … , z ) with the 2 r 2r 2 r -th directional derivative of f ∘ E x f\circ\mathcal{E}_x f ∘ E x at the origin in direction z z z . When we expand this derivative, all terms with an odd power of any z i z_i z i integrate to zero. The remaining terms are evaluated using
π − 1 / 2 ∫ R exp ( − s 2 ) s 2 k d s = ( 2 k ) ! 4 k k ! , k ≥ 0. \pi^{-1/2}\int_{\mathbb{R}}\exp(-s^2)s^{2k}\,ds=\frac{(2k)!}{4^k k!},\qquad k\ge0. π − 1/2 ∫ R exp ( − s 2 ) s 2 k d s = 4 k k ! ( 2 k )! , k ≥ 0.
Thus, summing over nonnegative integers k 1 , … , k ϱ k_1,\ldots,k_\varrho k 1 , … , k ϱ with total r r r and then applying the multinomial formula, we obtain
π − ϱ / 2 ( 2 r ) ! ∫ R ϱ exp ( − ∥ z ∥ 2 ) ( ∇ 2 r f ) x ( z , … , z ) d z = 1 4 r ∑ k 1 + ⋯ + k ϱ = r ∂ u 1 2 k 1 ⋯ ∂ u ϱ 2 k ϱ ( f ∘ E x ) ( 0 ) k 1 ! ⋯ k ϱ ! = 1 4 r r ! ( ∑ i = 1 ϱ ∂ u i 2 ) r ( f ∘ E x ) ( 0 ) . \begin{aligned}
\frac{\pi^{-\varrho/2}}{(2r)!}\int_{\mathbb{R}^{\varrho}}\exp(-\lVert z\rVert^2)(\nabla^{2r}f)_x(z,\ldots,z)\,dz
&=\frac{1}{4^r}\sum_{k_1+\cdots+k_\varrho=r}\frac{\partial_{u_1}^{2k_1}\cdots\partial_{u_\varrho}^{2k_\varrho}(f\circ\mathcal{E}_x)(0)}{k_1!\cdots k_\varrho!}\\
&=\frac{1}{4^r r!}\left(\sum_{i=1}^{\varrho}\partial_{u_i}^2\right)^r(f\circ\mathcal{E}_x)(0).
\end{aligned} ( 2 r )! π − ϱ /2 ∫ R ϱ exp ( − ∥ z ∥ 2 ) ( ∇ 2 r f ) x ( z , … , z ) d z = 4 r 1 k 1 + ⋯ + k ϱ = r ∑ k 1 ! ⋯ k ϱ ! ∂ u 1 2 k 1 ⋯ ∂ u ϱ 2 k ϱ ( f ∘ E x ) ( 0 ) = 4 r r ! 1 ( i = 1 ∑ ϱ ∂ u i 2 ) r ( f ∘ E x ) ( 0 ) .
By the coordinate formula for Δ X \Delta_X Δ X and the normal-coordinate identities in [12 Proposition 2.46 and Proposition 5.24] , its second-order coefficients agree with those of the Euclidean Laplacian to first order at the origin, while its first-order coefficients vanish there. Consequently, the terms of order 2 r 2r 2 r agree and those of order 2 r − 1 2r-1 2 r − 1 vanish at the origin when taking the r r r -th power, so for r ≥ 1 r\ge1 r ≥ 1 we obtain
1 4 r r ! ( ∑ i = 1 ϱ ∂ u i 2 ) r ( f ∘ E x ) ( 0 ) = 1 4 r r ! ( Δ X r f ) ( x ) + ( B ~ r f ) ( x ) , B ~ r ∈ Diff 2 r − 2 ( X ) . \frac{1}{4^r r!}\left(\sum_{i=1}^{\varrho}\partial_{u_i}^2\right)^r(f\circ\mathcal{E}_x)(0)=\frac{1}{4^r r!}(\Delta_X^r f)(x)+(\widetilde{B}_r f)(x),\qquad\widetilde{B}_r\in\operatorname{Diff}^{2r-2}(X). 4 r r ! 1 ( i = 1 ∑ ϱ ∂ u i 2 ) r ( f ∘ E x ) ( 0 ) = 4 r r ! 1 ( Δ X r f ) ( x ) + ( B r f ) ( x ) , B r ∈ Diff 2 r − 2 ( X ) .
We now consider the remaining Leibniz terms, in which at least one derivative falls on
exp ( − Ψ x ( t , z ) ) J x ( t z ) χ 0 ( t 2 ∥ z ∥ 2 ) . \exp(-\Psi_x(t,z))J_x(tz)\chi_0(t^2\lVert z\rVert^2). exp ( − Ψ x ( t , z )) J x ( t z ) χ 0 ( t 2 ∥ z ∥ 2 ) .
Since q x q_x q x and J x J_x J x are smooth, the estimates (A.9) imply, for every x ∈ X x\in X x ∈ X ,
D u J x ( 0 ) = 0 , D u 3 q x ( 0 ) = 0. D_uJ_x(0)=0,\qquad D_u^3q_x(0)=0. D u J x ( 0 ) = 0 , D u 3 q x ( 0 ) = 0.
Differentiating (A.10) with respect to t t t and evaluating at zero gives, for each fixed x , z x,z x , z ,
∂ t Ψ x ( 0 , z ) = ∫ 0 1 θ ( 1 − θ ) D u 3 q x ( 0 ) [ z , z , z ] d θ = 0. \partial_t\Psi_x(0,z)=\int_0^1\theta(1-\theta)D_u^3q_x(0)[z,z,z]\,d\theta=0. ∂ t Ψ x ( 0 , z ) = ∫ 0 1 θ ( 1 − θ ) D u 3 q x ( 0 ) [ z , z , z ] d θ = 0.
Together with the vanishing of D u J x ( 0 ) D_uJ_x(0) D u J x ( 0 ) and the constancy of the cutoff near zero, this gives
∂ t [ exp ( − Ψ x ( t , z ) ) J x ( t z ) χ 0 ( t 2 ∥ z ∥ 2 ) ] ∣ 0 = 0. \left.\partial_t\left[\exp(-\Psi_x(t,z))J_x(tz)\chi_0(t^2\lVert z\rVert^2)\right]\right|_0=0. ∂ t [ exp ( − Ψ x ( t , z )) J x ( t z ) χ 0 ( t 2 ∥ z ∥ 2 ) ] 0 = 0.
Thus the terms with exactly one derivative on this factor vanish. Every nonzero remaining term has at least two derivatives on this factor and hence at most 2 r − 2 2r-2 2 r − 2 derivatives on f f f . This proves that A r A_r A r has the stated form for r ≥ 1 r\ge1 r ≥ 1 .
Now we are ready to prove the final estimate (11) . For f ∈ C ∞ ( X ) f\in C^\infty(X) f ∈ C ∞ ( X ) and 0 < γ ≤ 1 0<\gamma\le1 0 < γ ≤ 1 , we start from the decomposition
T γ f = T γ f a r f + T γ l o c f = T γ f a r f + I γ f . T_\gamma f=T_\gamma^{\mathrm{far}}f+T_\gamma^{\mathrm{loc}}f=T_\gamma^{\mathrm{far}}f+I_\gamma f. T γ f = T γ far f + T γ loc f = T γ far f + I γ f .
The Taylor expansion of I t f I_t f I t f above retains terms through degree 2 N + 1 2N+1 2 N + 1 . Since its odd coefficients vanish, the definition of A r A_r A r gives
I γ f = ∑ r = 0 N γ 2 r A r f + R L , γ f . I_\gamma f=\sum_{r=0}^{N}\gamma^{2r}A_rf+R_{L,\gamma}f. I γ f = r = 0 ∑ N γ 2 r A r f + R L , γ f .
Substituting this identity into the decomposition yields
T γ f − ∑ r = 0 N γ 2 r A r f = T γ f a r f + R L , γ f . T_\gamma f-\sum_{r=0}^{N}\gamma^{2r}A_rf=T_\gamma^{\mathrm{far}}f+R_{L,\gamma}f. T γ f − r = 0 ∑ N γ 2 r A r f = T γ far f + R L , γ f .
By (A.8) and the fact that γ − ϱ − L exp ( − c δ 2 / γ 2 ) \gamma^{-\varrho-L}\exp(-c_\delta^2/\gamma^2) γ − ϱ − L exp ( − c δ 2 / γ 2 ) is bounded for 0 < γ ≤ 1 0<\gamma\le1 0 < γ ≤ 1 , we have
∥ T γ f a r f ∥ L 2 ( X ) ≤ C X , δ γ − ϱ exp ( − c δ 2 / γ 2 ) ∥ f ∥ L 2 ( X ) ≤ C γ L ∥ f ∥ L 2 ( X ) . \lVert T_\gamma^{\mathrm{far}}f\rVert_{L^2(X)}
\le C_{X,\delta}\gamma^{-\varrho}\exp(-c_\delta^2/\gamma^2)\lVert f\rVert_{L^2(X)}
\le C\gamma^L\lVert f\rVert_{L^2(X)}. ∥ T γ far f ∥ L 2 ( X ) ≤ C X , δ γ − ϱ exp ( − c δ 2 / γ 2 ) ∥ f ∥ L 2 ( X ) ≤ C γ L ∥ f ∥ L 2 ( X ) .
The local remainder estimate (A.16), evaluated at γ \gamma γ , gives
∥ R L , γ f ∥ L 2 ( X ) ≤ C γ L ∥ f ∥ H L ( X ) . \lVert R_{L,\gamma}f\rVert_{L^2(X)}\le C\gamma^L\lVert f\rVert_{H^L(X)}. ∥ R L , γ f ∥ L 2 ( X ) ≤ C γ L ∥ f ∥ H L ( X ) .
Thus,
∥ T γ f − ∑ r = 0 N γ 2 r A r f ∥ L 2 ( X ) ≤ ∥ T γ f a r f ∥ L 2 ( X ) + ∥ R L , γ f ∥ L 2 ( X ) ≤ C γ L ( ∥ f ∥ L 2 ( X ) + ∥ f ∥ H L ( X ) ) ≤ C N , X , d γ 2 N + 2 ∥ f ∥ H 2 N + 2 ( X ) . \left\lVert T_\gamma f-\sum_{r=0}^{N}\gamma^{2r}A_rf\right\rVert_{L^2(X)}
\le\lVert T_\gamma^{\mathrm{far}}f\rVert_{L^2(X)}+\lVert R_{L,\gamma}f\rVert_{L^2(X)}
\le C\gamma^L\left(\lVert f\rVert_{L^2(X)}+\lVert f\rVert_{H^L(X)}\right)
\le C_{N,X,d}\gamma^{2N+2}\lVert f\rVert_{H^{2N+2}(X)}. T γ f − r = 0 ∑ N γ 2 r A r f L 2 ( X ) ≤ ∥ T γ far f ∥ L 2 ( X ) + ∥ R L , γ f ∥ L 2 ( X ) ≤ C γ L ( ∥ f ∥ L 2 ( X ) + ∥ f ∥ H L ( X ) ) ≤ C N , X , d γ 2 N + 2 ∥ f ∥ H 2 N + 2 ( X ) .
This proves (11) for smooth f f f . For f ∈ H L ( X ) f\in H^L(X) f ∈ H L ( X ) , choose f m ∈ C ∞ ( X ) f_m\in C^\infty(X) f m ∈ C ∞ ( X ) with f m → f f_m\to f f m → f in H L ( X ) H^L(X) H L ( X ) . This convergence also holds in L 2 ( X ) L^2(X) L 2 ( X ) . For each fixed γ \gamma γ , the boundedness of T γ : L 2 ( X ) → L 2 ( X ) T_\gamma:L^2(X)\to L^2(X) T γ : L 2 ( X ) → L 2 ( X ) therefore gives T γ f m → T γ f T_\gamma f_m\to T_\gamma f T γ f m → T γ f in L 2 ( X ) L^2(X) L 2 ( X ) . Moreover, each A r A_r A r has order at most 2 r ≤ L 2r\le L 2 r ≤ L , so it is continuous from H L ( X ) H^L(X) H L ( X ) to L 2 ( X ) L^2(X) L 2 ( X ) and A r f m → A r f A_rf_m\to A_rf A r f m → A r f in L 2 ( X ) L^2(X) L 2 ( X ) . Applying the estimate just proved to f m f_m f m and passing to the limit proves (11) for every f ∈ H L ( X ) f\in H^L(X) f ∈ H L ( X ) . We thus complete the proof. □ \square □
Proof of Theorem 4.3. Set
N : = ⌈ α 2 ⌉ , R : = 2 N + 2. N:=\left\lceil\frac{\alpha}{2}\right\rceil,\qquad R:=2N+2. N := ⌈ 2 α ⌉ , R := 2 N + 2.
Then
R > α , 2 N ≤ ⌊ α ⌋ = s − 1. R>\alpha,\qquad2N\le\lfloor\alpha\rfloor=s-1. R > α , 2 N ≤ ⌊ α ⌋ = s − 1.
The definition (9) and the identities (10) give, for every γ > 0 \gamma>0 γ > 0 ,
S γ , s − I = ∑ j = 1 s c j [ T j γ − ∑ r = 0 N ( j γ ) 2 r A r ] . S_{\gamma,s}-I=\sum_{j=1}^{s}c_j\left[T_{j\gamma}-\sum_{r=0}^{N}(j\gamma)^{2r}A_r\right]. S γ , s − I = j = 1 ∑ s c j [ T j γ − r = 0 ∑ N ( j γ ) 2 r A r ] .
Suppose first that s γ ≤ 1 s\gamma\le1 s γ ≤ 1 . Then j γ ≤ 1 j\gamma\le1 j γ ≤ 1 for every 1 ≤ j ≤ s 1 \le j \le s 1 ≤ j ≤ s . Applying (11) to each term on the right and using the triangle inequality yields, for every f ∈ H R ( X ) f \in H^{R}(X) f ∈ H R ( X ) ,
∥ ( S γ , s − I ) f ∥ L 2 ( X ) ≤ C α γ R ∥ f ∥ H R ( X ) . \left\|(S_{\gamma,s}-I)f\right\|_{L_2(X)} \le C_{\alpha}\gamma^{R}\|f\|_{H^{R}(X)}. ∥ ( S γ , s − I ) f ∥ L 2 ( X ) ≤ C α γ R ∥ f ∥ H R ( X ) .
Since the operators T γ : L 2 ( X ) → L 2 ( X ) T_{\gamma}: L_2(X) \to L_2(X) T γ : L 2 ( X ) → L 2 ( X ) are uniformly bounded over γ > 0 \gamma> 0 γ > 0 , the definition (9) and the triangle inequality give
∥ S γ , s − I ∥ L 2 ( X ) → L 2 ( X ) ≤ C α . (41) \left\|S_{\gamma,s}-I\right\|_{L_2(X)\to L_2(X)} \le C_{\alpha}. \tag*{(41)} ∥ S γ , s − I ∥ L 2 ( X ) → L 2 ( X ) ≤ C α . (41)
Hence, for every decomposition f = f 0 + f 1 f=f_0+f_1 f = f 0 + f 1 with f 0 ∈ L 2 ( X ) f_0\in L_2(X) f 0 ∈ L 2 ( X ) and f 1 ∈ H R ( X ) f_1\in H^{R}(X) f 1 ∈ H R ( X ) ,
∥ ( S γ , s − I ) f ∥ L 2 ( X ) ≤ C α ( ∥ f 0 ∥ L 2 ( X ) + γ R ∥ f 1 ∥ H R ( X ) ) . \left\|(S_{\gamma,s}-I)f\right\|_{L_2(X)}
\le C_{\alpha}\left(\|f_0\|_{L_2(X)}+\gamma^{R}\|f_1\|_{H^{R}(X)}\right). ∥ ( S γ , s − I ) f ∥ L 2 ( X ) ≤ C α ( ∥ f 0 ∥ L 2 ( X ) + γ R ∥ f 1 ∥ H R ( X ) ) .
Taking the infimum over all decompositions yields
∥ ( S γ , s − I ) f ∥ L 2 ( X ) ≤ C α K ( γ R , f ; L 2 ( X ) , H R ( X ) ) . \left\|(S_{\gamma,s}-I)f\right\|_{L_2(X)}
\le C_{\alpha}K(\gamma^{R},f;L_2(X),H^{R}(X)). ∥ ( S γ , s − I ) f ∥ L 2 ( X ) ≤ C α K ( γ R , f ; L 2 ( X ) , H R ( X )) .
Since R > α R>\alpha R > α , use (3) with m = R m=R m = R and the corresponding norm equivalence to obtain
K ( γ R , f ; L 2 ( X ) , H R ( X ) ) ≤ C α , X γ α ∥ f ∥ B 2 , ∞ α ( X ) . K(\gamma^{R},f;L_2(X),H^{R}(X))
\le C_{\alpha,X}\gamma^{\alpha}\|f\|_{B^{\alpha}_{2,\infty}(X)}. K ( γ R , f ; L 2 ( X ) , H R ( X )) ≤ C α , X γ α ∥ f ∥ B 2 , ∞ α ( X ) .
This proves (12) for s γ ≤ 1 s\gamma\le1 s γ ≤ 1 . For 1 / s < γ ≤ 1 1/s<\gamma\le1 1/ s < γ ≤ 1 , the same conclusion follows from (A.17), the continuous embedding B 2 , ∞ α ( X ) ↪ L 2 ( X ) B^{\alpha}_{2,\infty}(X)\hookrightarrow L_2(X) B 2 , ∞ α ( X ) ↪ L 2 ( X ) , and the fact that γ α \gamma^{\alpha} γ α is bounded below by s − α s^{-\alpha} s − α .
Proof of Theorem 4.4. For P X P_X P X -almost every x x x , define the conditional logistic risk by
ℓ x ( t ) : = η ( x ) log ( 1 + e − t ) + ( 1 − η ( x ) ) log ( 1 + e t ) . \ell_x(t):=\eta(x)\log(1+e^{-t})+(1-\eta(x))\log(1+e^t). ℓ x ( t ) := η ( x ) log ( 1 + e − t ) + ( 1 − η ( x )) log ( 1 + e t ) .
Since f ∗ ( x ) f^*(x) f ∗ ( x ) is its minimizer and is finite, we have
ℓ x ′ ( f ∗ ( x ) ) = 0 , ℓ x ′ ′ ( t ) = e t ( 1 + e t ) 2 ≤ 1 4 . \ell_x'(f^*(x))=0,\qquad\ell_x''(t)=\frac{e^t}{(1+e^t)^2}\le\frac{1}{4}. ℓ x ′ ( f ∗ ( x )) = 0 , ℓ x ′′ ( t ) = ( 1 + e t ) 2 e t ≤ 4 1 .
Taylor’s formula therefore gives a quadratic bound on the conditional excess risk. Combining it with the 1-Lipschitz property of the logistic loss, for every t ∈ R t\in\mathbb{R} t ∈ R and 1 ≤ r ≤ 2 1\le r\le2 1 ≤ r ≤ 2 we obtain
0 ≤ ℓ x ( t ) − ℓ x ( f ∗ ( x ) ) ≤ min { ∣ t − f ∗ ( x ) ∣ , 1 8 ∣ t − f ∗ ( x ) ∣ 2 } ≤ 8 1 − r ∣ t − f ∗ ( x ) ∣ r . 0\le\ell_x(t)-\ell_x(f^*(x))
\le\min\left\{|t-f^*(x)|,\frac{1}{8}|t-f^*(x)|^2\right\}
\le8^{1-r}|t-f^*(x)|^r. 0 ≤ ℓ x ( t ) − ℓ x ( f ∗ ( x )) ≤ min { ∣ t − f ∗ ( x ) ∣ , 8 1 ∣ t − f ∗ ( x ) ∣ 2 } ≤ 8 1 − r ∣ t − f ∗ ( x ) ∣ r .
Let p ′ p' p ′ be the conjugate exponent of p p p and take r : = 2 / p ′ = 2 ( 1 − 1 / p ) r:=2/p'=2(1-1/p) r := 2/ p ′ = 2 ( 1 − 1/ p ) . Then 1 ≤ p ′ ≤ 2 1\le p'\le2 1 ≤ p ′ ≤ 2 , 1 ≤ r ≤ 2 1\le r\le2 1 ≤ r ≤ 2 , and r p ′ = 2 rp'=2 r p ′ = 2 . For any measurable f f f with f − f ∗ ∈ L 2 ( X ) f-f^*\in L_2(X) f − f ∗ ∈ L 2 ( X ) , Hölder’s inequality gives
R ( f ) − R ∗ ≤ 8 1 − r ∫ X ∣ f − f ∗ ∣ r g d V ≤ 8 1 − r ∥ g ∥ L p ( X ) ( ∫ X ∣ f − f ∗ ∣ r p ′ d V ) 1 / p ′ = 8 1 − r ∥ g ∥ L p ( X ) ∥ f − f ∗ ∥ L 2 ( X ) . (42) \begin{aligned}
\mathcal{R}(f)-\mathcal{R}^*
&\le8^{1-r}\int_X |f-f^*|^r g\,dV \\
&\le8^{1-r}\|g\|_{L_p(X)}
\left(\int_X |f-f^*|^{rp'}\,dV\right)^{1/p'} \\
&=8^{1-r}\|g\|_{L_p(X)}\|f-f^*\|_{L_2(X)}. \tag*{(42)}
\end{aligned} R ( f ) − R ∗ ≤ 8 1 − r ∫ X ∣ f − f ∗ ∣ r g d V ≤ 8 1 − r ∥ g ∥ L p ( X ) ( ∫ X ∣ f − f ∗ ∣ r p ′ d V ) 1/ p ′ = 8 1 − r ∥ g ∥ L p ( X ) ∥ f − f ∗ ∥ L 2 ( X ) . (42)
Take S γ , s f ∗ ∈ H γ ( X ) S_{\gamma,s}f^*\in H_{\gamma}(X) S γ , s f ∗ ∈ H γ ( X ) , with s s s chosen from the smoothness α \alpha α as in Theorem 4.3. Applying that theorem to (A.18) yields
R ( S γ , s f ∗ ) − R ∗ ≤ 8 1 − r ∥ g ∥ L p ( X ) ∥ S γ , s f ∗ − f ∗ ∥ L 2 ( X ) ≤ C α , p , X , d γ 2 α ( 1 − 1 / p ) ∥ g ∥ L p ( X ) ∥ f ∗ ∥ B 2 , ∞ 2 ( 1 − 1 / p ) ( X ) . \begin{aligned}
\mathcal{R}(S_{\gamma,s}f^*)-\mathcal{R}^*
&\le8^{1-r}\|g\|_{L_p(X)}\|S_{\gamma,s}f^*-f^*\|_{L_2(X)} \\
&\le C_{\alpha,p,X,d}\gamma^{2\alpha(1-1/p)}
\|g\|_{L_p(X)}\|f^*\|_{B^{2(1-1/p)}_{2,\infty}(X)}.
\end{aligned} R ( S γ , s f ∗ ) − R ∗ ≤ 8 1 − r ∥ g ∥ L p ( X ) ∥ S γ , s f ∗ − f ∗ ∥ L 2 ( X ) ≤ C α , p , X , d γ 2 α ( 1 − 1/ p ) ∥ g ∥ L p ( X ) ∥ f ∗ ∥ B 2 , ∞ 2 ( 1 − 1/ p ) ( X ) .
Moreover, Theorem 4.1 gives
λ ∥ S γ , s f ∗ ∥ H γ ( X ) 2 ≤ C α , X , d λ γ − ϱ ∥ f ∗ ∥ L 2 ( X ) 2 . \lambda\|S_{\gamma,s}f^*\|_{H_{\gamma}(X)}^2
\le C_{\alpha,X,d}\lambda\gamma^{-\varrho}\|f^*\|_{L_2(X)}^2. λ ∥ S γ , s f ∗ ∥ H γ ( X ) 2 ≤ C α , X , d λ γ − ϱ ∥ f ∗ ∥ L 2 ( X ) 2 .
Substitution into (13) proves (14) .
Proof of Proposition 4.5. (i) Set θ : = α / β ∈ ( 0 , 1 ] \theta:= \alpha/\beta\in(0,1] θ := α / β ∈ ( 0 , 1 ] . For a , b > 0 a,b>0 a , b > 0 , the inequality log ( 1 + s ) ≤ C θ s θ \log(1+s) \le C_{\theta}s^{\theta} log ( 1 + s ) ≤ C θ s θ for s ≥ 0 s \ge0 s ≥ 0 gives
∣ log a − log b ∣ ≤ C θ ∣ a − b ∣ θ ( a − θ + b − θ ) . (43) |\log a-\log b| \le C_{\theta}|a-b|^{\theta}\left(a^{-\theta}+b^{-\theta}\right). \tag*{(43)} ∣ log a − log b ∣ ≤ C θ ∣ a − b ∣ θ ( a − θ + b − θ ) . (43)
Condition (15) implies that η \eta η takes values in ( 0 , 1 ) (0,1) ( 0 , 1 ) outside a set of volume zero. At every x ∈ X x\in X x ∈ X where 0 < η ( x ) < 1 0<\eta(x)<1 0 < η ( x ) < 1 ,
∣ f ∗ ( x ) ∣ ≤ − log ( η ( x ) ( 1 − η ( x ) ) ) ≤ C θ [ η ( x ) ( 1 − η ( x ) ) ] − θ , |f^{*}(x)|\le-\log(\eta(x)(1-\eta(x)))\le C_{\theta}[\eta(x)(1-\eta(x))]^{-\theta}, ∣ f ∗ ( x ) ∣ ≤ − log ( η ( x ) ( 1 − η ( x ))) ≤ C θ [ η ( x ) ( 1 − η ( x )) ] − θ ,
so f ∗ ∈ L 2 ( X ) f^{*}\in L_{2}(X) f ∗ ∈ L 2 ( X ) . Applying (A.19) to η \eta η and 1 − η 1-\eta 1 − η , and using the Hölder regularity of η \eta η , yields
∣ f ∗ ( x ) − f ∗ ( y ) ∣ ≤ C ∥ η ∥ C 0 , β ( X ) d X ( x , y ) α ( [ η ( x ) ( 1 − η ( x ) ) ] − θ + [ η ( y ) ( 1 − η ( y ) ) ] − θ ) . (44) \begin{aligned}
|f^{*}(x)-f^{*}(y)|
&\le C\|\eta\|_{C^{0,\beta}(X)}d_{X}(x,y)^{\alpha}
\bigl([\eta(x)(1-\eta(x))]^{-\theta}\\
&\qquad+[\eta(y)(1-\eta(y))]^{-\theta}\bigr). \tag*{(44)}
\end{aligned} ∣ f ∗ ( x ) − f ∗ ( y ) ∣ ≤ C ∥ η ∥ C 0 , β ( X ) d X ( x , y ) α ( [ η ( x ) ( 1 − η ( x )) ] − θ + [ η ( y ) ( 1 − η ( y )) ] − θ ) . (44)
Use the fixed finite atlas { ( U ν , ϕ ν ) } ν = 1 J \{(U_{\nu},\phi_{\nu})\}_{\nu=1}^{J} {( U ν , ϕ ν ) } ν = 1 J and partition of unity { χ ν } ν = 1 J \{\chi_{\nu}\}_{\nu=1}^{J} { χ ν } ν = 1 J fixed above. Write
F ν : = ( χ ν f ∗ ) ∘ ϕ ν − 1 , F_{\nu}:=(\chi_{\nu}f^{*})\circ\phi_{\nu}^{-1}, F ν := ( χ ν f ∗ ) ∘ ϕ ν − 1 ,
extended by zero to R ϱ \mathbb{R}^{\varrho} R ϱ . Fix a compact neighborhood of supp χ ν \operatorname{supp}\chi_{\nu} supp χ ν inside U ν U_{\nu} U ν . For sufficiently small ∥ h ∥ \|h\| ∥ h ∥ , whenever either cutoff value at u u u or u + h u+h u + h is nonzero, both corresponding manifold points lie in this neighborhood. There, geodesic distances are bounded by a constant times coordinate distances, and the volume density is bounded above and below, see (A.9). For such u u u , write x : = ϕ ν − 1 ( u ) x:=\phi_{\nu}^{-1}(u) x := ϕ ν − 1 ( u ) and y : = ϕ ν − 1 ( u + h ) y:=\phi_{\nu}^{-1}(u+h) y := ϕ ν − 1 ( u + h ) . Then
F ν ( u + h ) − F ν ( u ) = χ ν ( y ) f ∗ ( y ) − χ ν ( x ) f ∗ ( x ) = χ ν ( y ) ( f ∗ ( y ) − f ∗ ( x ) ) + ( χ ν ( y ) − χ ν ( x ) ) f ∗ ( x ) . \begin{aligned}
F_{\nu}(u+h)-F_{\nu}(u)
&=\chi_{\nu}(y)f^{*}(y)-\chi_{\nu}(x)f^{*}(x)\\
&=\chi_{\nu}(y)(f^{*}(y)-f^{*}(x))+(\chi_{\nu}(y)-\chi_{\nu}(x))f^{*}(x).
\end{aligned} F ν ( u + h ) − F ν ( u ) = χ ν ( y ) f ∗ ( y ) − χ ν ( x ) f ∗ ( x ) = χ ν ( y ) ( f ∗ ( y ) − f ∗ ( x )) + ( χ ν ( y ) − χ ν ( x )) f ∗ ( x ) .
We apply (A.20) to the first term and use ∣ χ ν ( y ) − χ ν ( x ) ∣ ≤ C ∥ h ∥ |\chi_{\nu}(y)-\chi_{\nu}(x)|\le C\|h\| ∣ χ ν ( y ) − χ ν ( x ) ∣ ≤ C ∥ h ∥ for the second. Taking L 2 L_{2} L 2 norms and using the coordinate change then gives
∥ F ν ( ⋅ + h ) − F ν ∥ L 2 ( R ϱ ) ≤ C ∥ η ∥ C 0 , β ( X ) θ ∥ h ∥ α ( ∫ X [ η ( 1 − η ) ] − 2 θ d V ) 1 2 + C ∥ h ∥ ∥ f ∗ ∥ L 2 ( X ) . \|F_{\nu}(\cdot+h)-F_{\nu}\|_{L_{2}(\mathbb{R}^{\varrho})}
\le C\|\eta\|_{C^{0,\beta}(X)}^{\theta}\|h\|^{\alpha}
\left(\int_{X}[\eta(1-\eta)]^{-2\theta}\,dV\right)^{\frac{1}{2}}
+C\|h\|\|f^{*}\|_{L_{2}(X)}. ∥ F ν ( ⋅ + h ) − F ν ∥ L 2 ( R ϱ ) ≤ C ∥ η ∥ C 0 , β ( X ) θ ∥ h ∥ α ( ∫ X [ η ( 1 − η ) ] − 2 θ d V ) 2 1 + C ∥ h ∥∥ f ∗ ∥ L 2 ( X ) .
Since 0 < α < 1 0<\alpha<1 0 < α < 1 , the second term is also bounded by a constant times ∥ h ∥ α ∥ f ∗ ∥ L 2 ( X ) \|h\|^{\alpha}\|f^{*}\|_{L_{2}(X)} ∥ h ∥ α ∥ f ∗ ∥ L 2 ( X ) for ∥ h ∥ ≤ 1 \|h\|\le1 ∥ h ∥ ≤ 1 . Larger shifts are controlled by 2 ∥ F ν ∥ L 2 ( R ϱ ) 2\|F_{\nu}\|_{L_{2}(\mathbb{R}^{\varrho})} 2∥ F ν ∥ L 2 ( R ϱ ) . The Euclidean first-difference characterization in [22 ] , Theorem 2.6.1 applies to each F ν F_{\nu} F ν because 0 < α < 1 0<\alpha<1 0 < α < 1 . On the compact manifold X X X , the corresponding chartwise norm is equivalent to the interpolation norm in (3) by interpolating the localization and reconstruction maps over a finite atlas, see [22 ] , proof of Theorem 7.2.3, Step 2, and Section 7.3.1. The preceding estimates therefore give (16) .
(ii) Under the assumptions of the second case, we first establish the gradient estimates
∣ ∇ η ( x ) ∣ 2 ≤ C η ( x ) , ∣ ∇ η ( x ) ∣ 2 ≤ C ( 1 − η ( x ) ) . (45) |\nabla\eta(x)|^{2}\le C\eta(x),\qquad|\nabla\eta(x)|^{2}\le C(1-\eta(x)). \tag*{(45)} ∣∇ η ( x ) ∣ 2 ≤ C η ( x ) , ∣∇ η ( x ) ∣ 2 ≤ C ( 1 − η ( x )) . (45)
To prove these estimates, let M : = 1 + ∥ ∇ 2 η ∥ L ∞ ( X ) M:=1+\|\nabla^{2}\eta\|_{L_{\infty}(X)} M := 1 + ∥ ∇ 2 η ∥ L ∞ ( X ) . The estimates hold where the gradient vanishes. At a point where ∇ η ( x ) ≠ 0 \nabla\eta(x)\ne0 ∇ η ( x ) = 0 , take v : = ∇ η ( x ) / ∣ ∇ η ( x ) ∣ v:=\nabla\eta(x)/|\nabla\eta(x)| v := ∇ η ( x ) /∣∇ η ( x ) ∣ . Compactness makes X X X geodesically complete, so Taylor’s formula along the geodesic t ↦ E x ( − t v ) t\mapsto\mathcal{E}_{x}(-tv) t ↦ E x ( − t v ) gives, for every t ≥ 0 t\ge0 t ≥ 0 ,
0 ≤ η ( E x ( − t v ) ) ≤ η ( x ) − t ∣ ∇ η ( x ) ∣ + M 2 t 2 . 0\le\eta(\mathcal{E}_{x}(-tv))\le\eta(x)-t|\nabla\eta(x)|+\frac{M}{2}t^{2}. 0 ≤ η ( E x ( − t v )) ≤ η ( x ) − t ∣∇ η ( x ) ∣ + 2 M t 2 .
Choosing t = ∣ ∇ η ( x ) ∣ / M t=|\nabla\eta(x)|/M t = ∣∇ η ( x ) ∣/ M yields ∣ ∇ η ( x ) ∣ 2 ≤ 2 M η ( x ) |\nabla\eta(x)|^{2}\le2M\eta(x) ∣∇ η ( x ) ∣ 2 ≤ 2 M η ( x ) . Applying the same argument to 1 − η 1-\eta 1 − η proves the second estimate in (A.21).
To justify differentiation despite possible endpoint values of η \eta η , define
f ε : = log ( η + ε ) − log ( 1 − η + ε ) , 0 < ε ≤ 1. f_{\varepsilon}:=\log(\eta+\varepsilon)-\log(1-\eta+\varepsilon),\qquad0<\varepsilon\le1. f ε := log ( η + ε ) − log ( 1 − η + ε ) , 0 < ε ≤ 1.
These functions belong to C m ( X ) C^{m}(X) C m ( X ) . Use the fixed finite atlas { ( U ν , ϕ ν ) } ν = 1 J \{(U_{\nu},\phi_{\nu})\}_{\nu=1}^{J} {( U ν , ϕ ν ) } ν = 1 J and partition of unity { χ ν } ν = 1 J \{\chi_{\nu}\}_{\nu=1}^{J} { χ ν } ν = 1 J . For each ν \nu ν , choose an open neighborhood V ν V_{\nu} V ν of supp χ ν \operatorname{supp}\chi_{\nu} supp χ ν such that V ν ⋐ U ν V_{\nu}\Subset U_{\nu} V ν ⋐ U ν . For u = ( u 1 , … , u ϱ ) ∈ ϕ ν ( V ν ) u=(u_{1},\ldots,u_{\varrho})\in\phi_{\nu}(V_{\nu}) u = ( u 1 , … , u ϱ ) ∈ ϕ ν ( V ν ) and a multi-index a = ( a 1 , … , a ϱ ) ∈ N 0 ϱ a=(a_{1},\ldots,a_{\varrho})\in\mathbb{N}_{0}^{\varrho} a = ( a 1 , … , a ϱ ) ∈ N 0 ϱ with 1 ≤ ∣ a ∣ ≤ m 1\le|a|\le m 1 ≤ ∣ a ∣ ≤ m , the chain rule expresses
∂ u a [ log ( η ∘ ϕ ν − 1 + ε ) ] ( u ) \partial_{u}^{a}\left[\log(\eta\circ\phi_{\nu}^{-1}+\varepsilon)\right](u) ∂ u a [ log ( η ∘ ϕ ν − 1 + ε ) ] ( u )
as a finite linear combination of terms of the form
( η ( ϕ ν − 1 ( u ) ) + ε ) − k ∏ i = 1 k [ ∂ u b i ( η ∘ ϕ ν − 1 ) ] ( u ) , b 1 + ⋯ + b k = a , b i ∈ N 0 d , ∣ b i ∣ ≥ 1. \left(\eta(\phi_\nu^{-1}(u))+\varepsilon\right)^{-k}\prod_{i=1}^{k}\left[\partial_u^{b_i}\left(\eta\circ\phi_\nu^{-1}\right)\right](u),\qquad b_1+\cdots+b_k=a,\quad b_i\in\mathbb{N}_0^d,\quad|b_i|\geq1. ( η ( ϕ ν − 1 ( u )) + ε ) − k i = 1 ∏ k [ ∂ u b i ( η ∘ ϕ ν − 1 ) ] ( u ) , b 1 + ⋯ + b k = a , b i ∈ N 0 d , ∣ b i ∣ ≥ 1.
If ℓ \ell ℓ of the k k k factors are exactly the first derivatives, then
∣ a ∣ = ∑ i = 1 k ∣ b i ∣ ≥ ℓ + 2 ( k − ℓ ) = 2 k − ℓ . |a|=\sum_{i=1}^{k}|b_i|\geq\ell+2(k-\ell)=2k-\ell. ∣ a ∣ = i = 1 ∑ k ∣ b i ∣ ≥ ℓ + 2 ( k − ℓ ) = 2 k − ℓ .
By (A.21) and the boundedness of the coordinate vector fields on V ν V_\nu V ν , each first coordinate derivative is bounded by C η ( ϕ ν − 1 ( u ) ) + ε C\sqrt{\eta(\phi_\nu^{-1}(u))+\varepsilon} C η ( ϕ ν − 1 ( u )) + ε , while the higher coordinate derivatives are uniformly bounded by C ∥ η ∥ C m C\|\eta\|_{C^m} C ∥ η ∥ C m , with constants depending on the fixed chart. Multiplying these bounds, we find that the ℓ \ell ℓ first derivative factors contribute a power ℓ / 2 \ell/2 ℓ /2 , while the bounds for the remaining factors are absorbed into C C C . Hence
∣ ( η ( ϕ ν − 1 ( u ) ) + ε ) − k ∏ i = 1 k [ ∂ u b i ( η ∘ ϕ ν − 1 ) ] ( u ) ∣ ≤ C ( η ( ϕ ν − 1 ( u ) ) + ε ) − k ( η ( ϕ ν − 1 ( u ) ) + ε ) ℓ / 2 = C ( η ( ϕ ν − 1 ( u ) ) + ε ) − k + ℓ / 2 = C ( η ( ϕ ν − 1 ( u ) ) + ε ) − ∣ a ∣ / 2 ( η ( ϕ ν − 1 ( u ) ) + ε ) ∣ a ∣ / 2 − k + ℓ / 2 . \begin{aligned}
\left|\left(\eta(\phi_\nu^{-1}(u))+\varepsilon\right)^{-k}\prod_{i=1}^{k}\left[\partial_u^{b_i}\left(\eta\circ\phi_\nu^{-1}\right)\right](u)\right|
\leq C\left(\eta(\phi_\nu^{-1}(u))+\varepsilon\right)^{-k}\left(\eta(\phi_\nu^{-1}(u))+\varepsilon\right)^{\ell/2} \\
&=C\left(\eta(\phi_\nu^{-1}(u))+\varepsilon\right)^{-k+\ell/2} \\
&=C\left(\eta(\phi_\nu^{-1}(u))+\varepsilon\right)^{-|a|/2}\left(\eta(\phi_\nu^{-1}(u))+\varepsilon\right)^{|a|/2-k+\ell/2}.
\end{aligned} ( η ( ϕ ν − 1 ( u )) + ε ) − k i = 1 ∏ k [ ∂ u b i ( η ∘ ϕ ν − 1 ) ] ( u ) ≤ C ( η ( ϕ ν − 1 ( u )) + ε ) − k ( η ( ϕ ν − 1 ( u )) + ε ) ℓ /2 = C ( η ( ϕ ν − 1 ( u )) + ε ) − k + ℓ /2 = C ( η ( ϕ ν − 1 ( u )) + ε ) − ∣ a ∣/2 ( η ( ϕ ν − 1 ( u )) + ε ) ∣ a ∣/2 − k + ℓ /2 .
Since k − ℓ / 2 ≤ ∣ a ∣ / 2 k-\ell/2\leq|a|/2 k − ℓ /2 ≤ ∣ a ∣/2 and 0 ≤ ℓ ≤ k 0\leq\ell\leq k 0 ≤ ℓ ≤ k , the exponent in the last factor satisfies
0 ≤ ∣ a ∣ 2 − k + ℓ 2 ≤ ∣ a ∣ 2 ≤ m 2 . 0\leq\frac{|a|}{2}-k+\frac{\ell}{2}\leq\frac{|a|}{2}\leq\frac{m}{2}. 0 ≤ 2 ∣ a ∣ − k + 2 ℓ ≤ 2 ∣ a ∣ ≤ 2 m .
Together with 0 < η ( ϕ ν − 1 ( u ) ) + ε ≤ 2 0<\eta(\phi_\nu^{-1}(u))+\varepsilon\leq2 0 < η ( ϕ ν − 1 ( u )) + ε ≤ 2 , this gives
( η ( ϕ ν − 1 ( u ) ) + ε ) ∣ a ∣ / 2 − k + ℓ / 2 ≤ 2 m / 2 . \left(\eta(\phi_\nu^{-1}(u))+\varepsilon\right)^{|a|/2-k+\ell/2}\leq2^{m/2}. ( η ( ϕ ν − 1 ( u )) + ε ) ∣ a ∣/2 − k + ℓ /2 ≤ 2 m /2 .
Absorbing 2 m / 2 2^{m/2} 2 m /2 into C C C , we obtain the bound C ( η ( ϕ ν − 1 ( u ) ) + ε ) − ∣ a ∣ / 2 C(\eta(\phi_\nu^{-1}(u))+\varepsilon)^{-|a|/2} C ( η ( ϕ ν − 1 ( u )) + ε ) − ∣ a ∣/2 for each term, uniformly in ε \varepsilon ε . Applying the same reasoning to ( 1 − η ) ∘ ϕ ν − 1 (1-\eta)\circ\phi_\nu^{-1} ( 1 − η ) ∘ ϕ ν − 1 gives
∣ ∂ u a ( f ε ∘ ϕ ν − 1 ) ( u ) ∣ ≤ C [ ( η ( ϕ ν − 1 ( u ) ) + ε ) − ∣ a ∣ / 2 + ( 1 − η ( ϕ ν − 1 ( u ) ) + ε ) − ∣ a ∣ / 2 ] ≤ C [ η ( ϕ ν − 1 ( u ) ) ( 1 − η ( ϕ ν − 1 ( u ) ) ) ] − ∣ a ∣ / 2 . \begin{aligned}
\left|\partial_u^a(f_\varepsilon\circ\phi_\nu^{-1})(u)\right|
\leq C\left[\left(\eta(\phi_\nu^{-1}(u))+\varepsilon\right)^{-|a|/2}+\left(1-\eta(\phi_\nu^{-1}(u))+\varepsilon\right)^{-|a|/2}\right] \\
&\leq C\left[\eta(\phi_\nu^{-1}(u))(1-\eta(\phi_\nu^{-1}(u)))\right]^{-|a|/2}.
\end{aligned} ∂ u a ( f ε ∘ ϕ ν − 1 ) ( u ) ≤ C [ ( η ( ϕ ν − 1 ( u )) + ε ) − ∣ a ∣/2 + ( 1 − η ( ϕ ν − 1 ( u )) + ε ) − ∣ a ∣/2 ] ≤ C [ η ( ϕ ν − 1 ( u )) ( 1 − η ( ϕ ν − 1 ( u ))) ] − ∣ a ∣/2 .
The constants are independent of ε \varepsilon ε and can be chosen uniformly over the finitely many charts. Since the sets V ν V_\nu V ν cover X X X , the relation between coordinate and covariant derivatives therefore yields
∣ ∇ ∣ a ∣ f ε ∣ ≤ C [ η ( 1 − η ) ] − ∣ a ∣ / 2 , 1 ≤ ∣ a ∣ ≤ m . (46) |\nabla^{|a|}f_\varepsilon|\leq C[\eta(1-\eta)]^{-|a|/2},\qquad1\leq|a|\leq m. \tag*{(46)} ∣ ∇ ∣ a ∣ f ε ∣ ≤ C [ η ( 1 − η ) ] − ∣ a ∣/2 , 1 ≤ ∣ a ∣ ≤ m . (46)
Also, adding the same positive number to η \eta η and 1 − η 1-\eta 1 − η moves their ratio towards one, so
∣ f ε ∣ ≤ ∣ f ∗ ∣ ≤ − log ( η ( 1 − η ) ) ≤ C m [ η ( 1 − η ) ] − m / 2 . |f_\varepsilon|\leq|f^*|\leq-\log(\eta(1-\eta))\leq C_m[\eta(1-\eta)]^{-m/2}. ∣ f ε ∣ ≤ ∣ f ∗ ∣ ≤ − log ( η ( 1 − η )) ≤ C m [ η ( 1 − η ) ] − m /2 .
Since η ( 1 − η ) ≤ 1 / 4 < 1 \eta(1-\eta)\leq1/4<1 η ( 1 − η ) ≤ 1/4 < 1 , (A.22) and (17) give the uniform bound
∥ f ε ∥ H m ( X ) ≤ C ( ∫ X [ η ( 1 − η ) ] − m d V ) 1 / 2 . \|f_\varepsilon\|_{H^m(X)}\leq C\left(\int_X[\eta(1-\eta)]^{-m}\,dV\right)^{1/2}. ∥ f ε ∥ H m ( X ) ≤ C ( ∫ X [ η ( 1 − η ) ] − m d V ) 1/2 .
Note that the endpoint set has zero volume. As ε \varepsilon ε decreases to zero, f ε f_\varepsilon f ε converges to f ∗ f^* f ∗ almost everywhere. By the dominated convergence theorem, f ε → f ∗ f_\varepsilon\to f^* f ε → f ∗ in L 2 ( X ) L_2(X) L 2 ( X ) . By weak compactness in H m ( X ) H^m(X) H m ( X ) and the strong convergence in L 2 ( X ) L_2(X) L 2 ( X ) , we obtain f ∗ ∈ H m ( X ) f^*\in H^m(X) f ∗ ∈ H m ( X ) . Weak lower semicontinuity of the H m ( X ) H^m(X) H m ( X ) norm then yields (18) . By [22 Theorem 7.4.5 and Theorem 7.3.1] , we have the identity H m ( X ) = B 2 , 2 m ( X ) H^m(X)=B^m_{2,2}(X) H m ( X ) = B 2 , 2 m ( X ) . The embeddings
B 2 , 2 m ( X ) ↪ B 2 , ∞ m ( X ) ↪ B 2 , ∞ α ( X ) , 0 < α ≤ m , B^m_{2,2}(X)\hookrightarrow B^m_{2,\infty}(X)\hookrightarrow B^\alpha_{2,\infty}(X),\qquad0<\alpha\leq m, B 2 , 2 m ( X ) ↪ B 2 , ∞ m ( X ) ↪ B 2 , ∞ α ( X ) , 0 < α ≤ m ,
by [22 Theorem 7.4.2] give the Besov conclusion. Finally, we verify the two claims concerning (19) . Put w ( x ) : = min { η ( x ) , 1 − η ( x ) } w(x) := \min\{\eta(x), 1-\eta(x)\} w ( x ) := min { η ( x ) , 1 − η ( x )} and let r > 0 r>0 r > 0 . The small value bound implies that w > 0 w>0 w > 0 almost everywhere. Since η ( 1 − η ) = w ( 1 − w ) ≥ w / 2 \eta(1-\eta)=w(1-w)\ge w/2 η ( 1 − η ) = w ( 1 − w ) ≥ w /2 , we have
∫ X [ η ( 1 − η ) ] − r d V ≤ 2 r ∫ X w − r d V = 2 r [ 2 r V ( X ) + r ∫ 0 1 / 2 u − r − 1 V { w < u } d u ] ≤ 2 r [ 2 r V ( X ) + r C 0 ∫ 0 1 / 2 u κ − r − 1 d u ] < ∞ . \int_X [\eta(1-\eta)]^{-r}\,\mathrm{d}V \le2^r\int_X w^{-r}\,\mathrm{d}V
=2^r\left[2^rV(X)+r\int_0^{1/2}u^{-r-1}V\{w<u\}\,\mathrm{d}u\right]
\le2^r\left[2^rV(X)+rC_0\int_0^{1/2}u^{\kappa-r-1}\,\mathrm{d}u\right]<\infty. ∫ X [ η ( 1 − η ) ] − r d V ≤ 2 r ∫ X w − r d V = 2 r [ 2 r V ( X ) + r ∫ 0 1/2 u − r − 1 V { w < u } d u ] ≤ 2 r [ 2 r V ( X ) + r C 0 ∫ 0 1/2 u κ − r − 1 d u ] < ∞.
whenever κ > r \kappa>r κ > r . Taking r = 2 α / β r=2\alpha/\beta r = 2 α / β in (i) and r = m r=m r = m in (ii) proves the respective assertions.
Proof of Proposition 4.6. Applying [28 Lemma 5.4] at the fixed smoothness order m m m gives (20) for sufficiently small γ > 0 \gamma>0 γ > 0 and γ m − 1 < δ < 1 / 2 \gamma^{m-1}<\delta<1/2 γ m − 1 < δ < 1/2 . Choose γ m ∈ ( 0 , 1 ) \gamma_m\in(0,1) γ m ∈ ( 0 , 1 ) so that this estimate holds for every 0 < γ ≤ γ m 0<\gamma\le\gamma_m 0 < γ ≤ γ m and γ m m − 1 < 1 / 4 \gamma_m^{m-1}<1/4 γ m m − 1 < 1/4 . For 1 / 4 ≤ δ ≤ 1 1/4\le\delta\le1 1/4 ≤ δ ≤ 1 , monotonicity of covering numbers and the bound at δ = 1 / 4 \delta=1/4 δ = 1/4 give the same estimate after increasing C m C_m C m .
For the polynomial estimate, observe that, for 0 < q ≤ 1 / 2 0<q\le1/2 0 < q ≤ 1/2 ,
δ 2 q ( log 2 δ ) ϱ + 1 = 2 2 q exp ( − 2 q log 2 δ ) ( log 2 δ ) ϱ + 1 ≤ 2 sup t > 0 exp ( − 2 q t ) t ϱ + 1 = 2 ( ϱ + 1 2 e q ) ϱ + 1 . \delta^{2q}\left(\log\frac{2}{\delta}\right)^{\varrho+1}
=2^{2q}\exp\left(-2q\log\frac{2}{\delta}\right)\left(\log\frac{2}{\delta}\right)^{\varrho+1}
\le2\sup_{t>0}\exp(-2qt)t^{\varrho+1}
=2\left(\frac{\varrho+1}{2eq}\right)^{\varrho+1}. δ 2 q ( log δ 2 ) ϱ + 1 = 2 2 q exp ( − 2 q log δ 2 ) ( log δ 2 ) ϱ + 1 ≤ 2 t > 0 sup exp ( − 2 q t ) t ϱ + 1 = 2 ( 2 e q ϱ + 1 ) ϱ + 1 .
Substitution into (20) proves (21) with an enlarged constant C m C_m C m .
Proof of Theorem 4.7. Fix an integer m ≥ 2 m\ge2 m ≥ 2 such that
m − 1 > β + ϱ 2 . m-1>\beta+\frac{\varrho}{2}. m − 1 > β + 2 ϱ .
Let C m C_m C m and γ m \gamma_m γ m be as in Proposition 4.6. For all sufficiently large n n n , we have 0 < γ n ≤ γ m 0<\gamma_n\le\gamma_m 0 < γ n ≤ γ m , 0 < λ n ≤ 1 0<\lambda_n\le1 0 < λ n ≤ 1 , and
q n : = log 2 2 log n ∈ ( 0 , 1 / 2 ] , b n : = 2 γ n m − 1 ≤ 1 4 , A n : = C m q n − ( ϱ + 1 ) γ n − ϱ ≥ 1. q_n:=\frac{\log2}{2\log n}\in(0,1/2],\qquad b_n:=2\gamma_n^{m-1}\le\frac14,\qquad A_n:=C_mq_n^{-(\varrho+1)}\gamma_n^{-\varrho}\ge1. q n := 2 log n log 2 ∈ ( 0 , 1/2 ] , b n := 2 γ n m − 1 ≤ 4 1 , A n := C m q n − ( ϱ + 1 ) γ n − ϱ ≥ 1.
The estimate (21) verifies (7) for F = H γ n ( X ) \mathcal{F}=\mathcal{H}_{\gamma_n}(X) F = H γ n ( X ) with entropy exponent q n q_n q n , A = A n A=A_n A = A n , and b = b n b=b_n b = b n . Apply Corollary 3.3 and Theorem 4.4. Then with probability at least 1 − e − t 1-e^{-t} 1 − e − t , the excess risk is bounded, up to constants independent of n n n , by
γ n β + λ n γ n − ϱ + λ n − 1 ( A n n ) 2 / ( 1 + q n ) + λ n + λ n − 1 b n 2 + λ n − 1 t 2 n 2 . \gamma_n^\beta+\lambda_n\gamma_n^{-\varrho}+\lambda_n^{-1}\left(\frac{A_n}{n}\right)^{2/(1+q_n)}+\lambda_n+\lambda_n^{-1}b_n^2+\lambda_n^{-1}\frac{t^2}{n^2}. γ n β + λ n γ n − ϱ + λ n − 1 ( n A n ) 2/ ( 1 + q n ) + λ n + λ n − 1 b n 2 + λ n − 1 n 2 t 2 .
The choice of q n q_n q n gives
A n n ≍ ( log n ) ϱ + 1 n γ n − ϱ = ( ( log n ) ϱ + 1 n ) ( 2 β + ϱ ) / ( 2 β + 3 ϱ ) , \frac{A_n}{n}\asymp\frac{(\log n)^{\varrho+1}}{n}\gamma_n^{-\varrho}
=\left(\frac{(\log n)^{\varrho+1}}{n}\right)^{(2\beta+\varrho)/(2\beta+3\varrho)}, n A n ≍ n ( log n ) ϱ + 1 γ n − ϱ = ( n ( log n ) ϱ + 1 ) ( 2 β + ϱ ) / ( 2 β + 3 ϱ ) ,
( A n n ) 2 / ( 1 + q n ) = ( A n n ) 2 exp ( 2 q n 1 + q n log n A n ) ≍ ( A n n ) 2 . \left(\frac{A_n}{n}\right)^{2/(1+q_n)}
=\left(\frac{A_n}{n}\right)^2\exp\left(\frac{2q_n}{1+q_n}\log\frac{n}{A_n}\right)
\asymp\left(\frac{A_n}{n}\right)^2. ( n A n ) 2/ ( 1 + q n ) = ( n A n ) 2 exp ( 1 + q n 2 q n log A n n ) ≍ ( n A n ) 2 .
It follows that the first three terms have the same order:
γ n β ≍ λ n γ n − ϱ ≍ λ n − 1 ( A n n ) 2 / ( 1 + q n ) ≍ ( ( log n ) ϱ + 1 n ) 2 β / ( 2 β + 3 ϱ ) . \gamma_n^\beta\asymp\lambda_n\gamma_n^{-\varrho}\asymp\lambda_n^{-1}\left(\frac{A_n}{n}\right)^{2/(1+q_n)}
\asymp\left(\frac{(\log n)^{\varrho+1}}{n}\right)^{2\beta/(2\beta+3\varrho)}. γ n β ≍ λ n γ n − ϱ ≍ λ n − 1 ( n A n ) 2/ ( 1 + q n ) ≍ ( n ( log n ) ϱ + 1 ) 2 β / ( 2 β + 3 ϱ ) .
To control the error from the finite covering scale, note that λ n = γ n β + ϱ \lambda_n=\gamma_n^{\beta+\varrho} λ n = γ n β + ϱ and hence
λ n − 1 b n 2 γ n β = 4 γ n 2 ( m − 1 ) − 2 β − ϱ → 0 , \frac{\lambda_n^{-1}b_n^2}{\gamma_n^\beta}=4\gamma_n^{2(m-1)-2\beta-\varrho}\to0, γ n β λ n − 1 b n 2 = 4 γ n 2 ( m − 1 ) − 2 β − ϱ → 0 ,
by the choice of m m m . The remaining terms satisfy
λ n γ n β = γ n ϱ → 0 , \frac{\lambda_n}{\gamma_n^\beta}=\gamma_n^\varrho\to0, γ n β λ n = γ n ϱ → 0 ,
and
λ n − 1 t 2 n 2 ≤ t 2 ( ( log n ) ϱ + 1 n ) ( 2 β + 4 ϱ ) / ( 2 β + 3 ϱ ) = o ( ( ( log n ) ϱ + 1 n ) 2 β / ( 2 β + 3 ϱ ) ) \lambda_n^{-1}\frac{t^2}{n^2}\le t^2\left(\frac{(\log n)^{\varrho+1}}{n}\right)^{(2\beta+4\varrho)/(2\beta+3\varrho)}
=o\left(\left(\frac{(\log n)^{\varrho+1}}{n}\right)^{2\beta/(2\beta+3\varrho)}\right) λ n − 1 n 2 t 2 ≤ t 2 ( n ( log n ) ϱ + 1 ) ( 2 β + 4 ϱ ) / ( 2 β + 3 ϱ ) = o ( ( n ( log n ) ϱ + 1 ) 2 β / ( 2 β + 3 ϱ ) )
for fixed t t t . Since 2 β = 4 α ( 1 − 1 / p ) 2\beta=4\alpha(1-1/p) 2 β = 4 α ( 1 − 1/ p ) , this proves (22) .
Proof of Theorem 4.8 . We first establish a comparison bound for clipped predictions. For every u ∈ R u\in\mathbb{R} u ∈ R , y ∈ { − 1 , 1 } y\in\{-1,1\} y ∈ { − 1 , 1 } , and M > 0 M>0 M > 0 ,
ϕ ( y π M ( u ) ) ≤ ϕ ( y u ) + e − M . (47) \phi(y\pi_M(u))\le\phi(yu)+e^{-M}. \tag*{(47)} ϕ ( y π M ( u )) ≤ ϕ ( y u ) + e − M . (47)
Indeed, clipping decreases the loss when y u < − M yu<-M y u < − M and leaves it unchanged when ∣ u ∣ ≤ M |u|\le M ∣ u ∣ ≤ M . When y u > M yu>M y u > M , the increase is at most ϕ ( M ) ≤ e − M \phi(M)\le e^{-M} ϕ ( M ) ≤ e − M .
Fix 0 < λ ≤ 1 0<\lambda\le1 0 < λ ≤ 1 , γ > 0 \gamma>0 γ > 0 , M ≥ 1 M\ge1 M ≥ 1 , and a deterministic comparison function f 0 ∈ H γ ( X ) f_0\in\mathcal{H}_\gamma(X) f 0 ∈ H γ ( X ) . By the 1-Lipschitz property of ϕ \phi ϕ ,
0 ≤ ( ϕ ( Y f 0 ( X ) ) − ϕ ( Y π M f 0 ( X ) ) ) + ≤ ∥ f 0 ∥ ∞ . 0\le\left(\phi(Yf_0(X))-\phi(Y\pi_Mf_0(X))\right)_+\le\|f_0\|_\infty. 0 ≤ ( ϕ ( Y f 0 ( X )) − ϕ ( Y π M f 0 ( X )) ) + ≤ ∥ f 0 ∥ ∞ .
Moreover, (A.23) bounds the negative part of this loss difference by e − M e^{-M} e − M . Hence
E P ( ϕ ( Y f 0 ( X ) ) − ϕ ( Y π M f 0 ( X ) ) ) + ≤ R ( f 0 ) − R ∗ + e − M , \mathbb{E}_{P}\left(\phi(Yf_0(X))-\phi(Y\pi_Mf_0(X))\right)_+
\le\mathcal{R}(f_0)-\mathcal{R}^*+e^{-M}, E P ( ϕ ( Y f 0 ( X )) − ϕ ( Y π M f 0 ( X )) ) + ≤ R ( f 0 ) − R ∗ + e − M ,
E P ( ϕ ( Y f 0 ( X ) ) − ϕ ( Y π M f 0 ( X ) ) ) 2 ≤ ∥ f 0 ∥ ∞ ( R ( f 0 ) − R ∗ + e − M ) . \mathbb{E}_{P}\left(\phi(Yf_0(X))-\phi(Y\pi_Mf_0(X))\right)^2
\le\|f_0\|_\infty\left(\mathcal{R}(f_0)-\mathcal{R}^*+e^{-M}\right). E P ( ϕ ( Y f 0 ( X )) − ϕ ( Y π M f 0 ( X )) ) 2 ≤ ∥ f 0 ∥ ∞ ( R ( f 0 ) − R ∗ + e − M ) .
Bernstein’s inequality shows that, with probability at least 1 − e − t / 2 1-e^{-t/2} 1 − e − t /2 ,
E D ( ϕ ( Y f 0 ( X ) ) − ϕ ( Y π M f 0 ( X ) ) ) + ≤ 2 ( R ( f 0 ) − R ∗ + e − M ) + C ∥ f 0 ∥ ∞ ( t + log 2 ) n . (48) \mathbb{E}_{D}\left(\phi(Yf_0(X))-\phi(Y\pi_Mf_0(X))\right)_+
\le2\left(\mathcal{R}(f_0)-\mathcal{R}^*+e^{-M}\right)
+\frac{C\|f_0\|_\infty(t+\log2)}{n}. \tag*{(48)} E D ( ϕ ( Y f 0 ( X )) − ϕ ( Y π M f 0 ( X )) ) + ≤ 2 ( R ( f 0 ) − R ∗ + e − M ) + n C ∥ f 0 ∥ ∞ ( t + log 2 ) . (48)
(A.23) implies
E D ( ϕ ( Y π M f ^ n , λ , γ ( X ) ) − ϕ ( Y π M f 0 ( X ) ) ) ≤ E D ( ϕ ( Y f ^ n , λ , γ ( X ) ) − ϕ ( Y f 0 ( X ) ) ) + E D ( ϕ ( Y f 0 ( X ) ) − ϕ ( Y π M f 0 ( X ) ) ) + e − M ≤ λ ∥ f 0 ∥ H γ 2 + E D ( ϕ ( Y f 0 ( X ) ) − ϕ ( Y π M f 0 ( X ) ) ) + + e − M . (49) \begin{aligned}
\mathbb{E}_{D}\left(\phi(Y\pi_M\widehat{f}_{n,\lambda,\gamma}(X))-\phi(Y\pi_Mf_0(X))\right)
&\le\mathbb{E}_{D}\left(\phi(Y\widehat{f}_{n,\lambda,\gamma}(X))-\phi(Yf_0(X))\right)\\
&\quad+\mathbb{E}_{D}\left(\phi(Yf_0(X))-\phi(Y\pi_Mf_0(X))\right)+e^{-M}\\
&\le\lambda\|f_0\|_{\mathcal{H}_\gamma}^2+\mathbb{E}_{D}\left(\phi(Yf_0(X))-\phi(Y\pi_Mf_0(X))\right)_++e^{-M}. \tag*{(49)}
\end{aligned} E D ( ϕ ( Y π M f n , λ , γ ( X )) − ϕ ( Y π M f 0 ( X )) ) ≤ E D ( ϕ ( Y f n , λ , γ ( X )) − ϕ ( Y f 0 ( X )) ) + E D ( ϕ ( Y f 0 ( X )) − ϕ ( Y π M f 0 ( X )) ) + e − M ≤ λ ∥ f 0 ∥ H γ 2 + E D ( ϕ ( Y f 0 ( X )) − ϕ ( Y π M f 0 ( X )) ) + + e − M . (49)
Comparison of the empirical objective with zero gives ∥ f ^ n , λ , γ ∥ H γ ≤ λ − 1 / 2 \|\widehat{f}_{n,\lambda,\gamma}\|_{\mathcal{H}_\gamma}\le\lambda^{-1/2} ∥ f n , λ , γ ∥ H γ ≤ λ − 1/2 . Put
V = { π M f : f ∈ λ − 1 / 2 B H γ } . \mathcal{V}=\{\pi_Mf:f\in\lambda^{-1/2}B_{\mathcal{H}_\gamma}\}. V = { π M f : f ∈ λ − 1/2 B H γ } .
Suppose that V \mathcal{V} V has a finite δ \delta δ -net { v 1 , … , v N } ⊂ V \{v_1,\ldots,v_N\}\subset\mathcal{V} { v 1 , … , v N } ⊂ V in the supremum norm. For each net member, both v j v_j v j and π M f 0 \pi_Mf_0 π M f 0 are bounded by M M M , so their loss difference is bounded in absolute value by 2 M 2M 2 M . By Lemma 3.1,
E P ( ϕ ( Y v j ( X ) ) − ϕ ( Y π M f 0 ( X ) ) ) 2 ≤ 2 E P ( ϕ ( Y v j ( X ) ) − ψ ( X , Y ) ) 2 + 2 E P ( ϕ ( Y π M f 0 ( X ) ) − ψ ( X , Y ) ) 2 \mathbb{E}_{P}\left(\phi(Yv_j(X))-\phi(Y\pi_Mf_0(X))\right)^2
\le2\mathbb{E}_{P}\left(\phi(Yv_j(X))-\psi(X,Y)\right)^2+2\mathbb{E}_{P}\left(\phi(Y\pi_Mf_0(X))-\psi(X,Y)\right)^2 E P ( ϕ ( Y v j ( X )) − ϕ ( Y π M f 0 ( X )) ) 2 ≤ 2 E P ( ϕ ( Y v j ( X )) − ψ ( X , Y ) ) 2 + 2 E P ( ϕ ( Y π M f 0 ( X )) − ψ ( X , Y ) ) 2
≤ 2 ( M + 2 ) ( R ( v j ) − R ∗ + R ( π M f 0 ) − R ∗ ) . \le2(M+2)\left(\mathcal{R}(v_j)-\mathcal{R}^*+\mathcal{R}(\pi_Mf_0)-\mathcal{R}^*\right). ≤ 2 ( M + 2 ) ( R ( v j ) − R ∗ + R ( π M f 0 ) − R ∗ ) .
Apply Bernstein’s inequality to these N N N loss differences and use a union bound with parameter log N + t + log 2 \log N+t+\log2 log N + t + log 2 . Then, with probability at least 1 − e − t / 2 1-e^{-t/2} 1 − e − t /2 , simultaneously for every j j j ,
R ( v j ) − R ( π M f 0 ) ≤ E D ( ϕ ( Y v j ( X ) ) − ϕ ( Y π M f 0 ( X ) ) ) + 1 2 ( R ( v j ) − R ∗ + R ( π M f 0 ) − R ∗ ) + C ( M + 2 ) ( log N + t + log 2 ) n . \begin{aligned}
\mathcal{R}(v_j)-\mathcal{R}(\pi_Mf_0)
&\le\mathbb{E}_{D}\left(\phi(Yv_j(X))-\phi(Y\pi_Mf_0(X))\right)\\
&\quad+\frac{1}{2}\left(\mathcal{R}(v_j)-\mathcal{R}^*+\mathcal{R}(\pi_Mf_0)-\mathcal{R}^*\right)
+\frac{C(M+2)(\log N+t+\log2)}{n}.
\end{aligned} R ( v j ) − R ( π M f 0 ) ≤ E D ( ϕ ( Y v j ( X )) − ϕ ( Y π M f 0 ( X )) ) + 2 1 ( R ( v j ) − R ∗ + R ( π M f 0 ) − R ∗ ) + n C ( M + 2 ) ( log N + t + log 2 ) .
For each w ∈ V w \in\mathcal{V} w ∈ V , choose a net member with ∥ w − v j ∥ ∞ ≤ δ \lVert w-v_j\rVert_\infty\le\delta ∥ w − v j ∥ ∞ ≤ δ . The Lipschitz property of the loss controls both the population and empirical loss differences by δ \delta δ . Rearranging the preceding estimate therefore yields, on the same event,
R ( w ) − R ∗ ≤ 2 E D ( ϕ ( Y w ( X ) ) − ϕ ( Y π M f 0 ( X ) ) ) + 3 ( R ( π M f 0 ) − R ∗ ) + C ( M + 2 ) ( log N + t + log 2 ) n + 3 δ . \mathcal{R}(w)-\mathcal{R}^* \le2\mathbb{E}_D\bigl(\phi(Yw(X))-\phi(Y\pi_M f_0(X))\bigr)+3\bigl(\mathcal{R}(\pi_M f_0)-\mathcal{R}^*\bigr)
+\frac{C(M+2)(\log N+t+\log2)}{n}+3\delta. R ( w ) − R ∗ ≤ 2 E D ( ϕ ( Y w ( X )) − ϕ ( Y π M f 0 ( X )) ) + 3 ( R ( π M f 0 ) − R ∗ ) + n C ( M + 2 ) ( log N + t + log 2 ) + 3 δ .
In particular, this holds for w = π M f ^ n , λ , γ w=\pi_M\widehat{f}_{n,\lambda,\gamma} w = π M f n , λ , γ . Since R ( π M f 0 ) − R ∗ ≤ R ( f 0 ) − R ∗ + e − M \mathcal{R}(\pi_M f_0)-\mathcal{R}^* \le\mathcal{R}(f_0)-\mathcal{R}^*+e^{-M} R ( π M f 0 ) − R ∗ ≤ R ( f 0 ) − R ∗ + e − M by (A.23), combining the two events with (A.24) and (A.25) proves that, with probability at least 1 − e − t 1-e^{-t} 1 − e − t ,
R ( π M f ^ n , λ , γ ) − R ∗ ≤ C [ R ( f 0 ) − R ∗ + λ ∥ f 0 ∥ H γ 2 + ( M + 2 ) ( log N + t + log 2 ) n + ∥ f 0 ∥ ∞ ( t + log 2 ) n + e − M + δ ] . (50) \mathcal{R}(\pi_M\widehat{f}_{n,\lambda,\gamma})-\mathcal{R}^*
\le C\left[\mathcal{R}(f_0)-\mathcal{R}^*+\lambda\lVert f_0\rVert_{\mathcal{H}_\gamma}^2+\frac{(M+2)(\log N+t+\log2)}{n}
+\frac{\lVert f_0\rVert_\infty(t+\log2)}{n}+e^{-M}+\delta\right]. \tag*{(50)} R ( π M f n , λ , γ ) − R ∗ ≤ C [ R ( f 0 ) − R ∗ + λ ∥ f 0 ∥ H γ 2 + n ( M + 2 ) ( log N + t + log 2 ) + n ∥ f 0 ∥ ∞ ( t + log 2 ) + e − M + δ ] . (50)
We now apply this bound with the parameters in the theorem and choose
f 0 = S γ n , s f ∗ , s = ⌊ α ⌋ + 1 , δ = γ n β . f_0=S_{\gamma_n,s}f^*,\qquad s=\lfloor\alpha\rfloor+1,\qquad\delta=\gamma_n^\beta. f 0 = S γ n , s f ∗ , s = ⌊ α ⌋ + 1 , δ = γ n β .
Theorem 4.1 and the proof of Theorem 4.4 give
R ( f 0 ) − R ∗ ≤ C γ n β , ∥ f 0 ∥ H γ n 2 ≤ C γ n − ϱ , ∥ f 0 ∥ ∞ ≤ ∥ f 0 ∥ H γ n ≤ C γ n − ϱ / 2 . (51) \mathcal{R}(f_0)-\mathcal{R}^*\le C\gamma_n^\beta,\qquad\lVert f_0\rVert_{\mathcal{H}_{\gamma_n}}^2\le C\gamma_n^{-\varrho},\qquad\lVert f_0\rVert_\infty\le\lVert f_0\rVert_{\mathcal{H}_{\gamma_n}}\le C\gamma_n^{-\varrho/2}. \tag*{(51)} R ( f 0 ) − R ∗ ≤ C γ n β , ∥ f 0 ∥ H γ n 2 ≤ C γ n − ϱ , ∥ f 0 ∥ ∞ ≤ ∥ f 0 ∥ H γ n ≤ C γ n − ϱ /2 . (51)
Fix an integer m ≥ 2 m\ge2 m ≥ 2 such that
m − 1 > 3 β + ϱ 2 . m-1>\frac{3\beta+\varrho}{2}. m − 1 > 2 3 β + ϱ .
Since λ n = γ n β + ϱ \lambda_n=\gamma_n^{\beta+\varrho} λ n = γ n β + ϱ , the normalized covering radius satisfies
δ λ n 2 = 1 2 γ n ( 3 β + ϱ ) / 2 > γ n m − 1 \frac{\delta\sqrt{\lambda_n}}{2}=\frac{1}{2}\gamma_n^{(3\beta+\varrho)/2}>\gamma_n^{m-1} 2 δ λ n = 2 1 γ n ( 3 β + ϱ ) /2 > γ n m − 1
for all sufficiently large n n n . It is also at most one, and γ n ≤ γ m \gamma_n\le\gamma_m γ n ≤ γ m eventually. Thus Proposition 4.6 applies at this radius. Clipping does not increase supremum distances. Covering λ n − 1 / 2 B H γ n \lambda_n^{-1/2}B_{\mathcal{H}_{\gamma_n}} λ n − 1/2 B H γ n at radius δ / 2 \delta/2 δ /2 , clipping the centers, and choosing one member of V \mathcal{V} V from each nonempty covering ball produces a δ \delta δ -net in V \mathcal{V} V , with
log N ≤ C m γ n − ϱ ( log 4 δ λ n ) ϱ + 1 ≤ C γ n − ϱ ( log n ) ϱ + 1 . \log N\le C_m\gamma_n^{-\varrho}\left(\log\frac{4}{\delta\sqrt{\lambda_n}}\right)^{\varrho+1}\le C\gamma_n^{-\varrho}(\log n)^{\varrho+1}. log N ≤ C m γ n − ϱ ( log δ λ n 4 ) ϱ + 1 ≤ C γ n − ϱ ( log n ) ϱ + 1 .
The constants are independent of n n n .
Substituting this estimate and (A.27) into (A.26), and using M n = 3 log n M_n=3\log n M n = 3 log n , gives
R ( π M n f ^ n , λ n , γ n ) − R ∗ ≤ C [ γ n β + λ n γ n − ϱ + γ n − ϱ ( log n ) ϱ + 2 n + ( log n + γ n − ϱ / 2 ) ( 1 + t ) n + n − 3 ] . \mathcal{R}(\pi_{M_n}\widehat{f}_{n,\lambda_n,\gamma_n})-\mathcal{R}^*
\le C\left[\gamma_n^\beta+\lambda_n\gamma_n^{-\varrho}+\frac{\gamma_n^{-\varrho}(\log n)^{\varrho+2}}{n}+\frac{(\log n+\gamma_n^{-\varrho/2})(1+t)}{n}+n^{-3}\right]. R ( π M n f n , λ n , γ n ) − R ∗ ≤ C [ γ n β + λ n γ n − ϱ + n γ n − ϱ ( log n ) ϱ + 2 + n ( log n + γ n − ϱ /2 ) ( 1 + t ) + n − 3 ] .
The parameter choices imply that the main terms are
γ n β = λ n γ n − ϱ = γ n − ϱ ( log n ) ϱ + 2 n = ( ( log n ) ϱ + 2 n ) β / ( β + ϱ ) . \gamma_n^\beta=\lambda_n\gamma_n^{-\varrho}=\frac{\gamma_n^{-\varrho}(\log n)^{\varrho+2}}{n}=\left(\frac{(\log n)^{\varrho+2}}{n}\right)^{\beta/(\beta+\varrho)}. γ n β = λ n γ n − ϱ = n γ n − ϱ ( log n ) ϱ + 2 = ( n ( log n ) ϱ + 2 ) β / ( β + ϱ ) .
The remaining terms are of smaller order for fixed t t t . Absorbing the factor 1 + t 1+t 1 + t into C C C for fixed t t t proves the logistic bound in (24) . For the misclassification bound, note that
sgn ( π M n f ^ n , λ n , γ n ) = sgn ( f ^ n , λ n , γ n ) . \operatorname{sgn}(\pi_{M_n}\widehat{f}_{n,\lambda_n,\gamma_n})=\operatorname{sgn}(\widehat{f}_{n,\lambda_n,\gamma_n}). sgn ( π M n f n , λ n , γ n ) = sgn ( f n , λ n , γ n ) .
The calibration inequality (2) then gives the misclassification bound on the same event.