Title: CAPACITY AND TRAINABILITY IN RECURRENT NEURAL NETWORKS

Abstract: Two potential bottlenecks on the expressiveness of recurrent neural networks (RNNs) are their ability to store information about the task in their parameters, and to store information about the input history in their units. We show experimentally that all common RNN architectures achieve nearly the same per-task and per-unit capacity bounds with careful training, for a variety of tasks and stacking depths. They can store an amount of task information which is linear in the number of parameters, and is approximately 5 bits per parameter. They can additionally store approximately one real number from their input history per hidden unit. We further find that for several tasks it is the per-task parameter capacity bound that determines performance. These results suggest that many previous results comparing RNN architectures are driven primarily by differences in training effectiveness, rather than differences in capacity. Supporting this observation, we compare training difficulty for several architectures, and show that vanilla RNNs are far more difficult to train, yet have slightly higher capacity. Finally, we propose two novel RNN architectures, one of which is easier to train than the LSTM or GRU for deeply stacked architectures.

1. INTRODUCTION:

Research and application of recurrent neural networks (RNNs) have seen explosive growth over the last few years, (Martens & Sutskever, 2011; Graves et al., 2009) , and RNNs have become the central component for some very successful model classes and application domains in deep learning (speech recognition (Amodei et al., 2015) , seq2seq (Sutskever et al., 2014) , neural machine translation (Bahdanau et al., 2014) , the DRAW model (Gregor et al., 2015) , educational applications (Piech et al., 2015) , and scientific discovery (Mante et al., 2013) ). Despite these recent successes, it is widely acknowledged that designing and training the RNN components in complex models can be extremely tricky. Painfully acquired RNN expertise is still crucial to the success of most projects.

One of the main strategies involved in the deployment of RNN models is the use of the Long Short Term Memory (LSTM) networks (Hochreiter & Schmidhuber, 1997) , and more recently the Gated Recurrent Unit (GRU) proposed by Cho et al. (2014) ; Chung et al. (2014) (we refer to these as gated architectures). The resulting models are perceived as being more easily trained, and achieving lower error. While it is widely appreciated that RNNs are universal approximators (Doya, 1993) , an unresolved question is the degree to which gated models are more computationally powerful in practice, as opposed to simply being easier to train.

Here we provide evidence that the observed superiority of gated models over vanilla RNN models is almost exclusively driven by trainability. First we describe two types of capacity bottlenecks that various RNN architectures might be expected to suffer from: parameter efficiency related to learning the task, and the ability to remember input history. Next, we describe our experimental setup where we disentangle the effects of these two bottlenecks, including training with extremely thorough hyperparameter (HP) optimization. Finally, we describe our capacity experiment results (per-parameter and per-unit), as well as the results of trainability experiments (training on extremely hard tasks where gated models might reasonably be expected to perform better).

1.1. CAPACITY BOTTLENECKS:

There are several potential bottlenecks for RNNs, for example: How much information about the task can they store in their parameters? How much information about the input history can they store in their units? These first two bottlenecks can both be seen as memory capacities (one for the task, one for the inputs), for different types of memory.

Another, different kind of capacity stems from the set of computational primitives an RNN is able to perform. For example, maybe one wants to multiply two numbers. In terms of number of units and time steps, this task may be very straight-forward using some specific computational primitives and dynamics, but with others it may be extremely resource heavy. One might expect that differences in computational capacity due to different computational primitives would play a large role in performance. However, despite the fact that the gated architectures are outfitted with a multiplicative primitive between hidden units, while the vanilla RNN is not, we found no evidence of a computational bottleneck in our experiments. We therefore will focus only on the per-parameter capacity of an RNN to learn about its task during training, and on the per-unit memory capacity of an RNN to remember its inputs.

1.2. EXPERIMENTAL SETUP:

RNNs have many HPs, such as the scalings of matrices and biases, and the functional form of certain nonlinearities. There are additionally many HPs involved in training, such as the choice of optimizer, and the learning rate schedule. In order to train our models we employed a HP tuner that uses a Gaussian Process model similar to Spearmint (see Appendix, section on HP tuning and Desautels et al. (2014) ; Snoek et al. (2012) for related work). The basic idea is that one requests HP values from the tuner, runs the optimization to completion using those values, and then returns the validation loss. This loss is then used by the tuner, in combination with previously reported losses, to choose new HP values such that over many experiments, the validation loss is minimized with respect to the HPs. For our experiments, we report the evaluation loss (separate from the validation loss returned to the HP optimizer, except where otherwise noted) after the HP tuner has highly optimized the task (hundreds to many thousands of experiments for each architecture and task).

In our studies we used a variety of well-known RNN architectures: standard RNNs such as the vanilla RNN and the newer IRNN (Le et al., 2015) , as well as gated RNN architectures such as the GRU and LSTM. We rounded out our set of models by innovating two novel (to our knowledge) RNN architectures (see Section 1.4) we call the Update Gate RNN (UGRNN), and the Intersection RNN (+RNN). The UGRNN is a 'minimally gated' RNN architecture that has only a coupled gate between the recurrent hidden state, and the update to the hidden state. The +RNN uses coupled gates to gate both the recurrent and depth dimensions in a straightforward way.

