Year: 2019

Abstract:
Gating mechanisms are widely used in neural network models, where they allow gradients to backpropagate more easily through depth or time. However, their saturation property introduces problems of its own. For example, in recurrent models these gates need to have outputs near 1 to propagate information over long time-delays, which requires them to operate in their saturation regime and hinders gradient-based learning of the gate mechanism. We address this problem by deriving two synergistic modifications to the standard gating mechanism that are easy to implement, introduce no additional hyperparameters, and improve learnability of the gates when they are close to saturation. We show how these changes are related to and improve on alternative recently proposed gating mechanisms such as chrono-initialization and Ordered Neurons. Empirically, our simple gating mechanisms robustly improve the performance of recurrent models on a range of applications, including synthetic memorization tasks, sequential image classification, language modeling, and reinforcement learning, particularly when long-term dependencies are involved.

heading: 1 INTRODUCTION
text: Recurrent neural networks (RNNs) have become a standard machine learning tool for learning from sequential data. However, RNNs are prone to the vanishing gradient problem, which occurs when the gradients of the recurrent weights become vanishingly small as they get backpropagated through time (Hochreiter et al., 2001). A common approach to alleviate the vanishing gradient problem is to use gating mechanisms, leading to models such as the long short term memory (Hochreiter & Schmidhuber, 1997, LSTM) and gated recurrent units (Chung et al., 2014, GRUs). These gated RNNs have been very successful in several different application areas such as in reinforcement learning (Kapturowski et al., 2018; Espeholt et al., 2018) and natural language processing (Bahdanau et al., 2014; Kočiskỳ et al., 2018).
At every time step, gated recurrent models form a weighted combination of the history summarized by the previous state, and (a function of) the incoming inputs, to create the next state. The values of the gates, which are the coefficients of the combination, control the length of temporal dependencies that can be addressed. This weighted update can be seen as an additive or residual connection on the recurrent state, which helps signals propagate through time without vanishing. However, the gates themselves are prone to a saturating property which can also hamper gradient-based learning. This is particularly troublesome for RNNs, where carrying information for very long time delays requires gates to be very close to their saturated states.
We formulate and address two particular problems that arise with the standard gating mechanism of recurrent models. First, typical initialization of the gates is relatively concentrated. This restricts the range of timescales the model can address, as the timescale of a particular unit is dictated by its gates. Our first proposal, which we call uniform gate initialization (Section 2.2), addresses this by directly initializing the activations of these gates from a distribution that captures a wider spread of dependency lengths.
Second, learning when gates are in their saturation regime is difficult because of vanishing gradients through the gates. We derive a modification that uses an auxiliary refine gate to modulate a main gate, which allows it to have a wider range of activations without gradients vanishing as quickly.
Combining these two independent modifications yields our main proposal, which we call the URgating mechanism. These changes can be applied to any gate (i.e. bounded parametrized function) and have minimal to no overhead in terms of speed, memory, code complexity, and (hyper-)parameters.
We apply them to the forget gate of recurrent models, and evaluate on many benchmarks including synthetic long-term dependency tasks, sequential pixel-level image classification, language modeling, program execution, and reinforcement learning. Finally, we connect our methods to other proposed gating modifications, introduce a framework that allows each component to be replaced with similar ones, and perform theoretical analysis and extensive ablations of our method. Empirically, the UR- gating mechanism robustly improves on the standard forget and input gates of gated recurrent models. When applied to the LSTM, these simple modifications solve synthetic memory tasks that are pathologically difficult for the standard LSTM, achieve state-of-the-art results on sequential MNIST and CIFAR-10, and show consistent improvements in language modeling on the WikiText-103 dataset (Merity et al., 2016) and reinforcement learning tasks (Hung et al., 2018).

heading: 2 GATED RECURRENT NEURAL NETWORKS
text: Broadly speaking, RNNs are used to sweep over a sequence of input data xt to produce a sequence of recurrent states ht ∈Rd summarizing information seen so far. At a high level, an RNN is just a parametrized function in which each sequential application of the network computes a state update u : (xt,ht−1) 7→ht. Gating mechanisms were introduced to address the vanishing gradient problem (Bengio et al., 1994; Hochreiter et al., 2001), and have proven crucial to the success of RNNs. This mechanism essentially smooths out the update using the following equation,
ht=ft(xt,ht−1)◦ht−1+it(xt,ht−1)◦u(xt,ht−1), (1) where the forget gate ft and input gate it are [0,1]d-valued functions that control how fast information is forgotten or allowed into the memory state. When the gates are tied, i.e. ft+it=1 as in GRUs, they behave as a low-pass filter, deciding the time-scale on which the unit will respond (Tallec & Ollivier, 2018). For example, large forget gate activations close to ft=1 are necessary for recurrent models to address long-term dependencies.1
We will introduce our improvements to the gating mechanism primarily in the context of the LSTM, which is the most popular recurrent model. However, these techniques can be used in any model that makes similar use of gates.
ft=σ(Lf (xt, ht−1)) (2) it=σ(Li(xt, ht−1)) (3) ut=tanh(Lu(xt, ht−1)) (4) ct=ft◦ct−1+it◦ut (5) ot=σ(Lo(xt, ht−1)) (6) ht=ottanh(ct) (7) A typical LSTM (equations (2)-(7)) is an RNN whose state is represented by a tuple (ht,ct) consisting of a “hidden” state and “cell” state. The basic gate equation (1) is used to create the next cell state ct (5). Note that the gate and update activations are a function of the previous hidden state ht−1 instead of ct−1. Here,L? stands for a parameterized linear function of its inputs with bias b?, e.g.
Lf (xt,ht−1)=Wfxxt+Wfhht−1+bf , (8) and σ(·) refers to the standard sigmoid activation function which we will assume is used for defining [0,1]-valued activations in the rest of this paper. The gates of the LSTM were initially motivated as a binary mechanism, switching on or off to allow information and gradients to pass through. However, in reality this fails to happen due to a combination of initialization and saturation. This can be problematic, such as when very long dependencies are present.
2.1 THE UR-LSTM
bf ∼σ−1(U [0,1]) (9) ft=σ(Lf (xt, ht−1)) (10) rt=σ(Lr(xt, ht−1)) (11) gt=rt ·(1−(1−ft)2)
+(1−rt)·f2t (12) ct=gtct−1+(1−gt)ut (13) ot=σ(Lo(xt, ht−1)) (14) ht=ottanh(ct) (15)
We present two solutions which work in tandem to address the previously described issues. The first ensures a diverse range of gate values at the start of training by sampling the gate’s biases so that the activations will be approximately uniformly distributed at initialization. We call this Uniform Gate Initialization (UGI). The second allows better gradient flow by reparameterizing the gate using an auxiliary “refine” gate.
As our main application is for recurrent models, we present the full UR-LSTM model in equations (9)-(15). However, we note that
1In this work, we use “gate” to alternatively refer to a [0,1]-valued function or the value (“activation”) of that function.
these methods can be used to modify any gate (or more generally, bounded function) in any model. In this context the UR-LSTM is simply defined by applying UGI and a refine gate r on the original forget gate f to create an effective forget gate g (equation (12)). This effective gate is then used in the cell state update (13). Empirically, these small modifications to an LSTM are enough to allow it to achieve nearly binary activations and solve difficult memory problems (Figure 4). In the rest of Section 2, we provide theoretical justifications for UGI and refine gates.

