Source-linked AI summary

Random Walk Initialization for Training Very Deep Feedforward Networks

David Sussillo, L. F. Abbott

arXiv:1412.6558v3cs.NEcs.LGstat.ML

TL;DR

The paper addresses vanishing and exploding gradients in very deep networks, where back-propagated error norms can change exponentially. It analyzes successive random-matrix products and derives scaling rules for linear and ReLU networks. Correct scaling produces a log-norm random walk and supports training feedforward networks at depths above 200 layers, while wider layers reduce fluctuations.

  • Problem

    Very deep networks suffer from back-propagated error signals that can decrease or increase exponentially with depth, hindering the usefulness of additional layers.

    Source: §1 · p. 1
  • Method

    The paper combines random-matrix analysis with computational experiments on deep feedforward networks and derives scaling rules for linear and ReLU cases.

    Source: §1 · p. 2, §2.2 · p. 4, §3.2 · p. 6
  • Results

    Correctly initialized feedforward networks can be trained at depths upwards of 200 layers, with analytic scaling choices matching MNIST training-error results.

    Source: §3.3 · p. 7, §4 · p. 9
  • Takeaways & Limitations

    The log-norm random walk has fluctuations that can be reduced by increasing layer width, providing a practical way to mitigate vanishing gradients.

    Source: §1 · p. 2, §4 · p. 9

    For ReLU networks with width N below about 20, a layer may have no active units, so those cases were excluded from numerical studies.

    Source: §2.2 · p. 4

Abstract

from arXiv · show

Training very deep networks is an important open problem in machine learning. One of many difficulties is that the norm of the back-propagated error gradient can grow or decay exponentially. Here we show that training very deep feed-forward networks (FFNs) is not as difficult as previously thought. Unlike when back-propagation is applied to a recurrent network, application to an FFN amounts to multiplying the error gradient by a different random matrix at each layer. We show that the successive application of correctly scaled random matrices to an initial vector results in a random walk of the log of the norm of the resulting vectors, and we compute the scaling that makes this walk unbiased. The variance of the random walk grows only linearly with network depth and is inversely proportional to the size of each layer. Practically, this implies a gradient whose log-norm scales with the square root of the network depth and shows that the vanishing gradient problem can be mitigated by increasing the width of the layers. Mathematical analyses and experimental results using stochastic gradient descent to optimize tasks related to the MNIST and TIMIT datasets are provided to support these claims. Equations for the optimal matrix scaling are provided for the linear and ReLU cases.

1 INTRODUCTION

Vanishing gradients can grow or decay exponentially with depth, but feedforward networks apply different random matrices across layers, enabling a less drastic gradient behavior. The paper analyzes this effect mathematically and empirically and introduces Random Walk Initialization based on log-norm dynamics.

  • Feedforward back-propagated error signals typically decrease or increase exponentially with distance from the final layer, limiting the usefulness of adding depth.Source: §1 · p. 1
  • Recurrent-network analyses repeatedly apply similar matrices, so eigenvalue magnitudes above or below one produce exponential growth or decay.Source: §1 · p. 1
  • Appropriate initialization can constrain gradient behavior initially, but optimization almost surely destroys this fine-tuning.Source: §1 · p. 2
  • Feedforward networks differ because back-propagation applies a different random matrix at every layer, yielding a gradient magnitude that scales only with the square root of depth when correctly initialized.Source: §1 · p. 2
  • The paper studies successive random-matrix applications mathematically and experimentally, introducing Random Walk Initialization from the random walk of back-propagated log-norms.Source: §1 · p. 2

2 ANALYSIS AND PROPOSED INITIALIZATION

The analysis models back-propagated gradient magnitudes in randomly initialized FFNs as log-norm random walks and derives scaling choices that keep them near unity. Analytical results and simulations show how width, nonlinearity, and initialization scaling affect gradient stability.

2.1 THE MAGNITUDE OF THE ERROR GRADIENT IN FFNS

  • The network assumes width N, Gaussian weights with variance 1/N, zero biases, and a scale factor g normalized by f′(0) = 1.These assumptions define the initialization regime analyzed for linear, ReLU, and tanh nonlinearities.Source: §2.1 · p. 2
  • Different random matrices at successive FFN layers make the gradient-norm analysis a random-walk problem rather than repeated-matrix dynamics.The across-layer ratio Z is kept near order 1 by adjusting the matrix scale g before learning changes the weights.Source: §2.1 · p. 3, §2.2 · p. 3
  • Taking ln(Z) converts the product of layerwise gradient factors into a sum whose steps are ln(z_d), and Random Walk Initialization chooses g to make it unbiased.This logarithmic formulation avoids relying on the long-tailed extremes of the approximately log-normal Z distribution.Source: §2.1 · p. 3