To further explore the various strengths and weaknesses of each RNN architecture, we also used a variety of network depths: 1, 2, 4, 8, in our experiments. 1 In most experiments, we held the number of parameters fixed across different architectures and different depths. More precisely, for a given experiment, a maximum number of parameters was set, along with an input and output dimension. The number of hidden units per layer was then chosen such that the number of parameters, summed across all layers of the network, was as large as possible without exceeding the allowed maximum.

For each of our 6 tasks, 6 RNN variants, 4 depths, and 6+ model sizes, we ran the HP tuner in order to optimize the relevant loss function. Typically this resulted in many hundreds to several thousands of HP evaluations, each of which was a full training run up to millions of training steps. Taken together, this amounted to CPU-millennia worth of computation.

1.3. RELATED WORK:

While it is well known that RNNs are universal approximators of arbitrary dynamical systems (Doya, 1993) , there is little theoretical work on the task-capacity of RNNs. Koiran & Sontag (1998) studied the VC dimension of RNNs, which provides an upper bound on their task-capacity (defined in Section 2.1). These upper bounds are not a close match to our experimental results. For instance, we find that performance saturates rapidly in terms of the number of unrolling steps (Figure 2b ), while the relevant bound increases linearly with the number of unrolling steps. "Unrolling" refers to recurrent computation through time.

Empirically, Karpathy et al. (2015) have studied how LSTMs encode information in character-based text modeling tasks. Further, Sussillo & Barak (2013) have reverse-engineered the vanilla RNN trained on simple tasks, using the tools and language of nonlinear dynamical systems theory. In Foerster et al. (2016) the behavior of switched affine recurrent networks is carefully examined.

The ability of RNNs to store information about their input has been better studied, in both the context of machine learning and theoretical neuroscience. Previous work on short term memory traces explores the tradeoffs between memory fidelity and duration, for the case that a new input is presented to the RNN at every time step (Jaeger & Haas, 2004; Maass et al., 2002; White et al., 2004; Ganguli et al., 2008; Charles et al., 2014) . We use a simpler capacity measure consisting only of the ability of an RNN to store a single input vector. Our results suggest that, contrary to common belief, the capacity of RNNs to remember their input history is not a practical limiting factor on their performance.

The precise details of what makes an RNN architecture perform well is an extremely active research field (e.g. Jozefowicz et al. (2015) ). A highly related article is Greff et al. (2015) , in which the authors used random search of HPs, along with systematic removal of pieces of the LSTM architecture to determine which pieces of the LSTM were more important than the others. Our UGRNN architecture is directly inspired by the large impact of removing the forget gate from the LSTM (Gers et al., 1999) . Zhou et al. (2016) introduced an architecture with minimal gating that is similar to the UGRNN, but is directly inspired by the GRU. An in-depth comparison between RNNs and GRUs in the context of end-to-end speech recognition and a limited computational budget was conducted in Amodei et al. (2015) . Further, ideas from RNN architectures that improve ease of training, such as forget gates (Gers et al., 1999) , and copying recurrent state from one time step to another, are making their way into deep feed-forward networks as highway networks (Srivastava et al., 2015) and residual connections (He et al., 2015) , respectively. Indeed, the +RNN was inspired in part by the coupled depth gate of Srivastava et al. (2015) .

1.4. RECURRENT NEURAL NETWORK ARCHITECTURES:

Below we briefly define the RNN architectures used in this study. Unless otherwise stated W denotes a matrix, b denotes a vector of biases. The symbol x t is the input at time t, and h t is the hidden state at time t. Remaining vector variables represent intermediate values. The function σ(•) denotes the logistic sigmoid function and s(•) is either tanh or ReLU, set as a HP (see Appendix, Section RNN HPs for the complete list of HPs). Initial conditions for the networks were set to a learned bias. Finally, it is a well-known trick of the trade to initialize the gates of an LSTM or GRU with a large bias to induce better gradient flow. We included this parameter, denoted as b f g , and tuned it along with all other HPs. RNN, IRNN (LE ET AL., 2015)

h t = s W h h t-1 + W x x t + b h (1)

Note the IRNN is identical in structure to the vanilla RNN, but with an identity initialization for W h , zero initialization for the biases, and s = ReLU only.

UGRNN -UPDATE GATE RNN:

Based on Greff et al. (2015) , where they noticed the forget gate "was crucial" to LSTM performance, we tried an RNN variant where we began with a vanilla RNN and added a single gate. This gate determines whether the hidden state is carried over from the previous time step, or updated -hence, it is an update gate. An alternative way to view the UGRNN is a highway layer gated through time (Srivastava et al., 2015) .

c t = s W ch h t-1 + W cx x t + b c

(2)

g t = σ W gh h t-1 + W gx x t + b g + b f g (3) h t = g t • h t-1 + (1 -g t ) • c t (4)

GRU -GATED RECURRENT UNIT (CHO ET AL., 2014)

r t = σ W rh h t-1 + W rx x t + b r (5) u t = σ W uh h t-1 + W ux x t + b u + b f g (6) c t = s W ch (r t • h t-1 ) + W cx x t + b c (7) h t = u t • h t-1 + (1 -u t ) • c t (8)

LSTM -LONG SHORT TERM MEMORY(HOCHREITER & SCHMIDHUBER, 1997)

EQUATION

)

EQUATION

+RNN -INTERSECTION RNN:

Due to the success of the UGRNN for shallower architectures in this study (see later figures on trainability), as well as some of the observed trainability problems for both the LSTM and GRU for deeper architectures (e.g. Figure 4h ) we developed the Intersection RNN (denoted with a '+') architecture with a coupled depth gate in addition to a coupled recurrent gate. Additional influences for this architecture were the recurrent gating of the LSTM and GRU, and the depth gating from the highway network (Srivastava et al., 2015) . This architecture has recurrent input, h t-1 , and depth input, x t . It also has recurrent output, h t , and depth output, y t . Note that this architecture only applies between layers where x t and y t have the same dimension, and is not appropriate for networks with a depth of 1 (we exclude depth one +RNNs in our experiments).

EQUATION

EQUATION

EQUATION

EQUATION

)

h t = g h t • h t-1 + (1 -g h t ) • h in t (20)

In practice we used ReLU for s1 and tanh for s2.

2. CAPACITY EXPERIMENTS:



2.1. PER-PARAMETER CAPACITY:

A foundational result in machine learning is that a single-layer perceptron with N 2 parameters can store at least 2 bits of information per parameter (Cover, 1965; Gardner, 1988; Baldi & Venkatesh, 1987) . More precisely, a perceptron can implement a mapping from 2N , N -dimensional, input vectors to arbitrary N -dimensional binary output vectors, subject only to the extremely weak restriction that the input vectors be in general position. RNNs provide a far more complex input-output mapping, with hidden units, recurrent dynamics, and a diversity of nonlinearities. Nonetheless, we wondered if there were analogous capacity results for RNNs that we might be able to observe empirically.

As we will show in Section 3, tasks with complex temporal dynamics, such as language modeling, exhibit a per-parameter capacity bottleneck that explains the performance of RNNs far better than a per-unit bottleneck. To make the experimental design as simple as possible, and to remove potential confounds stemming from the choice of temporal dynamics, we study per-parameter capacity using a task inspired by Gardner (1988) . Specifically, to measure how much task-related information can be stored in the parameters of an RNN, we use a memorization task, where a random static input is injected into an RNN, and a random static output is read out some number of time steps later. We emphasize that the same per-parameter bottleneck that we find in this simplified task also arises in more temporally complex tasks, such as language modeling. At a high level, we draw a fixed set of random inputs and random labels, and train the RNN to map random inputs to randomly chosen labels via cross-entropy error. However, rather than returning the cross-entropy error to the HP tuner (as is normally done), we instead return the mutual information between the RNN outputs and the true labels. In this way, we can treat the number of input-output mappings as a HP, and the tuner will select for us the correct number of mappings so as to maximize the mutual information between the RNN outputs and the labels. From this mutual information we compute bits per parameter, which provides a normalized measurement of how much the RNN learned about the task.

More precisely, we draw datasets of binary inputs X and target binary labels Y at uniform from the set of all binary datasets, X ∼

EQUATION

, where b is the number of samples, and n in is the dimensionality of the inputs. Number of samples, b, is treated as a HP and in practice the optimal dataset size is very close to the bits of mutual information between true and predicted labels. This trend is demonstrated in Figure App.1 in the Appendix. For each value of b the RNN is trained to minimize the cross entropy of the network output with the true labels. We write the output of the RNN for all inputs as Ŷ = f (X), with corresponding random variable Ŷ. We are interested in the mutual information I Y; Ŷ between the true class labels and the class labels predicted by the RNN. This is the amount of (directly recoverable) information that the RNN has stored about the task. In this setting, it is calculated as

EQUATION

)

where p is the fraction of correctly classified samples. The number b is then adjusted, along with all the other HPs, so as to maximize the mutual information I Y; Ŷ . In practice p is computed using only a single draw of {X, Y} ∼ X × Y.

We performed this optimization of I Y; Ŷ for various RNN architectures, depths, and numbers of parameters. We plot the best value of I Y; Ŷ vs. number of parameters in Figure 1a . This captures the amount of information stored in the parameters about the mapping between X and Y. To get an estimate of bits per parameter, we divide by the number of parameters, as shown in Figure 1e .

2.1.2. RESULTS:

Five Bits per Parameter Examining the results of Figure 1 , we find the capacity of all architectures is roughly linear in the number of parameters, across several orders of magnitude of parameter count.

We further find that the capacity is between 3 and 6 bits per parameter, once again across all architectures, depths 1, 2 and 4, and across several orders of magnitude in terms of number of parameters. Given the possibility of small size effects, and a larger portion of weights used as biases at a small number of parameters, we believe our estimates for larger networks are more reliable. This leads us to a bits per parameter estimate of approximately 5, averaging over all architectures and all depths. Finally, we note that the per-parameter task capacity increases as a function of the number of unrollings, though with diminishing gains (Figure 2b ).

The finding that our results are consistent across diverse architectures and scales is even more surprising, since prior to these experiments it was not clear that capacity would even scale linearly with the number of parameters. For instance, previous results on model compression -by reducing the number of parameters (Yang et al., 2015) , or by reducing the bit depth of parameters (Hubara et al., 2016 ) -might lead one to predict that different architectures use parameters with vastly different efficiencies, and that task capacity increases only sublinearly with parameter count.