heading: 2.2 UNIFORM GATE INITIALIZATION
text: Standard initialization schemes for the gates can prevent the learning of long-term temporal correlations (Tallec & Ollivier, 2018). For example, supposing that a unit in the cell state has constant forget gate value ft, then the contribution of an input xt in k time steps will decay by (ft)k. This gives the unit an effective decay period or characteristic timescale ofO( 11−ft ).
2 Standard initialization of linear layersL sets the bias term to 0, which causes the forget gate values (2) to concentrate around 1/2. A common trick of setting the forget gate bias to bf =1.0 (Jozefowicz et al., 2015) does increase the value of the decay period to 11−σ(1.0)≈3.7. However, this is still relatively small, and moreover fixed, hindering the model from easily learning dependencies at varying timescales.
We instead propose to directly control the distribution of forget gates, and hence the corresponding distribution of decay periods. In particular, we propose to simply initialize the value of the forget gate activations ft according to a uniform distribution U(0,1), as described in Section 2.1. An important difference between UGI and standard or other (e.g. Tallec & Ollivier, 2018) initializations is that negative forget biases are allowed. The effect of UGI is that all timescales are covered, from units with very high forget activations remembering information (nearly) indefinitely, to those with low activations focusing solely on the incoming input. Additionally, it introduces no additional parameters; it even can have less hyperparameters than the standard gate initialization, which sometimes tunes the forget bias bf . Appendix B.2 and B.3 further discuss the theoretical effects of UGI on timescales.

heading: 2.3 THE REFINE GATE
text: Given a gate f=σ(Lf (x))∈ [0,1], the refine gate is an independent gate r=σ(Lr(x)), and modulates f to produce a value g∈ [0,1] that will be used in place of f downstream. It is motivated by considering how to modify the output of a gate f in a way that promotes gradient-based learning, derived below.
An additive modification The root of the saturation problem is that the gradient∇f of a gate, which can be written solely as a function of the activation value as f(1−f), decays rapidly as f approaches 0 or 1. Thus when the activation f is past a certain upper or lower threshold, learning effectively stops. This problem cannot be fully addressed only by modifying the input to the sigmoid, as in UGI and other techniques, as the gradient will still vanish by backpropagating through the activation function.
Therefore to better control activations near the saturating regime, instead of changing the input to the sigmoid in f = σ(L(x)), we consider modifying the output. In particular, we consider adjusting f with an input-dependent update φ(f,x) for some function φ, to create an effective gate g=f+φ(f,x) that will be used in place of f downstream such as in the main state update (1). This sort of additive (“residual”) connection is a common technique to increase gradient flow, and indeed was the motivation of the LSTM additive gated update (1) itself (Hochreiter & Schmidhuber, 1997).
Choosing the adjustment function Although many choices seem plausible for selecting the additive update φ, we reason backwards from necessary properties of the effective activation g to deduce a principled function φ. The refine gate will appear as a result.
First, note that ft might need to be increased or decreased, regardless of what its value is. For example, given a large activation ft near saturation, it may need to be even higher to address long-term dependencies in recurrent models; alternatively, if it is too high by initialization or needs to unlearn previous behavior, it may need to decrease. Therefore, the additive update to f should create an effective activation gt in the range ft±α for some α. Note that the allowed adjustment range α=α(ft) needs to be a function of f in order to keep g∈ [0,1].
2This corresponds to the number of timesteps it takes to decay by 1/e.
In particular, the additive adjustment range α(f) should satisfy the following natural properties:
Validity: α(f)≤min(f,1−f), to ensure g∈f±α(f)∈ [0,1]. Symmetry: Since 0 and 1 are completely symmetrical in the gating framework, α(f)=α(1−f). Differentiability: α(f) will be used in backpropagation, requiring α∈C1(R). Figure 2a illustrates the general appearance of α(f) based on these properties. In particular, Validity implies that that its derivative satisfies α′(0) ≤ 1 and α′(1) ≥ −1, Symmetry implies α′(f) = −α′(1−f), and Differentiability implies α′ is continuous. The simplest such function satisfying these is the linear α′(f)=1−2f , yielding α(f)=f−f2=f(1−f). Given such a α(f), recall that the goal is to produce an effective activation g=f+φ(f,x) such that g∈f±α(f) (Figure 2b). Our final observation is that the simplest such function φ satisfying this is φ(f,x)=α(f)ψ(f,x) for some ψ(·)∈ [−1,1]. Using the standard method for defining [−1,1]-valued functions via a tanh non-linearity leads to φ(f,x)=α(f)(2r−1) for another gate r=σ(L(x)). The full update is given in Equation (16),
g=f+α(f)(2r−1)=f+f(1−f)(2r−1)=(1−r)·f2+r·(1−(1−f)2) (16) Equation (16) has the elegant interpretation that the gate r linearly interpolates between the lower band f−α(f)=f2 and the symmetric upper band f+α(f)=1−(1−f)2 (Figure 2b). In other words, the original gate f is the coarse-grained determinant of the effective gate g, while the gate r “refines” it. This allows the effective gate g to reach much higher and lower activations than the constituent gates f and r (Figure 2c), bypassing the saturating gradient problem. For example, this allows the effective forget gate to reach g=0.99 when the forget gate is only f=0.9.

heading: 2.4 REFINING RECURRENT MODELS
text: Formally, the full mechanism of the refine gate as applied to gated recurrent models is defined in equations (11)-(13). Note that it is an isolated change where the forget gate (10) is modified before applying the the standard update (1). Figure 1 illustrates the refine gate in an LSTM cell. Figure 2 illustrates how the refine gate rt is defined and how it changes the effective gate ft to produce an effective gate gt.
Finally, to simplify comparisons and ensure that we always use the same number of parameters as the standard gates, when using the refine gate we tie the input gate to the effective forget gate, it=1−gt. However, we emphasize that these techniques are extremely simple and broad, and can be applied to any gate (or more broadly, any bounded function) to improve initialization distribution and help optimization. For example, our methods can be combined in different ways in recurrent models, e.g. an independent input gate can be modified with its own refine gate. Alternatively, the refine gate can also be initialized uniformly, which we do in our experiments whenever both UGI and refine gates are used.

