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.
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
- 2 Feed-forward neural networks
- 3 Popular models
- 4 Deep unsupervised learning
- 5 Representation power: approximation theory
- 6 Training deep neural nets
- 7 Generalization power
- 8 Discussion
Modern machine learning and statistics deal with the problem of learning from data: given a training dataset where is the input and is the output11 1 When the label 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 is not given., one seeks a function 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 ) corresponds to the raw image (reps. its category) and the goal is to find a mapping that can classify future images accurately. Decades of research efforts in statistical machine learning have been devoted to developing methods to find 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 is chosen.
Deep learning , in its simplest form, proposes the following compositional function class:
Here, for each , is some nonlinear function, and 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 , game playing , and machine translation . 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 , etc.
|Model||Year||# Layers||# Params||Top-5 error|
To get a better idea of the success of deep learning, let us take the ImageNet Challenge  (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 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 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?
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  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.
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  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  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 . 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 . 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 . 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 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|
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 , denote by the expected risk on fresh sample (a.k.a. test error, prediction error or generalization error), and by the empirical risk/training error averaged over a training dataset. Arguably, the key theoretical question in deep learning is
why is small, where is the classifier returned by the training algorithm?
We follow the conventional approximation-estimation decomposition (sometimes, also bias-variance tradeoff) to decompose the term into two parts. Let be the function space expressible by a family of neural nets. Define to be the best possible classifier and to be the best classifier in . Then, we can decompose the excess error into two parts:
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., ; see e.g. left panel of Figure 2). Moreover, its generalization error 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., 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 2–4, 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 where and are independent across , supervised learning aims at finding a (possibly random) function that predicts the outcome for a new input , assuming follows the same distribution as . In the terminology of machine learning, the input is often called the feature, the output called the label, and the pair is an example. The function is called the classifier, and estimation of is training or learning. The performance of is evaluated through the prediction error , which can be often estimated from a separate test dataset.
As with classical statistical estimation, for each , a classifier approximates the conditional probability using a function parametrized by . Then the category with the highest probability is predicted. Thus, learning is essentially estimating the parameters . In statistics, one of the most popular methods is (multinomial) logistic regression, which stipulates a specific form for the functions : let and where is a normalization factor to make a valid probability distribution. It is clear that logistic regression induces linear decision boundaries in , and hence it is restrictive in modeling nonlinear dependency between 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 is the number of hidden layers, and is usually called depth of a NN model. Letting , one can recursively define for all . The feed-forward neural networks, also called the multilayer perceptrons (MLPs), are neural nets with a specific choice of : for , define
where and are the weight matrix and the bias/intercept, respectively, associated with the -th layer, and is usually a simple given (known) nonlinear function called the activation function. In words, in each layer , the input vector 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:
Other choices of activation functions include leaky ReLU, function  and the classical sigmoid function , which is less used now.
Given an output from the final hidden layer and a label , 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 goes through an affine transformation and then the soft-max function:
Then the loss is defined to be the cross-entropy between the label (in the form of an indicator vector) and the score vector , which is exactly the negative log-likelihood of the multinomial logistic regression model:
where . As a final remark, the number of parameters scales with both the depth 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 and then iteratively updates the parameters by moving it in the direction of the negative gradient. The difference is that, in each update, a small subsample 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 . 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
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  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 , then , , and finally . 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
where denotes a diagonal matrix with elements given by the argument. Note that the calculation of 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
where if the -th element of is nonnegative, and otherwise. The step size , also called the learning rate, controls how much parameters are changed in a single update.
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 regularization:
where is the same as (6), and is a tuning parameter. A similar example is considered in . 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 compositions: first the input data multiplies the weight matrix resulting in , then it goes through the ReLU activation function relu resulting in , then it multiplies another weight matrix leading to , and finally it produces the cross-entropy with label 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, ) 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 (). 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.
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 . Here, the first two dimensions of indicate spatial coordinates of an image while the third indicates the number of channels. For instance, is for the raw inputs due to the red, green and blue channels, and 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.
Convolutional layer (CONV). A convolutional layer has the same functionality as described in (3), where the input feature 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 (), where is the size of the filter (typically 3 or 5) and denotes the total number of filters. Note that the third dimension of is equal to that of the input feature . For this reason, one usually says that the filter has size , while suppressing the third dimension . Each filter then convolves with the input feature to obtain one single feature map , where44 4 To simplify notation, we omit the bias/intercept term associated with each filter.
Here is a small “patch” of starting at location . See Figure 5 for an illustration of the convolution operation. If we view the 3D tensors and as vectors, then each filter essentially computes their inner product with a part of indexed by (which can be also viewed as convolution, as its name suggests). One then pack the resulted feature maps into a 3D tensor with size , where
The outputs of convolutional layers are then followed by nonlinear activation functions. In the ReLU case, we have
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 are shared across all locations . A patch of an input responds strongly (that is, producing a large value) to a filter if they are positively correlated. Therefore intuitively, each filter serves to extract features similar to .
As a side note, after the convolution (10), the spatial size of the input shrinks to 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 . In addition, a stride in the convolutional layer determines the gap and between two patches and : in (10) the stride is , and a larger stride would lead to feature maps with smaller sizes.
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 max-pooling filter. It computes , that is, the maximum of the neighborhood in the spatial coordinates; see Figure 6 for an illustration. Note that the pooling operation is done separately for each feature map . As a consequence, a max-pooling filter acting on will result in an output of size . 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 in a feature map is likely to contain the same information.
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 , and computes . 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 , which is composed of two sets of CONV-POOL layers and three fully-connected layers.
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 , machine translation , genome sequencing , 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|
3.2.1 Vanilla RNNs
Suppose we have general time series inputs . A vanilla RNN models the “hidden state” at time by a vector , which is subject to the recursive formula
Here, 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.
where are trainable weight matrices, are trainable bias vectors, and is the output at time . 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 is the number of categories for classification (e.g., size of the vocabulary in machine translation), and is the length of the output sequence. During the training, the gradients are computed in the reverse time order (from to ). 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 involves the product 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 , instead of the full sequence . Though effective in some cases, this technique alone does not fully address the issue of long-term dependency.
3.2.2 GRUs and LSTM
A GRU refines the recursive formula (13) by introducing gates, which are vectors of the same length as . The gates, which take values in elementwise, multiply with elementwise and determine how much they keep the old hidden states.
An LSTM similarly uses gates in the recursive formula. In addition to , 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):
where is a big weight matrix with appropriate dimensions. The cell state vector carries information of the sequence (e.g., singular/plural form in a sentence). The forget gate determines how much the values of are kept for time , the input gate controls the amount of update to the cell state, and the output gate gives how much reveals to . Ideally, the elements of these gates have nearly binary values. For example, an element of being close to may suggest the presence of a feature in the sequence data. Similar to the skip connections in residual nets, the cell state has an additive recursive formula, which helps back-propagation and thus captures long-range dependencies.
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 hidden layers now reads
Note that a multilayer RNN has two dimensions: the sequence length and depth . Two special cases are the feed-forward neural nets (where ) introduced in Section 2, and RNNs with one hidden layer (where ). Multilayer RNNs usually do not have very large depth (e.g., –), since 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.
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 . 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 filters is an tensor (where 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|
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 —for every two layers, an identity map is added:
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 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 .
4 Deep unsupervised learning
In supervised learning, given labelled training set , we focus on discriminative models, which essentially represents by a deep neural net with parameters . Unsupervised learning, in contrast, aims at extracting information from unlabeled data , where the labels are absent. In regard to this information, it can be a low-dimensional embedding of the data 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.
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 into a low-dimensional space via a linear function such that maximum variance can be explained. Equivalently, we want to find linear functions and () such that the difference between and is minimized. Formally, we let
Here, for simplicity, we assume that the intercept/bias terms for and are zero. Then, PCA amounts to minimizing the quadratic loss function
It is the same as minimizing subject to , where 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 to a hidden code/representation , and the decoder function , which maps the hidden representation to a point . Both functions can be multilayer neural networks as (3). See Figure 11 for an illustration of autoencoders. Let be a loss function that measures the difference between and in . 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
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 maps the data onto a space with a smaller dimension, i.e., . 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 of the hidden code is larger than the input dimension , and that admits a sparse representation. As with LASSO  or SCAD , one may add a regularization term to the reconstruction loss in (16) to encourage sparsity . A sparse autoencoder solves
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 are corrupted by small noise 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  achieves this robustness by explicitly building a noisy data as the new input, and then solves an optimization problem similar to (16) where is replaced by . 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 and 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 , density estimation aims to estimate the underlying probability density function from which the data is generated. Both parametric and nonparametric estimators  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)  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 , 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 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 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 , the probability of being a real sample from or not. Oftentimes, both the generating function and the discriminating function 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 , respectively. Then GAN tries to solve the following min-max problem:
Recall that 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 , one can then rewrite (17) as
Observe that the inner maximization problem is solved by the likelihood ratio, i.e.
As a result, (18) can be simplified as
where 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 , Wasserstein GAN (W-GAN) , MMD GAN , etc. We single out the Wasserstein GAN (W-GAN)  to introduce due to its popularity. As the name suggests, it minimizes the Wasserstein distance between and :
where 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 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 . 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 (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 .
5.1 Universal approximation theory for shallow NNs
The universal approximation theories study the approximation of in a space by a function represented by a one-hidden-layer neural net
where is certain activation function and 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 . See  for a comprehensive overview. Here we present representative results.
First, as , any continuous function can be approximated by some under mild conditions. Loosely speaking, this is because each component 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 ?
Let us restrict the domain of to a unit ball in . For and integer , consider the space and the Sobolev space with standard norms
where denotes partial derivatives indexed by . Let be the space of functions in the Sobolev space with . Note that functions in have bounded derivatives up to -th order, and that smoothness of functions is controlled by (larger means smoother). Denote by the space of functions with the form (21). The following general upper bound is due to .
Theorem 1 (Theorem 2.1 in ).
Assume is such that has arbitrary order derivatives in an open interval , and that is not a polynomial on . Then, for any , , and integer ,
where is independent of , the number of hidden units.
In the above theorem, the condition on is mainly technical. This upper bound is useful when the dimension 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 ; nor do we have control over the magnitude of the parameters (huge weights are impractical). While increasing the number of hidden units leads to better approximation, the exponent suggests the presence of the curse of dimensionality. The following (nearly) matching lower bound is stated in .
Theorem 2 (Theorem 5 in ).
Let , and . If the activation function is the standard sigmoid function , then
where is independent of .
For the natural space of smooth functions, the exponential dependence on in the upper and lower bounds may look unappealing. However,  showed that for a different function space, there is a good dimension-free approximation by the neural nets. Suppose that a function has a Fourier representation
where . Assume that and that the following quantity is finite
 uncovers the following dimension-free approximation guarantee.
Theorem 3 (Proposition 1 in ).
Fix a and an arbitrary probability measure on the unit ball in . For every function with and every , there exists some such that
Moreover, the coefficients of may be restricted to satisfy .
The upper bound is now independent of the dimension . However, may implicitly depend on , as the formula in (24) involves an integration over (so for some functions may depend exponentially on ). Nevertheless, this theorem does characterize an interesting function space with an improved upper bound. Details of the function space are discussed by . This theorem can be generalized; see  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 randomly from a given density in . Consider the approximation to (23) by
Then, is a one-hidden-layer neural network with units and the sinusoid activation function. Note that , where the expectation is taken with respect to randomness . Now, by independence, we have
if . Therefore, the rate is independent of the dimensionality , 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 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: . A unit or neuron refers to an element of vectors for any . For a multivariate polynomial , define to be the smallest integer such that, for any , there exists a neural net satisfying , with hidden layers (i.e., ) and no more than neurons in total. Essentially, is the minimum number of neurons required to approximate arbitrarily well.
Theorem 4 (Theorem 4.1 in ).
Let be a monomial with . Suppose that has derivatives of order at the origin, and that they are nonzero. Then,
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 , whereas linearly many neurons suffice for a deep network (with bounded ). The exponential dependence on , as shown in Theorem 4(i), is resonant with the curse of dimensionality widely seen in many fields; see . 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 neurons , so by composing multiple products, we can express monomials with neurons.
Recent advances in nonparametric regressions also support the idea that deep neural nets excel at representing composition of functions [15, 110]. In particular,  considered the nonparametric regression setting where we want to estimate a function from i.i.d. data . If the true regression function 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 children. See  for the precise definition. , then the error
has an optimal minimax convergence rate , rather than the usual rate that depends on the ambient dimension . Here 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
Here measures the discrepancy between the prediction of the neural network and the true label . Correspondingly, denote by the out-of-sample error, where is the joint distribution over . 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 and the number of parameters 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 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 is much larger than , 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
for , where is the step size (or learning rate), is an initial point and is chosen randomly from . It is easy to verify that is an unbiased estimate of . 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 ? 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 , 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 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: for some . See  for details., learning rates that satisfy
guarantee almost sure convergence to the unique minimizer , i.e., as [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  that for robust linear regression with fixed dimension , under the choice , is asymptotically normal under some regularity conditions (but is not asymptotically efficient in general). Moreover, by averaging the iterates of SGD,  proved that even with a larger step size , the averaged iterate 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 as long as the neural network is sufficiently over-parametrized. This phenomenon is formalized below.
Theorem 5 (Theorem 2 in ).
Let be a training set satisfying . Consider fitting the data using a feed-forward neural network (1) with ReLU activations. Denote by (resp. ) the depth (resp. width) of the network. Suppose that the neural network is sufficiently over-parametrized, i.e.,
where means a polynomial function. Then with high probability, running SGD (27) with certain random initialization and properly chosen step sizes yields in 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 , the loss function 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 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 forms the stochastic gradient through random samples:
where for each , is sampled uniformly from . Mini-batch SGD, which is an “interpolation” between gradient descent and stochastic gradient descent, achieves the best of both worlds: (1) using samples to estimate the gradient, one effectively reduces the variance and hence accelerates the convergence, and (2) by taking the batch size appropriately (say 64 or 128), the stochastic gradient 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 . 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 in the direction of the current stochastic gradient , smooth the past (stochastic) gradients to stabilize the update directions. Mathematically, let be the direction of update in the th iteration, i.e.,
Here and for
with . A typical choice of is 0.9. Notice that recovers the mini-batch SGD (30), where no past information of gradients is used. A simple unrolling of reveals that is actually an exponential averaging of the past gradients, i.e., 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:
with being a preconditioner at the -th step. Newton’s method can be viewed as one type of preconditioning where . 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 , builds a preconditioner based on information of the past gradients:
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 (say ) 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  is a popular remedy for this problem which incorporates the idea of exponential averaging:
Again, the decaying parameter is usually set to be . 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 or , and the derivative remains even for a large input. This is in sharp contrast with the standard sigmoid function 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.
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 instead. Now the gradient becomes
where is an identity matrix. By the chain rule, gradient update requires computing products of many components, e.g., , so it is desirable to keep the spectra (singular values) of each component close to . 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 , where . Let be the feature of the -th example in the -th layer ( corresponds to the input ). The batch normalization layer computes the normalized version of via the following steps:
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 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 .
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 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 where
Note that the bias (intercept) terms are not penalized. If is a least square loss, then regularization with weight decay gives precisely ridge regression. The penalty is a smooth function and thus it can be also implemented efficiently with back-propagation.
Dropout, introduced by , prevents overfitting by randomly dropping out subsets of features during training. Take the -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 . It is sometimes useful to rescale the features , 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 . It has been nicely shown by  that for generalized linear models, dropout serves as adaptive regularization. In the simplest case of linear regression, it is equivalent to regularization. Another possible way to understand the regularization effect of dropout is through the lens of bagging . 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 w.r.t. parameters , and we know a priori that certain transformation where (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 in the minimization formulation. Instead of solving a constrained optimization problem, data augmentation enlarges the training dataset by sampling and generating new data . 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  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 ) 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., , 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 is upper bounded by , 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 . [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 :
Here represents a certain constraint on the parameters. For instance, one can restrict the Frobenius norm of each parameter through the constraint , where 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 w.r.t. a dataset is defined as
where is composed of i.i.d. Rademacher random variables, i.e., .
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 .
Assume that for all and all we have . In addition, assume that for any fixed , the univariate function is Lipschitz with constant 1. Then with probability at least over the sample , one has for all
In English, the generalization gap of any function 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 ).
The upper bound of the empirical Rademacher complexity (38) is in a similar vein to that of linear functions with bounded norm, i.e., , where plays the role of in the latter case. Moreover, ignoring the term , the upper bound (38) does not depend on the size of the network in an explicit way if sharply concentrates around . 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., ) 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 be a function given by a one-hidden-layer neural net with hidden units, where is the ReLU activation function and parameters are suitably randomly initialized. Consider the regression setting where we want to minimize the population risk over parameters . A key observation is that this population risk depends on the parameters only through its empirical distribution, i.e., where is a point mass at . This motivates us to view express equivalently as , where is a functional that maps distributions to real numbers. Running SGD for —in a suitable scaling limit—results in a gradient flow on the space of distributions endowed with the Wasserstein metric that minimizes . It turns out that the empirical distribution of the parameters after steps of SGD is well approximated by the gradient follow, as long as the the neural net is over-parametrized (i.e., ) and the number of steps is not too large. In particular,  have shown that under certain regularity conditions,
where is an proxy for the step size of SGD and is the distribution of the gradient flow at time . In words, the out-of-sample error under generated by SGD is well-approximated by that of . 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,  further demonstrated that in some simple settings, the out-of-sample error of the distributional limit can be fully characterized. Nevertheless, how well does perform and how fast it converges remain largely open for general problems.
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 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 of size and returns an estimated parameter . Following , we have the following definition for stability.
An algorithm (possibly randomized) is -uniformly stable with respect to the loss function if for all datasets of size which differ in at most one example, one has
Here the expectation is taken w.r.t. the randomness in the algorithm and might depend on . The loss function takes an example (say ) and the estimated parameter (say ) 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 ).
Let be -uniformly stable. Then the expected generalization gap is no larger than , i.e.,
Theorem 8 (Theorem 3.12 in ).
Assume that for any fixed , the loss function , viewed as a function of , is -Lipschitz and -smooth. Consider running SGD on the empirical loss function with decaying step size , where is some small absolute constant. Then SGD is uniformly stable with
where we have ignored the dependency on and .
Theorem 8 reveals that SGD operating on nonconvex loss functions is indeed uniformly stable as long as the number of steps is not large compared with . 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 , we would like to have a high probability guarantee as in the convex case . 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 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., 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 and with much larger than . Running gradient descent on the loss from the origin (i.e., ) results in the solution with the minimum Euclidean norm, that is GD converges to
In words, without any regularization in the loss function, gradient descent automatically finds the solution with the least 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 with and , one aims to fit a logistic regression model by solving the following program:
Here, denotes the logistic loss. Further assume that the data is separable, i.e., there exists such that for all . Under this condition, the loss function (40) can be arbitrarily close to zero for certain with . What happens when we minimize (40) using gradient descent?  uncovers a striking phenomenon.
Theorem 9 (Theorem 3 in ).
Consider the logistic regression (40) with separable data. If we run GD
from any initialization with appropriate step size , then normalized converges to a solution with the maximum margin. That is,
where is the solution to the hard margin support vector machine:
The above theorem reveals that gradient descent, when solving logistic regression with separable data, implicitly regularizes the iterates towards the max margin vector (cf. (41)), without any explicit regularization as in (42). Similar results have been obtained by . In addition,  studied algorithms other than gradient descent and showed that coordinate descent produces a solution with the maximum 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  for an illustration in deep linear convolutional neural networks.
Due to space limitations, we have omitted several important deep learning models; notable examples include deep reinforcement learning , deep probabilistic graphical models , variational autoencoders , transfer learning , etc. Apart from the modeling aspect, interesting theories on generative adversarial networks [10, 11], recurrent neural networks , 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 . 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 .
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.
-  Martín Abadi and et. al. TensorFlow: Large-scale machine learning on heterogeneous systems, 2015. Software available from tensorflow.org.
-  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.
-  Zeyuan Allen-Zhu and Yuanzhi Li. Can SGD Learn Recurrent Neural Networks with Provable Generalization? ArXiv e-prints, abs/1902.01028, 2019.
-  Zeyuan Allen-Zhu, Yuanzhi Li, and Zhao Song. A convergence theory for deep learning via over-parameterization. arXiv preprint arXiv:1811.03962, 2018.
-  Martin Anthony and Peter L Bartlett. Neural network learning: Theoretical foundations. cambridge university press, 2009.
-  Martin Arjovsky, Soumith Chintala, and Léon Bottou. Wasserstein generative adversarial networks. 70:214–223, 06–11 Aug 2017.
-  Vladimir I Arnold. On functions of three variables. Collected Works: Representations of Functions, Celestial Mechanics and KAM Theory, 1957–1965, pages 5–8, 2009.
-  Sanjeev Arora and Boaz Barak. Computational complexity: a modern approach. Cambridge University Press, 2009.
-  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.
-  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.
-  Yu Bai, Tengyu Ma, and Andrej Risteski. Approximability of discriminators implies diversity in gans. arXiv preprint arXiv:1806.10586, 2018.
-  Andrew R Barron. Universal approximation bounds for superpositions of a sigmoidal function. IEEE Transactions on Information theory, 39(3):930–945, 1993.
-  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.
-  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.
-  Benedikt Bauer and Michael Kohler. On deep learning as a remedy for the curse of dimensionality in nonparametric regression. Technical report, Technical report, 2017.
-  Léon Bottou. Online learning and stochastic approximations. On-line learning in neural networks, 17(9):142, 1998.
-  Olivier Bousquet and André Elisseeff. Stability and generalization. Journal of machine learning research, 2(Mar):499–526, 2002.
-  Leo Breiman. Bagging predictors. Machine learning, 24(2):123–140, 1996.
-  Leo Breiman et al. Heuristics of instability and stabilization in model selection. The annals of statistics, 24(6):2350–2383, 1996.
-  Emmanuel J Candès and Terence Tao. The power of convex relaxation: Near-optimal matrix completion. arXiv preprint arXiv:0903.1476, 2009.
-  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.
-  Tian Qi Chen, Yulia Rubanova, Jesse Bettencourt, and David Duvenaud. Neural ordinary differential equations. arXiv preprint arXiv:1806.07366, 2018.
-  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.
-  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.
-  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.
-  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.
-  R Dennis Cook et al. Fisher lecture: Dimension reduction in regression. Statistical Science, 22(1):1–26, 2007.
-  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.
-  Luc Devroye and Terry Wagner. Distribution-free performance bounds for potential function rules. IEEE Transactions on Information Theory, 25(5):601–604, 1979.
-  David L Donoho. High-dimensional data analysis: The curses and blessings of dimensionality. AMS math challenges lecture, 1(2000):32, 2000.
-  David L Donoho and Jain M Johnstone. Ideal spatial adaptation by wavelet shrinkage. biometrika, 81(3):425–455, 1994.
-  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.
-  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.
-  Weinan E, Chao Ma, and Qingcan Wang. A priori estimates of the population risk for residual networks. arXiv preprint arXiv:1903.02154, 2019.
-  Ronen Eldan and Ohad Shamir. The power of depth for feedforward neural networks. In Conference on Learning Theory, pages 907–940, 2016.
-  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.
-  Vitaly Feldman and Jan Vondrak. High probability generalization bounds for uniformly stable algorithms with nearly optimal rate. arXiv preprint arXiv:1902.10710, 2019.
-  Jerome H Friedman and Werner Stuetzle. Projection pursuit regression. Journal of the American statistical Association, 76(376):817–823, 1981.
-  Haoyu Fu, Yuejie Chi, and Yingbin Liang. Local geometry of one-hidden-layer neural networks for logistic regression. arXiv preprint arXiv:1802.06463, 2018.
-  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.
-  Chao Gao, Jiyi Liu, Yuan Yao, and Weizhi Zhu. Robust estimation and generative adversarial nets. arXiv preprint arXiv:1810.02030, 2018.
-  Surbhi Goel, Adam Klivans, and Raghu Meka. Learning one convolutional layer with overlapping patches. arXiv preprint arXiv:1802.02547, 2018.
-  Noah Golowich, Alexander Rakhlin, and Ohad Shamir. Size-independent sample complexity of neural networks. arXiv preprint arXiv:1712.06541, 2017.
-  Gene H Golub and Charles F Van Loan. Matrix computations. JHU Press, 4 edition, 2013.
-  Ian Goodfellow, Yoshua Bengio, and Aaron Courville. Deep Learning. MIT Press, 2016.
-  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.
-  Suriya Gunasekar, Jason Lee, Daniel Soudry, and Nathan Srebro. Characterizing implicit bias in terms of optimization geometry. arXiv preprint arXiv:1802.08246, 2018.
-  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.
-  Moritz Hardt, Benjamin Recht, and Yoram Singer. Train faster, generalize better: Stability of stochastic gradient descent. arXiv preprint arXiv:1509.01240, 2015.
-  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.
-  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.
-  Geoffrey Hinton, Nitish Srivastava, and Kevin Swersky. Neural networks for machine learning lecture 6a overview of mini-batch gradient descent. 2012.
-  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.
-  Sepp Hochreiter and Jürgen Schmidhuber. Long short-term memory. Neural computation, 9(8):1735–1780, 1997.
-  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.
-  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.
-  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.
-  Sergey Ioffe and Christian Szegedy. Batch normalization: Accelerating deep network training by reducing internal covariate shift. arXiv preprint arXiv:1502.03167, 2015.
-  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.
-  Prateek Jain, Sham M Kakade, Rahul Kidambi, Praneeth Netrapalli, and Aaron Sidford. Accelerating stochastic gradient descent. arXiv preprint arXiv:1704.08227, 2017.
-  Adel Javanmard, Marco Mondelli, and Andrea Montanari. Analysis of a two-layer neural network via displacement convexity. arXiv preprint arXiv:1901.01375, 2019.
-  Ziwei Ji and Matus Telgarsky. Risk and parameter convergence of logistic regression. arXiv preprint arXiv:1803.07300, 2018.
-  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.
-  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.
-  Diederik P Kingma and Jimmy Ba. Adam: A method for stochastic optimization. arXiv preprint arXiv:1412.6980, 2014.
-  Diederik P Kingma and Max Welling. Auto-encoding variational bayes. arXiv preprint arXiv:1312.6114, 2013.
-  Jason M Klusowski and Andrew R Barron. Risk bounds for high-dimensional ridge function combinations including neural networks. arXiv preprint arXiv:1607.01434, 2016.
-  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.
-  Harold Kushner and G George Yin. Stochastic approximation and recursive algorithms and applications, volume 35. Springer Science & Business Media, 2003.
-  Yann LeCun, Yoshua Bengio, and Geoffrey Hinton. Deep learning. nature, 521(7553):436, 2015.
-  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.
-  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.
-  Ker-Chau Li. Sliced inverse regression for dimension reduction. Journal of the American Statistical Association, 86(414):316–327, 1991.
-  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.
-  Yujia Li, Kevin Swersky, and Rich Zemel. Generative moment matching networks. In International Conference on Machine Learning, pages 1718–1727, 2015.
-  Tengyuan Liang. How well can generative adversarial networks (gan) learn densities: A nonparametric view. arXiv preprint arXiv:1712.08244, 2017.
-  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.
-  Min Lin, Qiang Chen, and Shuicheng Yan. Network in network. arXiv preprint arXiv:1312.4400, 2013.
-  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.
-  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.
-  Yuly Makovoz. Random approximants and neural networks. Journal of Approximation Theory, 85(1):98–109, 1996.
-  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.
-  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.
-  Hrushikesh Mhaskar, Qianli Liao, and Tomaso Poggio. Learning functions: when is deep better than shallow. arXiv preprint arXiv:1603.00988, 2016.
-  Hrushikesh N Mhaskar. Neural networks for optimal approximation of smooth and analytic functions. Neural computation, 8(1):164–177, 1996.
-  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.
-  Marco Mondelli and Andrea Montanari. On the connection between learning two-layers neural networks and tensor decomposition. arXiv preprint arXiv:1802.07301, 2018.
-  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.
-  Behnam Neyshabur, Ryota Tomioka, and Nathan Srebro. Norm-based capacity control in neural networks. In Conference on Learning Theory, pages 1376–1401, 2015.
-  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.
-  Ian Parberry. Circuit complexity and neural networks. MIT press, 1994.
-  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.
-  Allan Pinkus. Approximation theory of the mlp model in neural networks. Acta numerica, 8:143–195, 1999.
-  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.
-  Boris T Polyak. Some methods of speeding up the convergence of iteration methods. USSR Computational Mathematics and Mathematical Physics, 4(5):1–17, 1964.
-  Boris T Polyak and Anatoli B Juditsky. Acceleration of stochastic approximation by averaging. SIAM Journal on Control and Optimization, 30(4):838–855, 1992.
-  Boris Teodorovich Polyak and Yakov Zalmanovich Tsypkin. Adaptive estimation algorithms: convergence, optimality, stability. Avtomatika i Telemekhanika, (3):71–84, 1979.
-  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.
-  Sashank J Reddi, Satyen Kale, and Sanjiv Kumar. On the convergence of adam and beyond. 2018.
-  Herbert Robbins and Sutton Monro. A stochastic approximation method. The Annals of Mathematical Statistics, 22(3):400–407, 1951.
-  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.
-  David Rolnick and Max Tegmark. The power of deeper networks for expressing natural functions. arXiv preprint arXiv:1705.05502, 2017.
-  Yaniv Romano, Matteo Sesia, and Emmanuel J Candès. Deep knockoffs. arXiv preprint arXiv:1811.06687, 2018.
-  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.
-  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.
-  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.
-  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.
-  Ruslan Salakhutdinov and Geoffrey Hinton. Deep boltzmann machines. In Artificial intelligence and statistics, pages 448–455, 2009.
-  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.
-  Johannes Schmidt-Hieber. Nonparametric regression using deep neural networks with relu activation function. arXiv preprint arXiv:1708.06633, 2017.
-  Shai Shalev-Shwartz and Shai Ben-David. Understanding machine learning: From theory to algorithms. Cambridge university press, 2014.
-  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.
-  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.
-  Bernard W Silverman. Density estimation for statistics and data analysis. Chapman & Hall, CRC, 1998.
-  Chandan Singh, W James Murdoch, and Bin Yu. Hierarchical interpretations for neural network predictions. arXiv preprint arXiv:1806.05337, 2018.
-  Justin Sirignano and Konstantinos Spiliopoulos. Mean field analysis of neural networks. arXiv preprint arXiv:1805.01053, 2018.
-  Mahdi Soltanolkotabi. Learning relus via gradient descent. In Advances in Neural Information Processing Systems, pages 2007–2017, 2017.
-  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.
-  David A Sprecher. On the structure of continuous functions of several variables. Transactions of the American Mathematical Society, 115:340–355, 1965.
-  Charles J Stone. Optimal global rates of convergence for nonparametric regression. The annals of statistics, pages 1040–1053, 1982.
-  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.
-  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.
-  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.
-  Matus Telgarsky. Benefits of depth in neural networks. arXiv preprint arXiv:1602.04485, 2016.
-  Robert Tibshirani. Regression shrinkage and selection via the lasso. Journal of the Royal Statistical Society: Series B (Methodological), 58(1):267–288, 1996.
-  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.
-  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.
-  Stefan Wager, Sida Wang, and Percy S Liang. Dropout training as adaptive regularization. In Advances in neural information processing systems, pages 351–359, 2013.
-  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.
-  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.
-  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.
-  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.
-  Jason Yosinski, Jeff Clune, Anh Nguyen, Thomas Fuchs, and Hod Lipson. Understanding neural networks through deep visualization. arXiv preprint arXiv:1506.06579, 2015.
-  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.