Smooth function approximation by deep neural networks with general activation functions

TL;DR

This paper analyzes the approximation capabilities of deep neural networks using a broad class of activation functions, establishing their efficiency in approximating smooth functions and their minimax optimality in regression and classification tasks.

Contribution

It extends the theoretical understanding of neural network expressivity to general activation functions beyond ReLU and sigmoid.

Findings

Derived depth, width, and sparsity requirements for approximation

Proved minimax optimality of estimators with general activation functions

Applicable to both regression and classification problems

Abstract

There has been a growing interest in expressivity of deep neural networks. However, most of the existing work about this topic focuses only on the specific activation function such as ReLU or sigmoid. In this paper, we investigate the approximation ability of deep neural networks with a broad class of activation functions. This class of activation functions includes most of frequently used activation functions. We derive the required depth, width and sparsity of a deep neural network to approximate any H\"older smooth function upto a given approximation error for the large class of activation functions. Based on our approximation error analysis, we derive the minimax optimality of the deep neural network estimators with the general activation functions in both regression and classification problems.

Equations

$$N_{\sigma }(\mathbf{x}|\mathbf{\theta })\text{mathchar}58={\mathsf{A}}_{L+1}\circ {\sigma }_L\circ {\mathsf{A}}_L\circ \cdots \circ {\sigma }_1\circ {\mathsf{A}}_1(\mathbf{x}),$$

$${\left|\mathbf{\theta }\right|}_0\text{mathchar}58=\sum _{l=1}^{L+1}\left({\left|\text{vec}({\mathbf{W}}_l)\right|}_0+{\left|{\mathbf{b}}_l\right|}_0\right),$$

$${\left|\mathbf{\theta }\right|}_{\infty }\text{mathchar}58=\max \left\{\underset{1\le l\le L+1}{\max }{\left|\text{vec}({\mathbf{W}}_l)\right|}_{\infty },\underset{1\le l\le L+1}{\max }{\left|{\mathbf{b}}_l\right|}_{\infty }\right\}.$$

$${\Theta }_{d,o}(L,N)\text{mathchar}58=\left\{\mathbf{\theta }\text{mathchar}58L(\mathbf{\theta })\le L,n_{\max }(\mathbf{\theta })\le N,\text{in}(\mathbf{\theta })=d,\text{out}(\mathbf{\theta })=o\right\}.$$

$${\Theta }_{d,o}(L,N,S,B)\text{mathchar}58=\left\{\mathbf{\theta }\in {\Theta }_{d,o}(L,N)\text{mathchar}58{\left|\mathbf{\theta }\right|}_0\le S,{\left|\mathbf{\theta }\right|}_{\infty }\le B\right\}.$$

$${\partial }^{\mathbf{m}}f\text{mathchar}58=\frac{{\partial }^{|\mathbf{m}|}f}{\partial {\mathbf{x}}^{\mathbf{m}}}=\frac{{\partial }^{|\mathbf{m}|}f}{\partial x_1^{m_1}\cdots \partial x_d^{m_d}},$$

$$[f{]}_s\text{mathchar}58=\underset{{\mathbf{x}}_1,{\mathbf{x}}_2\in \mathcal{X},{\mathbf{x}}_1\ne {\mathbf{x}}_2}{\sup }\frac{|f({\mathbf{x}}_1)-f({\mathbf{x}}_2)|}{|{\mathbf{x}}_1-{\mathbf{x}}_2{|}^s}.$$

$${\mathcal{H}}^{\alpha }(\mathcal{X})\text{mathchar}58=\left\{f\in {\mathcal{C}}^{\lfloor \alpha \rfloor }(\mathcal{X})\text{mathchar}58\Vert f{\Vert }_{{\mathcal{H}}^{\alpha }(\mathcal{X})}<\infty \right\},$$

$$\Vert f{\Vert }_{{\mathcal{H}}^{\alpha }(\mathcal{X})}\text{mathchar}58=\sum _{\mathbf{m}\in \text{mathbbm}{N}_0^d\text{mathchar}58|\mathbf{m}|\le \lfloor \alpha \rfloor }\Vert {\partial }^{\mathbf{m}}f{\Vert }_{\infty }+\sum _{\mathbf{m}\in \text{mathbbm}{N}_0^d\text{mathchar}58|\mathbf{m}|=\lfloor \alpha \rfloor }[{\partial }^{\mathbf{m}}f{]}_{\alpha -\lfloor \alpha \rfloor }.$$

$${\mathcal{H}}^{\alpha ,R}(\mathcal{X})\text{mathchar}58=\left\{f\in {\mathcal{H}}^{\alpha }(\mathcal{X})\text{mathchar}58\Vert f{\Vert }_{{\mathcal{H}}^{\alpha }(\mathcal{X})}\le R\right\}.$$

$${\mathbf{\theta }}_{ϵ}\in {\Theta }_{d,1}\left(L_0\log (1/ϵ),N_0{ϵ}^{-d/\alpha },S_0{ϵ}^{-d/\alpha }\log (1/ϵ),B_0{ϵ}^{-4(d/\alpha +1)}\right)$$

$$\underset{\mathbf{x}\in [0,1{]}^d}{\sup }\left|f(\mathbf{x})-N_{\sigma }(\mathbf{x}|{\mathbf{\theta }}_{ϵ})\right|\le ϵ.$$

$$\text{mathbbm}{G}_{d,M}\text{mathchar}58=\left\{\frac{1}{M}(m_1,\dots ,m_d)\text{mathchar}58m_j\in \{0,1,\dots ,M\},j=1,\dots ,d\right\}.$$

$$\underset{x\in [0,1{]}^d}{\sup }\left|f(\mathbf{x})-\sum _{\mathbf{z}\in \text{mathbbm}{G}_{d,M}}{g}_{\mathbf{z},M}(\mathbf{x})\right|\le C|\text{mathbbm}{G}_{d,M}{|}^{-\alpha /d},$$

$$P_M(\mathbf{x})\text{mathchar}58=\sum _{\mathbf{z}\in \text{mathbbm}{G}_{d,M}}{g}_{\mathbf{z},M}(\mathbf{x})$$

$$=\sum _{\mathbf{z}\in \text{mathbbm}{G}_{d,M}}\sum _{\mathbf{m}\text{mathchar}58|\mathbf{m}|\le \alpha }{\beta }_{\mathbf{z},\mathbf{m}}{\mathbf{x}}^{\mathbf{m}}{\varphi }_{\mathbf{z},M}(\mathbf{x}),$$

$$\mathcal{R}({\hat{f}}_n)-\underset{f\in \mathcal{F}}{\inf }\mathcal{R}(f)=\underset{\text{Estimation error}}{\underbrace{\left[\mathcal{R}({\hat{f}}_n)-\underset{f\in {\mathcal{F}}_n}{\inf }\mathcal{R}(f)\right]}}+\underset{\text{Approximation error}}{\underbrace{\left[\underset{f\in {\mathcal{F}}_n}{\inf }\mathcal{R}(f)-\underset{f\in \mathcal{F}}{\inf }\mathcal{R}(f)\right]}}.$$

$$|\sigma (x_1)-\sigma (x_2)|\le C_{\sigma }|x_1-x_2|,$$

$$\mathcal{N}(\delta ,\mathcal{F},\Vert \cdot {\Vert }_{\infty })\text{mathchar}58=\inf \left\{J\in \text{mathbbm}N\text{mathchar}58\exists f_1,\dots ,f_J\text{mbox}suchthat\mathcal{F}\subset \bigcup _{j=1}^J{B}_{\infty }(f_j,\delta )\right\},$$

$${\mathcal{F}}_{d,1}(L,N,S,B)\text{mathchar}58=\left\{N_{\sigma }(\cdot |\mathbf{\theta })\text{mathchar}58\mathbf{\theta }\in {\Theta }_{d,1}(L,N,S,B)\right\}.$$

$$\log \mathcal{N}$$

$${\mathcal{R}}_{2,f_0}(f)\text{mathchar}58={\mathsf{E}}_{f_0,{\mathsf{P}}_{\mathbf{x}}}(Y-f(\mathbf{X}){)}^2\text{mathchar}58={\mathsf{E}}_{Y|\mathbf{X}\sim \mathsf{N}(f_0(\mathbf{X}),1),\mathbf{X}\sim {\mathsf{P}}_{\mathbf{x}}}(Y-f(\mathbf{X}){)}^2,$$

$${\hat{f}}_n\in \underset{f\in {\mathcal{F}}_{\sigma ,n}}{\mathrm{argmin}}\sum _{i=1}^n{\left(y_i-f({\mathbf{x}}_i)\right)}^2,$$

$$\displaystyle\mathcal{F}_{\sigma,n}\mathrel{\mathop{\mathchar 58\relax}}=\Big{\{}N_{\sigma}(\cdot|\bm{\theta})\mathrel{\mathop{\mathchar 58\relax}}\mathinner{\!\left\lVert N_{\sigma}(\cdot|\bm{\theta})\right\rVert}_{\infty}\leq 2R,\bm{\theta}\in\Theta_{d,1}\mathinner{\left(L_{0}\log n,N_{0}n^{\frac{d}{2\alpha+d}},S_{0}n^{\frac{d}{2\alpha+d}}\log n,B_{0}n^{\kappa}\right)}\Big{\}}$$

$$\underset{f_0\in {\mathcal{H}}^{\alpha ,R}([0,1{]}^d)}{\sup }\mathsf{E}\left[{\mathcal{R}}_{2,f_0}({\hat{f}}_n)-\underset{f\in \mathcal{F}}{\inf }{\mathcal{R}}_{2,f_0}(f)\right]\le C{n}^{-\frac{2\alpha }{2\alpha +d}}{\log }^{3}n,$$

