A Selective Overview of Deep Learning

  • 2019-04-10 17:53:15
  • Jianqing Fan, Cong Ma, Yiqiao Zhong
  • 37


Deep learning has arguably achieved tremendous success in recent years. Insimple words, deep learning uses the composition of many nonlinear functions tomodel the complex dependency between input features and labels. While neuralnetworks have a long history, recent advances have greatly improved theirperformance in computer vision, natural language processing, etc. From thestatistical and scientific perspective, it is natural to ask: What is deeplearning? What are the new characteristics of deep learning, compared withclassical methods? What are the theoretical foundations of deep learning? Toanswer these questions, we introduce common neural network models (e.g.,convolutional neural nets, recurrent neural nets, generative adversarial nets)and training techniques (e.g., stochastic gradient descent, dropout, batchnormalization) from a statistical point of view. Along the way, we highlightnew characteristics of deep learning (including depth and over-parametrization)and explain their practical and theoretical benefits. We also sample recentresults on theories of deep learning, many of which are only suggestive. Whilea complete understanding of deep learning remains elusive, we hope that ourperspectives and discussions serve as a stimulus for new statistical research.


Quick Read (beta)

A Selective Overview of Deep Learning00footnotetext: Author names are sorted alphabetically.

Jianqing Fan Department of Operations Research and Financial Engineering, Princeton University, Princeton, NJ 08544, USA; Email: {jqfan, congm, yiqiaoz}@princeton.edu.    Cong Ma33footnotemark: 3    Yiqiao Zhong11footnotemark: 1

Deep learning has arguably achieved tremendous success in recent years. In simple words, deep learning uses the composition of many nonlinear functions to model the complex dependency between input features and labels. While neural networks have a long history, recent advances have greatly improved their performance in computer vision, natural language processing, etc. From the statistical and scientific perspective, it is natural to ask: What is deep learning? What are the new characteristics of deep learning, compared with classical methods? What are the theoretical foundations of deep learning?

To answer these questions, we introduce common neural network models (e.g., convolutional neural nets, recurrent neural nets, generative adversarial nets) and training techniques (e.g., stochastic gradient descent, dropout, batch normalization) from a statistical point of view. Along the way, we highlight new characteristics of deep learning (including depth and over-parametrization) and explain their practical and theoretical benefits. We also sample recent results on theories of deep learning, many of which are only suggestive. While a complete understanding of deep learning remains elusive, we hope that our perspectives and discussions serve as a stimulus for new statistical research.

Keywords: neural networks, over-parametrization, stochastic gradient descent, approximation theory, generalization error.

1 Introduction

Modern machine learning and statistics deal with the problem of learning from data: given a training dataset {(yi,𝒙i)}1in where 𝒙id is the input and yi is the output11 1 When the label y is given, this problem is often known as supervised learning. We mainly focus on this paradigm throughout this paper and remark sparingly on its counterpart, unsupervised learning, where y is not given., one seeks a function f:d from a certain function class that has good prediction performance on test data. This problem is of fundamental significance and finds applications in numerous scenarios. For instance, in image recognition, the input 𝒙 (reps. the output y) corresponds to the raw image (reps. its category) and the goal is to find a mapping f() that can classify future images accurately. Decades of research efforts in statistical machine learning have been devoted to developing methods to find f() efficiently with provable guarantees. Prominent examples include linear classifiers (e.g., linear/logistic regression, linear discriminant analysis), kernel methods (e.g., support vector machines), tree-based methods (e.g., decision trees, random forests), nonparametric regression (e.g., nearest neighbors, local kernel smoothing), etc. Roughly speaking, each aforementioned method corresponds to a different function class from which the final classifier f() is chosen.

Deep learning [70], in its simplest form, proposes the following compositional function class:

{f(𝒙;𝜽)=𝐖L𝝈L(𝐖L-1𝝈2(𝐖2𝝈1(𝐖1𝒙)))|𝜽={𝐖1,,𝐖L}}. (1)

Here, for each 1lL, 𝝈() is some nonlinear function, and 𝜽={𝐖1,,𝐖L} consists of matrices with appropriate sizes. Though simple, deep learning has made significant progress towards addressing the problem of learning from data over the past decade. Specifically, it has performed close to or better than humans in various important tasks in artificial intelligence, including image recognition [50], game playing [113], and machine translation [131]. Owing to its great promise, the impact of deep learning is also growing rapidly in areas beyond artificial intelligence; examples include statistics [15, 110, 76, 103, 41], applied mathematics [129, 22], clinical research [28], etc.

Table 1: Winning models for ILSVRC image classification challenge.
Model Year # Layers # Params Top-5 error
Shallow <2012 >25%
AlexNet 2012 8 61M 16.4%
VGG19 2014 19 144M 7.3%
GoogleNet 2014 22 7M 6.7%
ResNet-152 2015 152 60M 3.6%

To get a better idea of the success of deep learning, let us take the ImageNet Challenge [106] (also known as ILSVRC) as an example. In the classification task, one is given a training dataset consisting of 1.2 million color images with 1000 categories, and the goal is to classify images based on the input pixels. The performance of a classifier is then evaluated on a test dataset of 100 thousand images, and in the end the top-5 error22 2 The algorithm makes an error if the true label is not contained in the 5 predictions made by the algorithm. is reported. Table 1 highlights a few popular models and their corresponding performance. As can be seen, deep learning models (the second to the last rows) have a clear edge over shallow models (the first row) that fit linear models/tree-based models on handcrafted features. This significant improvement raises a foundational question:

  • Why is deep learning better than classical methods on tasks like image recognition?

Figure 1: Visualization of trained filters in the first layer of AlexNet. The model is pre-trained on ImageNet and is downloadable via PyTorch package torchvision.models. Each filter contains 11×11×3 parameters and is shown as an RGB color map of size 11×11.

1.1 Intriguing new characteristics of deep learning

It is widely acknowledged that two indispensable factors contribute to the success of deep learning, namely (1) huge datasets that often contain millions of samples and (2) immense computing power resulting from clusters of graphics processing units (GPUs). Admittedly, these resources are only recently available: the latter allows to train larger neural networks which reduces biases and the former enables variance reduction. However, these two alone are not sufficient to explain the mystery of deep learning due to some of its “dreadful” characteristics: (1) over-parametrization: the number of parameters in state-of-the-art deep learning models is often much larger than the sample size (see Table 1), which gives them the potential to overfit the training data, and (2) nonconvexity: even with the help of GPUs, training deep learning models is still NP-hard [8] in the worst case due to the highly nonconvex loss function to minimize. In reality, these characteristics are far from nightmares. This sharp difference motivates us to take a closer look at the salient features of deep learning, which we single out a few below.

1.1.1 Depth

Deep learning expresses complicated nonlinearity through composing many nonlinear functions; see (1). The rationale for this multilayer structure is that, in many real-world datasets such as images, there are different levels of features and lower-level features are building blocks of higher-level ones. See [133] for a visualization of trained features of convolutional neural nets; here in Figure 1, we sample and visualize weights from a pre-trained AlexNet model. This intuition is also supported by empirical results from physiology and neuroscience [56, 2]. The use of function composition marks a sharp difference from traditional statistical methods such as projection pursuit models [38] and multi-index models [73, 27]. It is often observed that depth helps efficiently extract features that are representative of a dataset. In comparison, increasing width (e.g., number of basis functions) in a shallow model leads to less improvement. This suggests that deep learning models excel at representing a very different function space that is suitable for complex datasets.

1.1.2 Algorithmic regularization

The statistical performance of neural networks (e.g., test accuracy) depends heavily on the particular optimization algorithms used for training [130]. This is very different from many classical statistical problems, where the related optimization problems are less complicated. For instance, when the associated optimization problem has a relatively simple structure (e.g., convex objective functions, linear constraints), the solution to the optimization problem can often be unambiguously computed and analyzed. However, in deep neural networks, due to over-parametrization, there are usually many local minima with different statistical performance [72]. Nevertheless, common practice runs stochastic gradient descent with random initialization and finds model parameters with very good prediction accuracy.

1.1.3 Implicit prior learning

It is well observed that deep neural networks trained with only the raw inputs (e.g., pixels of images) can provide a useful representation of the data. This means that after training, the units of deep neural networks can represent features such as edges, corners, wheels, eyes, etc.; see [133]. Importantly, the training process is automatic in the sense that no human knowledge is involved (other than hyper-parameter tuning). This is very different from traditional methods, where algorithms are designed after structural assumptions are posited. It is likely that training an over-parametrized model efficiently learns and incorporates the prior distribution p(𝒙) of the input, even though deep learning models are themselves discriminative models. With automatic representation of the prior distribution, deep learning typically performs well on similar datasets (but not very different ones) via transfer learning.

(a) MNIST images (b) training and test accuracies
Figure 2: (a) shows the images in the public dataset MNIST; and (b) depicts the training and test accuracies along the training dynamics. Note that the training accuracy is approaching 100% and the test accuracy is still high (no overfitting).

1.2 Towards theory of deep learning

Despite the empirical success, theoretical support for deep learning is still in its infancy. Setting the stage, for any classifier f, denote by 𝔼(f) the expected risk on fresh sample (a.k.a. test error, prediction error or generalization error), and by 𝔼n(f) the empirical risk/training error averaged over a training dataset. Arguably, the key theoretical question in deep learning is

why is E(f^n) small, where f^n is the classifier returned by the training algorithm?

We follow the conventional approximation-estimation decomposition (sometimes, also bias-variance tradeoff) to decompose the term 𝔼(f^n) into two parts. Let be the function space expressible by a family of neural nets. Define f*=argminf𝔼(f) to be the best possible classifier and f*=argminf𝔼(f) to be the best classifier in . Then, we can decompose the excess error 𝔼(f^n)-𝔼(f*) into two parts:

=𝔼(f*)-𝔼(f*)approximation error+𝔼(f^n)-𝔼(f*)estimation error. (2)

Both errors can be small for deep learning (cf. Figure 2), which we explain below.

  • The approximation error is determined by the function class . Intuitively, the larger the class, the smaller the approximation error. Deep learning models use many layers of nonlinear functions (Figure 3)that can drive this error small. Indeed, in Section 5, we provide recent theoretical progress of its representation power. For example, deep models allow efficient representation of interactions among variable while shallow models cannot.

  • The estimation error reflects the generalization power, which is influenced by both the complexity of the function class and the properties of the training algorithms. Interestingly, for over-parametrized deep neural nets, stochastic gradient descent typically results in a near-zero training error (i.e., 𝔼n(f^n)0; see e.g. left panel of Figure 2). Moreover, its generalization error 𝔼(f^n) remains small or moderate. This “counterintuitive” behavior suggests that for over-parametrized models, gradient-based algorithms enjoy benign statistical properties; we shall see in Section 7 that gradient descent enjoys implicit regularization in the over-parametrized regime even without explicit regularization (e.g., 2 regularization).

The above two points lead to the following heuristic explanation of the success of deep learning models. The large depth of deep neural nets and heavy over-parametrization lead to small or zero training errors, even when running simple algorithms with moderate number of iterations. In addition, these simple algorithms with moderate number of steps do not explore the entire function space and thus have limited complexities, which results in small generalization error with a large sample size. Thus, by combining the two aspects, it explains heuristically that the test error is also small.

1.3 Roadmap of the paper

We first introduce basic deep learning models in Sections 24, and then examine their representation power via the lens of approximation theory in Section 5. Section 6 is devoted to training algorithms and their ability of driving the training error small. Then we sample recent theoretical progress towards demystifying the generalization power of deep learning in Section 7. Along the way, we provide our own perspectives, and at the end we identify a few interesting questions for future research in Section 8. The goal of this paper is to present suggestive methods and results, rather than giving conclusive arguments (which is currently unlikely) or a comprehensive survey. We hope that our discussion serves as a stimulus for new statistics research.

2 Feed-forward neural networks

Before introducing the vanilla feed-forward neural nets, let us set up necessary notations for the rest of this section. We focus primarily on classification problems, as regression problems can be addressed similarly. Given the training dataset {(yi,𝒙i)}1in where yi[K]{1,2,,K} and 𝒙id are independent across i[n], supervised learning aims at finding a (possibly random) function f^(𝒙) that predicts the outcome y for a new input 𝒙, assuming (y,𝒙) follows the same distribution as (yi,𝒙i). In the terminology of machine learning, the input 𝒙i is often called the feature, the output yi called the label, and the pair (yi,𝒙i) is an example. The function f^ is called the classifier, and estimation of f^ is training or learning. The performance of f^ is evaluated through the prediction error (yf^(𝒙)), which can be often estimated from a separate test dataset.

As with classical statistical estimation, for each k[K], a classifier approximates the conditional probability (y=k|𝒙) using a function fk(𝒙;𝜽k) parametrized by 𝜽k. Then the category with the highest probability is predicted. Thus, learning is essentially estimating the parameters 𝜽k. In statistics, one of the most popular methods is (multinomial) logistic regression, which stipulates a specific form for the functions fk(𝒙;𝜽k): let zk=𝒙𝜷k+αk and fk(𝒙;𝜽k)=Z-1exp(zk) where Z=k=1Kexp(zk) is a normalization factor to make {fk(𝒙;𝜽k)}1kK a valid probability distribution. It is clear that logistic regression induces linear decision boundaries in d, and hence it is restrictive in modeling nonlinear dependency between y and 𝒙. The deep neural networks we introduce below provide a flexible framework for modeling nonlinearity in a fairly general way.

2.1 Model setup

From the high level, deep neural networks (DNNs) use composition of a series of simple nonlinear functions to model nonlinearity


where denotes composition of two functions and L is the number of hidden layers, and is usually called depth of a NN model. Letting 𝒉(0)𝒙, one can recursively define 𝒉(l)=𝐠(l)(𝒉(l-1)) for all =1,2,,L. The feed-forward neural networks, also called the multilayer perceptrons (MLPs), are neural nets with a specific choice of 𝐠(l): for =1,,L, define

𝒉()=𝐠(l)(𝒉(l-1))𝝈(𝐖()𝒉(-1)+𝒃()), (3)

