What's Hidden in a Randomly Weighted Neural Network?

  • 2019-11-29 18:56:53
  • Vivek Ramanujan, Mitchell Wortsman, Aniruddha Kembhavi, Ali Farhadi, Mohammad Rastegari
  • 13

Abstract

Training a neural network is synonymous with learning the values of theweights. In contrast, we demonstrate that randomly weighted neural networkscontain subnetworks which achieve impressive performance without ever trainingthe weight values. Hidden in a randomly weighted Wide ResNet-50 we show thatthere is a subnetwork (with random weights) that is smaller than, but matchesthe performance of a ResNet-34 trained on ImageNet. Not only do these"untrained subnetworks" exist, but we provide an algorithm to effectively findthem. We empirically show that as randomly weighted neural networks with fixedweights grow wider and deeper, an "untrained subnetwork" approaches a networkwith learned weights in accuracy.

 

Quick Read (beta)

What’s Hidden in a Randomly Weighted Neural Network?

            Vivek Ramanujan  22footnotemark: 2 equal contribution    Mitchell Wortsman 11footnotemark: 1  33footnotemark: 3    Aniruddha Kembhavi              Allen Institute for Artificial IntelligenceUniversity of Washington    Ali Farhadi 22footnotemark: 2  33footnotemark: 3    Mohammad Rastegari 22footnotemark: 2  33footnotemark: 3
Abstract

Training a neural network is synonymous with learning the values of the weights. In contrast, we demonstrate that randomly weighted neural networks contain subnetworks which achieve impressive performance without ever training the weight values. Hidden in a randomly weighted Wide ResNet-50 [28] we show that there is a subnetwork (with random weights) that is smaller than, but matches the performance of a ResNet-34 [8] trained on ImageNet [3]. Not only do these “untrained subnetworks” exist, but we provide an algorithm to effectively find them. We empirically show that as randomly weighted neural networks with fixed weights grow wider and deeper, an “untrained subnetwork” approaches a network with learned weights in accuracy.

1 Introduction

What lies hidden in an overparameterized neural network with random weights? If the distribution is properly scaled, then it contains a subnetwork which performs well without ever modifying the values of the weights (as illustrated by Figure 1).

The number of subnetworks is combinatorial in the size of the network, and modern neural networks contain millions or even billions of parameters [21]. We should expect that even a randomly weighted neural network contains a subnetwork that performs well on a given task. In this work, we provide an algorithm to find these subnetworks.

Figure 1: If a neural network with random weights (center) is sufficiently overparameterized, it will contain a subnetwork (right) that perform as well as a trained neural network (left) with the same number of parameters.

Finding subnetworks contrasts with the prevailing paradigm for neural network training – learning the values of the weights by stochastic gradient descent. Traditionally, the network structure is either fixed during training (e.g. ResNet [8] or MobileNet [9]), or optimized in conjunction with the weight values (e.g. Neural Architecture Search (NAS)). We instead optimize to find a good subnetwork within a fixed, randomly weighted network. We do not ever tune the value of any weights in the network, not even the batch norm [10] parameters or first or last layer.

In [4], Frankle and Carbin articulate The Lottery Ticket Hypothesis: neural networks contain sparse subnetworks that can be effectively trained from scratch when reset to their initialization. We offer a complimentary conjecture: within a sufficiently overparameterized neural network with random weights (e.g. at initialization), there exists a subnetwork that achieves competitive accuracy. Specifically, the test accuracy of the subnetwork is able to match the accuracy of a trained network with the same number of parameters.

This work is catalyzed by the recent advances of Zhou et al. [29]. By sampling subnetworks in the forward pass, they first demonstrate that subnetworks of randomly weighted neural networks can achieve impressive accuracy. However, we hypothesize that stochasticity may limit their performance. As the number of parameters in the network grows, they are likely to have a high variability in their sampled networks.

To this end we propose the edge-popup algorithm for finding effective subnetworks within randomly weighted neural networks. We show a signifigant boost in performance and scale to ImageNet. For each fixed random weight in the network, we consider a positive real-valued score. To choose a subnetwork we take the weights with the top-k% highest scores. With a gradient estimator we optimize the scores via SGD. We are therefore able to find a good neural network without ever changing the values of the weights. We empirically demonstrate the efficacy of our algorithm and formally show that under certain technical assumptions the loss decreases on the mini-batch with each modification of the subnetwork.

We experiment on small and large scale datasets for image recognition, namely CIFAR-10 [12] and Imagenet [3]. On CIFAR-10 we empirically demonstrate that as networks grow wider and deeper, untrained subnetworks perform just as well as the dense network with learned weights. On ImageNet, we find a subnetwork of a randomly weighted Wide ResNet50 which is smaller than, but matches the performance of a trained ResNet-34. Moreover, a randomly weighted ResNet-101 [8] with fixed weights contains a subnetwork that is much smaller, but surpasses the performance of VGG-16 [23]. In short, we validate the unreasonable effectiveness of randomly weighted neural networks for image recognition.