$${\mathcal{R}}_{01,\eta }(f)\text{mathchar}58={\mathsf{E}}_{\eta ,{\mathsf{P}}_{\mathbf{x}}}\text{mathds}1(Yf(\mathbf{X})<0)\text{mathchar}58={\mathsf{E}}_{Y|\mathbf{X}\sim 2\mathsf{Bern}(\eta (\mathbf{X}))-1,\mathbf{X}\sim {\mathsf{P}}_{\mathbf{x}}}\text{mathds}1(Yf(\mathbf{X})<0).$$

$${\mathsf{P}}_{\mathbf{x}}\left(|\eta (\mathbf{X})-1/2|<u\right)\le u^q.$$

$${\hat{f}}_n\in \underset{f\in {\mathcal{F}}_{\sigma ,n}}{\mathrm{argmin}}\sum _{i=1}^n{\ell }_{\mathsf{hinge}}(y_{i}f({\mathbf{x}}_i)),$$

$$\displaystyle\mathcal{F}_{\sigma,n}\mathrel{\mathop{\mathchar 58\relax}}=\Big{\{}N_{\sigma}(\cdot|\bm{\theta})\mathrel{\mathop{\mathchar 58\relax}}\mathinner{\!\left\lVert N_{\sigma}(\cdot|\bm{\theta})\right\rVert}_{\infty}\leq 1,\bm{\theta}\in\Theta_{d,1}\mathinner{\left(L_{0}\log n,N_{0}n^{\nu}\log^{-3\nu}n,S_{0}n^{\nu}\log^{-3\nu+1}n,B_{0}n^{\kappa}\right)}\Big{\}},$$

$$\underset{\eta \in {\mathcal{H}}^{\alpha ,R}([0,1{]}^d)}{\sup }\mathsf{E}\left[{\mathcal{R}}_{01,\eta }({\hat{f}}_n)-\underset{f\in \mathcal{F}}{\inf }{\mathcal{R}}_{01,\eta }(f)\right]\le C{\left(\frac{{\log }^{3}n}{n}\right)}^{\frac{\alpha (q+1)}{\alpha (q+2)+d}},$$

$$f^{\ast }(\mathbf{x})=2\text{mathds}1\left(2\eta (\mathbf{x})-1\ge 0\right)-1,$$

Full text

Keywords: function approximation; deep neural networks; activation functions; Hölder continuity; convergence rates

1 Introduction

Neural networks are learning machines motivated by the architecture of the human brain. Neural networks are comprised of multiple hidden layers, and each of the hidden layers has multiple hidden nodes which consist of an affine map of the outputs from the previous layer and a nonlinear map called an activation function. Deep neural networks have been leading tremendous success in various pattern recognition and machine learning tasks such as object recognition, image segmentation, machine translation and others. For an overview on the empirical success of deep neural networks, we refer to the review paper [19] and recent book [14].

Inspired by the success of deep neural networks, many researchers have tried to give theoretical supports for the success of deep neural networks. Much of the work upto date has focused on the expressivity of deep neural networks, i.e., ability to approximate a rich class of functions efficiently. The well-known classical result on this topic is the universal approximation theorem, which states that every continuous function can be approximated arbitrarily well by a neural network [11, 15, 12, 5, 20]. But these results do not specify the required numbers of layers and nodes of a neural network to achieve a given approximation accuracy.

Recently, several results about the effects of the numbers of layers and nodes of a deep neural network to its expressivity have been reported. They provide upper bounds of the numbers of layers and nodes required for neural networks to uniformly approximate all functions of interest. Examples of a class of functions include the space of rational functions of polynomials [30], the Hölder space [33, 27, 2, 21], Besov and mixed Besov spaces [29] and even a class of discontinuous functions [25, 16].

The nonlinear activation function is a central part that makes neural networks differ from the linear models, that is, a neural network becomes a linear function if the linear activation function is used. Therefore, the choice of an activation function substantially influences on the performance and computational efficiency. Numerous activation functions have been suggested to improve neural network learning [3, 6, 4, 26, 18, 32]. We refer to the papers [13, 26] for an overview of this topic.

There are also many recent theoretical studies about the approximation ability of deep neural networks. However, most of the studies focus on a specific type of the activation function such as ReLU [33, 27, 25, 16, 29], or small classes of activation functions such as sigmoidal functions with additional monotonicity, continuity, and/or boundedness conditions [24, 9, 8, 10, 7] and $m$-admissible functions which are sufficiently smooth and bounded [2]. For definitions of sigmoidal and $m$-admissible functions, see [9] and [2], respectively. Thus a unified theoretical framework still lacks.

In this paper, we investigate the approximation ability of deep neural networks with a quite general class of activation functions. We derive the required numbers of layers and nodes of a deep neural network to approximate any Hölder smooth function upto a given approximation error for the large class of activation functions. Our specified class of activation functions and the corresponding approximation ability of deep neural networks include most of previous results [33, 27, 24, 2] as special cases.

Our general theoretical results of the approximation ability of deep neural networks enables us to study statistical properties of deep neural networks. Schmidt-Hieber, [27] and Kim et al. [17] proved the minimax optimality of a deep neural network estimator with the ReLU activation function in regression and classification problems, respectively. In this paper, we derive similar results for general activation functions.

This paper is structured as follows. In Section 2, we introduce some notions about deep neural networks. In Section 3, we introduce two large classes of activation functions. In Section 4, we present our main result on the approximation ability of a deep neural network with the general activation function considered in Section 3. In Section 5, we apply the result in Section 4 to the supervised learning problems of regression and classification. Conclusions are given in Section 6. The proofs of all results are given in Appendix.

Notation

We denote by $\mathds{1}\left({\cdot}\right)$ the indicator function. Let $\mathbb{R}$ be the set of real numbers and $\mathbbm{N}$ be the set of natural numbers. For a real valued vector $\mathbf{x}\equiv(x_{1},\dots,x_{d})$, we let $\mathinner{\!\left\lvert\mathbf{x}\right\rvert}_{0}\mathrel{\mathop{\mathchar 58\relax}}=\sum_{j=1}^{d}\mathds{1}(x_{j}\neq 0)$, $\mathinner{\!\left\lvert\mathbf{x}\right\rvert}_{p}\mathrel{\mathop{\mathchar 58\relax}}=\mathinner{\left(\sum_{j=1}^{d}|x_{j}|^{p}\right)}^{1/p}$ for $p\in[1,\infty)$ and $\mathinner{\!\left\lvert\mathbf{x}\right\rvert}_{\infty}\mathrel{\mathop{\mathchar 58\relax}}=\max_{1\leq j\leq d}|x_{j}|$. For simplicity, we let $|\mathbf{x}|\mathrel{\mathop{\mathchar 58\relax}}=|\mathbf{x}|_{1}.$ For a real valued function $f(x)\mathrel{\mathop{\mathchar 58\relax}}\mathbb{R}\to\mathbb{R}$, we let $f^{\prime}(a),f^{\prime\prime}(a)$ and $f^{\prime\prime\prime}(a)$ are the first, second and third order derivatives of $f$ at $a$, respectively. We let $f^{\prime}(a+)\mathrel{\mathop{\mathchar 58\relax}}=\lim_{\epsilon\downarrow 0}(f(a+\epsilon)-f(a))/\epsilon$ and $f^{\prime}(a-)\mathrel{\mathop{\mathchar 58\relax}}=\lim_{\epsilon\downarrow 0}(f(a-\epsilon)-f(a))/\epsilon$. For $x\in\mathbb{R}$, we write $(x)_{+}\mathrel{\mathop{\mathchar 58\relax}}=\max\{x,0\}$.

2 Deep Neural Networks

In this section we provide a mathematical representation of deep neural networks. A neural network with $L\in\mathbbm{N}$ layers, $n_{l}\in\mathbbm{N}$ many nodes at the $l$-th hidden layer for $l=1,\dots,L$, input of dimension $n_{0}$, output of dimension $n_{L+1}$ and nonlinear activation function $\sigma\mathrel{\mathop{\mathchar 58\relax}}\mathbb{R}\to\mathbb{R}$ is expressed as

[TABLE]

where $\mathsf{A}_{l}\mathrel{\mathop{\mathchar 58\relax}}\mathbb{R}^{n_{l-1}}\to\mathbb{R}^{n_{l}}$ is an affine linear map defined by $\mathsf{A}_{l}(\mathbf{x})=\mathbf{W}_{l}\mathbf{x}+\mathbf{b}_{l}$ for given $n_{l}\times n_{l-1}$ dimensional weight matrix $\mathbf{W}_{l}$ and $n_{l}$ dimensional bias vector $\mathbf{b}_{l}$ and $\sigma_{l}\mathrel{\mathop{\mathchar 58\relax}}\mathbb{R}^{n_{l}}\to\mathbb{R}^{n_{l}}$ is an element-wise nonlinear activation map defined by $\sigma_{l}(\mathbf{z})\mathrel{\mathop{\mathchar 58\relax}}=(\sigma(z_{1}),\dots,\sigma(z_{n_{l}}))^{\top}$. Here, $\bm{\theta}$ denotes the set of all weight matrices and bias vectors \bm{\theta}\mathrel{\mathop{\mathchar 58\relax}}=\mathinner{\bigl{(}(\mathbf{W}_{1},\mathbf{b}_{1}),(\mathbf{W}_{2},\mathbf{b}_{2}),\dots,(\mathbf{W}_{L+1},\mathbf{b}_{L+1})\bigr{)}}, which we call $\bm{\theta}$ the parameter of the neural network, or simply, a network parameter.