Gating Slightly Reduces Capacity While overall, the different architectures performed very similarly, there are some capacity differences between architectures that appear to hold up across most depths and parameter counts. To quantify these differences we constructed a table showing the change in the number of parameters one would need to switch from one architecture to another, while maintaining equivalent capacity (Figure 1i ). One trend that emerged from our capacity experiments is a slightly reduced capacity as a function of "gatedness". Putting aside the IRNN, which performed the worst and is discussed below, we noticed that across all depths and all model sizes, the performance was on average RNN > UGRNN > GRU > LSTM > +RNN. The vanilla RNN has no gates, the UGRNN has one, while the remaining three have two or more.

(a) (b) (c) (d) (e) (f) (g) (h)

All depths 1 layer 2 layers 4 layers (i) ReLUs Reduce Capacity In our capacity tasks, the IRNN performed noticeably worse than all other architectures, reaching a maximum bits per parameter of roughly 3.5. To determine if this performance drop was due to the ReLU nonlinearity of the IRNN, or its identity initialization, we sorted through the RNN and UGRNN results (which both have ReLU and tanh as choices for the nonlinearity HP) and looked at the maximum bits per parameter when only optimizations using ReLU architectures of all depths on the memory throughput task. The curve shows the error plotted as a function of the number of units for a random input of dimension 64 (black vertical line). All networks with with less than 64 units have error in reconstruction, while all networks with number of units greater than 64 nearly perfectly reconstruct the random input.

are considered. Indeed, both the RNN and UGRNN bits per parameter dropped dramatically to the 3.5 range (Figure 2a ) when those architectures exclusively used ReLU, providing strong evidence that the ReLU activation function is problematic for this capacity task.

2.2. PER-UNIT CAPACITY TO REMEMBER INPUTS:

An additional capacity bottleneck in RNNs is their ability to store information about their inputs over time. It may be plainly obvious that an IRNN, which is essentially an integrator, can achieve perfect memory of its inputs if the number of inputs is less than or equal to the number of hidden units, but it is not so clear for some of the more complex architectures. So we measured the per-unit input memory empirically. Figure 2c shows the intuitive result that every RNN architecture (at every depth and number of parameters) we studied can reconstruct a random n in dimensional input at some time in the future, if and only if the number of hidden units per layer in the network, n h , is greater than or equal to n in Moreover, regardless of RNN architecture, the error in reconstructing the input follows the same curve as a function of the number of hidden units for all RNN variants, corresponding to reconstructing an n h dimensional subspace of the n in dimensional input.

We highlight this per-unit capacity to make the point that a per-parameter task capacity appears to be the limiting factor in our experiments (e.g. Figure 1 and Figure 3 ), and not a per-unit capacity, such as the per-unit capacity to remember previous inputs. Thus when comparing results between architectures, one should normalize different architectures by the number of parameters, and not the number of units, as is frequently done in the literature (e.g. when comparing vanilla RNNs to LSTMs). This makes further sense as, for all common RNN architectures, the computational cost of processing a single sample is linear in the number of parameters, and quadratic in the number of units per layer. As we show in Figure 3d , plotting the capacity results by numbers of units gives very misleading results.

3. ADDITIONAL TASKS WHERE ARCHITECTURES ACHIEVE VERY SIMILAR LOSS:

We studied additional tasks that we believed to be easy enough to train that the evaluation loss of different architectures would reveal variations in capacity rather than trainability. A critical aspect of these tasks is that they could not be learned perfectly by any of the model sizes in our experiments.

As we change model size, we therefore expect performance on the task to also change. The tasks are (see Appendix, section Task Definitions for further elaboration of these tasks):

• text8 -1-step ahead character-based prediction on the text8 Wikipedia dataset (100 million characters) (Mahoney, 2011) . • Random Continuous Functions (RCF) -A task similar to the per-parameter capacity task above, except the target outputs are real numbers (not categorical), and the number of training samples is held fixed.

The performance on these two tasks is shown in Figure 3 . The evaluation loss as a function of the number of parameters is plotted in panels a-c and e-g, for the text8 task, and RCF task, respectively. For all tasks in this section, the number of parameters rather than the number of units provided the bottleneck on performance, and all architectures performed extremely closely for the same number of parameters. By close performance we mean that, for one model to achieve the same loss as another the model, the number of parameters would have to be adjusted by only a small factor (exemplified in Figure 1i for the per-parameter capacity task). 

4. TASKS THAT ARE VERY HARD TO LEARN:

In practice it is widely appreciated that there is often a significant gap in performance between, for example, the LSTM and the vanilla RNN, with the LSTM nearly always outperforming the vanilla RNN. Our per-parameter capacity results provide evidence for a rough equivalence among a variety of RNN architectures, with slightly higher capacity in the vanilla RNN (Figure 1 ). To reconcile our per-parameter capacity results with widely held experience, we provide evidence that gated architectures, such as the LSTM, are far easier to train than the vanilla RNN (and often the IRNN).

We study two tasks that are difficult to learn: parallel parentheses counting of independent input streams, and mathematical addition of integers encoded in a character string (see Appendix, section Task Definitions). The parentheses task is moderately difficult to learn, while the arithmetic task is quite hard. The results of the HP optimizations are shown in Figure 4a -4h for the parentheses task, and in Figure 4i -4p for the arithmetic task. These tasks show that, while it is possible for a vanilla RNN to learn these tasks reasonably well, it is far more difficult than for a gated architecture. Note that the best achieved loss on the arithmetic task is still significantly decreasing, even after 2500 HP evaluations (2500 full complete optimizations over the training set), for the RNN and IRNN.