where 𝐖(l) and 𝒃(l) are the weight matrix and the bias/intercept, respectively, associated with the l-th layer, and 𝝈() is usually a simple given (known) nonlinear function called the activation function. In words, in each layer , the input vector 𝒉(-1) goes through an affine transformation first and then passes through a fixed nonlinear function 𝝈(). See Figure 3 for an illustration of a simple MLP with two hidden layers. The activation function 𝝈() is usually applied element-wise, and a popular choice is the ReLU (Rectified Linear Unit) function:

[𝝈(𝒛)]j=max{zj,0}. (4)

Other choices of activation functions include leaky ReLU, tanh function [79] and the classical sigmoid function (1+e-z)-1, which is less used now.

Figure 3: A feed-forward neural network with an input layer, two hidden layers and an output layer. The input layer represents raw features {𝒙i}1in. Both hidden layers compute an affine transform (a.k.s. indices) of the input and then apply an element-wise activation function 𝝈(). Finally, the output returns a linear transform followed by the softmax activation (resp. simply a linear transform) of the hidden layers for the classification (resp. regression) problem.

Given an output 𝒉(L) from the final hidden layer and a label y, we can define a loss function to minimize. A common loss function for classification problems is the multinomial logistic loss. Using the terminology of deep learning, we say that 𝒉(L) goes through an affine transformation and then the soft-max function:


Then the loss is defined to be the cross-entropy between the label y (in the form of an indicator vector) and the score vector (f1(𝒙;𝜽),,fK(𝒙;𝜽)), which is exactly the negative log-likelihood of the multinomial logistic regression model:

(𝒇(𝒙;𝜽),y)=-k=1K𝟙{y=k}logpk, (5)

where 𝜽{𝐖(),𝒃():1L+1}. As a final remark, the number of parameters scales with both the depth L and the width (i.e., the dimensionality of 𝐖()), and hence it can be quite large for deep neural nets.

2.2 Back-propagation in computational graphs

Training neural networks follows the empirical risk minimization paradigm that minimizes the loss (e.g., (5)) over all the training data. This minimization is usually done via stochastic gradient descent (SGD). In a way similar to gradient descent, SGD starts from a certain initial value 𝜽0 and then iteratively updates the parameters 𝜽t by moving it in the direction of the negative gradient. The difference is that, in each update, a small subsample [n] called a mini-batch—which is typically of size 32–512—is randomly drawn and the gradient calculation is only on instead of the full batch [n]. This saves considerably the computational cost in calculation of gradient. By the law of large numbers, this stochastic gradient should be close to the full sample one, albeit with some random fluctuations. A pass of the whole training set is called an epoch. Usually, after several or tens of epochs, the error on a validation set levels off and training is complete. See Section 6 for more details and variants on training algorithms.

The key to the above training procedure, namely SGD, is the calculation of the gradient (𝜽), where

(𝜽)||-1i(𝒇(𝒙i;𝜽),yi). (6)

Gradient computation, however, is in general nontrivial for complex models, and it is susceptible to numerical instability for a model with large depth. Here, we introduce an efficient approach, namely back-propagation, for computing gradients in neural networks.

Back-propagation [105] is a direct application of the chain rule in networks. As the name suggests, the calculation is performed in a backward fashion: one first computes /𝒉(L), then /𝒉(L-1), , and finally /𝒉(1). For example, in the case of the ReLU activation function33 3 The issue of non-differentiability at the origin is often ignored in implementation., we have the following recursive/backward relation

𝒉(-1)=𝒉()𝒉(-1)𝒉()=(𝐖())𝖽𝗂𝖺𝗀(𝟙{𝐖()𝒉(-1)+𝒃()𝟎})𝒉() (7)

where 𝖽𝗂𝖺𝗀() denotes a diagonal matrix with elements given by the argument. Note that the calculation of /𝒉(-1) depends on /𝒉(), which is the partial derivatives from the next layer. In this way, the derivatives are “back-propagated” from the last layer to the first layer. These derivatives {/𝒉()} are then used to update the parameters. For instance, the gradient update for 𝐖() is given by

𝐖()𝐖()-η𝐖(),whereWjm()=hj()σhm(-1), (8)

where σ=1 if the j-th element of 𝐖()𝒉(-1)+𝒃() is nonnegative, and σ=0 otherwise. The step size η>0, also called the learning rate, controls how much parameters are changed in a single update.

Figure 4: The computational graph illustrates the loss (9). For simplicity, we omit the bias terms. Symbols inside nodes represent functions, and symbols outside nodes represent function outputs (vectors/scalars). matmul is matrix multiplication, relu is the ReLU activation, cross entropy is the cross entropy loss, and SoS is the sum of squares.

A more general way to think about neural network models and training is to consider computational graphs. Computational graphs are directed acyclic graphs that represent functional relations between variables. They are very convenient and flexible to represent function composition, and moreover, they also allow an efficient way of computing gradients. Consider an MLP with a single hidden layer and an 2 regularization:

λ(𝜽)=(𝜽)+rλ(𝜽)=(𝜽)+λ(j,j(Wj,j(1))2+j,j(Wj,j(2))2), (9)

where (𝜽) is the same as (6), and λ0 is a tuning parameter. A similar example is considered in [45]. The corresponding computational graph is shown in Figure 4. Each node represents a function (inside a circle), which is associated with an output of that function (outside a circle). For example, we view the term (𝜽) as a result of 4 compositions: first the input data 𝒙 multiplies the weight matrix 𝐖(1) resulting in 𝒖(1), then it goes through the ReLU activation function relu resulting in 𝒉(1), then it multiplies another weight matrix 𝐖(2) leading to 𝒑, and finally it produces the cross-entropy with label y as in (5). The regularization term is incorporated in the graph similarly.

A forward pass is complete when all nodes are evaluated starting from the input 𝒙. A backward pass then calculates the gradients of λ with respect to all other nodes in the reverse direction. Due to the chain rule, the gradient calculation for a variable (say, /𝒖(1)) is simple: it only depends on the gradient value of the variables (/𝒉) the current node points to, and the function derivative evaluated at the current variable value (𝝈(𝒖(1))). Thus, in each iteration, a computation graph only needs to (1) calculate and store the function evaluations at each node in the forward pass, and then (2) calculate all derivatives in the backward pass.

Back-propagation in computational graphs forms the foundations of popular deep learning programming softwares, including TensorFlow [1] and PyTorch [92], which allows more efficient building and training of complex neural net models.

3 Popular models

Moving beyond vanilla feed-forward neural networks, we introduce two other popular deep learning models, namely, the convolutional neural networks (CNNs) and the recurrent neural networks (RNNs). One important characteristic shared by the two models is weight sharing, that is some model parameters are identical across locations in CNNs or across time in RNNs. This is related to the notion of translational invariance in CNNs and stationarity in RNNs. At the end of this section, we introduce a modular thinking for constructing more flexible neural nets.

3.1 Convolutional neural networks

The convolutional neural network (CNN) [71, 40] is a special type of feed-forward neural networks that is tailored for image processing. More generally, it is suitable for analyzing data with salient spatial structures. In this subsection, we focus on image classification using CNNs, where the raw input (image pixels) and features of each hidden layer are represented by a 3D tensor 𝑿d1×d2×d3. Here, the first two dimensions d1,d2 of 𝑿 indicate spatial coordinates of an image while the third d3 indicates the number of channels. For instance, d3 is 3 for the raw inputs due to the red, green and blue channels, and d3 can be much larger (say, 256) for hidden layers. Each channel is also called a feature map, because each feature map is specialized to detect the same feature at different locations of the input, which we will soon explain. We now introduce two building blocks of CNNs, namely the convolutional layer and the pooling layer.

  1. 1.

    Convolutional layer (CONV). A convolutional layer has the same functionality as described in (3), where the input feature 𝑿d1×d2×d3 goes through an affine transformation first and then an element-wise nonlinear activation. The difference lies in the specific form of the affine transformation. A convolutional layer uses a number of filters to extract local features from the previous input. More precisely, each filter is represented by a 3D tensor 𝑭kw×w×d3 (1kd~3), where w is the size of the filter (typically 3 or 5) and d~3 denotes the total number of filters. Note that the third dimension d3 of 𝑭k is equal to that of the input feature 𝑿. For this reason, one usually says that the filter has size w×w, while suppressing the third dimension d3. Each filter 𝑭k then convolves with the input feature 𝑿 to obtain one single feature map 𝑶k(d1-w+1)×(d1-w+1), where44 4 To simplify notation, we omit the bias/intercept term associated with each filter.

    Oijk=[𝑿]ij,𝑭k=i=1wj=1wl=1d3[𝑿]i+i-1,j+j-1,l[𝑭k]i,j,l. (10)

    Here [𝑿]ijw×w×d3 is a small “patch” of 𝑿 starting at location (i,j). See Figure 5 for an illustration of the convolution operation. If we view the 3D tensors [𝑿]ij and 𝑭k as vectors, then each filter essentially computes their inner product with a part of 𝑿 indexed by i,j (which can be also viewed as convolution, as its name suggests). One then pack the resulted feature maps {𝑶k} into a 3D tensor 𝑶 with size (d1-w+1)×(d1-w+1)×d~3, where

    [𝑶]ijk=[𝑶k]ij. (11)
    Figure 5: 𝑿28×28×3 represents the input feature consisting of 28×28 spatial coordinates in a total number of 3 channels / feature maps. 𝑭k5×5×3 denotes the k-th filter with size 5×5. The third dimension 3 of the filter automatically matches the number 3 of channels in the previous input. Every 3D patch of 𝑿 gets convolved with the filter 𝑭k and this as a whole results in a single output feature map X~:,:,k with size 24×24×1. Stacking the outputs of all the filters {𝑭k}1kK will lead to the output feature with size 24×24×K.

    The outputs of convolutional layers are then followed by nonlinear activation functions. In the ReLU case, we have

    X~ijk=σ(Oijk),i[d1-w+1],j[d2-w+1],k[d~3]. (12)

    The convolution operation (10) and the ReLU activation (12) work together to extract features 𝑿~ from the input 𝑿. Different from feed-forward neural nets, the filters 𝑭k are shared across all locations (i,j). A patch [𝑿]ij of an input responds strongly (that is, producing a large value) to a filter 𝑭k if they are positively correlated. Therefore intuitively, each filter 𝑭k serves to extract features similar to 𝑭k.

    As a side note, after the convolution (10), the spatial size d1×d2 of the input 𝑿 shrinks to (d1-w+1)×(d2-w+1) of 𝑿~. However one may want the spatial size unchanged. This can be achieved via padding, where one appends zeros to the margins of the input 𝑿 to enlarge the spatial size to (d1+w-1)×(d2+w-1). In addition, a stride in the convolutional layer determines the gap i-i and j-j between two patches 𝑿ij and 𝑿ij: in (10) the stride is 1, and a larger stride would lead to feature maps with smaller sizes.

  2. 2.

    Pooling layer (POOL). A pooling layer aggregates the information of nearby features into a single one. This downsampling operation reduces the size of the features for subsequent layers and saves computation. One common form of the pooling layer is composed of the 2×2 max-pooling filter. It computes max{Xi,j,k,Xi+1,j,k,Xi,j+1,k,Xi+1,j+1,k}, that is, the maximum of the 2×2 neighborhood in the spatial coordinates; see Figure 6 for an illustration. Note that the pooling operation is done separately for each feature map k. As a consequence, a 2×2 max-pooling filter acting on 𝑿d1×d2×d3 will result in an output of size d1/2×d2/2×d3. In addition, the pooling layer does not involve any parameters to optimize. Pooling layers serve to reduce redundancy since a small neighborhood around a location (i,j) in a feature map is likely to contain the same information.

Figure 6: A 2×2 max pooling layer extracts the maximum of 2 by 2 neighboring pixels/features across the spatial dimension.

In addition, we also use fully-connected layers as building blocks, which we have already seen in Section 2. Each fully-connected layer treats input tensor 𝑿 as a vector Vec(𝑿), and computes 𝑿~=𝝈(𝐖Vec(𝑿)). A fully-connected layer does not use weight sharing and is often used in the last few layers of a CNN. As an example, Figure 7 depicts the well-known LeNet 5 [71], which is composed of two sets of CONV-POOL layers and three fully-connected layers.

Figure 7: LeNet is composed of an input layer, two convolutional layers, two pooling layers and three fully-connected layers. Both convolutions are valid and use filters with size 5×5. In addition, the two pooling layers use 2×2 average pooling.

3.2 Recurrent neural networks

Recurrent neural nets (RNNs) are another family of powerful models, which are designed to process time series data and other sequence data. RNNs have successful applications in speech recognition [107], machine translation [131], genome sequencing [21], etc. The structure of an RNN naturally forms a computational graph, and can be easily combined with other structures such as CNNs to build large computational graph models for complex tasks. Here we introduce vanilla RNNs and improved variants such as long short-term memory (LSTM).

(a) One-to-many (b) Many-to-one (c) Many-to-many
Figure 8: Vanilla RNNs with different inputs/outputs settings. (a) has one input but multiple outputs; (b) has multiple inputs but one output; (c) has multiple inputs and outputs. Note that the parameters are shared across time steps.

3.2.1 Vanilla RNNs

Suppose we have general time series inputs 𝒙1,𝒙2,,𝒙T. A vanilla RNN models the “hidden state” at time t by a vector 𝒉t, which is subject to the recursive formula

𝒉t=𝒇𝜽(𝒉t-1,𝒙t). (13)

Here, f𝜽 is generally a nonlinear function parametrized by 𝜽. Concretely, a vanilla RNN with one hidden layer has the following form55 5 Similar to the activation function 𝝈(), the function 𝐭𝐚𝐧𝐡() means element-wise operations.

𝒉t =𝐭𝐚𝐧𝐡(𝐖hh𝒉t-1+𝐖xh𝒙t+𝒃𝒉),wheretanh(a)=e2a-1e2a+1,
𝒛t =𝝈(𝐖hy𝒉t+𝒃𝒛),