We introduce some notations related to the network parameter. For a network parameter $\bm{\theta}$, we write $L(\bm{\theta})$ for the number of hidden layers of the corresponding neural network, and write $n_{\max}(\bm{\theta})$ for the maximum of the numbers of hidden nodes at each layer. Following a standard convention, we say that $L(\bm{\theta})$ is the depth of the deep neural network and $n_{\max}(\bm{\theta})$ is the width of the deep neural network. We let $\mathinner{\!\left\lvert\bm{\theta}\right\rvert}_{0}$ be the number of nonzero elements of $\bm{\theta}$, i.e.,

[TABLE]

where $\text{vec}(\mathbf{W}_{l})$ transforms the matrix $\mathbf{W}_{l}$ into the corresponding vector by concatenating the column vectors. We call $\mathinner{\!\left\lvert\bm{\theta}\right\rvert}_{0}$ sparsity of the deep neural network. Let $\mathinner{\!\left\lvert\bm{\theta}\right\rvert}_{\infty}$ be the largest absolute value of elements of $\bm{\theta}$, i.e.,

[TABLE]

We call $\mathinner{\!\left\lvert\bm{\theta}\right\rvert}_{\infty}$ magnitude of the deep neural network. We let $\textsf{in}(\bm{\theta})$ and $\textsf{out}(\bm{\theta})$ be the input and output dimensions of the deep neural network, respectively. We denote by $\Theta_{d,o}(L,N)$ the set of network parameters with depth $L$, width $N$, input dimension $d$ and output dimension $o$, that is,

[TABLE]

We further define a subset of $\Theta_{d,o}(L,N)$ with restrictions on sparsity and magnitude as

[TABLE]

3 Classes of Activation Functions

In this section, we consider two classes of activation functions. These two classes include most of commonly used activation functions. Definitions and examples of each class of activation functions are provided in the consecutive two subsections.

3.1 Piecewise Linear Activation Functions

We first consider piecewise linear activation functions.

Definition 3.1**.**

A function $\sigma\mathrel{\mathop{\mathchar 58\relax}}\mathbb{R}\to\mathbb{R}$ is continuous piecewise linear if it is continuous and there exist a finite number of break points $a_{1}\leq a_{2}\leq\cdots\leq a_{K}\in\mathbb{R}$ with $K\in\mathbbm{N}$ such that $\sigma^{\prime}(a_{k}-)\neq\sigma^{\prime}(a_{k}+)$ for every $k=1,\dots,K$ and $\sigma(x)$ is linear on $(-\infty,a_{1}],[a_{1},a_{2}],\dots,[a_{K-1},a_{K}],[a_{K},\infty)$.

Throughout this paper, we write “picewise linear” instead of “continuous picewise linear” for notational simplicity unless there is a confusion. The representative examples of piecewise linear activation functions are as follows:

• •

ReLU: $\sigma(x)=\max\{x,0\}$.

• •

Leaky ReLU: : $\sigma(x)=\max\{x,ax\}$ for $a\in(0,1)$.

The ReLU activation function is the most popular choice in practical applications due to better gradient propagation and efficient computation [13]. In this reason, most of the recent results on the function approximation by deep neural networks are based on the ReLU activation function [33, 27, 25, 16, 29]. In Section 4, as Yarotsky, [33] did, we extend these results to any continuous piecewise linear activation function by showing that the ReLU activation function can be exactly represented by a linear combination of piecewise linear activation functions. A formal proof for this argument is presented in Appendix A.1.

3.2 Locally Quadratic Activation Functions

One of the basic building blocks in approximation by deep neural networks is the square function, which should be approximated precisely. Piecewise linear activation functions have zero curvature (i.e., constant first-order derivative) inside each interval divided by its break points, which makes it relatively difficult to approximate the square function efficiently. But if there is an interval on which the activation function has nonzero curvature, the square function can be approximated more efficiently, which is a main motivation of considering a new class of activation functions that we call locally quadratic.

Definition 3.2**.**

A function $\sigma\mathrel{\mathop{\mathchar 58\relax}}\mathbb{R}\to\mathbb{R}$ is locally quadratic if there exits an open interval $(a,b)\subset\mathbb{R}$ on which $\sigma$ is three times continuously differentiable with bounded derivatives and there exists $t\in(a,b)$ such that $\sigma^{\prime}(t)\neq 0$ and $\sigma^{\prime\prime}(t)\neq 0$.

We now give examples of locally quadratic activation functions. First of all, any nonlinear smooth activation function with nonzero second derivative, is locally quadratic. Examples are:

• •

Sigmoid: $\displaystyle\sigma(x)=\frac{1}{1+\textup{e}^{-x}}.$

• •

Tangent hyperbolic: $\displaystyle\sigma(x)=\frac{\textup{e}^{x}-\textup{e}^{-x}}{\textup{e}^{x}+\textup{e}^{-x}}.$

• •

Inverse square root unit (ISRU) [4]: $\displaystyle\sigma(x)=\frac{x}{\sqrt{1+ax^{2}}}$ for $a>0$.

• •

Soft clipping [18]: $\displaystyle\sigma(x)=\frac{1}{a}\log\mathinner{\left(\frac{1+\textup{e}^{ax}}{1+\textup{e}^{a(x-1)}}\right)}$ for $a>0$.

• •

SoftPlus [13]: $\sigma(x)=\log(1+\textup{e}^{x})$.

• •

Swish [26]: $\displaystyle\sigma(x)=\frac{x}{1+\textup{e}^{-x}}$.

In addition, piecewise smooth function having nonzero second derivative on at least one piece, is also locally quadratic. Examples are:

• •

Rectified power unit (RePU) [21]: $\displaystyle\sigma(x)=\max\{x^{k},0\}$ for $k\in\mathbbm{N}\setminus\{1\}$.

• •

Exponential linear unit (ELU) [6]: $\displaystyle\sigma(x)=a(\textup{e}^{x}-1)\mathds{1}\left({x\leq 0}\right)+x\mathds{1}\left({x>0}\right)$ for $a>0$.

• •

Inverse square root linear unit (ISRLU) [4]: $\displaystyle\sigma(x)=\frac{x}{\sqrt{1+ax^{2}}}\mathds{1}\left({x\leq 0}\right)+x\mathds{1}\left({x>0}\right)$ for $a>0$.

• •

Softsign [3]: $\displaystyle\sigma(x)=\frac{x}{1+|x|}.$

• •

Square nonlinearity [32]:

$\displaystyle\sigma(x)=\mathds{1}\left({x>2}\right)+(x-x^{2}/4)\mathds{1}\left({0\leq x\leq 2}\right)+(x+x^{2}/4)\mathds{1}\left({-2\leq x<0}\right)-\mathds{1}\left({x<-2}\right)$.

4 Approximation of Smooth Functions by Deep Neural Networks

In this section we introduce the function class we consider and show the approximation ability of the deep neural networks with a activation function considered in Section 3.

4.1 Hölder Smooth Functions

We recall the definition of Hölder smooth functions. For a $d$-dimensional multiple index $\mathbf{m}\equiv(m_{1},\dots,m_{d})\in\mathbbm{N}_{0}^{d}$ where $\mathbbm{N}_{0}\mathrel{\mathop{\mathchar 58\relax}}=\mathbbm{N}\cup\{0\}$, we let $\mathbf{x}^{\mathbf{m}}\mathrel{\mathop{\mathchar 58\relax}}=x_{1}^{m_{1}}\cdots x_{d}^{m_{d}}$ for $\mathbf{x}\in\mathbb{R}^{d}$. For a function $f\mathrel{\mathop{\mathchar 58\relax}}\mathcal{X}\to\mathbb{R}$, where $\mathcal{X}$ denotes the domain of the function, we let $\|f\|_{\infty}\mathrel{\mathop{\mathchar 58\relax}}=\sup_{\mathbf{x}\in\mathcal{X}}|f(\mathbf{x})|$. We use notation

[TABLE]

for $\mathbf{m}\in\mathbbm{N}_{0}^{d}$ to denote a derivative of $f$ of order $\mathbf{m}$. We denote by $\mathcal{C}^{m}(\mathcal{X})$, the space of $m$ times differentiable functions on $\mathcal{X}$ whose partial derivatives of order $\mathbf{m}$ with $|\mathbf{m}|\leq m$ are continuous. We define the Hölder coefficient of order $s\in(0,1]$ as

[TABLE]

For a positive real value $\alpha$, the Hölder space of order $\alpha$ is defined as

[TABLE]

where $\|f\|_{\mathcal{H}^{\alpha}(\mathcal{X})}$ denotes the Hölder norm defined by

[TABLE]

We denote by $\mathcal{H}^{\alpha,R}(\mathcal{X})$ the closed ball in the Hölder space of radius $R$ with respect to the Hölder norm, i.e.,

[TABLE]

4.2 Approximation of Hölder Smooth Functions

We present our main theorem in this section.

Theorem 4.1**.**

Let $d\in\mathbbm{N}$, $\alpha>0$ and $R>0$. Let the activation function $\sigma$ be either continuous piecewise linear or locally quadratic. Let $f\in\mathcal{H}^{\alpha,R}([0,1]^{d})$. Then there exist positive constants $L_{0}$, $N_{0}$, $S_{0}$ and $B_{0}$ depending only on $d$, $\alpha$, $R$ and $\sigma$ such that, for any $\epsilon>0$, there is a neural network

[TABLE]

satisfying

[TABLE]

The result of Theorem 4.1 is equivalent to the results on the approximation by ReLU neural networks [33, 27] in a sense that the upper bounds of the depth, width and sparsity are the same orders of those for ReLU, namely, depth $=O(\log(1/\epsilon))$, width $=O(\epsilon^{-d/\alpha})$ and sparsity $=O(\epsilon^{-d/\alpha}\log(1/\epsilon))$. We remark that each upper bound is equivalent to the corresponding lower bound established by [33] up to logarithmic factor.