There are three noteworthy trends in these trainability experiments. First, across both tasks, and all depths (1, 2, 4 and 8), the RNN and IRNN performed most poorly, and took the longest to learn the task. Note, however that both the RNN and IRNN always solved the tasks eventually, at least for depth 1. Second, as the stacking depth increased, the gated architectures became the only architectures that could solve the tasks. Third, the most trainable architecture for depth 1 was the GRU, and the most trainable architecture for depth 8 was the +RNN (which performed the best on both of our metrics for trainability, on both tasks).

rnn irnn ugrnn gru lstm +rnn (i-p) Same as (a-h), except for the arithmetic task. We note that the best loss for the vanilla RNN is still decreasing after 2400+ HP evaluations.

To achieve our results on capacity and trainability, we relied heavily on a HP tuner. Most practitioners do not have the time or resources to make use of such a tuner, typically only adjusting the HPs a few times themselves. So we wondered how the various architectures would perform if we set HPs randomly, within the ranges specified (see Appendix for ranges). We tried this 1000 times on the parentheses task, for all 200k parameter architectures at depths 1 and 8 (Figure 5 and Table 1 ). The noticeable trends are that the IRNN returned an infeasible error nearly half of the time, and the LSTM (depth 1) and GRU (depth 8) were infeasible the least number of times, where infeasibility means that the training loss diverged. For depth 1, the GRU gave the smallest error, and the smallest median error, and for depth 8, the +RNN delivered the smallest error and smallest median error. 

5. DISCUSSION:

Here we report that a number of RNN variants can hold between 3-6 bits per parameter about their task, and that these variants can remember a number of random inputs that is nearly equal to the number of hidden units in the RNN. The quantification of the number of bits per parameter an RNN can store about a task is particularly important, as it was not previously known whether the amount of information about a task that could be stored was even linear in the number of parameters.

While our results point to empirical capacity limits for both task memorization, and input memorization, apparently the requirement to remember features of the input through time is not a practical bottleneck. If it were, then the vanilla RNN and IRNN would perform better than the gated architectures in proportion to the ratio of the number of units, which they do not. Based on widespread results in the literature, and our own results on our difficult tasks, the loss of some memory capacity (and possibly a small amount of per-parameter storage capacity) for improved trainability seems a worthwhile trade off. Indeed, the input memory capacity did not obviously impact any task not explicitly designed to measure it, as the error curves -for instance for the language modeling taskoverlapped across architectures for the same number of parameters, but not the same number of units.

Our result on per-parameter task capacity, about 5 bits per parameter averaged over architectures, is in surprising agreement with recently published results on the capacity of synapses in biological neurons. This number was recently calculated to be about 4.7 bits per synapse, based on biological synapses in the hippocampus having roughly 26 measurable discrete sizes (Bartol et al., 2016) . Our capacity results have implications for compressed networks that employ quantization techniques. In particular, they provide an estimate of the number of bits which a weight may be compressed without loss in task performance. Coincidentally, in Han et al. (2015) , the authors used 5 bits per weight in the fully connected layers.

An additional observation about per-parameter task capacity in our experiments is that it increases for a few time steps beyond one (Figure 2b ), and then appears to saturate. We interpret this to suggest that recurrence endows additional capacity to a network with shared parameters, but that there are diminishing returns, and the total capacity remains bounded even as the number of time steps increases.

We also note that performance is nearly constant across RNN architectures if the number of parameters is held fixed. This may motivate the design and use of architectures with small compute per parameter ratios, such as mixture of experts RNNs (Shazeer et al., 2017) , and RNNs with large embedding dictionaries on input and output (Józefowicz et al., 2016) .

Despite our best efforts, we cannot claim that we perfectly trained any of the models. Potential problems in HP optimization could be local minima, as well as stochastic behavior in the HP optimization as a result of the stochasticity of batching or random draws for weight matrices. We tried to uncover these effects by running the best performing HPs 100 times, and did not observe any serious deviations from the best results (see Table App .1 in Appendix). Another form of validation comes from the fact that in our capacity task, essentially 3 independent experiments (one for each level of depth) yielded a clustering by architecture (Figure 1e ).

Do our results yield a framework for choosing a recurrent architecture? In total, we believe yes. As explored in Amodei et al. (2015) , a practical concern for recurrent models is speed of execution in a production environment. Our results suggest that if one has a large resource budget for training and confined resource budget for inference, one should choose the vanilla RNN. Conversely, if the training resource budget is small, but the inference budget large, one should choose a gated model. Another serious concern relates to task complexity. If the task is easy to learn, a vanilla RNN should yield good results. However if the task is even moderately difficult to learn, a gated architecture is the right choice. Our results point to the GRU as being the most learnable of gated RNNs for shallow architectures, followed by the UGRNN. The +RNN typically performed best for deeper architectures. Our results on trainability confirm the widely held view that the LSTM is an extremely reliable architecture, but it was almost never the best performer in our experiments. Of course further experiments will be required to fully vet the UGRNN and +RNN. All things considered, in an uncertain training environment, our results suggest using the GRU or +RNN.

• A HP determined whether the input vector X •j was presented to the RNN only at the first time step, or whether it was presented at every time step.