heading: 2.5 RELATED GATING MECHANISMS
text: We highlight a few recent works that also propose small gate changes to address problems of long-term or variable-length dependencies. Like ours, they can be applied to any gated update equation.
Tallec & Ollivier (2018) suggest an initialization strategy to capture long-term dependencies on the order of Tmax, by sampling the gate biases from bf ∼ log U(1,Tmax−1). Although similar to UGI in definition, chrono initialization (CI) has key differences in the timescales captured, for example by using an explicit timescale parameter and having no negative biases. Due to its relation to UGI, we provide a more detailed comparison in Appendix B.3. As mentioned in Section 2.3, techniques such as these that only modify the input to a sigmoid gate do not fully address the saturation problem.
The Ordered Neuron LSTM introduced by Shen et al. (2018) aims to induce an ordering over the units in the hidden states such that “higher-level” neurons retain information for longer and capture higher-level information. We highlight this work due to its recent success in NLP, and also because its novelties can be factored into introducing two mechanisms which only affect the forget and input gates, namely (i) the cumax:=cumsum◦softmax activation function which creates a monotonically increasing vector in [0,1], and (ii) a pair of “master gates” which are ordered by cumax and fine-tuned with another pair of gates.
In fact, we observe that these are related to our techniques in that one controls the distribution of a gate activation, and the other is an auxiliary gate with modulating behavior. Despite its important novelties, we find that the ON-LSTM has drawbacks including speed/stability issues and theoretical flaws in the scaling of its gates. We provide the formal definition and detailed analysis of the ON-LSTM in Appendix B.4. In particular we flesh out a deeper relation between the master and refine gates and show how they can be interchanged for each other.
We include a more thorough overview of other related works on RNNs in Appendix B.1. These methods are largely orthogonal to the isolated gate changes considered here and are not analyzed. We note that an important drawback common to all other approaches is the introduction of substantial hyperparameters in the form of constants, training protocol, and significant architectural changes. For example, even for chrono initialization, one of the less intrusive proposals, we experimentally find it to be particularly sensitive to the hyperparameter Tmax (Section 3).

heading: 2.6 GATE ABLATIONS
text: Our insights about previous work with related gate components allow us to perform extensive ablations of our contributions. We observe two independent axes of variation, namely, activation function/initialization (cumax, constant bias sigmoid, CI, UGI) and auxiliary modulating gates (master, refine), where different components can be replaced with each other. Therefore we propose several other gate combinations to isolate the effects of different gating mechanisms. We summarize a few ablations here; precise details are given in Appendix B.5. O-: Ordered gates. A natural simplification of the main idea of ON-LSTM, while keeping the hierarchical bias on the forget activations, is to simply drop the auxiliary master gates and define ft,it (2)-(3) using the cumax activation function. UM-: UGI master gates. This variant of the ON-LSTM’s gates ablates the cumax operation on the master gates,
replacing it with a sigmoid activation and UGI which maintains the same initial distribution on the activation values. OR-: Refine instead of master. A final variant in between the UR- gates and the ON-LSTM’s gates combines cumax with refine gates. In this formulation, as in UR- gates, the refine gate modifies the forget gate and the input gate is tied to the effective forget gate. The forget gate is ordered using cumax.
Table 1 summarizes the gating modifications we consider and their naming conventions. Note that we also denote the ON-LSTM method as “OM-LSTM” (M for master) for mnemonic ease. Finally, we remark that all methods here are controlled with the same number of parameters as the standard LSTM, aside from the OM-LSTM and UM-LSTM which use an additional 12C -fraction parameters whereC is the downsize factor on the master gates (Appendix B.4).

heading: 3 EXPERIMENTS
text: We first perform full ablations of the gating variants (Section 2.6) on benchmark synthetic memorization and pixel-by-pixel image classification tasks. We then evaluate our main method on important applications for recurrent models including language modeling and reinforcement learning, comparing against baseline methods where appropriate.
The vanilla LSTM uses forget bias 1.0 (Section 2.2). When chrono-initialization is used and not explicitly tuned, we set Tmax to be proportional to the hidden size. This heuristic uses the intuition that if dependencies of length T exist, then so should dependencies of all lengths≤T . Moreover, the amount of information that can be remembered is proportional to the number of hidden units.
All of our benchmarks have prior work with recurrent baselines, from which we used the same models, protocol, and hyperparameters whenever possible, changing only the gating mechanism. Since our simple gate changes are compatible with other recurrent cores, we evaluate them in tandem with recurrent models such as the GRU, Reconstructive Memory Agent (RMA; Hung et al., 2018), and Relational Memory Core (RMC; Santoro et al., 2018) whenever they were used on these tasks. Full protocols and details for all experiments are given in Appendix D.

heading: 3.1 SYNTHETIC TASKS
text: Our first set of experiments is on synthetic memory tasks (Hochreiter & Schmidhuber, 1997; Arjovsky et al., 2016) that are known to be hard for standard LSTMs to solve.
Copy task. The input is a sequence of N + 20 digits where the first 10 tokens (a0,a1,...,a9) are randomly chosen from {1,...,8}, the middle N tokens are set to 0, and the last ten tokens are 9. The goal of the recurrent model is to output (a0,...,a9) in order on the last 10 time steps, whenever the cue token 9 is presented. We trained our models using cross-entropy with baseline loss log(8) (Appendix D.1).
Adding task. The input consists of two sequences: 1. N numbers (a0,...,aN−1) sampled independently from U [0,1] 2. an index i0∈ [0,N/2) and i1∈ [N/2,N), together encoded as a two-hot sequence. The target output is ai0+ai1 and models are evaluated by the mean squared error with baseline loss 1/6.
Figure 3 shows the loss of various methods on the Copy and Adding tasks. The only gate combinations capable of solving Copy completely are OR-, UR-, O-, and C-LSTM. This confirms the mechanism of their gates: these are the only methods capable of producing high enough forget gate values either through the cumax non-linearity, the refine gate, or extremely high forget biases. The U-LSTM is the only other method able to make progress, but converges slower as it suffers from gate saturation
without the refine gate. The vanilla LSTM makes no progress. The OM-LSTM and UM-LSTM also get stuck at the baseline loss, despite the OM-LSTM’s cumax activation, which we hypothesize is due to the suboptimal magnitudes of the gates at initialization (Appendix B.4). On the Adding task, every method besides the basic LSTM is able to eventually solve it, with all refine gate variants fastest.
Figure 4 shows the distributions of forget gate activations of sigmoid-activation methods, before and after training on the Copy task. It shows that activations near 1.0 are important for a model’s ability to make progress or solve this task, and that adding the refine gate makes this significantly easier.