where 𝐖hh,𝐖xh,𝐖hy are trainable weight matrices, 𝒃𝒉,𝒃𝒛 are trainable bias vectors, and 𝒛t is the output at time t. Like many classical time series models, those parameters are shared across time. Note that in different applications, we may have different input/output settings (cf. Figure 8). Examples include

  • One-to-many: a single input with multiple outputs; see Figure 8(a). A typical application is image captioning, where the input is an image and outputs are a series of words.

  • Many-to-one: multiple inputs with a single output; see Figure 8(b). One application is text sentiment classification, where the input is a series of words in a sentence and the output is a label (e.g., positive vs. negative).

  • Many-to-many: multiple inputs and outputs; see Figure 8(c). This is adopted in machine translation, where inputs are words of a source language (say Chinese) and outputs are words of a target language (say English).

As the case with feed-forward neural nets, we minimize a loss function using back-propagation, where the loss is typically


where K is the number of categories for classification (e.g., size of the vocabulary in machine translation), and 𝒯[T] is the length of the output sequence. During the training, the gradients 𝒯/𝒉t are computed in the reverse time order (from T to t). For this reason, the training process is often called back-propagation through time.

One notable drawback of vanilla RNNs is that, they have difficulty in capturing long-range dependencies in sequence data when the length of the sequence is large. This is sometimes due to the phenomenon of exploding/vanishing gradients. Take Figure 8(c) as an example. Computing 𝒯/𝒉1 involves the product t=13(𝒉t+1/𝒉t) by the chain rule. However, if the sequence is long, the product will be the multiplication of many Jacobian matrices, which usually results in exponentially large or small singular values. To alleviate this issue, in practice, the forward pass and backward pass are implemented in a shorter sliding window {t1,t1+1,,t2}, instead of the full sequence {1,2,,T}. Though effective in some cases, this technique alone does not fully address the issue of long-term dependency.

3.2.2 GRUs and LSTM

There are two improved variants that alleviate the above issue: gated recurrent units (GRUs) [26] and long short-term memory (LSTM) [54].

  • A GRU refines the recursive formula (13) by introducing gates, which are vectors of the same length as 𝒉t. The gates, which take values in [0,1] elementwise, multiply with 𝒉t-1 elementwise and determine how much they keep the old hidden states.

  • An LSTM similarly uses gates in the recursive formula. In addition to 𝒉t, an LSTM maintains a cell state, which takes values in elementwise and are analogous to counters.

Here we only discuss LSTM in detail. Denote by the element-wise multiplication. We have a recursive formula in replace of (13):

(𝒊t𝒇t𝒐t𝒈t) =(𝝈𝝈𝝈𝐭𝐚𝐧𝐡)𝐖(𝒉t-1𝒙t1),
𝒄t =𝒇t𝒄t-1+𝒊t𝒈t,
𝒉t =𝒐t𝐭𝐚𝐧𝐡(𝒄t),

where 𝐖 is a big weight matrix with appropriate dimensions. The cell state vector 𝒄t carries information of the sequence (e.g., singular/plural form in a sentence). The forget gate 𝒇t determines how much the values of 𝒄t-1 are kept for time t, the input gate 𝒊t controls the amount of update to the cell state, and the output gate 𝒐t gives how much 𝒄t reveals to 𝒉t. Ideally, the elements of these gates have nearly binary values. For example, an element of 𝒇t being close to 1 may suggest the presence of a feature in the sequence data. Similar to the skip connections in residual nets, the cell state 𝒄t has an additive recursive formula, which helps back-propagation and thus captures long-range dependencies.

Figure 9: A vanilla RNN with two hidden layers. Higher-level hidden states 𝒉t are determined by the old states 𝒉t-1 and lower-level hidden states 𝒉t-1. Multilayer RNNs generalize both feed-forward neural nets and one-hidden-layer RNNs.

3.2.3 Multilayer RNNs

Multilayer RNNs are generalization of the one-hidden-layer RNN discussed above. Figure 9 shows a vanilla RNN with two hidden layers. In place of (13), the recursive formula for an RNN with L hidden layers now reads

𝒉t=𝐭𝐚𝐧𝐡[𝐖(𝒉t-1𝒉t-11)],for all[L],𝒉t0𝒙t.

Note that a multilayer RNN has two dimensions: the sequence length T and depth L. Two special cases are the feed-forward neural nets (where T=1) introduced in Section 2, and RNNs with one hidden layer (where L=1). Multilayer RNNs usually do not have very large depth (e.g., 25), since T is already very large.

Finally, we remark that CNNs, RNNs, and other neural nets can be easily combined to tackle tasks that involve different sources of input data. For example, in image captioning, the images are first processed through a CNN, and then the high-level features are fed into an RNN as inputs. Theses neural nets combined together form a large computational graph, so they can be trained using back-propagation. This generic training method provides much flexibility in various applications.

3.3 Modules

Deep neural nets are essentially composition of many nonlinear functions. A component function may be designed to have specific properties in a given task, and it can be itself resulted from composing a few simpler functions. In LSTM, we have seen that the building block consists of several intermediate variables, including cell states and forget gates that can capture long-term dependency and alleviate numerical issues.

This leads to the idea of designing modules for building more complex neural net models. Desirable modules usually have low computational costs, alleviate numerical issues in training, and lead to good statistical accuracy. Since modules and the resulting neural net models form computational graphs, training follows the same principle briefly described in Section 2.

Here, we use the examples of Inception and skip connections to illustrate the ideas behind modules. Figure 10(a) is an example of “Inception” modules used in GoogleNet [122]. As before, all the convolutional layers are followed by the ReLU activation function. The concatenation of information from filters with different sizes give the model great flexibility to capture spatial information. Note that 1×1 filters is an 1×1×d3 tensor (where d3 is the number of feature maps), so its convolutional operation does not interact with other spatial coordinates, only serving to aggregate information from different feature maps at the same coordinate. This reduces the number of parameters and speeds up the computation. Similar ideas appear in other work [78, 57].

(a) “Inception” module (b) Skip connections
Figure 10: (a) The “Inception” module from GoogleNet. Concat means combining all features maps into a tensor. (b) Skip connections are added every two layers in ResNets.

Another module, usually called skip connections, is widely used to alleviate numerical issues in very deep neural nets, with additional benefits in optimization efficiency and statistical accuracy. Training very deep neural nets are generally more difficult, but the introduction of skip connections in residual networks [50, 51] has greatly eased the task.

The high level idea of skip connections is to add an identity map to an existing nonlinear function. Let 𝐅(𝒙) be an arbitrary nonlinear function represented by a (fragment of) neural net, then the idea of skip connections is simply replacing 𝐅(𝒙) with 𝒙+𝐅(𝒙). Figure 10(b) shows a well-known structure from residual networks [50]—for every two layers, an identity map is added:

𝒙𝝈(𝒙+𝐅(𝒙))=𝝈(𝒙+𝐖𝝈(𝐖𝒙+𝒃)+𝒃), (14)

where 𝒙 can be hidden nodes from any layer and 𝐖,𝐖,𝒃,𝒃 are corresponding parameters. By repeating (namely composing) this structure throughout all layers, [50, 51] are able to train neural nets with hundreds of layers easily, which overcomes well-observed training difficulties in deep neural nets. Moreover, deep residual networks also improve statistical accuracy, as the classification error on ImageNet challenge was reduced by 46% from 2014 to 2015. As a side note, skip connections can be used flexibly. They are not restricted to the form in (14), and can be used between any pair of layers , [55].

4 Deep unsupervised learning

In supervised learning, given labelled training set {(yi,𝐱i)}, we focus on discriminative models, which essentially represents (y|𝒙) by a deep neural net f(𝒙;𝜽) with parameters 𝜽. Unsupervised learning, in contrast, aims at extracting information from unlabeled data {𝒙i}, where the labels {yi} are absent. In regard to this information, it can be a low-dimensional embedding of the data {𝒙i} or a generative model with latent variables to approximate the distribution 𝐗(𝒙). To achieve these goals, we introduce two popular unsupervised deep leaning models, namely, autoencoders and generative adversarial networks (GANs). The first one can be viewed as a dimension reduction technique, and the second as a density estimation method. DNNs are the key elements for both of these two models.

4.1 Autoencoders

Recall that in dimension reduction, the goal is to reduce the dimensionality of the data and at the same time preserve its salient features. In particular, in principal component analysis (PCA), the goal is to embed the data {𝒙i}1in into a low-dimensional space via a linear function 𝒇 such that maximum variance can be explained. Equivalently, we want to find linear functions 𝒇:dk and 𝒈:kd (kd) such that the difference between 𝒙i and 𝒈(𝒇(𝒙i)) is minimized. Formally, we let

𝒇(𝒙)=𝑾f𝒙𝒉and𝒈(𝒉)=𝑾g𝒉,where𝑾fk×d and 𝑾gd×k.

Here, for simplicity, we assume that the intercept/bias terms for 𝒇 and 𝒈 are zero. Then, PCA amounts to minimizing the quadratic loss function

minimize𝑾f,𝑾g  1ni=1n𝒙i-𝑾f𝑾g𝒙i22. (15)

It is the same as minimizing 𝐗-𝐖𝐗F2 subject to rank(𝐖)k, where 𝐗p×n is the design matrix. The solution is given by the singular value decomposition of 𝐗 [44, Thm. 2.4.8], which is exactly what PCA does. It turns out that PCA is a special case of autoencoders, which is often known as the undercomplete linear autoencoder.

More broadly, autoencoders are neural network models for (nonlinear) dimension reduction, which generalize PCA. An autoencoder has two key components, namely, the encoder function 𝒇(), which maps the input 𝒙d to a hidden code/representation 𝒉𝒇(𝒙)k, and the decoder function 𝒈(), which maps the hidden representation 𝒉 to a point 𝒈(𝒉)d. Both functions can be multilayer neural networks as (3). See Figure 11 for an illustration of autoencoders. Let (𝒙1,𝒙2) be a loss function that measures the difference between 𝒙1 and 𝒙2 in d. Similar to PCA, an autoencoder is used to find the encoder 𝒇 and decoder 𝒈 such that (𝒙,𝒈(𝒇(𝒙))) is as small as possible. Mathematically, this amounts to solving the following minimization problem

minimize𝒇,𝒈1ni=1n(𝒙i,𝒈(𝒉i))with𝒉i=𝒇(𝒙i),for all i[n]. (16)
Figure 11: First an input 𝒙 goes through the decoder 𝒇(), and we obtain its hidden representation 𝒉=𝒇(𝒙). Then, we use the decoder 𝒈() to get 𝒈(𝒉) as a reconstruction of 𝒙. Finally, the loss is determined from the difference between the original input 𝒙 and its reconstruction 𝒈(𝒇(𝒙)).

One needs to make structural assumptions on the functions 𝒇 and 𝐠 in order to find useful representations of the data, which leads to different types of autoencoders. Indeed, if no assumption is made, choosing 𝒇 and 𝒈 to be identity functions clearly minimizes the above optimization problem. To avoid this trivial solution, one natural way is to require that the encoder f maps the data onto a space with a smaller dimension, i.e., k<d. This is the undercomplete autoencoder that includes PCA as a special case. There are other structured autoencoders which add desired properties to the model such as sparsity or robustness, mainly through regularization terms. Below we present two other common types of autoencoders.

  • Sparse autoencoders. One may believe that the dimension k of the hidden code 𝒉i is larger than the input dimension d, and that 𝒉i admits a sparse representation. As with LASSO [125] or SCAD [36], one may add a regularization term to the reconstruction loss in (16) to encourage sparsity [98]. A sparse autoencoder solves

    min𝒇,𝒈1ni=1n(𝒙i,𝒈(𝒉i))loss+λ𝒉i1regularizerwith𝒉i=𝒇(𝒙i), for all i[n].

    This is similar to dictionary learning, where one aims at finding a sparse representation of input data on an overcomplete basis. Due to the imposed sparsity, the model can potentially learn useful features of the data.

  • Denoising autoencoders. One may hope that the model is robust to noise in the data: even if the input data 𝒙i are corrupted by small noise 𝝃i or miss some components (the noise level or the missing probability is typically small), an ideal autoencoder should faithfully recover the original data. A denoising autoencoder [127] achieves this robustness by explicitly building a noisy data 𝒙~i=𝒙i+𝝃i as the new input, and then solves an optimization problem similar to (16) where (𝒙i,𝒈(𝒉i)) is replaced by (𝒙i,𝒈(𝒇(𝒙~i))). A denoising autoencoder encourages the encoder/decoder to be stable in the neighborhood of an input, which is generally a good statistical property. An alternative way could be constraining f and g in the optimization problem, but that would be very difficult to optimize. Instead, sampling by adding small perturbations in the input provides a simple implementation. We shall see similar ideas in Section 6.3.3.

4.2 Generative adversarial networks

Given unlabeled data {𝒙i}1in, density estimation aims to estimate the underlying probability density function 𝑿 from which the data is generated. Both parametric and nonparametric estimators [114] have been proposed and studied under various assumptions on the underlying distribution. Different from these classical density estimators, where the density function is explicitly defined in relatively low dimension, generative adversarial networks (GANs) [46] can be categorized as an implicit density estimator in much higher dimension. The reasons are twofold: (1) GANs put more emphasis on sampling from the distribution 𝑿 than estimation; (2) GANs define the density estimation implicitly through a source distribution 𝒁 and a generator function g(), which is usually a deep neural network. We introduce GANs from the perspective of sampling from 𝑿 and later we will generalize the vanilla GANs using its relation to density estimators.

4.2.1 Sampling view of GANs

Suppose the data {𝒙i}1in at hand are all real images, and we want to generate new natural images. With this goal in mind, GAN models a zero-sum game between two players, namely, the generator 𝒢 and the discriminator 𝒟. The generator 𝒢 tries to generate fake images akin to the true images {𝒙i}1in while the discriminator 𝒟 aims at differentiating the fake ones from the true ones. Intuitively, one hopes to learn a generator 𝒢 to generate images where the best discriminator 𝒟 cannot distinguish. Therefore the payoff is higher for the generator 𝒢 if the probability of the discriminator 𝒟 getting wrong is higher, and correspondingly the payoff for the discriminator correlates positively with its ability to tell wrong from truth.