For piecewise linear activation functions, Yarotsky, [33] derived similar results to ours. For locally quadratic activation functions, only special classes of activation functions were considered in the previous work. Li et al. [21] considered the RePU activation function and Bauer and Kohler, [2] considered sufficiently smooth and bounded activation functions which include the sigmoid, tangent hyperbolic, ISRU and soft clipping activation functions. However, soft plus, swish, ELU, ISRLU, softsign and square nonlinearity activation functions are new ones only considered in our results.

Even if the orders of the depth, width and sparsity are the same for both both piecewise linear and locally quadratic activation functions, the ways of approximating a smooth function by use of these two activation function classes are quite different. To describe this point, let us provide an outline of the proof. We first consider equally spaced grid points with length $1/M$ inside the $d$-dimensional unit hypercube $[0,1]^{d}$. Let $\mathbbm{G}_{d,M}$ be the set of such grid points, namely,

[TABLE]

For a given Hölder smooth function $f$ of order $\alpha$, we first find a “local” function for each grid that approximates the target function near the grid point but vanishes at apart from the grid point. To be more specific, we construct the local functions $g_{\mathbf{z}}$, $\mathbf{z}\in\mathbbm{G}_{d,M}$ which satisfies:

[TABLE]

for some universal constant $C>0$. The inequality (4.3) implies that the more grid points we used, the more accurate approximation we get. Moreover, the quality of approximation is improved when the target function is more smooth (i.e., large $\alpha$) and low dimensional (i.e., small $d$ ). In fact, $g_{\mathbf{z},M}(\mathbf{x})$ is given by a product of the Taylor polynomial $P_{\mathbf{z},M}(\mathbf{x})\mathrel{\mathop{\mathchar 58\relax}}=\sum_{\mathbf{m}\in\mathbbm{N}_{0}^{d}\mathrel{\mathop{\mathchar 58\relax}}|\mathbf{m}|\leq\alpha}\mathinner{\left(\partial^{\mathbf{m}}f\right)}(\mathbf{z})\frac{(\mathbf{x}-\mathbf{z})^{\mathbf{m}}}{\mathbf{m}!}$ at $\mathbf{z}$ and the local basis function $\phi_{\mathbf{z},M}(\mathbf{x})\mathrel{\mathop{\mathchar 58\relax}}=\prod_{j=1}^{d}(1/M-|x_{j}-z_{j}|)_{+}$, where $\mathbf{m}!\mathrel{\mathop{\mathchar 58\relax}}=\prod_{j=1}^{d}m_{j}!$. By simple algebra, we have

[TABLE]

where $\beta_{\mathbf{z},\mathbf{m}}\mathrel{\mathop{\mathchar 58\relax}}=\sum_{\tilde{\mathbf{m}}\mathrel{\mathop{\mathchar 58\relax}}\tilde{\mathbf{m}}\geq\mathbf{m},|\tilde{\mathbf{m}}|\leq\alpha}\mathinner{\left(\partial^{\tilde{\mathbf{m}}}f\right)}(\mathbf{z})\frac{(-\mathbf{z})^{\tilde{\mathbf{m}}-\mathbf{m}}}{\mathbf{m}!(\tilde{\mathbf{m}}-\mathbf{m})!}$.

The second stage is to approximate each monomial $\mathbf{x}^{\mathbf{m}}$ and each local basis function $\phi_{\mathbf{z},M}(\mathbf{x})$ by deep neural networks. Each monomial can be approximated more efficiently by a deep neural network with a locally quadratic activation function than a piecewise linear activation function since each monomial has nonzero curvature. On the other hand, the local basis function can be approximated more efficiently by a deep neural network with a piecewise linear activation than a locally quadratic activation function since the local basis function is piecewise linear itself. That is, there is a trade-off in using either a piecewise linear or a locally quadratic activation function.

We close this section by giving a comparison of our result to the approximation error analysis of [2]. Bauer and Kohler, [2] studies approximation of the Hölder smooth function of order $\alpha$ by a two layer neural network with $m$-admissible activation functions with $m\geq\alpha$, where a function $\sigma$ is called $m$-admissible if (1) $\sigma$ is at least $m+1$ times continuously differentiable with bounded derivatives; (2) a point $t\in\mathbb{R}$ exists, where all derivatives up to the order $m$ of $\sigma$ are different from zero; and (3) $|\sigma(x)-1|\leq 1/x$ for $x>0$ and $|\sigma(x)|\leq 1/|x|$ for $x<0$. Our notion of locally quadratic activation functions is a generalized version of the $m$-admissibility.

In the proof of [2], the condition $m\geq\alpha$ is necessary because they approximate any monomial of order $\mathbf{m}$ with $|\mathbf{m}|\leq\alpha$ with a two layer neural network, which is impossible when $m<\alpha$. We drop the condition $m\geq\alpha$ by showing that any monomial of order $\mathbf{m}$ with $|\mathbf{m}|\leq\alpha$ can be approximated by deep neural network with a finite number of layers, which depends on $\alpha$.

5 Application to Statistical Learning Theory

In this section, we apply our results about the approximation error of neural networks to the supervised learning problems of regression and classification. Let $\mathcal{X}$ be the input space and $\mathcal{Y}$ the output space. Let $\mathcal{F}$ be a given class of measurable functions from $\mathcal{X}$ to $\mathcal{Y}$. Let $\mathsf{P}_{0}$ be the true but unknown data generating distribution on $\mathcal{X}\times\mathcal{Y}$. The aim of supervised learning is to find a predictive function that minimizes the population risk $\mathcal{R}(f)\mathrel{\mathop{\mathchar 58\relax}}=\mathsf{E}_{(\mathbf{X},Y)\sim\mathsf{P}_{0}}\ell(Y,f(\mathbf{X}))$ with respect to a given loss function $\ell$. Since $\mathsf{P}_{0}$ is unknown, we cannot directly minimize the population risk, and thus any estimator $\hat{f}$ inevitably has the excess risk which is defined as $\mathcal{R}(\hat{f})-\inf_{f\in\mathcal{F}}\mathcal{R}(f).$ For a given sample of size $n$, let $\mathcal{F}_{n}$ be a given subset of $\mathcal{F}$ called a sieve and let $(\mathbf{x}_{1},y_{1}),\ldots,(\mathbf{x}_{n},y_{n})$ be observed (training) data of input–output pairs assumed to be independent realizations of $(\mathbf{X},Y)$ following $\mathsf{P}_{0}.$ Let $\hat{f}_{n}$ be an estimated function among $\mathcal{F}_{n}$ based on the training data $(\mathbf{x}_{1},y_{1}),\ldots,(\mathbf{x}_{n},y_{n}).$ The excess risk of $\hat{f}_{n}$ is decomposed to approximation and estimation errors as

[TABLE]

There is a trade-off between approximation and estimation errors. If the function class $\mathcal{F}_{n}$ is sufficiently large to approximate the optimal estimator $f^{*}\mathrel{\mathop{\mathchar 58\relax}}=\operatorname*{\arg\!\min}_{f\in\mathcal{F}}\mathcal{R}(f)$ well, then the estimation error becomes large due to high variance. In contrast, if $\mathcal{F}_{n}$ is small, it leads to low estimation error but it suffers from large approximation error.

One of the advantages of deep neural networks is that we can construct a sieve which has good approximation ability as well as low complexity. Schmidt-Hieber, [27] and Kim et al. [17] proved that a neural network estimator can achieve the optimal balance between the approximation and estimation errors to obtain the minimax optimal convergence rates in regression and classification problems, respectively. But they only considered the ReLU activation function. Based on the results of Theorem 4.1, we can easily extend their results to general activation functions.

The main tool to derive the minimax optimal convergence rate is that the complexity of a class of functions generated by a deep neural network is not affected much by a choice of an activation function, provided that the activation function is Lipschitz continuous. The function $\sigma\mathrel{\mathop{\mathchar 58\relax}}\mathbb{R}\to\mathbb{R}$ is Lipschitz continuous if there is a constant $C_{\sigma}>0$ such that

[TABLE]

for any $x_{1},x_{2}\in\mathbb{R}$. Here, $C_{\sigma}$ is called the Lipschitz constant. We use the covering number with respect to the $L_{\infty}$ norm $\mathinner{\!\left\lVert\cdot\right\rVert}_{\infty}$ as a measure of complexity of function classes. We recall the definition of the covering number. Let $\mathcal{F}$ be a given class of real-valued functions defined on $\mathcal{X}$. Let $\delta>0$. A collection $\{f_{j}\in\mathcal{F}\mathrel{\mathop{\mathchar 58\relax}}j=1,\dots,J\}$ is called a $\delta$-covering set of $\mathcal{F}$ with respect to the $L_{\infty}$ norm if for all $f\in\mathcal{F}$, there exists $f_{j}$ in the collection such that $\|f-f_{j}\|_{\infty}\leq\delta$. The cardinality of the minimal $\delta$-covering set is called the $\delta$-covering number of $\mathcal{F}$ with respect to the $L_{\infty}$ norm which is denoted by $\mathcal{N}(\delta,\mathcal{F},\|\cdot\|_{\infty}).$ That is,

[TABLE]

where $B_{\infty}(f_{j},\delta)\mathrel{\mathop{\mathchar 58\relax}}=\{f\in\mathcal{F}\mathrel{\mathop{\mathchar 58\relax}}\|f-f_{j}\|_{\infty}\leq\delta\}$. The following proposition provides the covering number of a class of functions generated by neural networks.

Proposition 5.1**.**

Assume that the activation function $\sigma$ is Lipschitz continuous with the Lipschitz constant $C_{\sigma}$. Consider a class of functions generated by a deep neural network

[TABLE]

For any $\delta>0$,

[TABLE]

where $B\vee 1\mathrel{\mathop{\mathchar 58\relax}}=\max\{B,1\}$.