heading: 3.2 PIXEL-BY-PIXEL IMAGE CLASSIFICATION
text: These tasks involve feeding a recurrent model the pixels of an image in a scanline order before producing a classification label. We test on the sequential MNIST (sMNIST), permuted MNIST (pMNIST) (Le et al., 2015), and sequential CIFAR-10 (sCIFAR) tasks. Each LSTM method was ran with a learning rate sweep with 3 seeds each. The best validation score found over any run is reported in the first two rows of Table 2.3 We find in general that all methods are able to improve over the vanilla LSTM. However, the differences become even more pronounced when stability is considered. Although Table 2 reports the best validation accuracies found on any run, we found that many methods were quite unstable. Asterisks are marked next to a score denoting how many of the 3 seeds diverged, for the learning rate that score was found at.
3sMNIST is not included here as it is too easy, making it difficult to draw conclusions.
Conversely, Figure 5 shows the accuracy curves of each method at their best stable learning rate. The basic LSTM is noticeably worse than all of the others. This suggests that any of the gate modifications, whether better initialization, cumax non-linearity, or master or refine gates, are better than standard gates especially when long-term dependencies are present. Additionally, the uniform gate initialization methods are generally better than the ordered and chrono initialization, and the refine gate performs better than the master gate. We additionally consider applying other techniques developed for recurrent models that are independent of the gating mechanism. Table 2 also reports scores when the same gating mechanisms are applied to the GRU model instead of the LSTM, where similar trends hold across the gating variants. In particular, UR-GRU is the only method that is able to stably attain good performance. As another example, the addition of a generic regularization technique—we chose Zoneout (Krueger et al., 2016) with default hyperparameters (zc=0.5, zh=0.05)—continued improving the UR-LSTM/GRU, outperforming even non-recurrent models on sequential MNIST and CIFAR-10. Table 3 compares the test accuracy of our main model against other models.
From Sections 3.1 and 3.2, we draw a few conclusions about the comparative performance of different gate modifications. First, the refine gate is consistently better than comparable master gates. CI solves the synthetic memory tasks but is worse than any other variant outside of those. We find ordered (cumax) gates to be effective, but speed issues prevent us from using them in more complicated tasks. UR- gates are consistently among the best performing and most stable.
Model Valid. Test
LSTM 34.3 35.8 C-LSTM 35.0 36.4 C-LSTM (Tmax=8) 34.3 36.1 C-LSTM (Tmax=11) 34.6 35.8 OM-LSTM 34.0 34.7 U-LSTM 33.8 34.9 UR-LSTM 33.6 34.6
(a) Perplexities on the WikiText-103 dataset.
0 100000 200000 300000 400000 500000
Iteration
3.52
3.54
3.56
3.58
3.60
Lo g P
e rp
le x it
y
Method C-LSTM
C11-LSTM C8-LSTM LSTM OM-LSTM U-LSTM UR-LSTM
(b) Validation learning curves, illustrating training speed and generalization (i.e. overfitting) behavior.

heading: 3.3 LANGUAGE MODELING
text: We consider word-level language modeling on the WikiText-103 dataset, where (i) the dependency lengths are much shorter than in the synthetic tasks, (ii) language has an implicit hierarchical structure and timescales of varying lengths. We evaluate our gate modifications against the exact hyperparameters of a SOTA LSTM-based baseline (Rae et al., 2018) without additional tuning (Appendix D). Additionally, we compare against ON-LSTM, which was designed for this domain (Shen et al., 2018), and chrono initialization, which addresses dependencies of a particular timescale as opposed to timescaleagnostic UGI methods. In addition to our default hyperparameter-free initialization, we tested models with the chrono hyperparameter Tmax manually set to 8 and 11, values previously used for language modeling to mimic fixed biases of about 1.0 and 2.0 respectively (Tallec & Ollivier, 2018).
Table 6a shows Validation and Test set perplexities for various models. We find that the OM-LSTM, U-LSTM, and UR-LSTM improve over the standard LSTM with no additional tuning. However, although the OM-LSTM was designed to capture the hierarchical nature of language with the cumax activation, it does not perform better than the U-LSTM and UR-LSTM. The chrono initialization with our default initialization strategy is far too large. While manually tweaking the Tmax hyperparameter helps, it is still far from any UGI-based methods. We attribute these observations to the nature of language having dependencies on multiple widely-varying timescales, and that UGI is enough to capture these without resorting to strictly enforced hierarchies such as in OM-LSTM.

heading: 3.4 REINFORCEMENT LEARNING
text: In many partially observable reinforcement learning (RL) tasks, the agent can observe only part of the environment at a time and thus requires a memory model to summarize what it has seen previously. However, designing memory architectures for reinforcement learning problems has been a challenging task (Oh et al., 2016; Wayne et al., 2018). These are usually based on an LSTM core to summarize what an agent has seen into a state.
We investigated if changing the gates of these recurrent cores can improve the performance of RL agents, especially on difficult tasks involving memory and long-term credit assignment. We chose the Passive and Active Image Match tasks from Hung et al. (2018) using A3C agents (Mnih et al., 2016). In these tasks, agents are either initially shown a colored indicator (Passive) or must search for it (Active), before being teleported to a room in which they must press a switch with matching color to receive reward. In between these two phases is an intermediate phase where they can acquire distractor rewards, but the true objective reported is the final reward in the last phase. Episodes last 450−600 steps, so these tasks require memorization and credit assignment across long sequences.
Hung et al. (2018) evaluated agents with different recurrent cores: the basic LSTM, the DNC (an LSTM with memory), and the RMA (which also uses an LSTM core). We modified each of these with our gates. Figure 7 shows the results of different models on the Passive Matching and Active Matching tasks without distractors. These are the most similar to the synthetic tasks (Sec. 3.1), and we found that those trends largely transferred to the RL setting even with several additional confounders present such as agents learning via RL algorithms, being required to learn relevant features from pixels rather than being given the relevant tokens, and being required to explore in the Active Match case.
We found that the UR- gates substantially improved the performance of the basic LSTM core on both Passive Match and Active Match tasks, with or without distractor rewards. On the difficult Active Match task, it was the only method to achieve better than random behavior. Figure 8 shows performance of LSTM and RMA cores on the harder Active Match task with distractors. Here the UR- gates again learn the fastest and reach the highest reward. In particular, although the RMA is a memory architecture with an explicit memory bank designed for long-term credit assignment, its performance was also improved.

heading: 3.5 ADDITIONAL RESULTS AND EXPERIMENTAL CONCLUSIONS
text: Appendix (E.1) shows an additional synthetic experiment investigating the effect of refine gates on saturation. Appendix (E.3) has results on a program execution task, which is interesting for having explicit long and variable-length dependencies and hierarchical structure. It additionally shows another very different gated recurrent model where the UR- gates provide consistent improvement.
Finally, we would like to comment on the longevity of the LSTM, which for example was frequently found to outperform newer competitors when tuned better (Melis et al., 2017). Although many improvements have been suggested over the years, none have been proven to be as robust as the LSTM across an enormously diverse range of sequence modeling tasks. By experimentally starting from well-tuned LSTM baselines, we believe our simple isolated gate modifications to actually be robust improvements. In Appendix B.3 and B.4, we offer a few conclusions for the practitioner about the other gate components considered based on our experimental experience.

heading: 4 DISCUSSION
text: In this work, we introduce, analyze, and evaluate several modifications to the ubiquitous gating mechanism that appears in recurrent neural networks. We describe theoretically-justified methods that improve on the standard gating method by alleviating problems with initialization and optimization. The mechanisms considered include changes on independent axes, namely initialization method and auxiliary gates, and we perform extensive ablations on our improvements with previously considered modifications. Our main gate model robustly improves on standard gates across many different tasks and recurrent cores, while requiring less tuning Finally, we emphasize that these improvements are completely independent of the large body of research on neural network architectures that use gates, and hope that these insights can be applied to improve machine learning models at large.