Mathematically, the generator 𝒢 consists of two components, an source distribution 𝒁 (usually a standard multivariate Gaussian distribution with hundreds of dimensions) and a function 𝐠() which maps a sample 𝒛 from 𝒁 to a point 𝐠(𝒛) living in the same space as 𝒙. For generating images, 𝐠(𝒛) would be a 3D tensor. Here 𝐠(𝒛) is the fake sample generated from 𝒢. Similarly the discriminator 𝒟 is composed of one function which takes an image 𝒙 (real or fake) and return a number d(𝒙)[0,1], the probability of 𝒙 being a real sample from 𝑿 or not. Oftentimes, both the generating function 𝐠() and the discriminating function d() are realized by deep neural networks, e.g., CNNs introduced in Section 3.1. See Figure 12 for an illustration for GANs. Denote 𝜽𝒢 and 𝜽𝒟 the parameters in 𝐠() and d(), respectively. Then GAN tries to solve the following min-max problem:

Figure 12: GANs consist of two components, a generator 𝒢 which generates fake samples and a discriminator 𝒟 which differentiate the true ones from the fake ones.
min𝜽𝒢max𝜽𝒟  𝔼𝒙𝑿[log(d(𝒙))]+𝔼𝒛𝒁[log(1-d(𝐠(𝒛)))]. (17)

Recall that d(𝒙) models the belief / probability that the discriminator thinks that 𝒙 is a true sample. Fix the parameters 𝜽𝒢 and hence the generator 𝒢 and consider the inner maximization problem. We can see that the goal of the discriminator is to maximize its ability of differentiation. Similarly, if we fix 𝜽𝒟 (and hence the discriminator), the generator tries to generate more realistic images 𝐠(𝒛) to fool the discriminator.

4.2.2 Density estimation view of GANs

Let us now take a density-estimation view of GANs. Fixing the source distribution 𝒁, any generator 𝒢 induces a distribution 𝒢 over the space of images. Removing the restrictions on d(), one can then rewrite (17) as

min𝒢maxd()  𝔼𝒙𝑿[log(d(𝒙))]+𝔼𝒙𝒢[log(1-d(𝒙))]. (18)

Observe that the inner maximization problem is solved by the likelihood ratio, i.e.


As a result, (18) can be simplified as

min𝒢  JS(𝑿𝒢), (19)

where JS() denotes the Jensen–Shannon divergence between two distributions


In words, the vanilla GAN (17) seeks a density 𝒢 that is closest to 𝑿 in terms of the Jensen–Shannon divergence. This view allows to generalize GANs to other variants, by changing the distance metric. Examples include f-GAN [90], Wasserstein GAN (W-GAN) [6], MMD GAN [75], etc. We single out the Wasserstein GAN (W-GAN) [6] to introduce due to its popularity. As the name suggests, it minimizes the Wasserstein distance between 𝑿 and 𝒢:

min𝜽𝒢WS(𝑿𝒢)=min𝜽𝒢supf:f 1-Lipschitz𝔼𝒙𝑿[f(𝒙)]-𝔼𝒙𝒢[f(𝒙)], (20)

where f() is taken over all Lipschitz functions with coefficient 1. Comparing W-GAN (20) with the original formulation of GAN (17), one finds that the Lipschitz function f in (20) corresponds to the discriminator 𝒟 in (17) in the sense that they share similar objectives to differentiate the true distribution 𝑿 from the fake one 𝒢. In the end, we would like to mention that GANs are more difficult to train than supervised deep learning models such as CNNs [109]. Apart from the training difficulty, how to evaluate GANs objectively and effectively is an ongoing research.

5 Representation power: approximation theory

Having seen the building blocks of deep learning models in the previous sections, it is natural to ask: what is the benefits of composing multiple layers of nonlinear functions. In this section, we address this question from a approximation theoretical point of view. Mathematically, letting be the space of functions representable by neural nets (NNs), how well can a function f (with certain properties) be approximated by functions in . We first revisit universal approximation theories, which are mostly developed for shallow neural nets (neural nets with a single hidden layer), and then provide recent results that demonstrate the benefits of depth in neural nets. Other notable works include Kolmogorov-Arnold superposition theorem [7, 119], and circuit complexity for neural nets [91].

5.1 Universal approximation theory for shallow NNs

The universal approximation theories study the approximation of f in a space  by a function represented by a one-hidden-layer neural net

g(𝒙)=j=1Ncjσ*(𝒘j𝒙-bj), (21)

where σ*: is certain activation function and N is the number of hidden units in the neural net. For different space and activation function σ*, there are upper bounds and lower bounds on the approximation error f-g. See [93] for a comprehensive overview. Here we present representative results.

First, as N, any continuous function f can be approximated by some g under mild conditions. Loosely speaking, this is because each component σ*(𝒘j𝒙-bj) behaves like a basis function and functions in a suitable space admits a basis expansion. Given the above heuristics, the next natural question is: what is the rate of approximation for a finite N?

Let us restrict the domain of 𝒙 to a unit ball Bd in d. For p[1,) and integer m1, consider the Lp space and the Sobolev space with standard norms


where D𝒌f denotes partial derivatives indexed by 𝒌+d. Let pm be the space of functions f in the Sobolev space with fm,p1. Note that functions in  have bounded derivatives up to m-th order, and that smoothness of functions is controlled by m (larger m means smoother). Denote by N the space of functions with the form (21). The following general upper bound is due to [85].

Theorem 1 (Theorem 2.1 in [85]).

Assume σ*:RR is such that σ* has arbitrary order derivatives in an open interval I, and that σ* is not a polynomial on I. Then, for any p[1,), d2, and integer m1,


where Cd,m,p is independent of N, the number of hidden units.

In the above theorem, the condition on σ*() is mainly technical. This upper bound is useful when the dimension d is not large. It clearly implies that the one-hidden-layer neural net is able to approximate any smooth function with enough hidden units. However, it is unclear how to find a good approximator g; nor do we have control over the magnitude of the parameters (huge weights are impractical). While increasing the number of hidden units N leads to better approximation, the exponent -m/d suggests the presence of the curse of dimensionality. The following (nearly) matching lower bound is stated in [80].

Theorem 2 (Theorem 5 in [80]).

Let p1, m1 and N2. If the activation function is the standard sigmoid function σ(t)=(1+e-t)-1, then

supfpminfgNf-gpCd,m,p(NlogN)-m/d, (22)

where Cd,m,p is independent of N.

Results for other activation functions are also obtained by [80]. Moreover, the term logN can be removed if we assume an additional continuity condition [85].

For the natural space pm of smooth functions, the exponential dependence on d in the upper and lower bounds may look unappealing. However, [12] showed that for a different function space, there is a good dimension-free approximation by the neural nets. Suppose that a function f:d has a Fourier representation

f(𝒙)=dei𝝎,𝒙f~(𝝎)𝑑𝝎, (23)

where f~(𝝎). Assume that f(𝟎)=0 and that the following quantity is finite

Cf=d𝝎2|f~(𝝎)|𝑑𝝎. (24)

[12] uncovers the following dimension-free approximation guarantee.

Theorem 3 (Proposition 1 in [12]).

Fix a C>0 and an arbitrary probability measure μ on the unit ball Bd in Rd. For every function f with CfC and every N1, there exists some gHN such that


Moreover, the coefficients of g may be restricted to satisfy j=1N|cj|2C.

The upper bound is now independent of the dimension d. However, Cf may implicitly depend on d, as the formula in (24) involves an integration over d (so for some functions Cf may depend exponentially on d). Nevertheless, this theorem does characterize an interesting function space with an improved upper bound. Details of the function space are discussed by [12]. This theorem can be generalized; see [81] for an example.

To help understand why a dimensionality-free approximation holds, let us appeal to a heuristic argument given by Monte Carlo simulations. It is well-known that Monte Carlo approximation errors are independent of dimensionality in evaluation of high-dimensional integrals. Let us generate {𝝎j}1jN randomly from a given density p() in d. Consider the approximation to (23) by

gN(𝒙)=1Nj=1Ncjei𝝎j,𝒙,cj=f~(𝝎j)p(𝝎j). (25)

Then, gN(𝒙) is a one-hidden-layer neural network with N units and the sinusoid activation function. Note that 𝔼gN(𝒙)=f(𝒙), where the expectation is taken with respect to randomness {𝝎j}. Now, by independence, we have


if 𝔼cj2<. Therefore, the rate is independent of the dimensionality d, though the constant can be.

5.2 Approximation theory for multi-layer NNs

The approximation theory for multilayer neural nets is less understood compared with neural nets with one hidden layer. Driven by the success of deep learning, there are many recent papers focusing on expressivity of deep neural nets. As studied by [124, 35, 84, 94, 15, 110, 77, 102], deep neural nets excel at representing composition of functions. This is perhaps not surprising, since deep neural nets are themselves defined by composing layers of functions. Nevertheless, it points to a new territory rarely studied in statistics before. Below we present a result based on [77, 102].

Suppose that the inputs 𝒙 have a bounded domain [-1,1]d for simplicity. As before, let σ*: be a generic function, and 𝝈*=(σ*,,σ*) be element-wise application of σ*. Consider a neural net which is similar to (3) but with scaler output: g(𝒙)=𝐖𝝈*(𝝈*(𝐖2𝝈*(𝐖1𝒙))). A unit or neuron refers to an element of vectors 𝝈*(𝐖k𝝈*(𝐖2𝝈*(𝐖1𝒙))) for any k=1,,-1. For a multivariate polynomial p, define mk(p) to be the smallest integer such that, for any ϵ>0, there exists a neural net g(𝒙) satisfying sup𝒙|p(𝒙)-g(𝒙)|<ϵ, with k hidden layers (i.e., =k+1) and no more than mk(p) neurons in total. Essentially, mk(p) is the minimum number of neurons required to approximate p arbitrarily well.

Theorem 4 (Theorem 4.1 in [102]).

Let p(𝐱) be a monomial x1r1x2r2xdrd with q=j=1drj. Suppose that σ* has derivatives of order 2q at the origin, and that they are nonzero. Then,
(i) m1(p)=j=1d(rj+1);
(ii) minkmk(p)j=1d(7log2(rj)+4).

This theorem reveals a sharp distinction between shallow networks (one hidden layer) and deep networks. To represent a monomial function, a shallow network requires exponentially many neurons in terms of the dimension d, whereas linearly many neurons suffice for a deep network (with bounded rj). The exponential dependence on d, as shown in Theorem 4(i), is resonant with the curse of dimensionality widely seen in many fields; see [30]. One may ask: how does depth help? Depth circumvents this issue, at least for certain functions, by allowing us to represent function composition efficiently. Indeed, Theorem 4(ii) offers a nice result with clear intuitions: it is known that the product of two scalar inputs can be represented using 4 neurons [77], so by composing multiple products, we can express monomials with O(d) neurons.

Recent advances in nonparametric regressions also support the idea that deep neural nets excel at representing composition of functions [15, 110]. In particular, [15] considered the nonparametric regression setting where we want to estimate a function f^n(𝒙) from i.i.d. data 𝒟n={(yi,𝒙i)}1in. If the true regression function f(𝒙) has certain hierarchical structure with intrinsic dimensionality66 6 Roughly speaking, the true regression function can be represented by a tree where each node has at most d* children. See [15] for the precise definition. d*, then the error


has an optimal minimax convergence rate O(n-2q2q+d*), rather than the usual rate O(n-2q2q+d) that depends on the ambient dimension d. Here q is the smoothness parameter. This provides another justification for deep neural nets: if data are truly hierarchical, then the quality of approximators by deep neural nets depends on the intrinsic dimensionality, which avoids the curse of dimensionality.

We point out that the approximation theory for deep learning is far from complete. For example, in Theorem 4, the condition on σ* excludes the widely used ReLU activation function, there are no constraints on the magnitude of the weights (so they can be unreasonably large).

6 Training deep neural nets

The existence of a good function approximator in the NN function class does not explain why in practice we can easily find them. In this section, we introduce standard methods, namely stochastic gradient descent (SGD) and its variants, to train deep neural networks (or to find such a good approximator). As with many statistical machine learning tasks, training DNNs follows the empirical risk minimization (ERM) paradigm which solves the following optimization problem

minimize𝜽p  n(𝜽)1ni=1n(f(𝒙i;𝜽),yi). (26)

Here (f(𝒙i;𝜽),yi) measures the discrepancy between the prediction f(𝒙i;𝜽) of the neural network and the true label yi. Correspondingly, denote by (𝜽)𝔼(𝒙,y)𝒟[(f(𝒙;𝜽),𝒚)] the out-of-sample error, where 𝒟 is the joint distribution over (y,𝒙). Solving ERM (26) for deep neural nets faces various challenges that roughly fall into the following three categories.

  • Scalability and nonconvexity. Both the sample size n and the number of parameters p can be huge for modern deep learning applications, as we have seen in Table 1. Many optimization algorithms are not practical due to the computational costs and memory constraints. What is worse, the empirical loss function n(𝜽) in deep learning is often nonconvex. It is a priori not clear whether an optimization algorithm can drive the empirical loss (26) small.

  • Numerical stability. With a large number of layers in DNNs, the magnitudes of the hidden nodes can be drastically different, which may result in the “exploding gradients” or “vanishing gradients” issue during the training process. This is because the recursive relations across layers often lead to exponentially increasing/decreasing values in both forward passes and backward passes.

  • Generalization performance. Our ultimate goal is to find a parameter 𝜽^ such that the out-of-sample error (𝜽^) is small. However, in the over-parametrized regime where p is much larger than n, the underlying neural network has the potential to fit the training data perfectly while performing poorly on the test data. To avoid this overfitting issue, proper regularization, whether explicit or implicit, is needed in the training process for the neural nets to generalize.

In the following three subsections, we discuss practical solutions/proposals to address these challenges.

6.1 Stochastic gradient descent

Stochastic gradient descent (SGD) [100] is by far the most popular optimization algorithm to solve ERM (26) for large-scale problems. It has the following simple update rule:

𝜽t+1=𝜽t-ηtG(𝜽t)  with  G(𝜽t)=(f(𝒙it;𝜽t),yit) (27)

for t=0,1,2,, where ηt>0 is the step size (or learning rate), 𝜽0p is an initial point and it is chosen randomly from {1,2,,n}. It is easy to verify that G(𝜽t) is an unbiased estimate of n(𝜽t). The advantage of SGD is clear: compared with gradient descent, which goes over the entire dataset in every update, SGD uses a single example in each update and hence is considerably more efficient in terms of both computation and memory (especially in the first few iterations).