The result in Proposition 5.1 is very similar to the existing results in literature, e.g., Theorem 14.5 of [1], Lemma 5 of [27] and Lemma 3 of [29]. We employ similar techniques used in [1, 27, 29] to obtain the version presented here. We give the proof of this proposition in Appendix B.

All of the activation functions considered in Section 3 except RePU satisfy the Lipschitz condition (5.2) and hence Proposition 5.1 can be applied. An interesting implication of Proposition 5.1 is that the complexity of the function class generated by deep neural networks is not affected by the choice of an activation function. Hence, the remaining step to derive the convergence rate of a neural network estimator is that approximation accuracies by various activation functions are the same as that of the ReLU neural network.

5.1 Application to Regression

First we consider the regression problem. For simplicity, we let $\mathcal{X}=[0,1]^{d}.$ Suppose that the generated model is $Y|\mathbf{X}=\mathbf{x}\sim\mathsf{N}(f_{0}(\mathbf{x}),1)$ for some $f_{0}\mathrel{\mathop{\mathchar 58\relax}}[0,1]^{d}\to\mathbb{R}.$ The performance of an estimator is measured by the $L_{2}$ risk $\mathcal{R}_{2,f_{0}}(f)$, which is defined by

[TABLE]

where $\mathsf{P}_{\mathbf{x}}$ is the marginal distribution of $\mathbf{X}.$ The following theorem proves that the optimal convergence rate is obtained by the deep neural network estimator of the regression function $f_{0}$ for a general activation function.

Theorem 5.2**.**

Suppose that the activation function $\sigma$ is either piecewise linear or locally quadratic satisfying the Lipschitz condition (5.2). Then there are universal positive constants $L_{0}$, $N_{0}$, $S_{0}$ and $B_{0}$ such that the deep neural network estimator obtained by

[TABLE]

with

[TABLE]

for some $\kappa>0$ satisfies

[TABLE]

for some universal constant $C>0$, where the expectation is taken over the training data.

5.2 Application to Binary Classification

The aim of the binary classification is to find a classifier that predicts the label $y\in\{-1,1\}$ for any input $\mathbf{x}\in[0,1]^{d}$. An usual assumption on the data generating process is that $Y|\mathbf{X}=\mathbf{x}\sim 2\mathsf{Bern}(\eta(\mathbf{x}))-1$ for some $\eta\mathrel{\mathop{\mathchar 58\relax}}[0,1]^{d}\to[0,1]$, where $\mathsf{Bern}(p)$ denotes the Bernoulli distribution with parameter $p$. Note that $\eta(\mathbf{x})$ is the conditional probability function $\mathsf{P}_{0}(Y=1|\mathbf{X}=\mathbf{x}).$ A common approach is, instead of finding a classifier directly, to construct a real valued function $f$, a so-called classification function, and predict the label $y$ based on the sign of $f(\mathbf{x})$. The performance of a classification function is measured by the misclassification error $\mathcal{R}_{01,\eta}(f)$, which is defined by

[TABLE]

It is well known that the convergence rate of the excess risk for classification is faster than that of regression when the conditional probability function $\eta(\mathbf{x})$ satisfies the following condition: there is a constant $q\in[0,\infty]$ such that for any sufficiently small $u>0$, we have

[TABLE]

This condition is called the Tsybakov noise condition and $q$ is called the noise exponent [23, 31]. When $q$ is larger, the classification task is easier since the probability of generating vague samples become smaller. The following theorem proves that the optimal convergence rate can be obtained by the deep neural network estimator with an activation function considered in Section 3. As is done by [17], we consider the hinge loss $\ell_{\mathsf{hinge}}(z)\mathrel{\mathop{\mathchar 58\relax}}=\max\{1-z,0\}.$

Theorem 5.3**.**

Assume the Tsybakov noise condition (5.4) with the noise exponent $q\in[0,\infty]$. Suppose that the activation function $\sigma$, which is either piecewise linear or locally quadratic satisfying the Lipschitz condition (5.2), is used for all hidden layers except the last one and the ReLU activation function is used for the last hidden layer. Then there are universal positive constants $L_{0}$, $N_{0}$, $S_{0}$ and $B_{0}$ such that the deep neural network estimator obtained by

[TABLE]

with

[TABLE]

for $\nu\mathrel{\mathop{\mathchar 58\relax}}=d/\mathinner{\left\{\alpha(q+2)+d\right\}}$ and some $\kappa>0$ satisfies

[TABLE]

for some universal constant $C>0$, where the expectation is taken over the training data.

Note that the Bayes classifier $f^{*}\mathrel{\mathop{\mathchar 58\relax}}=\operatorname*{\arg\!\min}_{f\in\mathcal{F}}\mathcal{R}_{01,\eta}(f)$ is given by

[TABLE]

which is an indicator function. Since a neural network with the ReLU activation function can approximate indicator functions well [25, 16, 17], we use the ReLU activation function in the last layer in order to approximate the Bayes classifier more precisely and thus to achieve the optimal convergence rate.

6 Conclusions

In this study, we established the upper bounds of the required depth, width and sparsity of deep neural networks to approximate any Hölder continuous function for the general classes of activation functions. These classes of activation functions include most of the popularly used activation functions. The derived upper bounds of the depth, width and sparsity are optimal in a sense that they are equivalent to the lower bounds up to logarithmic factors. We used this generalization of the approximation error analysis to extend the statistical optimality of the deep neural network estimator in regression and classification problems, where the activation function is other than the ReLU.

Our construction of neural networks for approximation reveals that the piecewise linear activation functions are more efficient in approximating local basis functions while locally quadratic activation functions are more efficient in approximating polynomials. Hence if the activation function has both piecewise linear region and locally quadratic region, we could have a better approximation result. We leave the development of such activation functions as a future work.

Appendix A Proof of Theorem 4.1

A.1 Proof of Theorem 4.1 for Piecewise Linear Activation Functions

The main idea of the proof is that any deep neural network with the ReLU activation function can be exactly reconstructed by a neural network with a piecewise activation function whose proof is in the next lemma that is a slight modification of Proposition 1 (b) of [33].

Lemma A.1**.**

Let $\sigma$ be an any continuous peicewise linear activation function, and $\rho$ be the ReLU activation function. Let $\bm{\theta}\in\Theta_{d,1}(L,N,S,B)$. Then there exists $\bm{\theta}^{*}\in\Theta_{d,1}(L,2N,4S+2LN+1,C_{1}B)$ such that

[TABLE]

where $C_{1}>0$ is a constant depending on the activation function $\sigma$.

Proof.

Let $a$ be any break point of $\sigma$. Note that $\sigma(a-)\neq\sigma(a+)$. Let $r_{0}$ be the distance between $a$ and the closest other break point. Then $\sigma$ is linear on $[a-r_{0},a]$ and $[a,a+r_{0}]$. Then for any $r>0$, the ReLU activation function $\rho(x)\mathrel{\mathop{\mathchar 58\relax}}=(x)_{+}$ is expressed as

[TABLE]

for any $x\in[-r,r]$, where we define $u_{1}\mathrel{\mathop{\mathchar 58\relax}}=1/\mathinner{(\mathinner{(\sigma^{\prime}(a+)-\sigma^{\prime}(a-))}\frac{r_{0}}{2r})}$, $u_{2}\mathrel{\mathop{\mathchar 58\relax}}=-1/\mathinner{(\mathinner{(\sigma^{\prime}(a+)-\sigma^{\prime}(a-))}\frac{r_{0}}{2r})}$ and $v\mathrel{\mathop{\mathchar 58\relax}}=\mathinner{(-\sigma(a)+\sigma\mathinner{\left(a-r_{0}/2\right)})}/\mathinner{(\mathinner{(\sigma^{\prime}(a+)-\sigma^{\prime}(a-))}\frac{r_{0}}{2r})}$.

Let $\bm{\theta}\equiv((\mathbf{W}_{1},\mathbf{b}_{1}),\dots,(\mathbf{W}_{L+1},\mathbf{b}_{L+1}))\in\Theta_{d,1}(L,N,S,B)$ be given. Since both input $\mathbf{x}\in[0,1]^{d}$ and the network parameter $\bm{\theta}$ are bounded, we can take a sufficiently large $r$ so that Equation (A.1) holds for any hidden nodes of the network $\bm{\theta}$. We consider the deep neural network $\bm{\theta}^{*}\equiv((\mathbf{W}_{1}^{*},\mathbf{b}_{1}^{*})$, $\dots,(\mathbf{W}_{L+1}^{*},\mathbf{b}_{L+1}^{*}))\in\Theta_{d,1}(L,2N),$ where we set

[TABLE]

for $l=1,\dots,L$ and

[TABLE]

Here, $\mathbf{1}_{n}$ denotes the $n$-dimensional vector of $1^{\prime}s$. Then by Equation (A.1) and some algebra, we have that $N_{\sigma}(\mathbf{x}|\bm{\theta}^{*})=N_{\rho}(\mathbf{x}|\bm{\theta})$ for any $\mathbf{x}\in[0,1]^{d}$. For the sparsity of $\bm{\theta}^{*}$, we note that

[TABLE]

which implies that $\mathinner{\!\left\lvert\bm{\theta}^{*}\right\rvert}_{0}\leq 4\mathinner{\!\left\lvert\bm{\theta}\right\rvert}_{0}+2L(\bm{\theta})n_{\max}(\bm{\theta})+1$. ∎

Thanks to Lemma A.1, to prove Theorem 4.1 for piecewise linear activation functions, it suffices to show the approximation ability of the ReLU networks, which is already done by [27] as in the next lemma.

Lemma A.2** (Theorem 5 of [27]).**