heading: A PSEUDOCODE
text: We show how the gated update in a typical LSTM implementation can be easily replaced by UR- gates.
The following snippets show pseudocode for the the gated state updates for a vanilla LSTM model (top) and UR-LSTM (bottom).
forget_bias = 1.0 # hyperparameter ... f, i, u, o = Linear(x, prev_hidden) f_ = sigmoid(f + forget_bias) i_ = sigmoid(i) next_cell = f_ * prev_cell + i_ * tanh(u) next_hidden = o * tanh(next_cell)
Listing 1: LSTM
# Initialization u = np.random.uniform(low=1/hidden_size, high=1-1/hidden_size, hidden_size) forget_bias = -np.log(1/u-1) ... # Recurrent update f, r, u, o = Linear(x, prev_hidden) f_ = sigmoid(f + forget_bias) r_ = sigmoid(r - forget_bias) g = 2*r_*f_ + (1-2*r_)*f_**2 next_cell = g * prev_cell + (1-g) * tanh(u) next_hidden = o * tanh(next_cell)
Listing 2: UR-LSTM

heading: B FURTHER DISCUSSION ON RELATED METHODS
text: Section 2.5 briefly introduced chrono initialization (Tallec & Ollivier, 2018) and the ON-LSTM (Shen et al., 2018), closely related methods that modify the gating mechanism of LSTMs. We provide more detailed discussion on these in Sections B.3 and B.4 respectively. Section B.1 has a more thorough overview of related work on recurrent neural networks that address long-term dependencies or saturating gates.

heading: B.1 RELATED WORK
text: Several methods exist for addressing gate saturation or allowing more binary activations. Gulcehre et al. (2016) proposed to use piece-wise linear functions with noise in order to allow the gates to operate in saturated regimes. Li et al. (2018b) instead use the Gumbel trick (Maddison et al., 2016; Jang et al., 2016), a technique for learning discrete variables within a neural network, to train LSTM models with discrete gates. These stochastic approaches can suffer from issues such as gradient estimation bias, unstable training, and limited expressivity from discrete instead of continuous gates. Additionally they require more involved training protocols with an additional temperature hyperparameter that needs to be tuned explicitly.
Alternatively, gates can be removed entirely if strong constraints are imposed on other parts of the model. (Li et al., 2018a) use diagonal weight matrices and require stacked RNN layers to combine information between hidden units. A long line of work has investigated the use of identity or orthogonal initializations and constraints on the recurrent weights to control multiplicative gradients unrolled through time (Le et al., 2015; Arjovsky et al., 2016; Henaff et al., 2016). Chandar et al. (2019) proposed another RNN architecture using additive state updates and non-saturating activation functions instead of gates. However, although these gate-less techniques can be used to alleviate the vanishing gradient problem with RNNs, unbounded activation functions can cause less stable learning dynamics and exploding gradients.
As mentioned, a particular consequence of the inability of gates to approach extrema is that gated recurrent models struggle to capture very long dependencies. These problems have traditionally been addressed by introducing new components to the basic RNN setup. Some techniques include stacking layers in a hierarchy (Chung et al., 2016), adding skip connections and dilations (Koutnik et al., 2014; Chang et al., 2017), using an external memory (Graves et al., 2014; Weston et al., 2014; Wayne et al., 2018; Gulcehre et al., 2017), auxiliary semi-supervision (Trinh et al., 2018), and more. However, these approaches have not been widely adopted over the standard LSTM as they are often specialized for certain tasks, are not as robust, and introduce additional complexity. Recently the transformer model has been successful in many applications areas such as NLP (Radford et al., 2019; Dai et al., 2019). However, recurrent neural networks are still important and commonly used due their faster inference without the need to maintain the entire sequence in memory. We emphasize that the vast majority of proposed RNN changes are completely orthogonal to the simple gate improvements in this work, and we do not focus on them.
A few other recurrent cores that use the basic gated update (1) but use more sophisticated update functions u include the GRU, Reconstructive Memory Agent (RMA; Hung et al., 2018), and Relational Memory Core (RMC; Santoro et al., 2018), which we consider in our experiments.

heading: B.2 EFFECT OF PROPOSED METHODS ON TIMESCALES
text: We briefly review the connection between our methods and the effective timescales that gated RNNs capture. Recall that Section 2.2 defines the characteristic timescale of a neuron with forget activation ft as 1/(1−ft), which would be the number of timesteps it takes to decay that neuron by a constant. The fundamental principle of gated RNNs is that the activations of the gates affects the timescales that the model can address; for example, forget gate activations near 1.0 are necessary to capture long-term dependencies.
Thus, although our methods were defined in terms of activations gt, it is illustrative to reason with their characteristic timescales 1/(1−gt) instead, whence both UGI and refine gate also have clean interpretations.
First, UGI is equivalent to initializing the decay period from a particular heavy-tailed distribution, in contrast to standard initialization with a fixed decay period (1−σ(bf ))−1. Proposition 1. UGI is equivalent to to sampling the decay periodD=1/(1−ft) from a distribution with density proportional to P(D=x)∝ ddx (1−1/x)=x−2, i.e. a Pareto(α=2) distribution.
On the other hand, for any forget gate activation ft with timescale D = 1/(1−ft), the refine gate fine-tunes it betweenD=1/(1−f2t )=1/(1−ft)(1+ft) and 1/(1−ft)2. Proposition 2. Given a forget gate activation with timescaleD, the refine gate creates an effective forget gate with timescale in (D/2, D2).

heading: B.3 CHRONO INITIALIZATION
text: The chrono initialization bf ∼ log(U([1,Tmax−1])) (17) bi=−bf . (18)
was the first to explicitly attempt to initialize the activation of gates across a distributional range. It was motivated by matching the gate activations to the desired timescales.
They also elucidate the benefits of tying the input and forget gates, leading to the simple trick (18) for approximating tying the gates at initialization, which we borrow for UGI. (We remark that perfect tied initialization can be accomplished by fully tying the linear maps Lf ,Li, but (18) is a good approximation.)
However, the main drawback of CI is that the initialization distribution is too heavily biased toward large terms. This leads to empirical consequences such as difficult tuning (due to most units starting in the saturation regime, requiring different learning rates) and high sensitivity to the hyperparameter Tmax that represents the maximum potential length of dependencies. For example, Tallec & Ollivier (2018) set this parameter according to a different protocol for every task, with values ranging from 8 to 2000. Our experiments used a hyperparameter-free method to initialize Tmax (Section 3), and we found that chrono initialization generally severely over-emphasizes long-term dependencies if Tmax is not carefully controlled.
A different workaround suggested by Tallec & Ollivier (2018) is to sample from P(T =k)∝ 1 klog2(k+1) and setting bf =log(T ). Note that such an initialization would be almost equivalent to sampling the decay period from the distribution with density P(D= x)∝ (xlog2x)−1 (since the decay period is (1−f)−1 = 1+exp(bf )). This parameter-free initialization is thus similar in spirit to the uniform gate initialization (Proposition 1), but from a much heavier-tailed distribution that emphasizes very long-term dependencies.
These interpretations suggest that it is plausible to define a family of Pareto-like distributions from which to draw the initial decay periods from, with this distribution treated as a hyperparameter. However, with no additional prior information on the task, we believe the uniform gate initialization to be the best candidate, as it 1. is a simple distribution with easy implementation, 2. has characteristic timescale distributed as an intermediate balance between the heavy-tailed chrono initialization and sharply decaying standard initialization, and 3. is similar to the ON-LSTM’s cumax activation, in particular matching the initialization distribution of the cumax activation.
Table 4 summarizes the decay period distributions at initialization using different activations and initialization strategies.
In general, our experimental recommendation for CI is that it can be better than standard initialization or UGI when certain conditions are met (tasks with long dependencies and nearly fixed-length sequences as in Sections 3.1, 3.4) and/or when it can be explicitly tuned (both the hyperparameter Tmax, as well as the learning rate to compensate for almost all units starting in saturation). Otherwise, we recommend UGI or standard initialization. We found no scenarios where it outperformed UR- gates.