Apart from practical benefits of SGD, how well does SGD perform theoretically in terms of minimizing n(𝜽)? We begin with the convex case, i.e., the case where the loss function is convex w.r.t. 𝜽. It is well understood in literature that with proper choices of the step sizes {ηt}, SGD is guaranteed to achieve both consistency and asymptotic normality.

  • Consistency. If (𝜽) is a strongly convex function77 7 For results on consistency and asymptotic normality, we consider the case where in each step of SGD, the stochastic gradient is computed using a fresh sample (y,𝒙) from 𝒟. This allows to view SGD as an optimization algorithm to minimize the population loss (𝜽)., then under some mild conditions88 8 One example of such condition can be constraining the second moment of the gradients: 𝔼[(𝒙i,yi;𝜽t)22]C1+C2𝜽t-𝜽*22 for some C1,C2>0. See [16] for details., learning rates that satisfy

    t=0ηt=+  and  t=0ηt2<+ (28)

    guarantee almost sure convergence to the unique minimizer 𝜽*argmin𝜽(𝜽), i.e., 𝜽ta.s.𝜽* as t [100, 64, 16, 69]. The requirements in (28) can be viewed from the perspective of bias-variance tradeoff: the first condition ensures that the iterates can reach the minimizer (controlled bias), and the second ensures that stochasticity does not prevent convergence (controlled variance).

  • Asymptotic normality. It is proved by [97] that for robust linear regression with fixed dimension p, under the choice ηt=t-1, t(𝜽t-𝜽*) is asymptotically normal under some regularity conditions (but 𝜽t is not asymptotically efficient in general). Moreover, by averaging the iterates of SGD, [96] proved that even with a larger step size ηtt-α,α(1/2,1), the averaged iterate 𝜽¯t=t-1s=1t𝜽s is asymptotic efficient for robust linear regression. These strong results show that SGD with averaging performs as well as the MLE asymptotically, in addition to its computational efficiency.

These classical results, however, fail to explain the effectiveness of SGD when dealing with nonconvex loss functions in deep learning. Admittedly, finding global minima of nonconvex functions is computationally infeasible in the worst case. Nevertheless, recent work [4, 32] bypasses the worst case scenario by focusing on losses incurred by over-parametrized deep learning models. In particular, they show that (stochastic) gradient descent converges linearly towards the global minimizer of n(𝜽) as long as the neural network is sufficiently over-parametrized. This phenomenon is formalized below.

Theorem 5 (Theorem 2 in [4]).

Let {(yi,𝐱i)}1in be a training set satisfying mini,j:ij𝐱i-𝐱j2δ>0. Consider fitting the data using a feed-forward neural network (1) with ReLU activations. Denote by L (resp. W) the depth (resp. width) of the network. Suppose that the neural network is sufficiently over-parametrized, i.e.,

W𝗉𝗈𝗅𝗒(n,L,1δ), (29)

where poly means a polynomial function. Then with high probability, running SGD (27) with certain random initialization and properly chosen step sizes yields n(𝛉t)ε in tlog1ε iterations.

Two notable features are worth mentioning: (1) first, the network under consideration is sufficiently over-parametrized (cf. (29)) in which the number of parameters is much larger than the number of samples, and (2) one needs to initialize the weight matrices to be in near-isometry such that the magnitudes of the hidden nodes do not blow up or vanish. In a nutshell, over-parametrization and random initialization together ensure that the loss function (26) has a benign landscape99 9 In [4], the loss function n(𝜽) satisfies the PL condition. around the initial point, which in turn implies fast convergence of SGD iterates.

There are certainly other challenges for vanilla SGD to train deep neural nets: (1) training algorithms are often implemented in GPUs, and therefore it is important to tailor the algorithm to the infrastructure, (2) the vanilla SGD might converge very slowly for deep neural networks, albeit good theoretical guarantees for well-behaved problems, and (3) the learning rates {ηt} can be difficult to tune in practice. To address the aforementioned challenges, three important variants of SGD, namely mini-batch SGD, momentum-based SGD, and SGD with adaptive learning rates are introduced.

6.1.1 Mini-batch SGD

Modern computational infrastructures (e.g., GPUs) can evaluate the gradient on a number (say 64) of examples as efficiently as evaluating that on a single example. To utilize this advantage, mini-batch SGD with batch size K1 forms the stochastic gradient through K random samples:

𝜽t+1=𝜽t-ηtG(𝜽t)  with  G(𝜽t)=1Kk=1K(f(𝒙itk;𝜽t),yitk), (30)

where for each 1kK, itk is sampled uniformly from {1,2,,n}. Mini-batch SGD, which is an “interpolation” between gradient descent and stochastic gradient descent, achieves the best of both worlds: (1) using 1Kn samples to estimate the gradient, one effectively reduces the variance and hence accelerates the convergence, and (2) by taking the batch size K appropriately (say 64 or 128), the stochastic gradient G(𝜽t) can be efficiently computed using the matrix computation toolboxes on GPUs.

6.1.2 Momentum-based SGD

While mini-batch SGD forms the foundation of training neural networks, it can sometimes be slow to converge due to its oscillation behavior [121]. Optimization community has long investigated how to accelerate the convergence of gradient descent, which results in a beautiful technique called momentum methods [95, 88]. Similar to gradient descent with moment, momentum-based SGD, instead of moving the iterate 𝜽t in the direction of the current stochastic gradient G(𝜽t), smooth the past (stochastic) gradients {G(𝜽t)} to stabilize the update directions. Mathematically, let 𝒗tp be the direction of update in the tth iteration, i.e.,


Here 𝒗0=G(𝜽0) and for t=1,2,

𝒗t=ρ𝒗t-1+G(𝜽t) (31)

with 0<ρ<1. A typical choice of ρ is 0.9. Notice that ρ=0 recovers the mini-batch SGD (30), where no past information of gradients is used. A simple unrolling of 𝒗t reveals that 𝒗t is actually an exponential averaging of the past gradients, i.e., 𝒗t=j=0tρt-jG(𝜽j). Compared with vanilla mini-batch SGD, the inclusion of the momentum “smoothes” the oscillation direction and accumulates the persistent descent direction. We want to emphasize that theoretical justifications of momentum in the stochastic setting is not fully understood [63, 60].

6.1.3 SGD with adaptive learning rates

In optimization, preconditioning is often used to accelerate first-order optimization algorithms. In principle, one can apply this to SGD, which yields the following update rule:

𝜽t+1=𝜽t-ηt𝑷t-1G(𝜽t) (32)

with 𝑷tp×p being a preconditioner at the t-th step. Newton’s method can be viewed as one type of preconditioning where 𝑷t=2(𝜽t). The advantages of preconditioning are two-fold: first, a good preconditioner reduces the condition number by changing the local geometry to be more homogeneous, which is amenable to fast convergence; second, a good preconditioner frees practitioners from laboring tuning of the step sizes, as is the case with Newton’s method. AdaGrad, an adaptive gradient method proposed by [33], builds a preconditioner 𝑷t based on information of the past gradients:

𝑷t={𝖽𝗂𝖺𝗀(j=0tG(𝜽t)G(𝜽t))}1/2. (33)

Since we only require the diagonal part, this preconditioner (and its inverse) can be efficiently computed in practice. In addition, investigating (32) and (33), one can see that AdaGrad adapts to the importance of each coordinate of the parameters by setting smaller learning rates for frequent features, whereas larger learning rates for those infrequent ones. In practice, one adds a small quantity δ>0 (say 10-8) to the diagonal entries to avoid singularity (numerical underflow). A notable drawback of AdaGrad is that the effective learning rate vanishes quickly along the learning process. This is because the historical sum of the gradients can only increase with time. RMSProp [52] is a popular remedy for this problem which incorporates the idea of exponential averaging:

𝑷t={𝖽𝗂𝖺𝗀(ρ𝑷t-1+(1-ρ)G(𝜽t)G(𝜽t))}1/2. (34)

Again, the decaying parameter ρ is usually set to be 0.9. Later, Adam [65, 99] combines the momentum method and adaptive learning rate and becomes the default training algorithms in many deep learning applications.

6.2 Easing numerical instability

For very deep neural networks or RNNs with long dependencies, training difficulties often arise when the values of nodes have different magnitudes or when the gradients “vanish” or “explode” during back-propagation. Here we discuss three partial solutions to alleviate this problem.

6.2.1 ReLU activation function

One useful characteristic of the ReLU function is that its derivative is either 0 or 1, and the derivative remains 1 even for a large input. This is in sharp contrast with the standard sigmoid function (1+e-t)-1 which results in a very small derivative when inputs have large magnitude. The consequence of small derivatives across many layers is that gradients tend to be “killed”, which means that gradients become approximately zero in deep nets.

The popularity of the ReLU activation function and its variants (e.g., leaky ReLU) is largely attributable to the above reason. It has been well observed that the ReLU activation function has superior training performance over the sigmoid function [68, 79].

6.2.2 Skip connections

We have introduced skip connections in Section 3.3. Why are skip connections helpful for reducing numerical instability? This structure does not introduce a larger function space, since the identity map can be also represented with ReLU activations: 𝒙=𝝈(𝒙)-𝝈(-𝒙).

One explanation is that skip connections bring ease to the training/optimization process. Suppose that we have a general nonlinear function 𝐅(𝒙;𝜽). With a skip connection, we represent the map as 𝒙+1=𝒙+𝐅(𝒙;𝜽) instead. Now the gradient 𝒙+1/𝒙 becomes

𝒙+1𝒙=𝐈+𝐅(𝒙;𝜽)𝒙  instead of  𝐅(𝒙;𝜽)𝒙, (35)

where 𝐈 is an identity matrix. By the chain rule, gradient update requires computing products of many components, e.g., 𝒙L𝒙1==1L-1𝒙+1𝒙, so it is desirable to keep the spectra (singular values) of each component 𝒙+1𝒙 close to 1. In neural nets, with skip connections, this is easily achieved if the parameters have small values; otherwise, this may not be achievable even with careful initialization and tuning. Notably, training neural nets with hundreds of layers is possible with the help of skip connections.

6.2.3 Batch normalization

Recall that in regression analysis, one often standardizes the design matrix so that the features have zero mean and unit variance. Batch normalization extends this standardization procedure from the input layer to all the hidden layers. Mathematically, fix a mini-batch of input data {(𝒙i,yi)}i, where [n]. Let 𝒉i() be the feature of the i-th example in the -th layer (=0 corresponds to the input 𝒙i). The batch normalization layer computes the normalized version of 𝒉i() via the following steps:

𝝁 1||i𝒉i(),𝝈21||i(𝒉i()-𝝁)2  and  𝒉i,norm(l)𝒉i()-𝝁𝝈.

Here all the operations are element-wise. In words, batch normalization computes the z-score for each feature over the mini-batch and use that as inputs to subsequent layers. To make it more versatile, a typical batch normalization layer has two additional learnable parameters 𝜸() and 𝜷() such that


Again denotes the element-wise multiplication. As can be seen, 𝜸() and 𝜷() set the new feature 𝒉inew(l) to have mean 𝜷() and standard deviation 𝜸(). The introduction of batch normalization makes the training of neural networks much easier and smoother. More importantly, it allows the neural nets to perform well over a large family of hyper-parameters including the number of layers, the number of hidden units, etc. At test time, the batch normalization layer needs more care. For brevity we omit the details and refer to [58].

6.3 Regularization techniques

So far we have focused on training techniques to drive the empirical loss (26) small efficiently. Here we proceed to discuss common practice to improve the generalization power of trained neural nets.

6.3.1 Weight decay

One natural regularization idea is to add an 2 penalty to the loss function. This regularization technique is known as the weight decay in deep learning. We have seen one example in (9). For general deep neural nets, the loss to optimize is nλ(𝜽)=n(𝜽)+rλ(𝜽) where


Note that the bias (intercept) terms are not penalized. If n(𝜽) is a least square loss, then regularization with weight decay gives precisely ridge regression. The penalty rλ(𝜽) is a smooth function and thus it can be also implemented efficiently with back-propagation.

6.3.2 Dropout

Dropout, introduced by [53], prevents overfitting by randomly dropping out subsets of features during training. Take the l-th layer of the feed-forward neural network as an example. Instead of propagating all the features in 𝒉() for later computations, dropout randomly omits some of its entries by


where denotes element-wise multiplication as before, and 𝗆𝖺𝗌𝗄 is a vector of Bernoulli variables with success probability p. It is sometimes useful to rescale the features 𝒉inv drop()=𝒉drop()/p, which is called inverted dropout. During training, 𝗆𝖺𝗌𝗄 are i.i.d. vectors across mini-batches and layers. However, when testing on fresh samples, dropout is disabled and the original features 𝒉() are used to compute the output label y. It has been nicely shown by [128] that for generalized linear models, dropout serves as adaptive regularization. In the simplest case of linear regression, it is equivalent to 2 regularization. Another possible way to understand the regularization effect of dropout is through the lens of bagging [45]. Since different mini-batches has different masks, dropout can be viewed as training a large ensemble of classifiers at the same time, with a further constraint that the parameters are shared. Theoretical justification remains elusive.

6.3.3 Data augmentation

Data augmentation is a technique of enlarging the dataset when we have knowledge about invariance structure of data. It implicitly increases the sample size and usually regularizes the model effectively. For example, in image classification, we have strong prior knowledge about what invariance properties a good classifier should possess. The label of an image should not be affected by translation, rotation, flipping, and even crops of the image. Hence one can augment the dataset by randomly translating, rotating and cropping the images in the original dataset.

Formally, during training we want to minimize the loss n(𝜽)=i(f(𝒙i;𝜽),yi) w.r.t. parameters 𝜽, and we know a priori that certain transformation T𝒯 where T:dd (e.g., affine transformation) should not change the category/label of a training sample. In principle, if computation costs were not a consideration, we could convert this knowledge to a constraint f𝜽(T𝒙i)=f𝜽(𝒙i),T𝒯 in the minimization formulation. Instead of solving a constrained optimization problem, data augmentation enlarges the training dataset by sampling T𝒯 and generating new data {(T𝒙i,yi)}. In this sense, data augmentation induces invariance properties through sampling, which results in a much bigger dataset than the original one.