Let $\rho$ be the ReLU activation function. For any $f\in\mathcal{H}^{\alpha,R}([0,1]^{d})$ and any integers $m\geq 1$ and $M\geq\max\mathinner{\left\{(\alpha+1)^{d},(R+1)\textup{e}^{d}\right\}}$, there exists a network parameter $\bm{\theta}\in\Theta_{d,1}(L,N,S,1)$ such that

[TABLE]

where $L=8+(m+5)(1+\lceil\log_{2}(d\vee\alpha)\rceil)$, $N=6(d+\left\lceil{\alpha}\right\rceil)M$, and $S=141(d+\alpha+1)^{3+d}M(m+6)$.

Theorem 4.1 for piecewise linear activation functions is a direct consequence of Lemmas A.1 and A.2, which is summarized as follows.

Proof of Theorem 4.1 for piecewise linear activation functions.

Let $\rho$ be the ReLU activation function. By letting $M=3^{d}(2R)^{d/\alpha}\epsilon^{-d/\alpha}$ and $m=\log_{2}\left(2(2R+1)(1+d^{2}+\alpha^{2})18^{d}(2R)^{d/\alpha}\epsilon^{-d/\alpha-1}\right),$ Lemma A.2 implies that there exists a network parameter $\bm{\theta}^{\prime}$ such that $\sup_{\mathbf{x}\in[0,1]^{d}}|N_{\rho}(\mathbf{x}|\bm{\theta}^{\prime})-f(\mathbf{x})|\leq\epsilon$ with $L(\bm{\theta}^{\prime})\leq L_{0}^{\prime}\log\left(1/\epsilon\right)$, $n_{\max}(\bm{\theta}^{\prime})\leq N_{0}^{\prime}\epsilon^{-d/\alpha}$ and $|\bm{\theta}^{\prime}|_{0}\leq S_{0}^{\prime}\epsilon^{-d/\alpha}\log\left(1/\epsilon\right)$ for some positive constants $L_{0}^{\prime}$, $N_{0}^{\prime}$, and $S_{0}^{\prime}$ depending only on $\alpha$, $d$ and $R$. Hence by Lemma A.1, there is a network parameter $\bm{\theta}$ producing the same output of the ReLU neural network $N_{\rho}(\cdot|\bm{\theta})$ with $L(\bm{\theta})=L(\bm{\theta}^{\prime})$, $n_{\max}(\bm{\theta})=2n_{\max}(\bm{\theta}^{\prime})$, $|\bm{\theta}|_{0}\leq 4|\bm{\theta}^{\prime}|_{0}+2L(\bm{\theta}^{\prime})n_{\max}(\bm{\theta}^{\prime})+1\leq S_{0}\epsilon^{-d/\alpha}\log\left(1/\epsilon\right)$ and $|\bm{\theta}|_{\infty}\leq B_{0}|\bm{\theta}^{\prime}|_{\infty}$ for some $S_{0}>0$ depending only on $\alpha$, $d$, $R$ and $\sigma$, and some $B_{0}>0$ depending only on $\sigma$, which completes the proof. ∎

A.2 Proof of Theorem 4.1 for Locally Quadratic Activation Functions

Lemma A.3**.**

Assume that an activation function $\sigma$ is locally quadratic. There is a constant $K_{0}$ depending only on the activation function such that for any $K>K_{0}$ the following results hold.

1. (a)

There is a neural network $\bm{\theta}_{2}\in\Theta_{1,1}(1,3)$ with $\mathinner{\!\left\lvert\bm{\theta}_{2}\right\rvert}_{\infty}\leq K^{2}$ such that

[TABLE]

where $C_{1}>0$ is a constant depending only on $\sigma$. 2. (b)

Let $A>0$. There is a neural network parameter $\bm{\theta}_{\times,A}\in\Theta_{2,1}(1,9)$ with $\mathinner{\!\left\lvert\bm{\theta}_{\times,A}\right\rvert}_{\infty}\leq\max\{K^{2},2A^{2}\}$ such that

[TABLE] 3. (c)

Let $\alpha$ be a positive integer. For any multi-index $\mathbf{m}\in\mathbbm{N}_{0}^{d}$ with $|\mathbf{m}|\leq\alpha$, there is a network parameter $\bm{\theta}_{\mathbf{m}}\in\Theta_{d,1}(\left\lceil{\log_{2}\alpha}\right\rceil,9\alpha)$ with $\mathinner{\!\left\lvert\bm{\theta}_{\mathbf{m}}\right\rvert}_{\infty}\leq\max\{K^{2},C_{2}\}$ such that

[TABLE]

for some positive constants $C_{2}$ and $C_{3}$ depending only on $\sigma$ and $\alpha$. 4. (d)

There is a network parameter $\bm{\theta}_{1/2}\in\Theta_{1,1}(\left\lceil{\log K}\right\rceil,15)$ with $\mathinner{\!\left\lvert\bm{\theta}_{1/2}\right\rvert}_{\infty}\leq\max\{K^{2},C_{4}\}$ such that

[TABLE]

for some positive constants $C_{4}$ and $C_{5}$ depending only on $\sigma$. 5. (e)

There is a network parameter $\bm{\theta}_{\textup{abs}}\in\Theta_{1,1}(\left\lceil{\log K}\right\rceil,15)$ with $\mathinner{\!\left\lvert\bm{\theta}_{\textup{abs}}\right\rvert}_{\infty}\leq\max\{K^{2},C_{6}\}$ such that

[TABLE]

for some positive constants $C_{6}$ and $C_{7}$ depending only on $\sigma$.

Proof.

Recall that there is an interval $(a,b)$ on which $\sigma(x)$ is three times continuously differentiable with bounded derivatives and there is $t\in(a,b)$ such that $\sigma^{\prime}(t)\neq 0$ and $\sigma^{\prime\prime}(t)\neq 0$

Proof of (a). Take $K$ large so that $2/K<\min\{|t-b|,|t-a|\}$. Consider a neural network

[TABLE]

Since $\sigma$ is three times continuously differentiable on $(a,b)$ and $(k-1)x/K+t\in(a,b)$ if $x\in[0,1]$, it can be expanded in the Taylor series with Lagrange remainder around $t$ to have

[TABLE]

where $\xi_{k}\in[t-k|x|/K,t+k|x|/K]\subset(a,b).$ Since the third order derivative is bounded on $(a,b)$, we get the desired assertion by retaking $K\leftarrow\sqrt{2/\sigma^{\prime\prime}(t)}K$.

Proof of (b). The proof can be done straightforwardly by the polarization type identity:

[TABLE]

We construct the network as

[TABLE]

where $\bm{\theta}_{2}$ is defined in (A.3). Since $(x_{1}+x_{2})/2A,x_{1}/2A,x_{2}/2A\in[-1,1]$ for $\mathbf{x}\in[-A,A]^{2}$, the triangle inequality implies that $|N_{\sigma}(\mathbf{x}|\bm{\theta}_{\times,A})-x_{1}x_{2}|\leq 6A^{2}C_{1}/K$.

Proof of (c). Let $q\mathrel{\mathop{\mathchar 58\relax}}=\left\lceil{\log_{2}\alpha}\right\rceil$. We construct $\bm{\theta}_{\mathbf{m}}$ as follows. Fix $\mathbf{x}\equiv(x_{1},\dots,x_{d})\in[0,1]^{d}$. We first consider the affine map that transforms $(x_{1},\dots,x_{d})$ to $\mathbf{z}\in[0,1]^{2^{q}}$ which is given by

[TABLE]

The first hidden layer of $\bm{\theta}_{\mathbf{m}}$ pairs neighboring entries in $\mathbf{z}$ and applies the network $\bm{\theta}_{\times,A_{1}}$ defined in (b) with $A_{1}=1$ to each pair. That is, the first hidden layer of $\bm{\theta}_{\mathbf{m}}$ produces

[TABLE]

Note that $\sup_{1\leq j\leq 2^{q-1}}|g_{1,j}-z_{2j-1}z_{2j}|\leq 6C_{1}/K$ and $\sup_{1\leq j\leq 2^{q-1}}|g_{1,j}|\leq 6C_{1}/K+1$, where $6C_{1}/K+1$ can be bounded by some constant $A_{2}>1$ depending only on $C_{1}$ and $K_{0}$. Then the second hidden layer of $\bm{\theta}_{\mathbf{m}}$ pairs neighboring entries of $\mathinner{\left\{g_{1,j}\mathrel{\mathop{\mathchar 58\relax}}j=1,\dots,2^{q-1}\right\}}$ and applies $\bm{\theta}_{\times,A_{2}}$ to each pair to have

[TABLE]

Note that $\sup_{1\leq j\leq 2^{q-2}}|g_{2,j}-g_{1,2j-1}g_{1,2j}|\leq 6C_{1}A_{2}^{2}/K$ and $\sup_{1\leq j\leq 2^{q-2}}|g_{2,j}|\leq 6C_{1}A_{2}^{2}/K+1\leq A_{3}$ for some $A_{3}>1$ depending only on $C_{1}$ and $K_{0}$. We repeat this procedure to produce $\mathinner{\left\{g_{k,j}\mathrel{\mathop{\mathchar 58\relax}}j=1,\dots,2^{q-k}\right\}}$ for $k=3,\dots,q$ with

[TABLE]

for some $A_{k+1}>1$, and we set $N_{\sigma}(\mathbf{x}|\bm{\theta}_{\mathbf{m}})$ equal to $g_{q,1}$.

By applying the triangle inequality repeatedly, we have

[TABLE]

for some $C_{1}^{\prime}>0$ depending only on $C_{1}$, $K_{0}$ and $q$. Since we set $\mathbf{x}$ arbitrary, the proof is done.