heading: B.4 ON-LSTM
text: In this section we elaborate on the connection between the mechanism of Shen et al. (2018) and our methods. We define the full ON-LSTM and show how its gating mechanisms can be improved. For
example, there is a remarkable connection between its master gates and our refine gates – independently of the derivation of refine gates in Section 2.3, we show how a specific way of fixing the normalization of master gates becomes equivalent to a single refine gate.
First, we formally define the full ON-LSTM. The master gates are a cumax-activation gate
f̃t=cumax(Lf̃ (xt,ht−1)) (19) ĩt=1−cumax(Lĩ(xt,ht−1)). (20)
These combine with an independent pair of forget and input gates ft,it, meant to control fine-grained behavior, to create an effective forget/input gate f̂t ,̂it which are used to update the state (equation (1) or (5)).
ωt= f̃t◦ ĩt (21) f̂t=ft◦ωt+(f̃t−ωt) (22) ît= it◦ωt+(̃it−ωt). (23)
As mentioned in Section B.1, this model modifies the standard forget/input gates in two main ways, namely ordering the gates via the cumax activation, and supplying an auxiliary set of gates controlling fine-grained behavior. Both of these are important novelties and together allow recurrent models to better capture tree structures.
However, the UGI and refine gate can be viewed as improvements over each of these, respectively, demonstrated both theoretically (below) and empirically (Sections 3 and E.3), even on tasks involving hierarchical sequences.
Ordered gates Despite having the same parameter count and asymptotic efficiency as standard sigmoid gates, cumax gates seem noticeably slower and less stable in practice for large hidden sizes. Additionally, using auxiliary master gates creates additional parameters compared to the basic LSTM. Shen et al. (2018) alleviated both of these problems by defining a downsize operation, whereby neurons are grouped in chunks of sizeC, each of which share the same master gate values. However, this also creates an additional hyperparameter.
The speed and stability issues can be fixed by just using the sigmoid non-linearity instead of cumax. To recover the most important properties of the cumax—activations at multiple timescales—the equivalent sigmoid gate can be initialized so as to match the distribution of cumax gates at initialization. This is just uniform gate initialization (equation (9)).
However, we believe that the cumax activation is still valuable in many situations if speed and instability are not issues. These include when the hidden size is small, when extremal gate activations are desired, or when ordering needs to be strictly enforced to induce explicit hierarchical structure. For example, Section (3.1) shows that they can solve hard memory tasks by themselves.
Master gates We observe that the magnitudes of master gates are suboptimally normalized. A nice interpretation of gated recurrent models shows that they are a discretization of a continuous differential equation. This leads to the leaky RNN model ht+1=(1−α)ht+αut, where ut is the update to the model such as tanh(Wxxt+Whht+b). Learning α as a function of the current time step leads to the
simplest gated recurrent model4
ft=σ(Lf (xt,ht−1)) ut=tanh(Lu(xt,ht−1)) ht=ftht−1+(1−ft)ut.
Tallec & Ollivier (2018) show that this exactly corresponds to the discretization of a differential equation that is invariant to time warpings and time rescalings. In the context of the LSTM, this interpretation requires the values of the forget and input gates to be tied so that ft+ it = 1. This weight-tying is often enforced, for example in the most popular LSTM variant, the GRU (Cho et al., 2014), or our URgates. In a large-scale LSTM architecture search, it was found that removing the input gate was not significantly detrimental (Greff et al., 2016).
However, the ON-LSTM does not satisfy this conventional wisdom that the input and forget gates should sum to close to 1. Proposition 3. At initialization, the expected value of the average effective forget gate activation f̂t is 5/6.
Let us consider the sum of the effective forget and input gates at initialization. Adding equations (22) and (23) yields
f̂t+ ît=(ft+it)◦ωt+(f̃t+ ĩt−2ωt) = f̃t+ ĩt+(ft+it−2)◦ωt.
Note that the master gates (19), (20) sum 1 in expectation at initialization, as do the original forget and input gates. Looking at individual units in the ordered master gates, we have Ef̂ (j)= jn ,Eî
(j)=1− jn . Thus the above simplifies to
E[f̂t+ ît]=1−Eωt E[f̂ (j)t + î (j) t ]=1− j
n (1− j n )
E [ Ej∈[n]f̂ (j) t + î (j) t ] ≈1− ∫ 1 0 xdx+ ∫ 1 0 x2dx
= 5
6 .
The gate normalization can be fixed by re-scaling equations (22) and (23). It turns out that tying the master gates and re-scaling is exactly equivalent to the mechanism of a refine gate. In this equivalence, the role of the master and forget gates of the ON-LSTM are played by our forget and refine gate respectively.
Proposition 4. Suppose the master gates f̃t ,̃it are tied and the equations (22)-(23) defining the effective gates f̂t ,̂it are rescaled such as to ensure E[f̂t+ ît]=1 at initialization. The resulting gate mechanism is exactly equivalent to that of the refine gate.
Consider the following set of equations where the master gates are tied (f̃t+ ĩt=1,ft+it=1) and (22)-(23) are modified with an extra coefficient (rescaling in bold):
ĩt=1−f̃t (24) ωt= f̃t · ĩt (25) f̂t=2·ft ·ωt+(f̃t−ωt) (26) ît=2·it ·ωt+(̃it−ωt) (27)
Now we have f̂t+ ît= f̃t+ ĩt+2(ft+it−1)·ωt
=1+2(ft+it−1)·ωt 4In the literature, this is called the JANET (van der Westhuizen & Lasenby, 2018), which is also equivalent to the GRU without a reset gate (Chung et al., 2014), or a recurrent highway network with depth L=1 (Zilly et al., 2017).
which has the correct scaling, i.e. E[f̂t + ît] = 1 at initialization assuming that E[ft + it] = 1 at initialization.
But (26) can be rewritten as follows:
f̂=2·f ·ω+(f̃−ω) =2·f ·f̃ ·(1−f̃)+(f̃−f̃ ·(1−f̃)) =2f ·f̃−2f ·f̃2+f̃2
=f ·2f̃−f ·f̃2−f ·f̃2+f̃2
=f ·(1−(1−f̃))2+(1−f)·f̃2.
This is equivalent to the refine gate, where the master gate plays the role of the forget gate and the forget gate plays the role of the refine gate. It can be shown that in this case, the effective input gate ît (27) is also defined through a refine gate mechanism, where ĩt=1−f̃t is refined by it:
î= i·(1−(1− ĩ))2+(1−i)· ĩ2.
Based on our experimental findings, in general we would recommend the refine gate in place of the master gate.