Some optimization algorithms had additional parameters such as ADAM's second order decay rate, or epsilon parameter. These were set to their default values and not optimized. The batch size was set individually by hand for all experiments. The same seed was used to initialize the random number generator for all task parameters, whereas the generator was randomly seeded for network parameters (e.g. initializations). Note that for each network, the initial condition was set to a learned vector.

Figure App.1: In the capacity task, the optimal dataset size found by the HP tuner was only slightly larger than the mutual information in bits reported in Figure 1a , for all architectures at all sizes and depths.

B TASK DEFINITIONS PERCEPTRON CAPACITY:

While at a high-level, for the perceptron capacity task, we wanted to optimize the amount of information the RNN carried about true random labels, in practice, the training objective was standard cross-entropy. However, when returning a validation loss to the HP tuner, we returned the mutual information I Y; Ŷ|X . Conceptually, this is as if there is one nested optimization inside another.

The inner loop optimizes the RNN for the set of HPs, training cross entropy, but returning mutual information. The outer loop then chooses the HPs, in particular, the number of samples b, in equation ( 21), so as to maximize the amount of mutual information. This implementation is necessitated because there is no straightforward way to differentiate mutual information with respect to number of samples. During training, cross entropy error is evaluated beginning after 5 time steps.

MEMORY CAPACITY:

In the Memory Capacity task, we wanted to know how much information an RNN can reconstruct about its inputs at some later time point. We picked an input dimension, 64, and varied the number of parameters in the networks such that the number of hidden units was roughly centered around 64. After 12 time steps the target of the network was exact reconstruction of the input, with a square error loss. The inputs were random values drawn from a uniform distribution between -√ 3 and √ 3 (corresponding to a variance of 1).

RANDOM CONTINUOUS FUNCTION:

A dataset was constructed consisting of N = 10 6 random unit norm Gaussian input vectors x, with size d = 50. Target scalar outputs y were generated for each input vector, and were also drawn from a unit norm Gaussian. Each sample i was assigned a power law weighting β i = (i+τ ) -1

Z:

, where Z was a normalization constant such that the weightings summed to 1, and the characteristic time constant τ = 5000. The loss function for training was calculated after 50 time steps and was weighted square error on the y i , with the β i acting as the weighting terms.

TEXT8:

In the text8 task, the task was to predict one character ahead in the text8 dataset (1e8 characters of Wikipedia) (Mahoney, 2011) . Input was a hot-one encoded sequence, as was the output. The loss was cross-entropy loss on a softmax output layer. Rather than use partial unrolling as is common in language modeling, we generated random pointers into the text. The first 13 time steps (where T = 50) were used to initialize the RNN into a normal operating mode, and remaining steps were used for training or inference.

PARENTHESES COUNTING TASK:

The parentheses counting task independently counts the number of opened 'parens', e.g. '(', without the closing ')'. Here parens is used to mean any of 10 parens type pairs, e.g. '<>' or '[]'. Additionally, there were 10 noise characters, 'a' to 'j'. For each paren type, there was a 20D + 10D = 30D hot-one encoding of all paren and noise symbols, for a total of 300 inputs. The output for each paren type was a hot-one encoding of the digits 0-9, which represented the count of the opened parens of that type. If the count exceeded 9, the the network kept the count at 9, if the paren was closed, the count decreased. The loss was the sum of cross-entropy losses, one for each paren type. Finally, for each paren input stream, 50% random noise characters were drawn, and 50% random paren characters were drawn, e.g. 10 streams like '(a<a<bcb>[[[)'. Parens of other types were treated as noise for the current type, e.g. for the above string if the paren type was '<>', the answer is '1' at the end. The loss was defined only at the final time point, T , and T = 175.

ARITHMETIC TASK:

In the arithmetic task, a hot-one encoded character sequence of an addition problem was presented as input to the network, e.g., '-343243+93851= ', and the output was the hot-one encoded answer, including the correct amount of left padded spaces, '-249392'. An additional HP for this task was the number of compute steps (1-6) between the input of the '=' and the first non-space character in the target output sequence. The two numbers in the input were randomly, uniformly selected in [-1e7, 1e7] . After 36 time steps, cross-entropy loss was calculated. We found this task to be extremely difficult for the networks to learn, but when the task was learned, certain of the network architectures could perform the task nearly perfectly.

C HP ROBUSTNESS:

We wondered how robust the HPs are to the variability of both random batching of data, and random initialization of parameters. So we identified the best HPs from the parentheses experiments of 100k parameter, 1 layer architectures, and reran the parameter optimization 100 times. We measured the number of infeasible experiments, as well as a number of statistics of the loss for the reruns (Table App 

:

Not all experiments used a depth of 8, due to limits on computational resources.

Figures and Tables:

Figure 1: All neural network architectures can store approximately five bits per parameter about a task, with only small variations across architectures. (a) Stored bits as a function of network size. These numbers represent the maximum stored bits across 1000+ HP optimizations with 5 time steps unrolled at each network size for all levels of depth. (b-d) Same as (a), but each level of depth shown separately. (e-h) Same as (a-d) but showing bits per parameter as a function of network size. (i) The value in cell (x, y) is the multiplier for the number of parameters needed to give the architecture on the x-axis the same capacity as the architecture on the y-axis. Capacities are measured by averaging the maximum stored bits per parameter for each architecture across all sizes and levels of depth.

Figure 2: Additional RNN capacity analysis. (a) The effect of the ReLU nonlinearity on capacity.Solid lines indicate bits per parameter for 1-layer architectures (same as Figure1b), where both tanh and ReLU are nonlinearity choices for the HP tuner. Dashed lines show the maximum bits per parameter for each architecture when only results achieved by the ReLU nonlinearity are considered. (b) Bits per parameter as a function of the number of time steps unrolled. (c) L2 error curve for all architectures of all depths on the memory throughput task. The curve shows the error plotted as a function of the number of units for a random input of dimension 64 (black vertical line). All networks with with less than 64 units have error in reconstruction, while all networks with number of units greater than 64 nearly perfectly reconstruct the random input.

Figure3: All RNN architectures achieved near identical performance given the same number of parameters, on a language modeling and random function fitting task. (a-c) text8 Wikipedia number of parameters vs bits per character for all RNN architectures. From left to right: 1 layer, 2 layer, 4 layer models. (d) text8 number of hidden units vs bits per character for 1 layer architectures. We note that this is almost always a misleading way to compare architectures as the more heavily gated architectures appear to do better when compared per-unit. (e-g) Same as (a-c), except showing square error for different model sizes trained on RCFs.

Figure 4: Some RNN architectures are far easier to train than others. Results of HP searches on extremely difficult tasks. (a) Median evaluation error as a function of HP optimization iteration for 1 layer architectures on the parentheses task. Dots indicate evaluation loss achieved on that HP iteration. (b-d) Same as (a), but for 2, 4 and 8 layer architectures. (e-h) Minimum evaluation error as a function of HP optimization iteration for parentheses task. Same depth order as (a-d).(i-p) Same as (a-h), except for the arithmetic task. We note that the best loss for the vanilla RNN is still decreasing after 2400+ HP evaluations.

For randomly generated hyperparameters, GRU and +RNN are the most easily trainable architectures. Evaluation losses from 1000 iterations of randomly chosen HP sets for 1 and 8 layer, 200k parameter models on the parentheses task. Statistics from a Welch's t-test for equality of means on all pairs of architectures are presented in TableApp.2. (a) Box and whisker plot of evaluation losses for the 1 layer model. (b) Same as (a) but for 8 layers. Fraction infeasible trials as a result of 1000 iterations of randomly chosen HP sets for 1 and 8 layer, 200k parameter models trained on the parentheses task.

Table content: <table><tr><td>(a)</td><td>1 layer</td><td>(b)</td><td>8 layers</td></tr><tr><td>Parentheses</td><td/><td/><td/></tr><tr><td colspan="4">Figure 5: Architecture % Infeasible (1 layer) % Infeasible (8 layer)</td></tr><tr><td>+RNN</td><td>-</td><td/><td>8.8 %</td></tr><tr><td>GRU</td><td>15.5 %</td><td/><td>3.2 %</td></tr><tr><td>IRNN</td><td>56.7 %</td><td/><td>44.6 %</td></tr><tr><td>LSTM</td><td>12.0 %</td><td/><td>4.0 %</td></tr><tr><td>RNN</td><td>21.5 %</td><td/><td>18.7 %</td></tr><tr><td>UGRNN</td><td>20.2 %</td><td/><td>11.5 %</td></tr></table>

.1). These results show that the best HPs yielded a distribution of losses very close to the originally reported loss value.Table App.1: Results of 100 runs on the parentheses task using the best HPs for each architecture, at depth 1. HPs were chosen to be the set which achieved the minimum loss. Table shows original loss achieved by the HP tuner, amount of infeasible trials, minimum loss from running 100 iterations of the same HPs, mean loss, maximum loss, standard deviation, and standard deviation divided by the mean. Table App.2: Results of Welch's t-test for equality of means on evaluation losses of architecture pairs trained on the parentheses task with randomly sampled HPs. 8 layer GRU and UGRNN, IRNN and RNN, and LSTM and RNN pairs have loss distributions that are different with statistical significance (p > 0.05). Negative t-statistic indicates that the mean of the second architecture in the pair is larger than the first.