7 Generalization power

Section 6 has focused on the in-sample/ training error obtained via SGD, but this alone does not guarantee good performance with respect to the out-of-sample/test error. The gap between the in-sample error and the out-of-sample error, namely the generalization gap, has been the focus of statistical learning theory since its birth; see [111] for an excellent introduction to this topic.

While understanding the generalization power of deep neural nets is difficult, we sample recent endeavors in this section. From a high level point of view, these approaches can be divided into two categories, namely algorithm-independent controls and algorithm-dependent controls. More specifically, algorithm-independent controls focus solely on bounding the complexity of the function class represented by certain deep neural networks. In contrast, algorithm-dependent controls take into account the algorithm (e.g., SGD) used to train the neural network.

7.1 Algorithm-independent controls: uniform convergence

The key to algorithm-independent controls is the notion of complexity of the function class parametrized by certain neural networks. Informally, as long as the complexity is not too large, the generalization gap of any function in the function class is well-controlled. However, the standard complexity measure (e.g., VC dimension [126]) is at least proportional to the number of weights in a neural network [5, 111], which fails to explain the practical success of deep learning. The caveat here is that the function class under consideration is all the functions realized by certain neural networks, with no restrictions on the size of the weights at all. On the other hand, for the class of linear functions with bounded norm, i.e., {𝒙𝒘𝒙|𝒘2M}, it is well understood that the complexity of this function class (measured in terms of the empirical Rademacher complexity) with respect to a random sample {𝒙i}1in is upper bounded by maxi𝒙i2M/n, which is independent of the number of parameters in 𝒘. This motivates researchers to investigate the complexity of norm-controlled deep neural networks1010 10 Such attempts have been made in the seminal work [13]. [89, 14, 43, 74]. Setting the stage, we introduce a few necessary notations and facts. The key object under study is the function class parametrized by the following fully-connected neural network with depth L:

L{𝒙𝑾L𝝈(𝑾L-1𝝈(𝑾2𝝈(𝑾1𝒙)))|(𝑾1,,𝑾L)𝒲}. (36)

Here (𝑾1,𝑾2,,𝑾L)𝒲 represents a certain constraint on the parameters. For instance, one can restrict the Frobenius norm of each parameter 𝑾l through the constraint 𝑾lFMF(l), where MF(l) is some positive quantity. With regard to the complexity measure, it is standard to use Rademacher complexity to control the capacity of the function class of interest.

Definition 1 (Empirical Rademacher complexity).

The empirical Rademacher complexity of a function class F w.r.t. a dataset S{𝐱i}1in is defined as

S()=𝔼𝜺[supf1ni=1nεif(𝒙i)], (37)

where 𝛆(ε1,ε2,,εn) is composed of i.i.d. Rademacher random variables, i.e., P(εi=1)=P(εi=-1)=1/2.

In words, Rademacher complexity measures the ability of the function class to fit the random noise represented by 𝜺. Intuitively, a function class with a larger Rademacher complexity is more prone to overfitting. We now formalize the connection between the empirical Rademacher complexity and the out-of-sample error; see Chapter 24 in [111].

Theorem 6.

Assume that for all fF and all (y,𝐱) we have |L(f(𝐱),y)|1. In addition, assume that for any fixed y, the univariate function L(,y) is Lipschitz with constant 1. Then with probability at least 1-δ over the sample S{(yi,𝐱i)}1ini.i.d.D, one has for all fF


In English, the generalization gap of any function f that lies in is well-controlled as long as the Rademacher complexity of is not too large. With this connection in place, we single out the following complexity bound.

Theorem 7 (Theorem 1 in [43]).

Consider the function class FL in (36), where each parameter 𝐖l has Frobenius norm at most MF(l). Further suppose that the element-wise activation function σ() is 1-Lipschitz and positive-homogeneous (i.e., σ(cx)=cσ(x) for all c0). Then the empirical Rademacher complexity (37) w.r.t. S{𝐱i}1in satisfies

S(L)maxi𝒙i24Ll=1LMF(l)n. (38)

The upper bound of the empirical Rademacher complexity (38) is in a similar vein to that of linear functions with bounded norm, i.e., maxi𝒙i2M/n, where Ll=1LMF(l) plays the role of M in the latter case. Moreover, ignoring the term L, the upper bound (38) does not depend on the size of the network in an explicit way if MF(l) sharply concentrates around 1. This reveals that the capacity of the neural network is well-controlled, regardless of the number of parameters, as long as the Frobenius norm of the parameters is bounded. Extensions to other norm constraints, e.g., spectral norm constraints, path norm constraints have been considered by [89, 14, 74, 67, 34]. This line of work improves upon traditional capacity analysis of neural networks in the over-parametrized setting, because the upper bounds derived are often size-independent. Having said this, two important remarks are in order: (1) the upper bounds (e.g., l=1LMF(l)) involve implicit dependence on the size of the weight matrix and the depth of the neural network, which is hard to characterize; (2) the upper bound on the Rademacher complexity offers a uniform bound over all functions in the function class, which is a pure statistical result. However, it stays silent about how and why standard training algorithms like SGD can obtain a function whose parameters have small norms.

7.2 Algorithm-dependent controls

In this subsection, we bring computational thinking into statistics and investigate the role of algorithms in the generalization power of deep learning. The consideration of algorithms is quite natural and well motivated: (1) local/global minima reached by different algorithms can exhibit totally different generalization behaviors due to extreme nonconvexity, which marks a huge difference from traditional models, (2) the effective capacity of neural nets is possibly not large, since a particular algorithm does not explore the entire parameter space.

These demonstrate the fact that on top of the complexity of the function class, the inherent property of the algorithm we use plays an important role in the generalization ability of deep learning. In what follows, we survey three different ways to obtain upper bounds on the generalization errors by exploiting properties of the algorithms.

7.2.1 Mean field view of neural nets

As we have emphasized, modern deep learning models are highly over-parametrized. A line of work [83, 116, 104, 25, 82, 61] approximates the ensemble of weights by an asymptotic limit as the number of hidden units tends to infinity, so that the dynamics of SGD can be studied via certain partial different equations.

More specifically, let f^(𝒙;𝜽)=N-1i=1Nσ(𝜽i𝒙) be a function given by a one-hidden-layer neural net with N hidden units, where σ() is the ReLU activation function and parameters 𝜽[𝜽1,,𝜽N]N×d are suitably randomly initialized. Consider the regression setting where we want to minimize the population risk RN(𝜽)=𝔼[(y-f^(𝒙;𝜽))2] over parameters 𝜽. A key observation is that this population risk depends on the parameters 𝜽 only through its empirical distribution, i.e., ρ^(N)=N-1i=1Nδ𝜽i where δ𝜽i is a point mass at 𝜽i. This motivates us to view express RN(𝜽) equivalently as R(ρ^(N)), where R() is a functional that maps distributions to real numbers. Running SGD for RN()—in a suitable scaling limit—results in a gradient flow on the space of distributions endowed with the Wasserstein metric that minimizes R(). It turns out that the empirical distribution ρ^k(N) of the parameters after k steps of SGD is well approximated by the gradient follow, as long as the the neural net is over-parametrized (i.e., Nd) and the number of steps is not too large. In particular, [83] have shown that under certain regularity conditions,


where ε>0 is an proxy for the step size of SGD and ρkε is the distribution of the gradient flow at time kε. In words, the out-of-sample error under 𝜽k generated by SGD is well-approximated by that of ρkε. Viewing the optimization problem from the distributional aspect greatly simplifies the problem conceptually, as the complicated optimization problem is now passed into its limit version—for this reason, this analytical approach is called the mean field perspective. In particular, [83] further demonstrated that in some simple settings, the out-of-sample error R(ρkε) of the distributional limit can be fully characterized. Nevertheless, how well does R(ρkε) perform and how fast it converges remain largely open for general problems.

7.2.2 Stability

A second way to understand the generalization ability of deep learning is through the stability of SGD. An algorithm is considered stable if a slight change of the input does not alter the output much. It has long been observed that a stable algorithm has a small generalization gap; examples include k nearest neighbors [101, 29], bagging [18, 19], etc. The precise connection between stability and generalization gap is stated by [17, 112]. In what follows, we formalize the idea of stability and its connection with the generalization gap. Let 𝒜 denote an algorithm (possibly randomized) which takes a sample S{(yi,𝒙i)}1in of size n and returns an estimated parameter 𝜽^𝒜(S). Following [49], we have the following definition for stability.

Definition 2.

An algorithm (possibly randomized) A is ε-uniformly stable with respect to the loss function L(,) if for all datasets S,S of size n which differ in at most one example, one has


Here the expectation is taken w.r.t. the randomness in the algorithm A and ε might depend on n. The loss function L(,) takes an example (say (𝐱,y)) and the estimated parameter (say A(S)) as inputs and outputs a real value.

Surprisingly, an ε-uniformly stable algorithm incurs small generalization gap in expectation, which is stated in the following lemma.

Lemma 1 (Theorem 2.2 in [49]).

Let A be ε-uniformly stable. Then the expected generalization gap is no larger than ε, i.e.,

|𝔼𝒜,S[1ni=1n(f(𝒙i;𝒜(S)),yi)-𝔼(𝒙,y)𝒟[(f(𝒙;𝒜(S)),y)]]|ε. (39)

With Lemma 1 in hand, it suffices to prove stability bound on specific algorithms. It turns out that SGD introduced in Section 6 is uniformly stable when solving smooth nonconvex functions.

Theorem 8 (Theorem 3.12 in [49]).

Assume that for any fixed (y,𝐱), the loss function L(f(x;𝛉),y), viewed as a function of 𝛉, is L-Lipschitz and β-smooth. Consider running SGD on the empirical loss function with decaying step size αtc/t, where c is some small absolute constant. Then SGD is uniformly stable with


where we have ignored the dependency on β,c and L.

Theorem 8 reveals that SGD operating on nonconvex loss functions is indeed uniformly stable as long as the number of steps T is not large compared with n. This together with Lemma 1 demonstrates the generalization ability of SGD in expectation. Nevertheless, two important limitations are worth mentioning. First, Lemma 1 provides an upper bound on the out-of-sample error in expectation, but ideally, instead of an on-average guarantee under 𝔼𝒜,S, we would like to have a high probability guarantee as in the convex case [37]. Second, controlling the generalization gap alone is not enough to achieve a small out-of-sample error, since it is unclear whether SGD can achieve a small training error within T steps.

7.2.3 Implicit regularization

In the presence of over-parametrization (number of parameters larger than the sample size), conventional wisdom informs us that we should apply some regularization techniques (e.g., 1/2 regularization) so that the model will not overfit the data. However, in practice, neural networks without explicit regularization generalize well. This phenomenon motivates researchers to look at the regularization effects introduced by training algorithms (e.g., SGD) in this over-parametrized regime. While there might exits multiple, if not infinite global minima of the empirical loss (26), it is possible that practical algorithms tend to converge to solutions with better generalization powers.

Take the underdetermined linear system 𝑿𝜽=𝒚 as a starting point. Here 𝑿n×p and 𝜽p with p much larger than n. Running gradient descent on the loss 12𝑿𝜽-𝒚22 from the origin (i.e., 𝜽0=𝟎) results in the solution with the minimum Euclidean norm, that is GD converges to

min𝜽p 𝜽2  subject to𝑿𝜽=𝒚.

In words, without any 2 regularization in the loss function, gradient descent automatically finds the solution with the least 2 norm. This phenomenon, often called as implicit regularization, not only has been empirically observed in training neural networks, but also has been theoretically understood in some simplified cases, e.g., logistic regression with separable data. In logistic regression, given a training set {(yi,𝒙i)}1in with 𝒙ip and yi{1,-1}, one aims to fit a logistic regression model by solving the following program:

min𝜽p  1ni=1n(yi𝒙i𝜽t). (40)

Here, (u)log(1+e-u) denotes the logistic loss. Further assume that the data is separable, i.e., there exists 𝜽*p such that yi𝜽*𝒙i>0 for all i. Under this condition, the loss function (40) can be arbitrarily close to zero for certain 𝜽 with 𝜽2. What happens when we minimize (40) using gradient descent? [118] uncovers a striking phenomenon.

Theorem 9 (Theorem 3 in [118]).

Consider the logistic regression (40) with separable data. If we run GD


from any initialization 𝛉0 with appropriate step size η>0, then normalized 𝛉t converges to a solution with the maximum 2 margin. That is,

limt𝜽t𝜽t2=𝜽^, (41)

where 𝛉^ is the solution to the hard margin support vector machine:

𝜽^argmin𝜽p𝜽2,subject toyi𝒙i𝜽1for all 1in. (42)

The above theorem reveals that gradient descent, when solving logistic regression with separable data, implicitly regularizes the iterates towards the 2 max margin vector (cf. (41)), without any explicit regularization as in (42). Similar results have been obtained by [62]. In addition, [47] studied algorithms other than gradient descent and showed that coordinate descent produces a solution with the maximum 1 margin.

Moving beyond logistic regression, which can be viewed as a one-layer neural net, the theoretical understanding of implicit regularization in deeper neural networks is still limited; see [48] for an illustration in deep linear convolutional neural networks.

8 Discussion

Due to space limitations, we have omitted several important deep learning models; notable examples include deep reinforcement learning [86], deep probabilistic graphical models [108], variational autoencoders [66], transfer learning [132], etc. Apart from the modeling aspect, interesting theories on generative adversarial networks [10, 11], recurrent neural networks [3], connections with kernel methods [59, 9] are also emerging. We have also omitted the inverse-problem view of deep learning where the data are assumed to be generated from a certain neural net and the goal is to recover the weights in the NN with as few examples as possible. Various algorithms (e.g., GD with spectral initialization) have been shown to recover the weights successfully in some simplified settings [134, 117, 42, 87, 23, 39].