Figure 2: In the edge-popup Algorithm, we associate a score with each edge. On the forward pass we choose the top edges by score. On the backward pass we update the scores of all the edges with the straight-through estimator, allowing helpful edges that are “dead” to re-enter the subnetwork. We never update the value of any weight in the network, only the score associated with each weight.

2 Related Work

Lottery Tickets and Supermasks

In [4], Frankle and Carbin offer an intriguing hypothesis: neural networks contain sparse subnetworks that can be effectively trained from scratch when reset to their initialization. These so-called winning tickets have won the “initialization lottery”. Frankle and Carbin find winning tickets by iteratively shrinking the size of the network, masking out weights which have the lowest magnitude at the end of each training run.

Follow up work by Zhou et al. [29] demonstrates that winning tickets achieve better than random performance without training. Motivated by this result they propose an algorithm to identify a “supermask” – a subnetwork of a randomly initialized neural network that achieves high accuracy without training. On CIFAR-10, they are able to find subnetworks of randomly initialized neural networks that achieve 65.4% accuracy.

The algorithm presented by Zhou et al. is as follows: for each weight w in the network they learn an associated probability p. On the forward pass they include weight w with probability p and otherwise zero it out. Equivalently, they use weight w~=wX where X is a Bernoulli(p) random variable (X is 1 with probability p and 0 otherwise). The probabilities p are the output of a sigmoid, and are learned using stochastic gradient descent. The terminology supermask” arises as finding a subnetwork is equivalent to learning a binary mask for the weights.

Our work builds upon Zhou et al., though we recognize that the stochasticity of their algorithm may limit performance. In section 3.1 we provide more intuition for this claim. We show a significant boost in performance with an algorithm that does not sample supermasks on the forward pass. For the first time we are able to match the performance of a dense network with a supermask.

Neural Architecture Search (NAS)

The advent of modern neural networks has shifted the focus from feature engineering to feature learning. However, researchers may now find themselves manually engineering the architecture of the network. Methods of Neural Architecture Search (NAS) [30, 2, 16, 24] instead provide a mechanism for learning the architecture of neural network jointly with the weights. Models powered by NAS have recently obtained state of the art classification performance on ImageNet [25].

As highlighted by Xie et al. [27], the connectivity patterns in methods of NAS remain largely constrained. Surprisingly, Xie et al. establish that randomly wired neural networks can achieve competitive performance. Accordingly, Wortsman et al. [26] propose a method of Discovering Neural Wirings (DNW) – where the weights and structure are jointly optimized free from the typical constraints of NAS. We highlight DNW as we use a similar method of analysis and gradient estimator to optimize our supermasks. In DNW, however, the subnetwork is chosen by taking the weights with the highest magnitude. There is therefore no way to learn supermasks with DNW as the weights and connectivity are inextricably linked – there is no way to separate the weights and the structure.

Weight Agnostic Neural Networks

In Weight Agnostic Neural Networks (WANNs) [5], Gaier and Ha question if an architecture alone may encode the solution to a problem. They present a mechanism for building neural networks that achieve high performance when each weight in the network has the same shared value. Importantly, the performance of the network is agnostic to the value itself. They are able to obtain 92% accuracy on MNIST [15].

We are quite inspired by WANNs, though we would like to highlight some important distinctions. Instead of each weight having the same value, we explore the setting where each weight has a random value. In Section A.2.2 of their appendix, Gaier and Ha mention that they were not successful in this setting. However, we find a good subnetwork for a given random initialization – the supermasks we find are not agnostic to the weights. Finally, Gaier and Ha construct their network architectures, while we look for supermasks within standard architectures.

Linear Classifiers

Linear classifiers on top of randomly weighted neural networks are often used as baselines in unsupervised learning [18]. Our work is different in motivation, we explicitly find untrained subnetworks which achieve high performance without changing any weight values, including the final layer.

3 Method

In this section we present our optimization method for finding effective subnetworks within randomly weighted neural networks. We begin by building intuition in an unusual setting – the infinite width limit. Next we motivate and present our algorithm for finding effective subnetworks.

3.1 Intuition

The Existence of Good Subnetworks

Modern neural networks have a staggering number of possible subnetworks. Consequently, even at initialization, a neural network should contain a subnetwork which performs well.

To build intuition we will consider an extreme case – a neural network N in the infinite width limit (for a convolutional neural networks, the width of the network is the number of channels). As in Figure 1, let τ be a network with the same structure of N that achieves achieves good accuracy. If the weights of N are initialized using any standard scaling of a normal distribution, e.g. xavier [6] or kaiming [7], then we may show there exists a subnetwork of N that achieves the same performance as τ without training. Let q be the probability that a given subnetwork of N has weights that are close enough to τ to obtain the same accuracy. This probability q is extremely small, but it is still nonzero. Therefore, the probability that no subnetwork of N is close enough to τ is effectively (1-q)S where S is the number of subnetworks. S grows very quickly with the width of the network, and this probability becomes arbitrarily small.