heading: B.5 GATE ABLATION DETAILS
text: For clarity, we formally define the gate ablations considered which mix and match different gate components.
We remark that other combinations are possible, for example combining CI with either auxiliary gate type, which would lead to CR- or CM- gates. Alternatively, the master or refine gates could be defined using different activation and initialization strategies. We chose not to consider these methods due to lack of interpretation and theoretical soundness.
O- This ablation uses the cumax activation to order the forget/input gates and has no auxiliary gates. ft=cumax(Lf (xt,ht−1)) (28) it=1−cumax(Li(xt,ht−1)). (29)
We note that one difficulty with this in practice is the reliance on the expensive cumax, and hypothesize that this is perhaps the ON-LSTM’s original motivation for the second set of gates combined with downsizing.
UM- This variant of the ON-LSTM ablates the cumax operation on the master gates, replacing it with a sigmoid activation initialized with UGI. Equations (19), (20) are replaced with
u=U(0,1) (30) bf =σ −1(u) (31)
f̃t=σ(Lf̃ (xt,ht−1)+bf ) (32) ĩt=σ(Lĩ(xt,ht−1)−bf ) (33)
Equations (21)-(23) are then used to define effective gates f̂t ,̂it which are used in the gated update (1) or (5).
OR- This ablation combines ordered main gates with an auxilliary refine gate.
f̃t=cumax(Lf̃ (xt,ht−1)+bf ) (34) rt=σ(Lr(xt, ht−1)+br) (35) gt=rt ·(1−(1−ft)2)+(1−rt)·f2t (36) it=1−gt (37)
gt,it are used as the effective forget and input gates.

heading: C ANALYSIS DETAILS
text: The gradient analysis in Figure 2 was constructed as follows. Let f,r,g be the forget, refine, and effective gates g=2rf+(1−2r)f2. Then
∇xg=2rf(1−f)+(1−2r)(2f)(f(1−f))=2f(1−f)[r+(1−2r)f ] ∇yg=2fr(1−r)+(−2f2)r(1−r)=2fr(1−r)(1−f) ‖∇g‖2=[2f(1−f)]2 [ (r+f−2fr)2+r2(1−r)2 ] .
Substituting the relation
r= g−f2
2f(1−f) ,
this reduces to ‖∇g‖2=[2f(1−f)]2 [(
g−f2+2f2(1−f)−2f(g−f2) 2f(1−f)
)2 + (g−f2)2(2f(1−f)−g+f2)2
(2f(1−f))4
]
=((g−f2)(1−2f)+2f2(1−f))2+(g−f2)2 ( 1− g−f 2
2f(1−f)
)2 .
Given the constraint f2≤g≤1−(1−f)2, this function can be minimized and maximized in terms of g to produce the upper and lower bounds in Figure 2d. This was performed numerically.

heading: D EXPERIMENTAL DETAILS
text: To normalize the number of parameters used for models using master gates, i.e. the OM- and UMgating mechanisms, we used a downsize factor on the main gates (see Section B.4). This was set to C =16 for the synthetic and image classification tasks, and C =32 for the language modeling and program execution tasks which used larger hidden sizes.

heading: D.1 SYNTHETIC TASKS
text: All models consisted of single layer LSTMs with 256 hidden units, trained with the Adam optimizer Kingma & Ba (2014) with learning rate 1e-3. Gradients were clipped at 1.0.
The training data consisted of randomly generated sequences for every minibatch rather than iterating through a fixed dataset. Each method ran 3 seeds, with the same training data for every method.
Our version of the Copy task is a very minor variant of other versions reported in the literature, with the main difference being that the loss is considered only over the last 10 output tokens which need to be memorized. This normalizes the loss so that losses approaching 0 indicate true progress. In contrast, this task is usually defined with the model being required to output a dummy token at the firstN+10 steps, meaning it can be hard to evaluate performance since low average losses simply indicate that the model learns to output the dummy token.
For Figure 3, the log loss curves show the median of 3 seeds, and the error bars indicate 60% confidence.
For Figure 4, each histogram represents the distribution of forget gate values of the hidden units (of which there are 256). The values are created by averaging units over time and samples, i.e., reducing a minibatch of forget gate activations of shape (batch size, sequence length, hidden size) over the first two diensions, to produce the average activation value for every unit.
D.2 IMAGE CLASSIFICATION
All models used a single hidden layer recurrent network (LSTM or GRU). Inputs x to the model were given in batches as a sequence of shape (sequence length, num channels), (e.g. (1024,3)
for CIFAR-10), by flattening the input image left-to-right, top-to-bottom. The outputs of the model of shape (sequence length, hidden size)were processed independently with a single ReLU hidden layer of size 256 before the final fully-connected layer outputting softmax logits. All training was performed with the Adam optimizer, batch size 50, and gradients clipped at 1.0. MNIST trained for 150 epochs, CIFAR-10 used 100 epochs over the training set.
Table 2 All models (LSTM and GRU) used hidden state size 512. Learning rate swept in {2e− 4,5e−4,1e−3,2e−3}with three seeds each.
Table 2 reports the highest validation score found. The GRU model swept over learning rates {2e− 4,5e−4}; all methods were unstable at higher learning rates. Figure 5 shows the median validation accuracy with quartiles (25/75% confidence intervals) over the seeds, for the best-performing stable learning rate (i.e. the one with highest average validation score on the final epoch).
Table 3 The UR-LSTM and UR-GRU used 1024 hidden units for the sequential and permuted MNIST task, and 2048 hidden units for the sequential CIFAR task. The vanilla LSTM baseline used 512 hidden units for MNIST and 1024 for CIFAR. Larger hidden sizes were found to be unstable.
Zoneout parameters were fixed to reasonable default settings based on Krueger et al. (2016), which are zc=0.5,zh=0.05 for LSTM and z=0.1 for GRU. When zoneout was used, standard Dropout (Srivastava et al., 2014) with probability 0.5 was also applied to the output classification hidden layer.

heading: D.3 LANGUAGE MODELING
text: Hyperparameters are taken from Rae et al. (2018) tuned for the vanilla LSTM, which consist of (chosen parameter bolded out of sweep): {1,2} LSTM layer, {0.0,0.1,0.2,0.3} embedding dropout, {yes,no} layer norm, and {shared,not shared} input/output embedding parameters. Our only divergence is using a hidden size of 3072 instead of 2048, which we found improved the performance of the vanilla LSTM. Training was performed with Adam at learning rate 1e-3, gradients clipped to 0.1, sequence length 128, and batch size 128 on TPU. The LSTM state was reset between article boundaries. Figure 6b shows smoothed validation perplexity curves showing the 95% confidence intervals over the last 1% of data.
Reinforcement Learning The Active Match and Passive Match tasks were borrowed from Hung et al. (2018) with the same settings. For Figures 7 and 10, the discount factor in the environment was set to γ= .96. For Figure 8, the discount factor was γ= .998. Figure 10 corresponds to the full Active Match task in Hung et al. (2018), while Figure 8 is their version with small distractor rewards where the apples in the distractor phase give 1 instead of 5 reward.