Proof of (d). By (b), it is easy to verify that there is a network $\bm{\theta}_{1}\in\Theta_{1,1}(1,6)$ with $|\bm{\theta}_{1}|_{\infty}\leq\max\{K^{2},2\}$ such that $|\sigma(x)-x|\leq C_{1}^{\prime}/K$ for any $x\in[-1,1]$ and some constant $C_{1}^{\prime}>0$. The Taylor series with Lagrange remainder around 1 of $\sqrt{x}$ is given by

[TABLE]

where $\xi\in[0,2],$ and thus

[TABLE]

for some $C_{1}^{\prime}>0,$ where the last inequality is because $n!\geq(n/\textup{e})^{n}\textup{e}$.

Now, we will construct a neural network $\bm{\theta}_{p,J}$ that approximates the polynomial $\sum_{k=0}^{J}\frac{(x-1)^{k}}{k!}$ as follows. The first hidden layer computes $(N_{\sigma}(x-1|\bm{\theta}_{2})/2,N_{\sigma}(x-1|\bm{\theta}_{1}))$ from the input $x$. Then

[TABLE]

for any $x\in[0,1]$ and some constant $C_{2}^{\prime}>0$. The next hidden layer computes $(N_{\sigma}((u,v)|\bm{\theta}_{\times,1+C_{2}^{\prime}/K})/3,N_{\sigma}(u+v|\bm{\theta}_{1}))$ from the input $(u,v)$ from the first hidden layer. Using the triangle inequality, we have that the second hidden layer approximates the vector $((x-1)^{3}/3!,(x-1)^{2}/2+(x-1))$ by error $\leq 2C_{3}^{\prime}/K$ for some $C_{3}^{\prime}>0$. Repeating this procedure, we construct the network $\bm{\theta}_{p,J}\in\Theta_{1,1}(J,15)$ which approximates $\sum_{k=0}^{J}\frac{(x-1)^{k}}{k!}$ by error $\leq C_{4}^{\prime}J/K$ for some $C_{4}^{\prime}>0$. Taking $J=\left\lceil{\log K}\right\rceil$, we observe that $(\textup{e}/J+1)^{J+1}\leq(\textup{e}/\log K)^{\log K+1}\leq\textup{e}K/(\log K)^{\log K}\leq 1/K$ for all sufficiently large $K$, which implies the desired result.

Proof of (e). Let $\zeta\in(0,1)$. Since for any $x\in\mathbb{R}$,

[TABLE]

the function $\sqrt{x^{2}+\zeta^{2}}$ approximates the absolute value function $|x|$ by error $\zeta$. For $\bm{\theta}_{2}$ in (a) and $\bm{\theta}_{1/2}$ in (d), we have that

[TABLE]

for some constant $C_{1}^{\prime}>0.$ We now set $\zeta=1/\sqrt{K}$ and $N_{\sigma}(x|\bm{\theta}_{\textup{abs}})\mathrel{\mathop{\mathchar 58\relax}}=N_{\sigma}(N_{\sigma}(x|\bm{\theta}_{2})+K^{-1}|\bm{\theta}_{1/2})$. Since $(\log K)/K=o(1/\sqrt{K})$, the proof is done. ∎

Proof of Theorem 4.1 for locally quadratic activation functions.

Recall that

[TABLE]

Then by Lemma B.1 of [27],

[TABLE]

From the equivalent representation of the ReLU function $(x)_{+}=(x+|x|)/2$, we can easily check that the neural network $N_{\sigma}(x|\bm{\theta}_{\textup{relu}})\mathrel{\mathop{\mathchar 58\relax}}=\mathinner{\left(N(x|\bm{\theta}_{\textup{abs}})+N_{\sigma}(x|\bm{\theta}_{1})\right)}/2$ with $\bm{\theta}_{\textup{relu}}\in\Theta_{1,1}(\left\lceil{\log K}\right\rceil,21)$ approximates the ReLU function by error $\leq C_{1}^{\prime}/\sqrt{K}$ for some $C_{1}^{\prime}>0$, where $\bm{\theta}_{1}\in\Theta_{1,1}(1,6)$ is defined in the proof of (d) of Lemma A.3 and $\bm{\theta}_{\textup{abs}}\in\Theta_{1,1}(\left\lceil{\log K}\right\rceil,15)$ is defined in (e) of Lemma A.3. For $z\in(0,1)$ and $M\in\mathbbm{N}$, we define

[TABLE]

Then it approximates the function $(1/M-|x-z|)_{+}$ by error $\leq C_{2}^{\prime}/\sqrt{K}$ for some $C_{2}^{\prime}>0.$ In turn, for $\mathbf{z}\in\mathbbm{G}_{d,M}$, by invoking the similar construction used in (c) of Lemma A.3 to approximates the product of $d$ components, we can construct the network $\bm{\theta}_{\phi,\mathbf{z},M}\in\Theta_{1,1}(\left\lceil{\log K}\right\rceil+\left\lceil{\log_{2}d}\right\rceil,21d)$ with $\mathinner{\!\left\lvert\bm{\theta}_{\phi,\mathbf{z},M}\right\rvert}_{\infty}\leq C_{3}^{\prime}K^{2}$ for some $C_{3}^{\prime}>0$ such that

[TABLE]

for some $C_{4}^{\prime}>0$. For each $\mathbf{m}\in\mathbbm{N}_{0}^{d}$ with $|\mathbf{m}|\leq\alpha$, we have the neural network $\bm{\theta}_{\mathbf{m}}$ in (c) of Lemma A.3 that approximates $\mathbf{x}^{\mathbf{m}}$. The number of these networks is $\binom{d+\alpha}{\alpha}$, which is denoted by $A_{\alpha}$. Also there are $|\mathbbm{G}_{d,M}|=(M+1)^{d}$ networks $\bm{\theta}_{\phi,\mathbf{z},M}$ for $\mathbf{z}\in\mathbbm{G}_{d,M}$. We need approximation of each product $\mathbf{x}^{\mathbf{m}}\phi_{\mathbf{z},M}$, which requires additional $A_{\alpha}(M+1)^{d}$ many networks $\bm{\theta}_{\times,A}\in\Theta_{2,1}(1,9)$, where $\bm{\theta}_{\times,A}$ is defined as in (A.4) for some $A>1$ not depending on $M$ and $K$. Finally we construct the output layer which computes the weighted sum of $\mathinner{\left\{N_{\sigma}\mathinner{\left((N_{\sigma}(\mathbf{x}|\bm{\theta}_{\mathbf{m}}),N_{\sigma}(\mathbf{x}|\bm{\theta}_{\phi,\mathbf{z},M}))|\bm{\theta}_{\times,A}\right)}\mathrel{\mathop{\mathchar 58\relax}}\mathbf{m}\in\mathbbm{N}_{0}^{d},|\mathbf{m}|\leq\alpha,\mathbf{z}\in\mathbbm{G}_{d,M}\right\}}$. Letting $\bm{\theta}_{f,K,M}$ be the network constructed above, we can check that

[TABLE]