How Should We Find A Good Subnetwork

Even if there are good subnetworks in randomly weighted neural networks, how should we find them?

Zhou et al. learn an associated probability p with each weight w in the network. On the forward pass they include weight w with probability p (where p is the output of a sigmoid) and otherwise zero it out. The infinite width limit provides intuition for a possible shortcoming of the algorithm presented by Zhou et al. [29]. Even if the parameters p are fixed, the algorithm will likely never observe the same subnetwork twice. As such, the gradient estimate becomes more unstable, and this in turn may make training difficult.

Our algorithm for finding a good subnetwork is illustrated by Figure 2. With each weight w in the neural network we learn a positive, real valued popup score s. The subnetwork is then chosen by selecting the weights in each layer corresponding to the top-k% highest scores. For simplicity we use the same value of k for all layers.

How should we update the score suv? Consider a single edge in a fully connected layer which connects neuron u to neuron v. Let wuv be the weight of this edge, and suv the associated score. If this score is initially low then wuv is not selected in the forward pass. But we would still like a way to update its score to allow it to pop back up. Informally, with backprop [22] we compute how the loss “wants” node v’s input to change (i.e. the negative gradient). We then examine the weighted output of node u. If this weighted output is aligned with the negative gradient, then node u can take node v’s output where the loss “wants” it to go. Accordingly, we should increase the score. If this alignment happens consistently, then the score will continue to increase and the edge will re-enter the chosen subnetwork (i.e. popup).

More formally, if wuv𝒵u denotes the weighted output of neuron u, and v denotes the input of neuron v, then we update suv as

suvsuv-αv𝒵uwuv. (1)

This argument and the analysis that follows is motivated and guided by the work of [26]. In their work, however, they do not consider a score and are instead directly updating the weights. In the forward pass they use the top k% of edges by magnitude, and therefore there is no way of learning a subnetwork without learning the weights. Their goal is to train sparse neural networks, while we aim to showcase the efficacy of randomly weighted neural networks.

3.2 The edge-popup Algorithm and Analysis

We now formally detail the edge-popup algorithm.

For clarity, we first describe our algorithm for a fully connected neural network. In Section B.2 we provide the straightforward extension to convolutions along with code in PyTorch [20].

A fully connected neural network consists of layers 1,,L where layer has n nodes 𝒱()={v1(),,vn()}. We let v denote the input to node v and let 𝒵v denote the output, where 𝒵v=σ(v) for some non-linear activation function σ (e.g. ReLU [13]). The input to neuron v in layer is a weighted sum of all neurons in the preceding layer. Accordingly, we write v as

v=u𝒱(-1)wuv𝒵u (2)

where wuv are the network parameters for layer . The output of the network is taken from the final layer while the input data is given to the very first layer. Before training, the weights wuv for layer are initialized by independently sampling from distribution 𝒟. For example, if we are using kaiming normal initialization [7] with ReLU activations, then 𝒟=𝒩(0,2/n-1) where 𝒩 denotes the normal distribution.

Normally, the weights wuv are optimized via stochastic gradient descent. In our edge-popup algorithm, we instead keep the weights at their random initialization, and optimize to find a subnetwork 𝒢=(𝒱,). We then compute the input of node v in layer as

v=(u,v)wuv𝒵u (3)

where 𝒢 is a sub-graph of the original fully connected network11 1 The original network has edges fc==1L-1(𝒱×𝒱+1) where × denotes the cross-product..

As mentioned above, for each weight wuv in the original network we learn a popup score suv. We choose the subnetwork 𝒢 by selecting the weights in each layer which have the top-k% highest scores. Equation 3 may therefore be written equivalently as

v=u𝒱(-1)wuv𝒵uh(suv) (4)

where h(suv)=1 if suv is among the top k% highest scores in layer and h(suv)=0 otherwise. Since the gradient of h is 0 everywhere it is not possible to directly compute the gradient of the loss with respect to suv. We instead use the straight-through gradient estimator [1], in which h is treated as the identity in the backwards pass – the gradient goes “straight-through” h.

Consequently, we approximate the gradient to suv as

g^suv=\pdvv\pdvvsuv=\pdvvwuv𝒵u (5)

where is the loss we are trying to minimize. The scores suv are then updated via stochastic gradient descent with learning rate α. If we ignore momentum and weight decay [14] then we update suv as

s~uv=suv-α\pdvvwuv𝒵u (6)

where s~uv denotes the score after the gradient step22 2 To ensure that the scores are positive we take the absolute value..

As the scores change certain edges in the subnetwork will be replaced with others. Motivated by the analysis of [26] we show that when swapping does occur, the loss decreases for the mini-batch.

Theorem 1: When edge (i,ρ) replaces (j,ρ) and the rest of the subnetwork remains fixed then the loss decreases for the mini-batch (provided the learning rate α is sufficiently small, and the loss is smooth).