heading: D.4 PROGRAM EVALUATION
text: Protocol was taken from Santoro et al. (2018) with minor changes to the hyperparameter search. All models were trained with the Adam optimizer, the Mix curriculum strategy from Zaremba & Sutskever (2014), and batch size 128.
RMC: The RMC models used a fixed memory slot size of 512 and swept over {2,4}memories and {2,4} attention heads for a total memory size of 1024 or 2048. They were trained for 2e5 iterations. LSTM: Instead of two-layer LSTMs with sweeps over skip connections and output concatenation, single-layer LSTMs of size 1024 or 2048 were used. Learning rate was swept in {5e-4, 1e-3}, and models were trained for 5e5 iterations. Note that training was still faster than the RMC models despite the greater number of iterations.

heading: D.5 ADDITIONAL DETAILS
text: Implementation Details The inverse sigmoid function (9) can be unstable if the input is too close to {0,1}. Uniform gate initialization was instead implemented by sampling from the distribution U [1/d,1−1/d] instead of U [0,1], where d is the hidden size, to avoid any potential numerical edge cases. This choice is justified by the fact that with perfect uniform sampling, the expected smallest and largest samples would be 1/(d+1) and 1−1/(d+1).
For distributional initialization strategies, a trainable bias vector was sampled independently from the chosen distribution (i.e. equation (17) or (9)) and added/subtracted to the forget and input gate ((2)-(3)) before the non-linearity. Additionally, each linear model such asWxfxt+Whfht−1 had its own trainable bias vector, effectively doubling the learning rate on the pre-activation bias terms on the forget and input gates. This was an artifact of implementation and not intended to affect performance.
The refine gate update equation (12) can instead be implemented as
gt=rt ·(1−(1−ft)2)+(1−rt)·f2t =2rt ·ft+(1−2rt)·f2t
Permuted image classification In an effort to standardize the permutation used in the Permuted MNIST benchmark, we use a particular deterministic permutation rather than a random one. After flattening the input image into a one-dimensional sequence, we apply the bit reversal permutation. This permutation sends the index i to the index j such that j’s binary representation is the reverse of i’s binary representation. The intuition is that if two indices i,i′ are close, they must differ in their lower-order bits. Then the bit-reversed indices will be far apart. Therefore the bit-reversal permutation destroys spatial and temporal locality, which is desirable for these sequence classification tasks meant to test long-range dependencies rather than local structure.

heading: E ADDITIONAL EXPERIMENTS

heading: E.1 SYNTHETIC FORGETTING
text: This raises the question: what happens if the initialization distribution does not match the task at hand; could the gates learn back to a more moderate regime? We point out that such a phenomenon could occur non-pathologically on more complex setups, such as a scenario where a model trains to remember on a Copy-like task and then needs to “unlearn” as part of a meta-learning or continual learning setup.
Here, we consider such a synthetic scenario and experimentally show that the addition of a refine gate helps models train much faster while in a saturated regime with extremal activations. We also point to the poor performance of C- outside of synthetic memory tasks when using our high hyperparameter-free initialization as more evidence that it is very difficult for standard gates to unlearn undesired saturated behavior.
For this experiment, we initialize the biases of the gates extremely high (effective forget activation ≈σ(6). We then consider the Adding task (Section 3.1 of length 500, hidden size 64, learning rate 1e-4. The R-LSTM is able to solve the task, while the LSTM is stuck after 1e4 iterations.

heading: E.2 REINFORCEMENT LEARNING
text: Figures 7 and 8 evaluated our gating methods with the LSTM and RMA models on the Passive Match and Active Match tasks, with and without distractors. We additionally ran the agents on an even harder
version of the Active Match task with larger distractor rewards (the full Active Match from Hung et al. (2018)). Learning curves are shown in Figure 10. Similarly to the other results, the UR- gated core is noticeably better than the others. For the DNC model, it is the only one that performs better than random chance.

heading: E.3 PROGRAM EXECUTION
text: The Learning to Execute (Zaremba & Sutskever, 2014) dataset consists of algorithmic snippets from a programming language of pseudo-code. An input is a program from this language presented one character at a time, and the target output is a numeric sequence of characters representing the execution output of the program. There are three categories of tasks: Addition, Control, and Program, with distinctive types of input programs. We use the most difficult setting from Zaremba & Sutskever (2014), which uses the parameters nesting=4, length=9, referring to the nesting depth of control structure and base length of numeric literals, respectively. Examples of input programs are shown in previous works (Zaremba & Sutskever, 2014; Santoro et al., 2018).
We are interested in this task for several reasons. First, we are interested in comparing against the Cand OM- gate methods, because
• The maximum sequence length is fairly long (several hundred tokens), meaning our Tmax heuristic for C- gates is within the right order of magnitude of dependency lengths.
• The task has highly variable sequence lengths, wherein the standard training procedure randomly samples inputs of varying lengths (called the ”Mix” curriculum in Zaremba & Sutskever (2014)). Additionally, the Control and Program tasks contain complex control flow and nested structure. They are thus a measure of a sequence model’s ability to model dependencies of differing lengths, as well as hierarchical information. Thus we are interested in comparing the effects of UGI methods, as well as the full OM- gates which are designed for hierarchical structures (Shen et al., 2018).
Finally, this task has prior work using a different type of recurrent core, the Relational Memory Core (RMC), that we also use as a baseline to evaluate our gates on different models Santoro et al. (2018). Both the LSTM and RMC were found to outperform other recurrent baselines such as the Differential Neural Computer (DNC) and EntNet.
Training curves are shown in Figure 11, which plots the median accuracy with confidence intervals. We point out a few observations. First, despite having a Tmax value on the right order of magnitude, the Cgated methods have very poor performance across the board, reaffirming the chrono initialization’s high sensitivity to this hyperparameter.
Second, the U-LSTM and U-RMC are the best methods on the Addition task. Additionally, the UR-RMC vs. RMC on Addition is one of the very few tasks we have found where a generic substitution of the UR- gate does not improve on the basic gate. We have not investigated what property of this task caused these phenomena.
Aside from the U-LSTM on addition, the UR-LSTM outperforms all other LSTM cores. The UR-RMC is also the best core on both Control and Program, the tasks involving hierarchical inputs and longer dependencies. For the most part, the improved mechanisms of the UR- gates seem to transfer to this recurrent core as well. We highlight that this is not true of similar gating mechanisms. In particular, the OM-LSTM, which is supposed to model hierarchies, has good performance on Control and Program as expected (although not better than the UR-LSTM). However, the OM- gates’ performance plummets when transferred to the RMC core.
Interestingly, the -LSTM cores are consistently better than the -RMC versions, contrary to previous findings on easier versions of this task using similar protocol and hyperparameters Santoro et al. (2018). We did not explore different hyperparameter regimes on this more difficult setting.

Review summaries:
This paper introduces two novel techniques to help long term signal propagation in RNNs. One is an initialization strategy which uses inverse sigmoid function to avoid the decay of the contribution of the input earlier in time and another is a new design of a refine gate which pushes the value of the gate closer to 0 or 1.  The authors conduct exhaustive ablation and empirical studies on copy task, sequential MNIST, language modeling and reinforcement learning. 

This paper proposes to improve the learnability of the gating mechanism in RNN by two modifications on the standard RNN structure, uniform gate initialization and refine gate. The authors give some propositions to show that the refine gate can maintain an effective forget effect within a larger range of timescale. The authors conduct experiments on four different tasks and compare the proposed modification with baseline methods.