In the end, we identify a few important directions for future research.

  • New characterization of data distributions. The success of deep learning relies on its power of efficiently representing complex functions relevant to real data. Comparatively, classical methods often have optimal guarantee if a problem has a certain known structure, such as smoothness, sparsity, and low-rankness [120, 31, 20, 24], but they are insufficient for complex data such as images. How to characterize the high-dimensional real data that can free us from known barriers, such as the curse of dimensionality is an interesting open question?

  • Understanding various computational algorithms for deep learning. As we have emphasized throughout this survey, computational algorithms (e.g., variants of SGD) play a vital role in the success of deep learning. They allow fast training of deep neural nets and probably contribute towards the good generalization behavior of deep learning in practice. Understanding these computational algorithms and devising better ones are crucial components in understanding deep learning.

  • Robustness. It has been well documented that DNNs are sensitive to small adversarial perturbations that are indistinguishable to humans [123]. This raises serious safety issues once if deploy deep learning models in applications such as self-driving cars, healthcare, etc. It is therefore crucial to refine current training practice to enhance robustness in a principled way [115].

  • Low SNRs. Arguably, for image data and audio data where the signal-to-noise ratio (SNR) is high, deep learning has achieved great success. In many other statistical problems, the SNR may be very low. For example, in financial applications, the firm characteristic and covariates may only explain a small part of the financial returns; in healthcare systems, the uncertainty of an illness may not be predicted well from a patient’s medical history. How to adapt deep learning models to excel at such tasks is an interesting direction to pursue?