Proof. Let s~uv denote the score of weight suv after the gradient update. If edge (i,ρ) replaces (j,ρ) then our algorithm dictates that siρ<sjρ but s~iρ>s~jρ. Accordingly,

s~iρ-siρ>s~jρ-sjρ (7)

which implies that

-α\pdvρwiρ𝒵i>-α\pdvρwjρ𝒵j (8)

by the update rule given in Equation 6. Let ~ρ denote the input to node k after the swap is made and ρ denote the original input. Note that ~ρ-ρ=wiρZi-wjρZj by Equation 3. We now wish to show that (~ρ)<(ρ).

When the learning rate is sufficiently small (and the loss is smooth) we may assume that ~ρ is close to ρ and ignore second-order terms in a Taylor expansion:

(~ρ) =(ρ+(~ρ-ρ)) (9)
(ρ)+ρ(~ρ-ρ) (10)
=(ρ)+ρ(wiρ𝒵i-wjρ𝒵j) (11)

and from equation 8 we have that ρ(wiρ𝒵i-wjρ𝒵j)<0 and so (~ρ)<(ρ) as needed.

We examine a more general case of Theorem 1 in Section B.1 of the supplementary material.

Figure 3: Going Deeper: Experimenting with shallow to deep neural networks on CIFAR-10 [12]. As the network becomes deeper, we are able to find subnetworks at initialization that perform as well as the dense original network when trained. The baselines are drawn as a horizontal line as we are not varying the % of weights. When we write Weights 𝒟 we mean that the weights are randomly drawn from distribution 𝒟 and are never tuned. Instead we find subnetworks with size (% of Weights)/100 * (Total # of Weights).

4 Experiments

We demonstrate the unreasonable effectiveness of randomly weighted neural networks image recognition on standard benchmark datasets CIFAR-10 [12] and ImageNet [3]. This section is organized as follows: in Section 4.1 we discuss the experimental setup and hyperparameters. We perform a series of ablations at small scale: we examine the effect of k, the % of Weights which remain in the subnetwork, and the effect of width. In Section 4.4 we compare against the algorithm of Zhou et al., followed by Section 4.5 in which we study the effect of the distribution used to sample the weights. We conclude with Section 4.6, where we otpimize to find subnetworks of randomly weighted neural networks which achieve good performance on ImageNet [3].

4.1 Experimental Setup

We use two different distributions for the weights in our network:

  • Kaiming Normal [7], which we denote 𝒩k. Following the notation in section 3.2 the Kaiming Normal distribution is defined as 𝒩k=𝒩(0,2/n-1) where 𝒩 denotes the normal distribution.

  • Signed Kaiming Constant which we denote Uk. Here we set each weight to be a constant and randomly choose its sign to be + or -. The constant we choose is the standard deviation of Kaiming Normal, and as a result the variance is the same. We use the notation Uk as we are sampling uniformly from the set {-σk,σk} where σk is the standard deviation for Kaiming Normal (i.e. 2/n-1).

In Section 4.5 we reflect on the importance of the random distribution and experiment with alternatives.

On CIFAR-10 [12] we experiment with simple VGG-like architectures of varying depth. These architectures are also used by Frankle and Carbin [4] and Zhou et al. [29] and are provided in Table 1. On ImageNet we experiment with ResNet-50 and ResNet-101 [8], as well as their wide variants [28].

Model Conv2 Conv4 Conv6 Conv8
Conv Layers 64, 64, pool 64, 64, pool 128, 128, pool 64, 64, pool 128, 128, pool 256, 256, pool 64, 64, pool 128, 128, pool 256, 256, pool 512, 512, pool
FC 256, 256, 10 256, 256, 10 256, 256, 10 256, 256, 10
Table 1: For completeness we provide the architecture of the simple VGG-like [23] architectures used for CIFAR-10 [12], which are identical to those used by Frankle and Carbin [4] and Zhou et al. [29]. However, the slightly deeper Conv8 does not appear in the previous work. Each model first performs convolutions followed by the fully connected (FC) layers, and pool denotes max-pooling.

In every experiment we train for 100 epochs and report the last epoch accuracy on the validation set. When we optimize with Adam [11] we do not decay the learning rate. When we optimize with SGD we use cosine learning rate decay [17]. On CIFAR-10 [12] we train our models with weight decay 1e-4, momentum 0.9, batch size 128, and learning rate 0.1. We also often run both an Adam and SGD baseline where the weights are learned. The Adam baseline uses the same learning rate and batch size as in [4, 29]33 3 Batch size 60, learning rate 2e-4, 3e-4 and 3e-4 for Conv2, Conv4, and Conv6 respectively Conv8 is not tested in [4], though we use find that learning rate 3e-4 still performs well.. For the SGD baseline we find that training does not converge with learning rate 0.1, and so we use 0.01. As standard we also use weight decay 1e-4, momentum 0.9, and batch size 128.

Figure 4: Going Wider: Varying the width (i.e. number of channels) of Conv4 and Conv6 for CIFAR-10 [12]. When Conv6 is wide enough, a subnetwork of the randomly weighted model (with %𝑊𝑒𝑖𝑔ℎ𝑡𝑠=50) performs just as well as the full model when it is trained.

For the ImageNet experiments we use the hyperparameters found on NVIDIA’s public github example repository for training ResNet [19]. For simplicity, our edge-popup algorithm does not modify batch norm parameters, they are frozen at their default initialization in PyTorch (i.e. bias 0, scale 1).

This discussion has encompassed the extent of the hyperparameter tuning for our models. We do, however, perform hyperparameter tuning for the Zhou et al. [29] baseline and improve accuracy significantly. We include further discussion of this in Section 4.4.

In all experiments on CIFAR-10 [12] we use 5 different random seeds and plot the mean accuracy ± one standard deviation. Moreover, on all figures, Learned Dense Weights denotes the standard training the full model (all weights remaining).

4.2 Varying the % of Weights

Our algorithm has one associated parameter: the % of weights which remain in the subnetwork, which we refer to as k. Figure 3 illustrates how the accuracy of the subnetwork we find varies with k, a trend which we will now dissect. We consider k[10,30,50,70,90] and plot the dense model when it is trained as a horizontal line (as it has 100% of the weights).

We recieve the worst accuracy when k approaches 0 or 100. When k approaches 0, we are not able to perform well as our subnetwork has very few weights. On the other hand, when k approaches 100, our network outputs are random.

The best accuracy occurs when k[30,70], and we make a combinatorial argument for this trend. We are choosing kn weights out of n, and there are (nkn) ways of doing so. The number of possible subnetworks is therefore maximized when k0.5, and at this value our search space is at its largest.

4.3 Varying the Width

Our intuition from Section 3.1 suggests that as the network gets wider, a subnetwork of a randomly weighted model should approach the trained model in accuracy. How wide is wide enough?

Figure 5: Varying the width of Conv4 on CIFAR-10 [12] while modifying k so that the # of Parameters is fixed along each curve. c1,c2,c3 are constants which coincide with # of Parameters for k=[30,50,70] for width multiplier 1.
Figure 6: Comparing the performance of edge-popup with the algorithm presented by Zhou et al. [29] on CIFAR-10 [12].

In Figure 4 we vary the width of Conv4 and Conv6. The width of a linear layer is the number of “neurons”, and the width of a convolution layer is the number of channels. The width multiplier is the factor by which the width of all layers is scaled. A width multiplier of 1 corresponds to the models tested in Figure 3.

As the width multiplier increases, the gap shrinks between the accuracy a subnetwork found with edge-popup and the dense model when it is trained. Notably, when Conv6 is wide enough, a subnetwork of the randomly weighted model (with %𝑊𝑒𝑖𝑔ℎ𝑡𝑠=50) performs just as well as the dense model when it is trained.

Moreover, this boost in performance is not solely from the subnetwork having more parameters. Even when the # of parameters is fixed, increasing the width and therefore the search space leads to better performance. In Figure 5 we fix the number of parameters and while modifying k and the width multiplier. Specifically, we test k[30,50,70] for subnetworks of constant size c1,c2 and c3. On Figure 5 we use || denote the size of the subnetwork.

4.4 Comparing with Zhou et al. [29]

In Figure 6 we compare the performance of edge-popup with Zhou et al. Their work considers distributions 𝒩x and Ux, which are identical to those presented in Section 4.1 but with xavier normal [6] instead of kaiming normal [7] – the factor of 2 is omitted from the standard deviation. By running their algorithm with 𝒩k and Uk we witness a significant improvement. However, even the Nx and Ux results exceed those in the paper as we perform some hyperparameter tuning. As in our experiments on CIFAR-10, we use SGD with weight decay 1e-4, momentum 0.9, batch size 128, and a cosine scheduler [17]. We double the learning rate until we see their performance become worse, and settle on 20044 4 An absurdly high learning rate is required as mentioned in their work..

Figure 7: Testing different weight distributions on CIFAR-10 [12].

4.5 Effect of The Distribution

The distribution that the random weights are sampled from is very important. As illustrated by Figure 7, the performance of our algorithm vastly decreases when we switch to using xavier normal [6] or kaiming uniform [7].

Following the derivation in [7], the variance of the forward pass is not exactly 1 when we consider a subnetwork with only k% of the weights. To reconcile for this we could scale standard deviation by 1/k. This distribution is referred to as “Scaled Kaiming Normal” on Figure 7. We may also consider this scaling for the Signed Kaiming Constant distribution which is described in Section 4.1.

Figure 8: Testing our Algorithm on ImageNet [3]. We use a fixed k=30%, and find subnetworks within a randomly weighted ResNet-50 [8], Wide ResNet-50 [28], and ResNet-101. Notably, a randomly weighted Wide ResNet-50 contains a subnetwork which is smaller than, but matches the performance of ResNet-34. Note that for the non-dense models, # of Parameters denotes the size of the subnetwork.

4.6 ImageNet [3] Experiments

On ImageNet we observe similar trends to CIFAR-10. As ImageNet is much harder, computationally feasible models are not overparameterized to the same degree. As a consequence, the performance a randomly weighted subnetwork does not match the full model with learned weights. However, we still witness a very encouraging trend – the performance increases with the width and depth of the network.

As illustrated by Figure 8, a randomly weighted Wide ResNet-50 contains a subnetwork that is smaller than, but matches the accuracy of ResNet-34 when trained on ImageNet [3]. As strongly suggested by our trends, better and larger “parent” networks would result in even stronger performance on ImageNet [3]. A table which reports the numbers in Figure 8 may be found in Section A of the supplementary material.

Figure 9 illustrates the effect of k, which follows an almost identical trend: k[30,70] performs best though 30 now provides the best performance. Figure 9 also demonstrates that we significantly outperform Zhou et al. at scale (they do not test on ImageNet in their paper). Their algorithm does not allow an explicit choice of the % of weights remaining in the subnetwork, and we found the algorithm unstable outside of the range reported.

The choice of the random distribution matters more for ImageNet. The “Scaled” distribution we discuss in Section 4.5 did not show any discernable difference on CIFAR-10. However, Figure 10 illustrates that on ImageNet it is much better. Recall that the “Scaled” distribution adds a factor of 1/k, which has less of an effect when k approaches 100%=1. This result highlights the possibility of finding better distributions which works best for this algorithm.

Figure 9: Examining the effect of % weights on ImageNet for edge-popup and the method of Zhou et al.
Figure 10: Examining the effect of using the “Scaled” initialization detailed in Section 4.5 on ImageNet.

5 Conclusion

Hidden within randomly weighted neural networks we find subnetworks with compelling accuracy. This work provides an avenue for many areas of exploration. For example, we anticipate the development of faster algorithms, or the alternating optimization of the structure and the weights.

Finally, we hope that our findings serve as a useful step in the pursuit the understanding of the optimization of neural networks.

Acknowledgments

We thank Jesse Dodge and Nicholas Lourie for many helpful discussions, Gabriel Ilharco for valuable feedback, and Sarah Pratt for (n2).

References

  • [1] Yoshua Bengio, Nicholas Léonard, and Aaron Courville. Estimating or propagating gradients through stochastic neurons for conditional computation, 2013.
  • [2] Han Cai, Ligeng Zhu, and Song Han. ProxylessNAS: Direct neural architecture search on target task and hardware. In ICLR, 2019.
  • [3] Jia Deng, Wei Dong, Richard Socher, Li-Jia Li, Kai Li, and Li Fei-Fei. Imagenet: A large-scale hierarchical image database. In CVPR, 2009.
  • [4] Jonathan Frankle and Michael Carbin. The lottery ticket hypothesis: Finding sparse, trainable neural networks. In ICLR 2019, 2019.
  • [5] Adam Gaier and David Ha. Weight agnostic neural networks. 2019. https://weightagnostic.github.io.
  • [6] Xavier Glorot and Yoshua Bengio. Understanding the difficulty of training deep feedforward neural networks. In Yee Whye Teh and Mike Titterington, editors, Proceedings of the Thirteenth International Conference on Artificial Intelligence and Statistics, volume 9 of Proceedings of Machine Learning Research, pages 249–256, Chia Laguna Resort, Sardinia, Italy, 13–15 May 2010. PMLR.
  • [7] Kaiming He, Xiangyu Zhang, Shaoqing Ren, and Jian Sun. Delving deep into rectifiers: Surpassing human-level performance on imagenet classification. In Proceedings of the 2015 IEEE International Conference on Computer Vision (ICCV), ICCV ’15, pages 1026–1034, Washington, DC, USA, 2015. IEEE Computer Society.
  • [8] Kaiming He, Xiangyu Zhang, Shaoqing Ren, and Jian Sun. Deep residual learning for image recognition. In CVPR, 2016.
  • [9] Andrew G. Howard, Menglong Zhu, Bo Chen, Dmitry Kalenichenko, Weijun Wang, Tobias Weyand, Marco Andreetto, and Hartwig Adam. Mobilenets: Efficient convolutional neural networks for mobile vision applications. ArXiv, abs/1704.04861, 2017.
  • [10] Sergey Ioffe and Christian Szegedy. Batch normalization: Accelerating deep network training by reducing internal covariate shift. In ICML, 2015.
  • [11] Diederik P. Kingma and Jimmy Ba. Adam: A method for stochastic optimization, 2014.
  • [12] Alex Krizhevsky. Learning multiple layers of features from tiny images. Technical report, University of Toronto, 2009.
  • [13] Alex Krizhevsky, Ilya Sutskever, and Geoffrey E Hinton. Imagenet classification with deep convolutional neural networks. In F. Pereira, C. J. C. Burges, L. Bottou, and K. Q. Weinberger, editors, Advances in Neural Information Processing Systems 25, pages 1097–1105. Curran Associates, Inc., 2012.
  • [14] Anders Krogh and John A. Hertz. A simple weight decay can improve generalization. In NIPS, 1991.
  • [15] Yann LeCun and Corinna Cortes. MNIST handwritten digit database. 2010.
  • [16] Chenxi Liu, Barret Zoph, Maxim Neumann, Jonathon Shlens, Wei Hua, Li-Jia Li, Li Fei-Fei, Alan Yuille, Jonathan Huang, and Kevin Murphy. Progressive neural architecture search. In Proceedings of the European Conference on Computer Vision (ECCV), pages 19–34, 2018.
  • [17] Ilya Loshchilov and Frank Hutter. Sgdr: Stochastic gradient descent with warm restarts, 2016.
  • [18] Mehdi Noroozi and Paolo Favaro. Unsupervised learning of visual representations by solving jigsaw puzzles, 2016.
  • [19] NVIDIA. Deeplearningexamples/pytorch/classification. https://github.com/NVIDIA/DeepLearningExamples/tree/master/PyTorch/Classification/RN50v1.5, 2019.
  • [20] 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.
  • [21] Colin Raffel, Noam Shazeer, Adam Roberts, Katherine Lee, Sharan Narang, Michael Matena, Yanqi Zhou, Wei Li, and Peter J. Liu. Exploring the limits of transfer learning with a unified text-to-text transformer. arXiv e-prints, 2019.
  • [22] David E. Rumelhart, Geoffrey E. Hinton, and Ronald J. Williams. Learning representations by back-propagating errors. Nature, 323:533–536, 1986.
  • [23] Karen Simonyan and Andrew Zisserman. Very deep convolutional networks for large-scale image recognition. 2015.
  • [24] Mingxing Tan, Bo Chen, Ruoming Pang, Vijay Vasudevan, and Quoc V. Le. Mnasnet: Platform-aware neural architecture search for mobile. CoRR, abs/1807.11626, 2018.
  • [25] Mingxing Tan and Quoc V. Le. Efficientnet: Rethinking model scaling for convolutional neural networks. In ICML, 2019.
  • [26] Mitchell Wortsman, Ali Farhadi, and Mohammad Rastegari. Discovering neural wirings. ArXiv, abs/1906.00586, 2019.
  • [27] Saining Xie, Alexander Kirillov, Ross Girshick, and Kaiming He. Exploring randomly wired neural networks for image recognition. arXiv preprint arXiv:1904.01569, 2019.
  • [28] Sergey Zagoruyko and Nikos Komodakis. Wide residual networks, 2016.
  • [29] Hattie Zhou, Janice Lan, Rosanne Liu, and Jason Yosinski. Deconstructing lottery tickets: Zeros, signs, and the supermask, 2019.
  • [30] Barret Zoph and Quoc V Le. Neural architecture search with reinforcement learning. arXiv preprint arXiv:1611.01578, 2016.
Method Model Initialization % of Weights # of Parameters Accuracy
ResNet-34 [8] - - 21.8M 73.3%
Learned Dense Weights (SGD) ResNet-50 [8] - - 25M 76.1%
Wide ResNet-50 [28] - - 69M 78.1%
ResNet-50 Kaiming Normal 30% 7.6M 61.71%
edge-popup ResNet-101 Kaiming Normal 30% 13M 66.15%
Wide ResNet-50 Kaiming Normal 30% 20.6M 67.95%
ResNet-50 Signed Constant 30% 7.6M 68.6%
edge-popup ResNet-101 Signed Constant 30% 13M 72.3%
Wide ResNet-50 Signed Constant 30% 20.6M 73.3%
Table 2: ImageNet [3] classification results corresponding to Figure 8. Note that for the non-dense models, # of Parameters denotes the size of the subnetwork.

Appendix A Table of ImageNet Results

In Table 2 we provide a table of the results for image classification with ImageNet [3]. These results correspond exactly to Figure 8.

Appendix B Additional Technical Details

In this section we first prove a more general case of Theorem 1 then provide an extension of edge-popup for convolutions along with code in PyTorch [20], found in Algorithm B.2.

B.1 A More General Case of Theorem 1

We now examine a more general case of Theorem 1, where the two swapped edges are not connected to the same node. Again we are motivated by the analysis of [26], though we tackle a more general case.

Theorem 1 (more general): When a nonzero number of edges are swapped in one layer and the rest of the network remains fixed then the loss decreases for the mini-batch (provided the learning rate α is sufficiently small, and the loss is smooth).

Proof. As before, we let s~uv denote the score of weight suv after the gradient update. Additionally, let ~v denote the input to node v after the gradient update whereas v is the input to node v before the update. Finally, let i1,,in denote the n nodes in layer -1 and j1,,jm denote the m notes in layer . Our goal is to show that

(~j1,,~jm)<(j1,,jm) (12)

where the loss is written as a function of layer ’s input for brevity. Since α is small and the loss is smooth we may assume that each ~jk is close to jk and ignore second-order terms in a Taylor expansion:

(~j1,,~jm) (13)
=(j1+(~j1-j1),,jm+(~jm-jm)) (14)
=(j1,,jm)+k=1mjk(~jk-jk) (15)

And so, in order to show Equation 12 it suffices to show that

k=1mjk(~jk-jk)<0. (16)

It is helpful to rewrite the sum to be over edges. Specifically, we will consider the sets old and new where new contains all edges that entered the network after the gradient update and old consists of edges which were previously in the subnetwork, but have now exited. As the total number of edges is conserved we know that |new|=|old|, and by assumption |new|>0.

Using the definition of k and ~k from Equation 3 we may rewrite Equation 16 as

(ia,jb)newjbwiajbZia-(ic,jd)oldjdwicjdZic<0 (17)

which, by Equation 6 and factoring out 1/α becomes

(ia,jb)new(siajb-s~iajb)-(ic,jd)old(sicjd-s~icjd)<0. (18)

We now show that

(siajb-s~iajb)-(sicjd-s~icjd)<0 (19)

for any pair of edges (ia,jb)new and (ic,jd)old. Since |new|=|old|>0 we are then able to conclude that Equation 18 holds.

As (ia,jb) was not in the edge set before the gradient update, but (ic,jd) was, we can conclude

siajb-sicjd<0. (20)

Likewise, since (ia,jb) is in the edge set after the gradient update, but (ic,jd) isn’t, we can conclude

s~icjd-s~iajb<0. (21)

By adding Equation 21 and Equation 20 we find that Equation 19 is satisfied as needed.

B.2 Extension to Convolutional Neural Networks

In order to show that our method extends to convolutional layers we recall that convolutions may be written in a form that resembles Equation 2. Let κ be the kernel size which we assume is odd for simplicity, then for w{1,,W} and h{1,,H} we have

vw,h=u𝒱(-1)κ1=1κκ2=1κwuv(κ1,κ2)𝒵u(w+κ1-κ2,h+κ2-κ2) (22)

where instead of “neurons”, we now have “channels”. The input v and output 𝒵v are now two dimensional and so 𝒵v(w,h) is a scalar. As before, 𝒵v=σ(v) where σ is a nonlinear function. However, in the convolutional case σ is often batch norm [10] followed by ReLU (and then implicitly followed by zero padding).

Instead of simply having weights wuv we now have weights wuv(κ1,κ2) for κ1{1,,κ}, κ2{1,,κ}. Likewise, in our edge-popup Algorithm we now consider scores suv(κ1,κ2) and again use the top k% in the forwards pass. As before, let h(suv(κ1,κ2))=1 if suv(κ1,κ2) is among the top k% highest scores in the layer and h(suv(κ1,κ2))=0 otherwise. Then in edge-popup we are performing a convolution as

vw,h=u𝒱(-1)κ1=1κκ2=1κwuv(κ1,κ2)𝒵u(w+κ1-κ2,h+κ2-κ2)h(suv(κ1,κ2)) (23)

which mirrors the formulation of edge-popup in Equation 4. In fact, when κ=W=H=1 (i.e. a 1x1 convolution on a 1x1 feature map) then Equation 23 and Equation 4 are equivalent.

The update for the scores is quite similar, though we must now sum over all spatial (i.e. w and h) locations as given below:

suv(κ1,κ2)suv(κ1,κ2)-αw=1Wh=1H\pdvvw,hwuv(κ1,κ2)𝒵u(w+κ1-κ2,h+κ2-κ2) (24)

In summary, we now have κ2 edges between each u and v. The PyTorch [20] code is given by Algorithm B.2, where h is GetSubnet. The gradient goes straight through h in the backward pass, and PyTorch handles the implementation of these equations.

{algorithm}

[t] PyTorch code for an edge-popup Conv.

class GetSubnet(autograd.Function):
    @staticmethod
    def forward(ctx, scores, k):
         # Get the subnetwork by sorting the scores and using the top k%
        out = scores.clone()
        _, idx = scores.flatten().sort()
        j = int((1-k) * scores.numel())
        # flat_out and out access the same memory.
        flat_out = out.flatten()
        flat_out[idx[:j]] = 0
        flat_out[idx[j:]] = 1
        return out
    @staticmethod
    def backward(ctx, g):
        # send the gradient g straight-through on the backward pass.
        return g, None
class SubnetConv(nn.Conv2d):
    # self.k is the % of weights remaining, a real number in [0,1]
    # self.popup_scores is a Parameter which has the same shape as self.weight
    # Gradients to self.weight, self.bias have been turned off.
    def forward(self, x):
        # Get the subnetwork by sorting the scores.
        adj = GetSubnet.apply(
                    self.popup_scores.abs(), self.k)
        # Use only the subnetwork in the forward pass.
        w = self.weight * adj
        x = F.conv2d(
            x, w, self.bias, self.stride, self.padding, self.dilation, self.groups
        )
        return x