Table content: <table><tr><td/><td/><td/><td>1 layer</td><td/><td/><td/><td>8 layer</td></tr><tr><td/><td/><td>t-stat</td><td>df</td><td>p-value</td><td colspan="2">t-stat</td><td>df</td><td>p-value</td></tr><tr><td/><td>+RNN/GRU</td><td>-</td><td>-</td><td>-</td><td colspan="3">-23.6 1080 &lt; 0.001</td></tr><tr><td/><td>+RNN/IRNN</td><td>-</td><td>-</td><td>-</td><td colspan="2">-25.7</td><td>954</td><td>&lt; 0.001</td></tr><tr><td/><td>+RNN/LSTM</td><td>-</td><td>-</td><td>-</td><td colspan="2">-26.1</td><td>941</td><td>&lt; 0.001</td></tr><tr><td/><td>+RNN/RNN</td><td>-</td><td>-</td><td>-</td><td colspan="2">-25.8</td><td>946</td><td>&lt; 0.001</td></tr><tr><td colspan="2">+RNN/UGRNN</td><td>-</td><td>-</td><td>-</td><td colspan="3">-24.3 1050 &lt; 0.001</td></tr><tr><td/><td>GRU/IRNN</td><td colspan="6">-7.74 696 &lt; 0.001 -3.51 1360 &lt; 0.001</td></tr><tr><td/><td>GRU/LSTM</td><td colspan="6">-6.65 1750 &lt; 0.001 -4.84 1290 &lt; 0.001</td></tr><tr><td/><td>GRU/RNN</td><td colspan="6">-26.5 1340 &lt; 0.001 -3.93 1330 &lt; 0.001</td></tr><tr><td/><td>GRU/UGRNN</td><td colspan="6">-4.11 1620 &lt; 0.001 -1.13 1840</td><td>0.261</td></tr><tr><td/><td>IRNN/LSTM</td><td>2.23</td><td>652</td><td>0.0264</td><td colspan="3">-2.04 1250 0.0420</td></tr><tr><td/><td>IRNN/RNN</td><td colspan="6">-12.7 426 &lt; 0.001 -0.571 1250</td><td>0.568</td></tr><tr><td/><td>IRNN/UGRNN</td><td>4.03</td><td colspan="2">719 &lt; 0.001</td><td colspan="2">2.37</td><td>1320 0.0178</td></tr><tr><td/><td>LSTM/RNN</td><td colspan="3">-19.6 1500 &lt; 0.001</td><td colspan="2">1.53</td><td>1730</td><td>0.125</td></tr><tr><td colspan="5">LSTM/UGRNN 2.25 1640 0.0247</td><td colspan="2">3.81</td><td>1260 &lt; 0.001</td></tr><tr><td/><td>RNN/UGRNN</td><td colspan="3">20.7 1210 &lt; 0.001</td><td colspan="2">2.81</td><td>1300 0.00498</td></tr><tr><td colspan="3">Architecture Original Infeasible</td><td>Min</td><td>Mean</td><td/><td colspan="2">Max</td><td>S.D.</td><td>S.D./Mean</td></tr><tr><td>RNN</td><td>1.16e-2</td><td>0 %</td><td colspan="3">1.41e-2 8.21e-2</td><td colspan="2">0.294</td><td>5.22e-2</td><td>0.636</td></tr><tr><td>IRNN</td><td>4.20e-4</td><td>48 %</td><td colspan="5">2.24e-4 5.02e-4 8.69e-4 1.35e-4</td><td>0.269</td></tr><tr><td>UGRNN</td><td>1.02e-4</td><td>0 %</td><td colspan="5">3.66e-5 2.71e-4 6.06e-3 7.12e-4</td><td>2.63</td></tr><tr><td>GRU</td><td>2.80e-4</td><td>1 %</td><td colspan="5">7.66e-5 1.89e-4 5.48e-4 9.08e-5</td><td>0.480</td></tr><tr><td>LSTM</td><td>7.96e-4</td><td>0 %</td><td colspan="5">8.10e-4 2.02e-3 0.0145 2.31e-3</td><td>1.14</td></tr></table>

Peer review (iclr2017):
The authors investigate a variety of existing and two new RNN architectures to obtain more insight about the effectiveness at which these models can store task information in their parameters and activations. The experimental setups look sound. To generalize comparisons between different architectures it’s necessary to consider multiple tasks and control for the effect of the hyperparameters. This work uses multiple tasks of varying complexities, principled hyperparameter tuning methodology and a number of tuning iterations that can currently only be achieved by the computational resources of some of the larger industrial research groups. The descriptions of the models and the objective where very clear to me. The descriptions of the experiments and presentation of the results were not always clear to me at times, even with the additional details in the appendix available. Most of these issues can easily be resolved by editing the text. For example, in the memory task the scaling of the inputs (and hence also outputs) is not provided so it’s hard to interpret the squared error scores in Figure 2c. It’s not clear to me what the term ‘unrollings’ refers to in Figure 2b. Is this a time lag with additional hidden state updates between the presentation of the input sequence and the generation of the output? Since the perceptron capacity task is somewhat central to the paper, I think a slightly more precise description of how and when the predictions are computed would be helpful. Due to the large number of graphs, it can be somewhat hard to find the most relevant results. Perhaps some of the more obvious findings (like Figure 1(b-d) given Figure 1a) could move to the appendix to make space for more detailed task descriptions. Novelty is not really the aim of this paper since it mostly investigates existing architectures. To use the mutual information to obtain bits per parameter scores in highly non-linear parameterized functions is new to me. The paper also proposed to new architectures that seem to have practical value. The paper adds to the currently still somewhat neglected research effort to employ the larger computational resources we currently have towards a better understanding of architectures which were designed when such resources were not yet available. I’d argue that the paper is original enough for that reason alone. The paper provides some interesting new insights into the properties of RNNs. While observed before, it is interesting to see the importance of gated units for maintaining trainability of networks with many layers. It is also interesting to see a potential new use for vanilla RNNs for simpler tasks where a high capacity per parameter may be required due to hardware constraints. The proposed +RNN may turn out to have practical value as well and the hyperparameter robustness results shed some light on the popularity of certain architectures when limited time for HP tuning is available. The large body of results and hyperparameter analysis should be useful to many researchers who want to use RNNs in the future. All in all, I think this paper would make a valuable addition to the ICLR conference but would benefit from some improvements to the text. Pros: * Thorough analysis. * Seemingly proper experiments. * The way of quantifying capacity in neural networks adds to the novelty of the paper. * The results have some practical value and suggest similar analysis of other architectures. * The results provide useful insights into the relative merits of different RNN architectures. Cons: * It’s hard to isolate the most important findings (some plots seem redundant). * Some relevant experimental details are missing.