2.2 CALCULATION OF THE OPTIMAL g VALUES

  • Gaussian matrix projections make each normalized back-propagated vector Gaussian and z chi-squared distributed, enabling layer-independent calculations of the optimal g.For linear and ReLU networks, the analysis uses the Gaussian structure of the matrix projections; tanh is treated numerically.Source: §2.2 · pp. 3–4
  • The derived optimal scaling and random-walk variance are accurate across tested widths, with variance slope inversely proportional to N.The formulas are lowest-order 1/N approximations, but numerical studies found them more accurate than expressions retaining 1/N^2 terms across the tested range.Source: §2.2 · p. 4
  • For ReLU, the derivative randomly removes about half the matrix rows, while tanh lacks an analytic calculation and is expected to require an intermediate g.ReLU analysis averages over binomially distributed active-row counts and excludes zero-activity cases in small layers; tanh uses numerical results.Source: §2.2 · pp. 4–5

2.3 COMPUTATIONAL VERIFICATION

  • Simulations confirm unbiased linear-network walks and linear variance growth, while predicted optimal g values match numerical data across linear, ReLU, and tanh networks.For tanh, saturation produces a wider serviceable range of g than in the linear or ReLU cases.Source: §2.3 · p. 5, Figure 1 · p. 5, Figure 2 · p. 6

3 RESULTS OF TRAINING DEEP NETWORKS WITH RANDOM WALK INITIALIZATION

The section specifies Random Walk Initialization and evaluates it on deep feedforward networks using controlled experiments on MNIST and TIMIT. Correct scaling of g supports training across substantial depths, although greater depth does not improve training error on the examined tasks.

3.1 RANDOM WALK INITIALIZATION

  • g is set using analytical values for linear and ReLU networks, while tanh works well in practice for g between 1.1 and 1.3; input and output scaling may also require adjustment.The procedure tunes input scaling, final-layer scaling, and the network-wide g, with g identified as the most important parameter.Source: §3.1 · p. 6

3.2 EXPERIMENTAL METHODS

  • Experiments evaluate MNIST and TIMIT classification plus MNIST auto-encoder reconstruction, focusing on training error rather than generalization.Networks use a standard feedforward architecture, with parameter counts controlled across depths.Source: §3.2 · pp. 6–7
  • Learning rates are varied by depth through an exponential schedule so shallow and deep networks receive a fair comparison.The schedule is anchored at λout in the output layer and uses λin at the first layer only when D equals Dmax.Source: §3.2 · p. 7

3.3 PERFORMANCE RESULTS

  • For MNIST, minimum training error occurs near g = 1.1–1.4 with tanh and g = 1.4–1.55 with ReLU, matching analytical predictions.Figure 3 reports these results across depths and learning-rate combinations under a 4e6 parameter limit.Source: §3.3 · p. 7, Figure 3 · p. 8
  • Depth 2 gives the best MNIST classification training error, while depths 16 and 32 tie for second; the shallowest auto-encoder is also best, but depth 128 is only roughly 2x worse.A 1000-layer MNIST network achieved about 50 training mistakes, and TIMIT performed best at depth 16 with depth 32 essentially tied.Source: §3.3 · p. 8, Figure 4 · p. 9
  • Across the examined tasks, increasing depth did not improve training error, but the results provide evidence that Random Walk Initialization is reasonable for very deep nonlinear feedforward networks.The experiments compare networks under equal parameter limits and varied learning-rate and g settings.Source: §3.3 · p. 8, Figure 4 · p. 9

4 DISCUSSION

Random Walk Initialization enables correctly initialized feedforward networks to train at depths above 200 layers, while extremely deep first-order optimization remains constrained by curvature and sensitivity. The experiments support using the prescribed g values rather than g = 1, although greater depth did not clearly improve training error.

Experimental results

  • The initialization scheme allowed training to proceed for networks deeper than 200 layers, although depth did not clearly improve training error.Source: §4 · pp. 9–10

Practical implications

  • Increasing layer width N decreases fluctuations in the norm of back-propagated errors.Source: §4 · p. 9

Optimization limits

  • For extremely deep networks around 1000 layers, curvature can make effective optimization with a first-order method impossible.The authors set g = 1.05 in Figure 4C to address this issue.Source: §4 · pp. 9–10

Experimental results

  • Correctly chosen g values enabled nearly all non-broken models to reach zero or near-zero training classification error, or extremely low reconstruction error, regardless of depth.Models could fail when learning-rate hyperparameters mismatched the architecture.Source: §4 · p. 10

Practical implications

  • The authors recommend Random Walk Initialization with g from Figure 2 or the stated equations instead of g = 1 as an easily implemented default.Source: §4 · pp. 9–10
Loading 1412.6558v3…