for some positive constants $C_{5}^{\prime}$ and $C_{6}^{\prime}.$ In addition, we have $L(\bm{\theta}_{f,K,M})\leq 1+(\left\lceil{\log K}\right\rceil+\left\lceil{\log_{2}(\alpha\vee d)}\right\rceil\leq C_{7}^{\prime}\left\lceil{\log K}\right\rceil$ and $n_{\max}(\bm{\theta}_{f,K,M})\leq C_{8}^{\prime}A_{\alpha}(M+1)^{d}$ for some positive constants $C_{7}^{\prime}$ and $C_{8}^{\prime}$. For sparsity of the network, we have

[TABLE]

for some $C_{9}^{\prime}>0.$ Taking $M+1=\epsilon^{-1/\alpha}$ and $K=\epsilon^{-2d/\alpha-2}$, we have

[TABLE]

so that $\mathinner{\!\left\lVert P_{M}-N_{\sigma}(\cdot|\bm{\theta}_{f,K,M})\right\rVert}_{\infty}\leq C_{10}^{\prime}\epsilon$ for some $C_{10}^{\prime}>0.$ Since $\mathinner{\!\left\lVert f-P_{M}\right\rVert}_{\infty}\leq RM^{-\alpha}\leq C_{11}^{\prime}\epsilon$ for some $C_{11}^{\prime}>0,$ the proof is done. ∎

Appendix B Proofs of Proposition 5.1

Proof.

Given a deep neural network $\bm{\theta}=((\mathbf{W}_{1},\mathbf{b}_{1}),\dots,(\mathbf{W}_{L+1},\mathbf{b}_{L+1}))\in\Theta_{d,1}(L,N,S,B)$, we define $\check{N}_{l,\sigma,\bm{\theta}}\mathrel{\mathop{\mathchar 58\relax}}\mathbb{R}^{d}\to\mathbb{R}^{n_{l-1}}$ and $\hat{N}_{l,\sigma,\bm{\theta}}\mathrel{\mathop{\mathchar 58\relax}}\mathbb{R}^{n_{l}}\to\mathbb{R}$ as

[TABLE]

for $l\in 2,\dots,L$, where $\mathsf{A}_{l}\mathbf{x}=\mathbf{W}_{l}\mathbf{x}+\mathbf{b}_{l}$. Corresponding to the last and first layer, we define $\check{N}_{1,\sigma,\bm{\theta}}(\mathbf{x})=\mathbf{x}$ and $\hat{N}_{L+1,\sigma,\bm{\theta}}(\mathbf{x})=\mathbf{x}$. Note that $N_{\sigma}(\mathbf{x}|\bm{\theta})=\hat{N}_{l+1,\sigma,\bm{\theta}}\circ\mathsf{A}_{l}\circ\check{N}_{l,\sigma,\bm{\theta}}(\mathbf{x})$. For given $\delta>0,$ let $\bm{\theta}=((\mathbf{W}_{1},\mathbf{b}_{1}),\dots,(\mathbf{W}_{L+1},\mathbf{b}_{L+1}))\in\Theta_{d,1}(L,N,S,B)$ and $\bm{\theta}^{*}=((\mathbf{W}_{1}^{*},\mathbf{b}_{1}^{*}),\dots,(\mathbf{W}_{L+1}^{*},\mathbf{b}_{L+1}^{*}))\in\Theta_{d,1}(L,N,S,B)$ be two neural network parameter such that $\mathinner{\!\left\lvert\text{vec}(\mathbf{W}_{l}-\mathbf{W}_{l}^{*})\right\rvert}_{\infty}\leq\delta$ and $\mathinner{\!\left\lvert\mathbf{b}_{l}-\mathbf{b}_{l}^{*}\right\rvert}_{\infty}\leq\delta$ for $l=1,\dots,L+1$. Let $C_{\sigma}$ be the Lipschitz constant of $\sigma$. We observe that

[TABLE]

and similarly, $\mathinner{\!\left\lVert\hat{N}_{l,\sigma,\bm{\theta}}\right\rVert}_{\infty}\leq(C_{\sigma}BN)^{L-l+1}$. Letting $\mathsf{A}_{l}^{*}\mathbf{x}=\mathbf{W}_{l}^{*}\mathbf{x}+\mathbf{b}_{l}^{*},$ we have

[TABLE]

Thus, for a fixed sparsity pattern (i.e., the location of nonzero elements in $\bm{\theta}$), the covering number is bounded by $\mathinner{\left[\delta/L\mathinner{\left\{C_{\sigma}(B\vee 1)(N+1)\right\}}^{L}\right]}^{-S}$. Since the number of the sparsity patterns is bounded by $\binom{(N+1)^{L}}{S}\leq(N+1)^{LS}$, the log of covering number is bounded above by

[TABLE]

which completes the proof. ∎

Appendix C Proof of Theorem 5.2

The proof Theorem 5.2 is based on the following oracle inequality.

Lemma C.1** (Lemma 4 of [27]).**

Assume that $Y|\mathbf{X}=\mathbf{x}\sim\mathsf{N}(f_{0}(\mathbf{x}),1)$ for some $f_{0}$ with $\mathinner{\!\left\lVert f_{0}\right\rVert}_{\infty}\leq R$. Let $\mathcal{F}^{\dagger}$ be a given function class from $[0,1]^{d}$ to $[-2R,2R]$, and let $\hat{f}$ be any estimator in $\mathcal{F}^{\dagger}$. Then for any $\delta\in(0,1]$, we have

[TABLE]

with

[TABLE]

where the expectations are taken over the training data.

Proof of Theorem 5.2.

We apply Lemma C.1 to $\mathcal{F}^{\dagger}=\mathcal{F}_{\sigma,n}$ and $\hat{f}=\hat{f}_{n}\in\operatorname*{\arg\!\min}_{f\in\mathcal{F}_{\sigma,n}}\sum_{i=1}^{n}\mathinner{\left(y_{i}-f(\mathbf{x}_{i})\right)}^{2}$. By definition of $\hat{f}_{n}$, we have $\Delta_{n}=0$. Also it can be easily verified that $f_{0}=\operatorname*{\arg\!\min}_{f\in\mathcal{F}}\mathcal{R}_{2,f_{0}}(f)$ and $\mathsf{E}_{f_{0},\mathsf{P}_{\mathbf{x}}}\mathinner{\left(\hat{f}_{n}(\mathbf{X})-f_{0}(\mathbf{X})\right)}^{2}=\mathcal{R}_{2,f_{0}}(\hat{f}_{n})-\mathcal{R}_{2,f_{0}}(f_{0}).$ Set $\delta=1/n$. By Proposition 5.1,

[TABLE]

for some $C_{1}^{\prime}>0$. If a function $f_{n}$ is approximates $f_{0}$ by error $\epsilon$ which is sufficeintly small, then $\mathinner{\!\left\lVert f_{n}\right\rVert}_{\infty}\leq 2R$ since $\mathinner{\!\left\lVert f_{0}\right\rVert}_{\infty}\leq R$. Now, Theorem 4.1 implies that there is $f_{n}\in\mathcal{F}_{\sigma,n}$ such that

[TABLE]

which completes the proof. ∎

Appendix D Proof of Theorem 5.3

For a given real-valued function $f$, let $\mathcal{R}_{\mathsf{hinge},\eta}(f)\mathrel{\mathop{\mathchar 58\relax}}=\mathsf{E}_{Y|\mathbf{X}\sim 2\mathsf{Bern}(\eta(\mathbf{X}))-1,\mathbf{X}\sim\mathsf{P}_{\mathbf{x}}}\ell_{\mathsf{hinge}}(Yf(\mathbf{X}))$, which we call the hinge risk. The proof of Theorem 5.3 is based on the following theorem, which is given in [17].

Lemma D.1** (Theorem 6 of [17]).**

Assume that $\eta(\mathbf{x})$ satisfies the Tsybakov noise condition (5.4) with the noise exponent $q\in[0,\infty]$. Assume that there exists a sequence $(\delta_{n})_{n\in\mathbbm{N}}$ such that

• •

there exists a sequence of classes of functions $\{\mathcal{F}_{n}\}_{n\in\mathbbm{N}}$ with $\sup_{n\in\mathbbm{N}}\sup_{f\in\mathcal{F}_{n}}\mathinner{\!\left\lVert f\right\rVert}_{\infty}\leq F$ for some $F>0$ such that there is $f_{n}\in\mathcal{F}_{n}$ with $\mathcal{R}_{\mathsf{hinge},\eta}(f_{n})-\min_{f\in\mathcal{F}}\mathcal{R}_{\mathsf{hinge},\eta}(f)\leq C_{1}\delta_{n}$ for some universal constant $C_{1}>0$;

• •

$\log\mathcal{N}(\delta_{n},\mathcal{F}_{n},\mathinner{\!\left\lVert\cdot\right\rVert}_{\infty})\leq C_{2}n\delta_{n}^{(q+2)/(q+1)}$* for some universal constant $C_{2}>0$.*

Then the estimator $\hat{f}_{n}$ obtained by

[TABLE]

satisfies

[TABLE]

for some universal constant $C_{3}>0$, where the expectation is taken over the training data.

Proof of Theorem 5.3.

It is well known that $f^{*}=2\mathds{1}\left({\eta(\cdot)\geq 1/2}\right)-1=\operatorname*{\arg\!\min}_{f\in\mathcal{F}}\mathcal{R}_{\mathsf{hinge},\eta}(f)$, i.e., the hinge risk minimizer is equal to the Bayes classifier [22]. The first step is to find a function $f_{n}\in\mathcal{F}_{\sigma,n}$ which approximates the Bayes classifier $f^{*}$ well. Let $(\xi_{n})_{n\in\mathbbm{N}}$ be a given sequence of positive integers. Since $\eta\in\mathcal{H}^{\alpha,R}([0,1]^{d})$, by Theorem 6, for each $\xi_{n}$ there exists $\bm{\theta}_{n}$ such that $\|N_{\sigma}(\cdot|\bm{\theta}_{n})-\eta(\cdot)\|_{\infty}\leq\xi_{n}$ with at most $O(\log(1/\xi_{n}))$ layers, $O(\xi_{n}^{-d/\alpha})$ nodes at each layer and $O(\xi_{n}^{-d/\alpha}\log(1/\xi_{n}))$ nonzero parameters. We construct the neural network $f_{n}$ by adding one ReLU layer to $N_{\sigma}(\cdot|\bm{\theta}_{n})$ to have

[TABLE]

where $\rho$ is the ReLU activation function. Note that $f_{n}(\mathbf{x})$ is equal to $1$ if $N_{\sigma}(\mathbf{x}|\bm{\theta}_{n})\geq 1/2+\xi_{n}$, $(N_{\sigma}(\mathbf{x}|\bm{\theta}_{n})-1/2)/\xi_{n}$ if $1/2\leq(N_{\sigma}(\mathbf{x}|\bm{\theta}_{n})<1/2+\xi_{n}$ and $-1$ otherwise. Let

[TABLE]

Then on $B(4\xi_{n})$, $|f_{n}(\mathbf{x})-f^{*}(\mathbf{x})|=0$, since $N_{\sigma}(\mathbf{x}|\bm{\theta}_{n})-1/2=(\eta(\mathbf{x})-1/2)-((N_{\sigma}(\mathbf{x}|\bm{\theta}_{n})-\eta(\mathbf{x}))\geq\xi_{n}$ when $2\eta(\mathbf{x})-1>4\xi_{n}.$ Similarly we can show that $N_{\sigma}(\mathbf{x}|\bm{\theta}_{n})-1/2<-\xi_{n}$ when $2\eta(\mathbf{x})-1<-4\xi_{n}$. Therefore the Tsybakov noise condition (5.4) implies

[TABLE]

for some $C_{1}^{\prime}>0,$ where the first equality follows from Theorem 2.31 of [28].

We take $\delta_{n}=C_{1}^{\prime}\xi_{n}^{q+1}$. Then there are positive constants $L_{0}$, $N_{0}$, $S_{0}$ and $B_{0}$ such that $f_{n}\in\mathcal{F}_{\sigma,n}$ where

[TABLE]

for some $\kappa^{\prime}>0$. Propostion 5.1 implies that the log covering number of $\mathcal{F}_{\sigma,n}$ is bounded above by

[TABLE]

Note that to satisfy the entropy condition of Lemma D.1, $\delta_{n}$ should satisfy

[TABLE]

for some $C_{2}^{\prime}>0$. If we let $\delta_{n}=(\log^{3}n/n)^{\alpha(q+1)/(\alpha(q+2)+d)}$, the condition (D.1) holds and so the proof is done. ∎