We thank Ruying Bao, Yuxin Chen, Chenxi Liu, Weijie Su, Qingcan Wang and Pengkun Yang for helpful comments and discussions.


  • [1] Martín Abadi and et. al. TensorFlow: Large-scale machine learning on heterogeneous systems, 2015. Software available from tensorflow.org.
  • [2] Reza Abbasi-Asl, Yuansi Chen, Adam Bloniarz, Michael Oliver, Ben DB Willmore, Jack L Gallant, and Bin Yu. The deeptune framework for modeling and characterizing neurons in visual cortex area v4. bioRxiv, page 465534, 2018.
  • [3] Zeyuan Allen-Zhu and Yuanzhi Li. Can SGD Learn Recurrent Neural Networks with Provable Generalization? ArXiv e-prints, abs/1902.01028, 2019.
  • [4] Zeyuan Allen-Zhu, Yuanzhi Li, and Zhao Song. A convergence theory for deep learning via over-parameterization. arXiv preprint arXiv:1811.03962, 2018.
  • [5] Martin Anthony and Peter L Bartlett. Neural network learning: Theoretical foundations. cambridge university press, 2009.
  • [6] Martin Arjovsky, Soumith Chintala, and Léon Bottou. Wasserstein generative adversarial networks. 70:214–223, 06–11 Aug 2017.
  • [7] Vladimir I Arnold. On functions of three variables. Collected Works: Representations of Functions, Celestial Mechanics and KAM Theory, 1957–1965, pages 5–8, 2009.
  • [8] Sanjeev Arora and Boaz Barak. Computational complexity: a modern approach. Cambridge University Press, 2009.
  • [9] Sanjeev Arora, Simon S Du, Wei Hu, Zhiyuan Li, and Ruosong Wang. Fine-grained analysis of optimization and generalization for overparameterized two-layer neural networks. arXiv preprint arXiv:1901.08584, 2019.
  • [10] Sanjeev Arora, Rong Ge, Yingyu Liang, Tengyu Ma, and Yi Zhang. Generalization and equilibrium in generative adversarial nets (gans). In Proceedings of the 34th International Conference on Machine Learning-Volume 70, pages 224–232. JMLR. org, 2017.
  • [11] Yu Bai, Tengyu Ma, and Andrej Risteski. Approximability of discriminators implies diversity in gans. arXiv preprint arXiv:1806.10586, 2018.
  • [12] Andrew R Barron. Universal approximation bounds for superpositions of a sigmoidal function. IEEE Transactions on Information theory, 39(3):930–945, 1993.
  • [13] Peter L Bartlett. The sample complexity of pattern classification with neural networks: the size of the weights is more important than the size of the network. IEEE transactions on Information Theory, 44(2):525–536, 1998.
  • [14] Peter L Bartlett, Dylan J Foster, and Matus J Telgarsky. Spectrally-normalized margin bounds for neural networks. In I. Guyon, U. V. Luxburg, S. Bengio, H. Wallach, R. Fergus, S. Vishwanathan, and R. Garnett, editors, Advances in Neural Information Processing Systems 30, pages 6240–6249. Curran Associates, Inc., 2017.
  • [15] Benedikt Bauer and Michael Kohler. On deep learning as a remedy for the curse of dimensionality in nonparametric regression. Technical report, Technical report, 2017.
  • [16] Léon Bottou. Online learning and stochastic approximations. On-line learning in neural networks, 17(9):142, 1998.
  • [17] Olivier Bousquet and André Elisseeff. Stability and generalization. Journal of machine learning research, 2(Mar):499–526, 2002.
  • [18] Leo Breiman. Bagging predictors. Machine learning, 24(2):123–140, 1996.
  • [19] Leo Breiman et al. Heuristics of instability and stabilization in model selection. The annals of statistics, 24(6):2350–2383, 1996.
  • [20] Emmanuel J Candès and Terence Tao. The power of convex relaxation: Near-optimal matrix completion. arXiv preprint arXiv:0903.1476, 2009.
  • [21] Chensi Cao, Feng Liu, Hai Tan, Deshou Song, Wenjie Shu, Weizhong Li, Yiming Zhou, Xiaochen Bo, and Zhi Xie. Deep learning and its applications in biomedicine. Genomics, proteomics & bioinformatics, 16(1):17–32, 2018.
  • [22] Tian Qi Chen, Yulia Rubanova, Jesse Bettencourt, and David Duvenaud. Neural ordinary differential equations. arXiv preprint arXiv:1806.07366, 2018.
  • [23] Yuxin Chen, Yuejie Chi, Jianqing Fan, and Cong Ma. Gradient descent with random initialization: Fast global convergence for nonconvex phase retrieval. Mathematical Programming, pages 1–33, 2019.
  • [24] Yuxin Chen, Yuejie Chi, Jianqing Fan, Cong Ma, and Yuling Yan. Noisy matrix completion: Understanding statistical guarantees for convex relaxation via nonconvex optimization. arXiv preprint arXiv:1902.07698, 2019.
  • [25] Lenaic Chizat and Francis Bach. On the global convergence of gradient descent for over-parameterized models using optimal transport. In Advances in neural information processing systems, pages 3040–3050, 2018.
  • [26] Kyunghyun Cho, Bart Van Merriënboer, Caglar Gulcehre, Dzmitry Bahdanau, Fethi Bougares, Holger Schwenk, and Yoshua Bengio. Learning phrase representations using rnn encoder-decoder for statistical machine translation. arXiv preprint arXiv:1406.1078, 2014.
  • [27] R Dennis Cook et al. Fisher lecture: Dimension reduction in regression. Statistical Science, 22(1):1–26, 2007.
  • [28] Jeffrey De Fauw, Joseph R Ledsam, Bernardino Romera-Paredes, Stanislav Nikolov, Nenad Tomasev, Sam Blackwell, Harry Askham, Xavier Glorot, Brendan O’Donoghue, Daniel Visentin, et al. Clinically applicable deep learning for diagnosis and referral in retinal disease. Nature medicine, 24(9):1342, 2018.
  • [29] Luc Devroye and Terry Wagner. Distribution-free performance bounds for potential function rules. IEEE Transactions on Information Theory, 25(5):601–604, 1979.
  • [30] David L Donoho. High-dimensional data analysis: The curses and blessings of dimensionality. AMS math challenges lecture, 1(2000):32, 2000.
  • [31] David L Donoho and Jain M Johnstone. Ideal spatial adaptation by wavelet shrinkage. biometrika, 81(3):425–455, 1994.
  • [32] Simon S Du, Jason D Lee, Haochuan Li, Liwei Wang, and Xiyu Zhai. Gradient descent finds global minima of deep neural networks. arXiv preprint arXiv:1811.03804, 2018.
  • [33] John Duchi, Elad Hazan, and Yoram Singer. Adaptive subgradient methods for online learning and stochastic optimization. Journal of Machine Learning Research, 12(Jul):2121–2159, 2011.
  • [34] Weinan E, Chao Ma, and Qingcan Wang. A priori estimates of the population risk for residual networks. arXiv preprint arXiv:1903.02154, 2019.
  • [35] Ronen Eldan and Ohad Shamir. The power of depth for feedforward neural networks. In Conference on Learning Theory, pages 907–940, 2016.
  • [36] Jianqing Fan and Runze Li. Variable selection via nonconcave penalized likelihood and its oracle properties. Journal of the American statistical Association, 96(456):1348–1360, 2001.
  • [37] Vitaly Feldman and Jan Vondrak. High probability generalization bounds for uniformly stable algorithms with nearly optimal rate. arXiv preprint arXiv:1902.10710, 2019.
  • [38] Jerome H Friedman and Werner Stuetzle. Projection pursuit regression. Journal of the American statistical Association, 76(376):817–823, 1981.
  • [39] Haoyu Fu, Yuejie Chi, and Yingbin Liang. Local geometry of one-hidden-layer neural networks for logistic regression. arXiv preprint arXiv:1802.06463, 2018.
  • [40] Kunihiko Fukushima and Sei Miyake. Neocognitron: A self-organizing neural network model for a mechanism of visual pattern recognition. In Competition and cooperation in neural nets, pages 267–285. Springer, 1982.
  • [41] Chao Gao, Jiyi Liu, Yuan Yao, and Weizhi Zhu. Robust estimation and generative adversarial nets. arXiv preprint arXiv:1810.02030, 2018.
  • [42] Surbhi Goel, Adam Klivans, and Raghu Meka. Learning one convolutional layer with overlapping patches. arXiv preprint arXiv:1802.02547, 2018.
  • [43] Noah Golowich, Alexander Rakhlin, and Ohad Shamir. Size-independent sample complexity of neural networks. arXiv preprint arXiv:1712.06541, 2017.
  • [44] Gene H Golub and Charles F Van Loan. Matrix computations. JHU Press, 4 edition, 2013.
  • [45] Ian Goodfellow, Yoshua Bengio, and Aaron Courville. Deep Learning. MIT Press, 2016.
  • [46] Ian Goodfellow, Jean Pouget-Abadie, Mehdi Mirza, Bing Xu, David Warde-Farley, Sherjil Ozair, Aaron Courville, and Yoshua Bengio. Generative adversarial nets. In Advances in neural information processing systems, pages 2672–2680, 2014.
  • [47] Suriya Gunasekar, Jason Lee, Daniel Soudry, and Nathan Srebro. Characterizing implicit bias in terms of optimization geometry. arXiv preprint arXiv:1802.08246, 2018.
  • [48] Suriya Gunasekar, Jason D Lee, Daniel Soudry, and Nati Srebro. Implicit bias of gradient descent on linear convolutional networks. In Advances in Neural Information Processing Systems, pages 9482–9491, 2018.
  • [49] Moritz Hardt, Benjamin Recht, and Yoram Singer. Train faster, generalize better: Stability of stochastic gradient descent. arXiv preprint arXiv:1509.01240, 2015.
  • [50] Kaiming He, Xiangyu Zhang, Shaoqing Ren, and Jian Sun. Deep residual learning for image recognition. In Proceedings of the IEEE conference on computer vision and pattern recognition, pages 770–778, 2016.
  • [51] Kaiming He, Xiangyu Zhang, Shaoqing Ren, and Jian Sun. Identity mappings in deep residual networks. In European conference on computer vision, pages 630–645. Springer, 2016.
  • [52] Geoffrey Hinton, Nitish Srivastava, and Kevin Swersky. Neural networks for machine learning lecture 6a overview of mini-batch gradient descent. 2012.
  • [53] Geoffrey E Hinton, Nitish Srivastava, Alex Krizhevsky, Ilya Sutskever, and Ruslan R Salakhutdinov. Improving neural networks by preventing co-adaptation of feature detectors. arXiv preprint arXiv:1207.0580, 2012.
  • [54] Sepp Hochreiter and Jürgen Schmidhuber. Long short-term memory. Neural computation, 9(8):1735–1780, 1997.
  • [55] Gao Huang, Zhuang Liu, Laurens Van Der Maaten, and Kilian Q Weinberger. Densely connected convolutional networks. In Proceedings of the IEEE conference on computer vision and pattern recognition, pages 4700–4708, 2017.
  • [56] David H Hubel and Torsten N Wiesel. Receptive fields, binocular interaction and functional architecture in the cat’s visual cortex. The Journal of physiology, 160(1):106–154, 1962.
  • [57] Forrest N Iandola, Song Han, Matthew W Moskewicz, Khalid Ashraf, William J Dally, and Kurt Keutzer. Squeezenet: Alexnet-level accuracy with 50x fewer parameters and< 0.5 mb model size. arXiv preprint arXiv:1602.07360, 2016.
  • [58] Sergey Ioffe and Christian Szegedy. Batch normalization: Accelerating deep network training by reducing internal covariate shift. arXiv preprint arXiv:1502.03167, 2015.
  • [59] Arthur Jacot, Franck Gabriel, and Clément Hongler. Neural tangent kernel: Convergence and generalization in neural networks. In Advances in neural information processing systems, pages 8580–8589, 2018.
  • [60] Prateek Jain, Sham M Kakade, Rahul Kidambi, Praneeth Netrapalli, and Aaron Sidford. Accelerating stochastic gradient descent. arXiv preprint arXiv:1704.08227, 2017.
  • [61] Adel Javanmard, Marco Mondelli, and Andrea Montanari. Analysis of a two-layer neural network via displacement convexity. arXiv preprint arXiv:1901.01375, 2019.
  • [62] Ziwei Ji and Matus Telgarsky. Risk and parameter convergence of logistic regression. arXiv preprint arXiv:1803.07300, 2018.
  • [63] Rahul Kidambi, Praneeth Netrapalli, Prateek Jain, and Sham Kakade. On the insufficiency of existing momentum schemes for stochastic optimization. In 2018 Information Theory and Applications Workshop (ITA), pages 1–9. IEEE, 2018.
  • [64] Jack Kiefer, Jacob Wolfowitz, et al. Stochastic estimation of the maximum of a regression function. The Annals of Mathematical Statistics, 23(3):462–466, 1952.
  • [65] Diederik P Kingma and Jimmy Ba. Adam: A method for stochastic optimization. arXiv preprint arXiv:1412.6980, 2014.
  • [66] Diederik P Kingma and Max Welling. Auto-encoding variational bayes. arXiv preprint arXiv:1312.6114, 2013.
  • [67] Jason M Klusowski and Andrew R Barron. Risk bounds for high-dimensional ridge function combinations including neural networks. arXiv preprint arXiv:1607.01434, 2016.
  • [68] Alex Krizhevsky, Ilya Sutskever, and Geoffrey E Hinton. Imagenet classification with deep convolutional neural networks. In Advances in neural information processing systems, pages 1097–1105, 2012.
  • [69] Harold Kushner and G George Yin. Stochastic approximation and recursive algorithms and applications, volume 35. Springer Science & Business Media, 2003.
  • [70] Yann LeCun, Yoshua Bengio, and Geoffrey Hinton. Deep learning. nature, 521(7553):436, 2015.
  • [71] Yann LeCun, Léon Bottou, Yoshua Bengio, and Patrick Haffner. Gradient-based learning applied to document recognition. Proceedings of the IEEE, 86(11):2278–2324, 1998.
  • [72] Hao Li, Zheng Xu, Gavin Taylor, Christoph Studer, and Tom Goldstein. Visualizing the loss landscape of neural nets. In Advances in Neural Information Processing Systems, pages 6391–6401, 2018.
  • [73] Ker-Chau Li. Sliced inverse regression for dimension reduction. Journal of the American Statistical Association, 86(414):316–327, 1991.
  • [74] Xingguo Li, Junwei Lu, Zhaoran Wang, Jarvis Haupt, and Tuo Zhao. On tighter generalization bound for deep neural networks: Cnns, resnets, and beyond. arXiv preprint arXiv:1806.05159, 2018.
  • [75] Yujia Li, Kevin Swersky, and Rich Zemel. Generative moment matching networks. In International Conference on Machine Learning, pages 1718–1727, 2015.
  • [76] Tengyuan Liang. How well can generative adversarial networks (gan) learn densities: A nonparametric view. arXiv preprint arXiv:1712.08244, 2017.
  • [77] Henry W Lin, Max Tegmark, and David Rolnick. Why does deep and cheap learning work so well? Journal of Statistical Physics, 168(6):1223–1247, 2017.
  • [78] Min Lin, Qiang Chen, and Shuicheng Yan. Network in network. arXiv preprint arXiv:1312.4400, 2013.
  • [79] Andrew L Maas, Awni Y Hannun, and Andrew Y Ng. Rectifier nonlinearities improve neural network acoustic models. In Proc. icml, volume 30, page 3, 2013.
  • [80] VE Maiorov and Ron Meir. On the near optimality of the stochastic approximation of smooth functions by neural networks. Advances in Computational Mathematics, 13(1):79–103, 2000.
  • [81] Yuly Makovoz. Random approximants and neural networks. Journal of Approximation Theory, 85(1):98–109, 1996.
  • [82] Song Mei, Theodor Misiakiewicz, and Andrea Montanari. Mean-field theory of two-layers neural networks: dimension-free bounds and kernel limit. arXiv preprint arXiv:1902.06015, 2019.
  • [83] Song Mei, Andrea Montanari, and Phan-Minh Nguyen. A mean field view of the landscape of two-layer neural networks. Proceedings of the National Academy of Sciences, 115(33):E7665–E7671, 2018.
  • [84] Hrushikesh Mhaskar, Qianli Liao, and Tomaso Poggio. Learning functions: when is deep better than shallow. arXiv preprint arXiv:1603.00988, 2016.
  • [85] Hrushikesh N Mhaskar. Neural networks for optimal approximation of smooth and analytic functions. Neural computation, 8(1):164–177, 1996.
  • [86] Volodymyr Mnih, Koray Kavukcuoglu, David Silver, Andrei A Rusu, Joel Veness, Marc G Bellemare, Alex Graves, Martin Riedmiller, Andreas K Fidjeland, Georg Ostrovski, et al. Human-level control through deep reinforcement learning. Nature, 518(7540):529, 2015.
  • [87] Marco Mondelli and Andrea Montanari. On the connection between learning two-layers neural networks and tensor decomposition. arXiv preprint arXiv:1802.07301, 2018.
  • [88] Yurii E Nesterov. A method for solving the convex programming problem with convergence rate o (1/k^ 2). In Dokl. Akad. Nauk SSSR, volume 269, pages 543–547, 1983.
  • [89] Behnam Neyshabur, Ryota Tomioka, and Nathan Srebro. Norm-based capacity control in neural networks. In Conference on Learning Theory, pages 1376–1401, 2015.
  • [90] Sebastian Nowozin, Botond Cseke, and Ryota Tomioka. f-gan: Training generative neural samplers using variational divergence minimization. In Advances in Neural Information Processing Systems, pages 271–279, 2016.
  • [91] Ian Parberry. Circuit complexity and neural networks. MIT press, 1994.
  • [92] Adam Paszke, Sam Gross, Soumith Chintala, Gregory Chanan, Edward Yang, Zachary DeVito, Zeming Lin, Alban Desmaison, Luca Antiga, and Adam Lerer. Automatic differentiation in pytorch. 2017.
  • [93] Allan Pinkus. Approximation theory of the mlp model in neural networks. Acta numerica, 8:143–195, 1999.
  • [94] Tomaso Poggio, Hrushikesh Mhaskar, Lorenzo Rosasco, Brando Miranda, and Qianli Liao. Why and when can deep-but not shallow-networks avoid the curse of dimensionality: a review. International Journal of Automation and Computing, 14(5):503–519, 2017.
  • [95] Boris T Polyak. Some methods of speeding up the convergence of iteration methods. USSR Computational Mathematics and Mathematical Physics, 4(5):1–17, 1964.
  • [96] Boris T Polyak and Anatoli B Juditsky. Acceleration of stochastic approximation by averaging. SIAM Journal on Control and Optimization, 30(4):838–855, 1992.
  • [97] Boris Teodorovich Polyak and Yakov Zalmanovich Tsypkin. Adaptive estimation algorithms: convergence, optimality, stability. Avtomatika i Telemekhanika, (3):71–84, 1979.
  • [98] Christopher Poultney, Sumit Chopra, Yann L Cun, et al. Efficient learning of sparse representations with an energy-based model. In Advances in neural information processing systems, pages 1137–1144, 2007.
  • [99] Sashank J Reddi, Satyen Kale, and Sanjiv Kumar. On the convergence of adam and beyond. 2018.
  • [100] Herbert Robbins and Sutton Monro. A stochastic approximation method. The Annals of Mathematical Statistics, 22(3):400–407, 1951.
  • [101] William H Rogers and Terry J Wagner. A finite sample distribution-free performance bound for local discrimination rules. The Annals of Statistics, pages 506–514, 1978.
  • [102] David Rolnick and Max Tegmark. The power of deeper networks for expressing natural functions. arXiv preprint arXiv:1705.05502, 2017.
  • [103] Yaniv Romano, Matteo Sesia, and Emmanuel J Candès. Deep knockoffs. arXiv preprint arXiv:1811.06687, 2018.
  • [104] Grant M Rotskoff and Eric Vanden-Eijnden. Neural networks as interacting particle systems: Asymptotic convexity of the loss landscape and universal scaling of the approximation error. arXiv preprint arXiv:1805.00915, 2018.
  • [105] David E Rumelhart, Geoffrey E Hinton, and Ronald J Williams. Learning internal representations by error propagation. Technical report, California Univ San Diego La Jolla Inst for Cognitive Science, 1985.
  • [106] Olga Russakovsky, Jia Deng, Hao Su, Jonathan Krause, Sanjeev Satheesh, Sean Ma, Zhiheng Huang, Andrej Karpathy, Aditya Khosla, Michael Bernstein, Alexander C. Berg, and Li Fei-Fei. ImageNet Large Scale Visual Recognition Challenge. International Journal of Computer Vision (IJCV), 115(3):211–252, 2015.
  • [107] Haşim Sak, Andrew Senior, and Françoise Beaufays. Long short-term memory recurrent neural network architectures for large scale acoustic modeling. In Fifteenth annual conference of the international speech communication association, 2014.
  • [108] Ruslan Salakhutdinov and Geoffrey Hinton. Deep boltzmann machines. In Artificial intelligence and statistics, pages 448–455, 2009.
  • [109] Tim Salimans, Ian Goodfellow, Wojciech Zaremba, Vicki Cheung, Alec Radford, and Xi Chen. Improved techniques for training gans. In Advances in Neural Information Processing Systems, pages 2234–2242, 2016.
  • [110] Johannes Schmidt-Hieber. Nonparametric regression using deep neural networks with relu activation function. arXiv preprint arXiv:1708.06633, 2017.
  • [111] Shai Shalev-Shwartz and Shai Ben-David. Understanding machine learning: From theory to algorithms. Cambridge university press, 2014.
  • [112] Shai Shalev-Shwartz, Ohad Shamir, Nathan Srebro, and Karthik Sridharan. Learnability, stability and uniform convergence. Journal of Machine Learning Research, 11(Oct):2635–2670, 2010.
  • [113] David Silver, Julian Schrittwieser, Karen Simonyan, Ioannis Antonoglou, Aja Huang, Arthur Guez, Thomas Hubert, Lucas Baker, Matthew Lai, Adrian Bolton, et al. Mastering the game of go without human knowledge. Nature, 550(7676):354, 2017.
  • [114] Bernard W Silverman. Density estimation for statistics and data analysis. Chapman & Hall, CRC, 1998.
  • [115] Chandan Singh, W James Murdoch, and Bin Yu. Hierarchical interpretations for neural network predictions. arXiv preprint arXiv:1806.05337, 2018.
  • [116] Justin Sirignano and Konstantinos Spiliopoulos. Mean field analysis of neural networks. arXiv preprint arXiv:1805.01053, 2018.
  • [117] Mahdi Soltanolkotabi. Learning relus via gradient descent. In Advances in Neural Information Processing Systems, pages 2007–2017, 2017.
  • [118] Daniel Soudry, Elad Hoffer, Mor Shpigel Nacson, Suriya Gunasekar, and Nathan Srebro. The implicit bias of gradient descent on separable data. The Journal of Machine Learning Research, 19(1):2822–2878, 2018.
  • [119] David A Sprecher. On the structure of continuous functions of several variables. Transactions of the American Mathematical Society, 115:340–355, 1965.
  • [120] Charles J Stone. Optimal global rates of convergence for nonparametric regression. The annals of statistics, pages 1040–1053, 1982.
  • [121] Ilya Sutskever, James Martens, George Dahl, and Geoffrey Hinton. On the importance of initialization and momentum in deep learning. In International conference on machine learning, pages 1139–1147, 2013.
  • [122] Christian Szegedy, Wei Liu, Yangqing Jia, Pierre Sermanet, Scott Reed, Dragomir Anguelov, Dumitru Erhan, Vincent Vanhoucke, and Andrew Rabinovich. Going deeper with convolutions. In Proceedings of the IEEE conference on computer vision and pattern recognition, pages 1–9, 2015.
  • [123] Christian Szegedy, Wojciech Zaremba, Ilya Sutskever, Joan Bruna, Dumitru Erhan, Ian Goodfellow, and Rob Fergus. Intriguing properties of neural networks. arXiv preprint arXiv:1312.6199, 2013.
  • [124] Matus Telgarsky. Benefits of depth in neural networks. arXiv preprint arXiv:1602.04485, 2016.
  • [125] Robert Tibshirani. Regression shrinkage and selection via the lasso. Journal of the Royal Statistical Society: Series B (Methodological), 58(1):267–288, 1996.
  • [126] VN Vapnik and A Ya Chervonenkis. On the uniform convergence of relative frequencies of events to their probabilities. Theory of Probability & Its Applications, 16(2):264–280, 1971.
  • [127] Pascal Vincent, Hugo Larochelle, Yoshua Bengio, and Pierre-Antoine Manzagol. Extracting and composing robust features with denoising autoencoders. In Proceedings of the 25th international conference on Machine learning, pages 1096–1103. ACM, 2008.
  • [128] Stefan Wager, Sida Wang, and Percy S Liang. Dropout training as adaptive regularization. In Advances in neural information processing systems, pages 351–359, 2013.
  • [129] E Weinan, Jiequn Han, and Arnulf Jentzen. Deep learning-based numerical methods for high-dimensional parabolic partial differential equations and backward stochastic differential equations. Communications in Mathematics and Statistics, 5(4):349–380, 2017.
  • [130] Ashia C Wilson, Rebecca Roelofs, Mitchell Stern, Nati Srebro, and Benjamin Recht. The marginal value of adaptive gradient methods in machine learning. In I. Guyon, U. V. Luxburg, S. Bengio, H. Wallach, R. Fergus, S. Vishwanathan, and R. Garnett, editors, Advances in Neural Information Processing Systems 30, pages 4148–4158. Curran Associates, Inc., 2017.
  • [131] Yonghui Wu, Mike Schuster, Zhifeng Chen, Quoc V Le, Mohammad Norouzi, Wolfgang Macherey, Maxim Krikun, Yuan Cao, Qin Gao, Klaus Macherey, et al. Google’s neural machine translation system: Bridging the gap between human and machine translation. arXiv preprint arXiv:1609.08144, 2016.
  • [132] Jason Yosinski, Jeff Clune, Yoshua Bengio, and Hod Lipson. How transferable are features in deep neural networks? In Advances in neural information processing systems, pages 3320–3328, 2014.
  • [133] Jason Yosinski, Jeff Clune, Anh Nguyen, Thomas Fuchs, and Hod Lipson. Understanding neural networks through deep visualization. arXiv preprint arXiv:1506.06579, 2015.
  • [134] Kai Zhong, Zhao Song, Prateek Jain, Peter L Bartlett, and Inderjit S Dhillon. Recovery guarantees for one-hidden-layer neural networks. In Proceedings of the 34th International Conference on Machine Learning-Volume 70, pages 4140–4149. JMLR. org, 2017.