Title: SELF-TUNING NETWORKS: BILEVEL OPTIMIZATION OF HYPERPARAMETERS US-ING STRUCTURED BEST-RESPONSE FUNCTIONS

Abstract: Hyperparameter optimization can be formulated as a bilevel optimization problem, where the optimal parameters on the training set depend on the hyperparameters. We aim to adapt regularization hyperparameters for neural networks by fitting compact approximations to the best-response function, which maps hyperparameters to optimal weights and biases. We show how to construct scalable best-response approximations for neural networks by modeling the best-response as a single network whose hidden units are gated conditionally on the regularizer. We justify this approximation by showing the exact best-response for a shallow linear network with L 2 -regularized Jacobian can be represented by a similar gating mechanism. We fit this model using a gradient-based hyperparameter optimization algorithm which alternates between approximating the best-response around the current hyperparameters and optimizing the hyperparameters using the approximate best-response function. Unlike other gradient-based approaches, we do not require differentiating the training loss with respect to the hyperparameters, allowing us to tune discrete hyperparameters, data augmentation hyperparameters, and dropout probabilities. Because the hyperparameters are adapted online, our approach discovers hyperparameter schedules that can outperform fixed hyperparameter values. Empirically, our approach outperforms competing hyperparameter optimization methods on large-scale deep learning problems. We call our networks, which update their own hyperparameters online during training, Self-Tuning Networks (STNs). * Equal contribution. 1 The uniqueness of the arg min is assumed.

1. INTRODUCTION:

Regularization hyperparameters such as weight decay, data augmentation, and dropout (Srivastava et al., 2014) are crucial to the generalization of neural networks, but are difficult to tune. Popular approaches to hyperparameter optimization include grid search, random search (Bergstra & Bengio, 2012) , and Bayesian optimization (Snoek et al., 2012) . These approaches work well with low-dimensional hyperparameter spaces and ample computational resources; however, they pose hyperparameter optimization as a black-box optimization problem, ignoring structure which can be exploited for faster convergence, and require many training runs.

We can formulate hyperparameter optimization as a bilevel optimization problem. Let w denote parameters (e.g. weights and biases) and λ denote hyperparameters (e.g. dropout probability). Let L T and L V be functions mapping parameters and hyperparameters to training and validation losses, respectively. We aim to solve 1 :

EQUATION

Substituting the best-response function w * (λ) = arg min w L T (λ, w) gives a single-level problem:

EQUATION

If the best-response w * is known, the validation loss can be minimized directly by gradient descent using Equation 2, offering dramatic speed-ups over black-box methods. However, as the solution to a high-dimensional optimization problem, it is difficult to compute w * even approximately.

Following Lorraine & Duvenaud (2018) , we propose to approximate the best-response w * directly with a parametric function ŵφ . We jointly optimize φ and λ, first updating φ so that ŵφ ≈ w * in a neighborhood around the current hyperparameters, then updating λ by using ŵφ as a proxy for w * in Eq. 2:

EQUATION

Finding a scalable approximation ŵφ when w represents the weights of a neural network is a significant challenge, as even simple implementations entail significant memory overhead. We show how to construct a compact approximation by modelling the best-response of each row in a layer's weight matrix/bias as a rank-one affine transformation of the hyperparameters. We show that this can be interpreted as computing the activations of a base network in the usual fashion, plus a correction term dependent on the hyperparameters. We justify this approximation by showing the exact best-response for a shallow linear network with L 2 -regularized Jacobian follows a similar structure. We call our proposed networks Self-Tuning Networks (STNs) since they update their own hyperparameters online during training.

STNs enjoy many advantages over other hyperparameter optimization methods. First, they are easy to implement by replacing existing modules in deep learning libraries with "hyper" counterparts which accept an additional vector of hyperparameters as input2 . Second, because the hyperparameters are adapted online, we ensure that computational effort expended to fit φ around previous hyperparameters is not wasted. In addition, this online adaption yields hyperparameter schedules which we find empirically to outperform fixed hyperparameter settings. Finally, the STN training algorithm does not require differentiating the training loss with respect to the hyperparameters, unlike other gradient-based approaches (Maclaurin et al., 2015; Larsen et al., 1996) , allowing us to tune discrete hyperparameters, such as the number of holes to cut out of an image (DeVries & Taylor, 2017) , data-augmentation hyperparameters, and discrete-noise dropout parameters. Empirically, we evaluate the performance of STNs on large-scale deep-learning problems with the Penn Treebank (Marcus et al., 1993) and CIFAR-10 datasets (Krizhevsky & Hinton, 2009) , and find that they substantially outperform baseline methods.

2. BILEVEL OPTIMIZATION:

A bilevel optimization problem consists of two sub-problems called the upper-level and lower-level problems, where the upper-level problem must be solved subject to optimality of the lower-level problem. Minimax problems are an example of bilevel programs where the upper-level objective equals the negative lower-level objective. Bilevel programs were first studied in economics to model leader/follower firm dynamics (Von Stackelberg, 2010) and have since found uses in various fields (see Colson et al. (2007) for an overview). In machine learning, many problems can be formulated as bilevel programs, including hyperparameter optimization, GAN training (Goodfellow et al., 2014) , meta-learning, and neural architecture search (Zoph & Le, 2016) .

Even if all objectives and constraints are linear, bilevel problems are strongly NP-hard (Hansen et al., 1992; Vicente et al., 1994) . Due to the difficulty of obtaining exact solutions, most work has focused on restricted settings, considering linear, quadratic, and convex functions. In contrast, we focus on obtaining local solutions in the nonconvex, differentiable, and unconstrained setting. Let F, f : R n × R m → R denote the upper-and lower-level objectives (e.g., L V and L T ) and λ ∈ R n , w ∈ R m denote the upper-and lower-level parameters. We aim to solve:

EQUATION

subject to w ∈ arg min

EQUATION

It is desirable to design a gradient-based algorithm for solving Problem 4, since using gradient information provides drastic speed-ups over black-box optimization methods (Nesterov, 2013) . The simplest method is simultaneous gradient descent, which updates λ using ∂F /∂λ and w using ∂f /∂w. However, simultaneous gradient descent often gives incorrect solutions as it fails to account for the dependence of w on λ. Consider the relatively common situation where F doesn't depend directly on λ , so that ∂F /∂λ ≡ 0 and hence λ is never updated.

2.1. GRADIENT DESCENT VIA THE BEST-RESPONSE FUNCTION:

A more principled approach to solving Problem 4 is to use the best-response function (Gibbons, 1992) . Assume the lower-level Problem 4b has a unique optimum w * (λ) for each λ. Substituting the best-response function w * converts Problem 4 into a single-level problem:

EQUATION

If w * is differentiable, we can minimize Eq. 5 using gradient descent on F * with respect to λ. This method requires a unique optimum w * (λ) for Problem 4b for each λ and differentiability of w * . In general, these conditions are difficult to verify. We give sufficient conditions for them to hold in a neighborhood of a point (λ 0 , w 0 ) where w 0 solves Problem 4b given λ 0 . Lemma 1. (Fiacco & Ishizuka, 1990 ) Let w 0 solve Problem 4b for λ 0 . Suppose f is C 2 in a neighborhood of (λ 0 , w 0 ) and the Hessian ∂ 2 f /∂w 2 (λ 0 , w 0 ) is positive definite. Then for some neighborhood U of λ 0 , there exists a continuously differentiable function w * : U → R m such that w * (λ) is the unique solution to Problem 4b for each λ ∈ U and w * (λ 0 ) = w 0 .

Proof. See Appendix B.1.

The gradient of F * decomposes into two terms, which we term the direct gradient and the response gradient. The direct gradient captures the direct reliance of the upper-level objective on λ, while the response gradient captures how the lower-level parameter responds to changes in the upper-level parameter:

∂F * ∂λ (λ 0 ) = ∂F ∂λ (λ 0 , w * (λ 0 )) Direct gradient + ∂F ∂w (λ 0 , w * (λ 0 )) ∂w * ∂λ (λ 0 ) Response gradient (6)

Even if ∂F /∂λ ≡ 0 and simultaneous gradient descent is possible, including the response gradient can stabilize optimization by converting the bilevel problem into a single-level one, as noted by Metz et al. (2016) for GAN optimization. Conversion to a single-level problem ensures that the gradient vector field is conservative, avoiding pathological issues described by Mescheder et al. (2017) .

2.2. APPROXIMATING THE BEST-RESPONSE FUNCTION:

In general, the solution to Problem 4b is a set, but assuming uniqueness of a solution and differentiability of w * can yield fruitful algorithms in practice. In fact, gradient-based hyperparameter optimization methods can often be interpreted as approximating either the best-response w * or its Jacobian ∂w * /∂λ, as detailed in Section 5. However, these approaches can be computationally expensive and often struggle with discrete hyperparameters and stochastic hyperparameters like dropout probabilities, since they require differentiating the training loss with respect to the hyperparameters. Promising approaches to approximate w * directly were proposed by Lorraine & Duvenaud (2018) , and are detailed below.

1. Global Approximation. The first algorithm proposed by Lorraine & Duvenaud (2018) approximates w * as a differentiable function ŵφ with parameters φ. If w represents neural net weights, then the mapping ŵφ is a hypernetwork (Schmidhuber, 1992; Ha et al., 2016) . If the distribution p(λ) is fixed, then gradient descent with respect to φ minimizes:

EQUATION

If support(p) is broad and ŵφ is sufficiently flexible, then ŵφ can be used as a proxy for w * in Problem 5, resulting in the following objective:

EQUATION

2. Local Approximation. In practice, ŵφ is usually insufficiently flexible to model w * on support(p). The second algorithm of Lorraine & Duvenaud (2018) locally approximates w * in a neighborhood around the current upper-level parameter λ. They set p( |σ) to a factorized Gaussian noise distribution with a fixed scale parameter σ ∈ R n + , and found φ by minimizing the objective:

EQUATION

Intuitively, the upper-level parameter λ is perturbed by a small amount, so the lower-level parameter learns how to respond. An alternating gradient descent scheme is used, where φ is updated to minimize equation 9 and λ is updated to minimize equation 8. This approach worked for problems using L 2 regularization on MNIST (LeCun et al., 1998) . However, it is unclear if the approach works with different regularizers or scales to larger problems. It requires ŵφ , which is a priori unwieldy for high dimensional w. It is also unclear how to set σ, which defines the size of the neighborhood on which φ is trained, or if the approach can be adapted to discrete and stochastic hyperparameters.

3. SELF-TUNING NETWORKS:

In this section, we first construct a best-response approximation ŵφ that is memory efficient and scales to large neural networks. We justify this approximation through analysis of simpler situations. Then, we describe a method to automatically adjust the scale of the neighborhood φ is trained on. Finally, we formally describe our algorithm and discuss how it easily handles discrete and stochastic hyperparameters. We call the resulting networks, which update their own hyperparameters online during training, Self-Tuning Networks (STNs).

3.1. AN EFFICIENT BEST-RESPONSE APPROXIMATION FOR NEURAL NETWORKS:

We propose to approximate the best-response for a given layer's weight matrix W ∈ R Dout×Din and bias b ∈ R Dout as an affine transformation of the hyperparameters λ3 :

Ŵφ (λ) = W elem + (V λ) row W hyper , bφ (λ) = b elem + (Cλ) b hyper (10)

Here, indicates elementwise multiplication and row indicates row-wise rescaling. This architecture computes the usual elementary weight/bias, plus an additional weight/bias which has been scaled by a linear transformation of the hyperparameters. Alternatively, it can be interpreted as directly operating on the pre-activations of the layer, adding a correction to the usual pre-activation to account for the hyperparameters:

EQUATION

This best-response architecture is tractable to compute and memory-efficient: it requires D out (2D in + n) parameters to represent Ŵφ and D out (2 + n) parameters to represent bφ , where n is the number of hyperparameters. Furthermore, it enables parallelism: since the predictions can be computed by transforming the pre-activations (Equation 11), the hyperparameters for different examples in a batch can be perturbed independently, improving sample efficiency. In practice, the approximation can be implemented by simply replacing existing modules in deep learning libraries with "hyper" counterparts which accept an additional vector of hyperparameters as input4 .

3.2. EXACT BEST-RESPONSE FOR TWO-LAYER LINEAR NETWORKS:

Given that the best-response function is a mapping from R n to the high-dimensional weight space R m , why should we expect to be able to represent it compactly? And why in particular would equation 10 be a reasonable approximation? In this section, we exhibit a model whose best-response function can be represented exactly using a minor variant of equation 10: a linear network with Jacobian norm regularization. In particular, the best-response takes the form of a network whose hidden units are modulated conditionally on the hyperparameters.

Consider using a 2-layer linear network with weights w Suppose we use a squared-error loss regularized with an L 2 penalty on the Jacobian ∂y /∂x, where the penalty weight λ lies in R and is mapped using exp to lie R + :

= (Q, s) ∈ R D×D × R D to predict targets t ∈ R from inputs x ∈ R D :

L T (λ, w) = (x,t)∈D (y(x; w) -t) 2 + 1 |D| exp(λ) ∂y ∂x (x; w) 2 (13)

Theorem 2. Let w 0 = (Q 0 , s 0 ), where Q 0 is the change-of-basis matrix to the principal components of the data matrix and s 0 solves the unregularized version of Problem 13 given Q 0 . Then there exist v, c ∈ R D such that the best-response function 5 w * (λ) = (Q * (λ), s * (λ)) is:

Q * (λ) = σ(λv + c) row Q 0 , s * (λ) = s 0 ,

where σ is the sigmoid function.

Proof. See Appendix B.2.

Observe that y(x; w * (λ)) can be implemented as a regular network with weights w 0 = (Q 0 , s 0 ) with an additional sigmoidal gating of its hidden units a(x; w * (λ)):

EQUATION

This architecture is shown in Figure 1 . Inspired by this example, we use a similar gating of the hidden units to approximate the best-response for deep, nonlinear networks.

3.3. LINEAR BEST-RESPONSE APPROXIMATIONS:

The sigmoidal gating architecture of the preceding section can be further simplified if one only needs to approximate the best-response function for a small range of hyperparameter values. In particular, for a narrow enough hyperparameter distribution, a smooth best-response function can be approximated by an affine function (i.e. its first-order Taylor approximation). Hence, we replace the sigmoidal gating with linear gating, in order that the weights be affine in the hyperparameters. The following theorem shows that, for quadratic lower-level objectives, using an affine approximation to the best-response function and minimizing E ∼p( |σ) [f (λ + , ŵφ (λ + ))] yields the correct best-response Jacobian, thus ensuring gradient descent on the approximate objective F (λ, ŵφ (λ)) converges to a local optimum: The effect of the sampled neighborhood. Left: If the sampled neighborhood is too small (e.g., a point mass) the approximation learned will only match the exact best-response at the current hyperparameter, with no guarantee that its gradient matches that of the best-response. Middle: If the sampled neighborhood is not too small or too wide, the gradient of the approximation will match that of the best-response. Right: If the sampled neighborhood is too wide, the approximation will be insufficiently flexible to model the best-response, and again the gradients will not match.

EQUATION

3.4. ADAPTING THE HYPERPARAMETER DISTRIBUTION:

The entries of σ control the scale of the hyperparameter distribution on which φ is trained. If the entries are too large, then ŵφ will not be flexible enough to capture the best-response over the samples. However, the entries must remain large enough to force ŵφ to capture the shape locally around the current hyperparameter values. We illustrate this in Figure 2 . As the smoothness of the loss landscape changes during training, it may be beneficial to vary σ.

To address these issues, we propose adjusting σ during training based on the sensitivity of the upperlevel objective to the sampled hyperparameters. We include an entropy term weighted by τ ∈ R + which acts to enlarge the entries of σ. The resulting objective is:

EQUATION

This is similar to a variational inference objective, where the first term is analogous to the negative log-likelihood, but τ = 1. As τ ranges from 0 to 1, our objective interpolates between variational optimization (Staines & Barber, 2012) and variational inference, as noted by Khan et al. (2018) . Similar objectives have been used in the variational inference literature for better training (Blundell et al., 2015) and representation learning (Higgins et al., 2017) .

Minimizing the first term on its own eventually moves all probability mass towards an optimum λ * , resulting in σ = 0 if λ * is an isolated local minimum. This compels σ to balance between shrinking to decrease the first term while remaining sufficiently large to avoid a heavy entropy penalty. When benchmarking our algorithm's performance, we evaluate F (λ, ŵφ (λ)) at the deterministic current hyperparameter λ 0 . (This is a common practice when using stochastic operations during training, such as batch normalization or dropout.)

3.5. TRAINING ALGORITHM:

We now describe the complete STN training algorithm and discuss how it can tune hyperparameters that other gradient-based algorithms cannot, such as discrete or stochastic hyperparameters. We use an unconstrained parametrization λ ∈ R n of the hyperparameters. Let r denote the element-wise function which maps λ to the appropriate constrained space, which will involve a non-differentiable discretization for discrete hyperparameters.

Let L T and L V denote training and validation losses which are (possibly stochastic, e.g., if using dropout) functions of the hyperparameters and parameters. Define functions f, F by f (λ, w) = L T (r(λ), w) and F (λ, w) = L V (r(λ), w). STNs are trained by a gradient descent scheme which alternates between updating φ for T train steps to minimize E ∼p( |σ) [f (λ + , ŵφ (λ + ))] (Eq. 9) and updating λ and σ for T valid steps to minimize . 15) . We give our complete algorithm as Algorithm 1 and show how it can be implemented in code in Appendix G. The possible non-differentiability of r due to discrete hyperparameters poses no problem. To estimate the derivative of E ∼p( |σ) [f (λ + , ŵφ (λ + ))] with respect to φ, we can use the reparametrization trick and compute ∂f /∂w and ∂ ŵφ/∂φ, neither of whose computation paths involve the discretization r.

E ∼p( |σ) [F (λ + , ŵφ (λ + ))] -τ H[p( |σ)] (Eq

To differentiate E ∼p( |σ) [F (λ + , ŵφ (λ + ))] -τ H[p( |σ)]

with respect to a discrete hyperparameter λ i , there are two cases we must consider:

Algorithm 1 STN Training Algorithm:

Initialize: Best-response approximation parameters φ, hyperparameters λ, learning rates

{α i } 3 i=1 while not converged do for t = 1, . . . , T train do ∼ p( |σ) φ ← φ -α 1 ∂ ∂φ f (λ + , ŵφ (λ + )) for t = 1, . . . , T valid do ∼ p( |σ) λ ← λ-α 2 ∂ ∂λ (F (λ + , ŵφ (λ + )) -τ H[p( |σ)]) σ ← σ-α 3 ∂ ∂σ (F (λ + , ŵφ (λ + )) -τ H[p( |σ)])

Case 1: For most regularization schemes, L V and hence F does not depend on λ i directly and thus the only gradient is through ŵφ . Thus, the reparametrization gradient can be used.

Case 2: If L V relies explicitly on λ i , then we can use the REINFORCE gradient estimator (Williams, 1992) to estimate the derivative of the expectation with respect to λ i . The number of hidden units in a layer is an example of a hyperparameter that requires this approach since it directly affects the validation loss. We do not show this in Algorithm 1, since we do not tune any hyperparameters which fall into this case.

4. EXPERIMENTS:

We applied our method to convolutional networks and LSTMs (Hochreiter & Schmidhuber, 1997) , yielding self-tuning CNNs (ST-CNNs) and self-tuning LSTMs (ST-LSTMs). We first investigated the behavior of STNs in a simple setting where we tuned a single hyperparameter, and found that STNs discovered hyperparameter schedules that outperformed fixed hyperparameter values. Next, we compared the performance of STNs to commonly-used hyperparameter optimization methods on the CIFAR-10 ( Krizhevsky & Hinton, 2009) and PTB (Marcus et al., 1993) datasets.

4.1. HYPERPARAMETER SCHEDULES:

Due to the joint optimization of the hypernetwork weights and hyperparameters, STNs do not use a single, fixed hyperparameter during training. Instead, STNs discover schedules for adapting the hyperparameters online, which can outperform any fixed hyperparameter. We examined this behavior in detail on the PTB corpus (Marcus et al., 1993) using an ST-LSTM to tune the output dropout rate applied to the hidden units.

The schedule discovered by an ST-LSTM for output dropout, shown in Figure 3 , outperforms the best, fixed output dropout rate (0.68) found by a fine-grained grid search, achieving 82.58 vs 85.83 validation perplexity. We claim that this is a consequence of the schedule, and not of regularizing effects from sampling hyperparameters or the limited capacity of ŵφ .

To rule out the possibility that the improved performance is due to stochasticity introduced by sampling hyperparameters during STN training, we trained a standard LSTM while perturbing its dropout rate around the best value found by grid search. We used (1) random Gaussian perturbations, and (2) sinusoid perturbations for a cyclic regularization schedule. STNs outperformed both perturbation methods (Table 1 ), showing that the improvement is not merely due to hyperparameter stochasticity. Details and plots of each perturbation method are provided in Appendix F. To determine whether the limited capacity of ŵφ acts as a regularizer, we trained a standard LSTM from scratch using the schedule for output dropout discovered by the ST-LSTM. Using this schedule, the standard LSTM performed nearly as well as the STN, providing evidence that the schedule itself (rather than some other aspect of the STN) was responsible for the improvement over a fixed dropout rate. To further demonstrate the importance of the hyperparameter schedule, we also trained a standard LSTM from scratch using the final dropout value found by the STN (0.78), and found that it did not perform as well as when following the schedule. The final validation and test perplexities of each variant are shown in Table 1 .

Next, we show in Figure 3 that the STN discovers the same schedule regardless of the initial hyperparameter values. Because hyperparameters adapt over a shorter timescale than the weights, we find that at any given point in training, the hyperparameter adaptation has already equilibrated. As shown empirically in Appendix F, low regularization is best early in training, while higher regularization is better later on. We found that the STN schedule implements a curriculum by using a low dropout rate early in training, aiding optimization, and then gradually increasing the dropout rate, leading to better generalization.

4.2. LANGUAGE MODELING:

We evaluated an ST-LSTM on the PTB corpus (Marcus et al., 1993) , which is widely used as a benchmark for RNN regularization due to its small size (Gal & Ghahramani, 2016; Merity et al., 2018; Wen et al., 2018) . We used a 2-layer LSTM with 650 hidden units per layer and 650dimensional word embeddings. We tuned 7 hyperparameters: variational dropout rates for the input, hidden state, and output; embedding dropout (that sets rows of the embedding matrix to 0); Drop-Connect (Wan et al., 2013) on the hidden-to-hidden weight matrix; and coefficients α and β that control the strength of activation regularization and temporal activation regularization, respectively.

For LSTM tuning, we obtained the best results when using a fixed perturbation scale of 1 for the hyperparameters. Additional details about the experimental setup and the role of these hyperparameters can be found in Appendix D. We compared STNs to grid search, random search, and Bayesian optimization. 6 Figure 4a shows the best validation perplexity achieved by each method over time. STNs outperform other meth- We evaluated ST-CNNs on the CIFAR-10 ( Krizhevsky & Hinton, 2009) dataset, where it is easy to overfit with high-capacity networks. We used the AlexNet architecture (Krizhevsky et al., 2012) , and tuned: (1) continuous hyperparameters controlling per-layer activation dropout, input dropout, and scaling noise applied to the input, (2) discrete data augmentation hyperparameters controlling the length and number of cut-out holes (DeVries & Taylor, 2017), and (3) continuous data augmentation hyperparameters controlling the amount of noise to apply to the hue, saturation, brightness, and contrast of an image. In total, we considered 15 hyperparameters.

We compared STNs to grid search, random search, and Bayesian optimization. Figure 5 shows the lowest validation loss achieved by each method over time, and Table 2 shows the final validation and test losses for each method. Details of the experimental setup are provided in Appendix E. Again, STNs find better hyperparameter configurations in less time than other methods. The hyperparameter schedules found by the STN are shown in Figure 6 .

5. RELATED WORK:

Bilevel Optimization. Colson et al. (2007) provide an overview of bilevel problems, and a comprehensive textbook was written by Bard (2013) . When the objectives/constraints are restricted to be linear, quadratic, or convex, a common approach replaces the lower-level problem with its KKT conditions added as constraints for the upper-level problem (Hansen et al., 1992; Vicente et al., 1994) .

In the unrestricted setting, our work loosely resembles trust-region methods (Colson et al., 2005) , which repeatedly approximate the problem locally using a simpler bilevel program. In closely related work, Sinha et al. (2013) used evolutionary techniques to estimate the best-response function iteratively.

Hypernetworks. First considered by Schmidhuber (1993; 1992) , hypernetworks are functions mapping to the weights of a neural net. Predicting weights in CNNs has been developed in various forms (Denil et al., 2013; Yang et al., 2015) . Ha et al. (2016) used hypernetworks to generate weights for modern CNNs and RNNs. Brock et al. (2017) used hypernetworks to globally approximate a bestresponse for architecture search. Because the architecture is not optimized during training, they require a large hypernetwork, unlike ours which locally approximates the best-response.

Gradient-Based Hyperparameter Optimization. There are two main approaches. The first approach approximates w * (λ 0 ) using w T (λ 0 , w 0 ), the value of w after T steps of gradient descent on f with respect to w starting at (λ 0 , w 0 ). The descent steps are differentiated through to approximate ∂w * /∂λ(λ 0 ) ≈ ∂w T /∂λ(λ 0 , w 0 ). This approach was proposed by Domke (2012) and used by Maclaurin et al. (2015) , Luketina et al. (2016) and Franceschi et al. (2018) . The second approach uses the Implicit Function Theorem to derive ∂w * /∂λ(λ 0 ) under certain conditions. This was first developed for hyperparameter optimization in neural networks (Larsen et al., 1996) and developed further by Pedregosa (2016) . Similar approaches have been used for hyperparameter optimization in log-linear models (Foo et al., 2008) , kernel selection (Chapelle et al., 2002; Seeger, 2007) , and image reconstruction (Kunisch & Pock, 2013; Calatroni et al., 2015) . Both approaches struggle with certain hyperparameters, since they differentiate gradient descent or the training loss with respect to the hyperparameters. In addition, differentiating gradient descent becomes prohibitively expensive as the number of descent steps increases, while implicitly deriving ∂w * /∂λ requires using Hessian-vector products with conjugate gradient solvers to avoid directly computing the Hessian.

Model-Based Hyperparameter Optimization. A common model-based approach is Bayesian optimization, which models p(r|λ, D), the conditional probability of the performance on some metric r given hyperparameters λ and a dataset D = {(λ i , r i )}. We can model p(r|λ, D) with various methods (Hutter et al., 2011; Bergstra et al., 2011; Snoek et al., 2012; 2015) . D is constructed iteratively, where the next λ to train on is chosen by maximizing an acquisition function C(λ; p(r|λ, D)) which balances exploration and exploitation. Training each model to completion can be avoided if assumptions are made on learning curve behavior (Swersky et al., 2014; Klein et al., 2017) . These approaches require building inductive biases into p(r|λ, D) which may not hold in practice, do not take advantage of the network structure when used for hyperparameter optimization, and do not scale well with the number of hyperparameters. However, these approaches have consistency guarantees in the limit, unlike ours.

Model-Free Hyperparameter Optimization. Model-free approaches include grid search and random search. Bergstra & Bengio (2012) advocated using random search over grid search. Successive Halving (Jamieson & Talwalkar, 2016) and Hyperband (Li et al., 2017) extend random search by adaptively allocating resources to promising configurations using multi-armed bandit techniques. These methods ignore structure in the problem, unlike ours which uses rich gradient information. However, it is trivial to parallelize model-free methods over computing resources and they tend to perform well in practice.

Hyperparameter Scheduling. Population Based Training (PBT) (Jaderberg et al., 2017) considers schedules for hyperparameters. In PBT, a population of networks is trained in parallel. The performance of each network is evaluated periodically, and the weights of under-performing networks are replaced by the weights of better-performing ones; the hyperparameters of the better network are also copied and randomly perturbed for training the new network clone. In this way, a single model can experience different hyperparameter settings over the course of training, implementing a schedule. STNs replace the population of networks by a single best-response approximation and use gradients to tune hyperparameters during a single training run.

6. CONCLUSION:

We introduced Self-Tuning Networks (STNs), which efficiently approximate the best-response of parameters to hyperparameters by scaling and shifting their hidden units. This allowed us to use gradient-based optimization to tune various regularization hyperparameters, including discrete hyperparameters. We showed that STNs discover hyperparameter schedules that can outperform fixed hyperparameters. We validated the approach on large-scale problems and showed that STNs achieve better generalization performance than competing approaches, in less time. We believe STNs offer a compelling path towards large-scale, automated hyperparameter tuning for neural networks. The (validation loss) direct (hyperparameter) gradient ∂w * (λ) ∂λ

A TABLE OF NOTATION:

The (elementary parameter) response gradient

∂L T (λ,w * (λ)) ∂w * (λ) ∂w * (λ) ∂λ

The (validation loss) response gradient dL T (λ,w) dλ

The hyperparameter gradient: a sum of the validation losses direct and response gradients

B PROOFS B.1 LEMMA 1:

Because w 0 solves Problem 4b given λ 0 , by the first-order optimality condition we must have:

∂f ∂w (λ 0 , w 0 ) = 0 (16)

The Jacobian of ∂f /∂w decomposes as a block matrix with sub-blocks given by:

∂ 2 f ∂λ∂w ∂ 2 f ∂w 2 (17)

We know that f is C 2 in some neighborhood of (λ 0 , w 0 ), so ∂f /∂w is continuously differentiable in this neighborhood. By assumption, the Hessian ∂ 2 f /∂w 2 is positive definite and hence invertible at (λ 0 , w 0 ). By the Implicit Function Theorem, there exists a neighborhood V of λ 0 and a unique continuously differentiable function w * : V → R m such that ∂f /∂w(λ, w * (λ)) = 0 for λ ∈ V and w * (λ 0 ) = w 0 .

Furthermore, by continuity we know that there is a neighborhood

W 1 × W 2 of (λ 0 , w 0 ) such that ∂ 2 f /∂w 2 is positive definite on this neighborhood. Setting U = V ∩ W 1 ∩ (w * ) -1 (W 2 ), we can conclude that ∂ 2 f /∂w 2 (λ, w * (λ)

) 0 for all λ ∈ U . Combining this with ∂f /∂w(λ, w * (λ)) = 0 and using second-order sufficient optimality conditions, we conclude that w * (λ) is the unique solution to Problem 4b for all λ ∈ U .

B.2 LEMMA 2:

This discussion mostly follows from Hastie et al. (2001) . We let X ∈ R N ×D denote the data matrix where N is the number of training examples and D is the dimensionality of the data. We let t ∈ R N denote the associated targets. We can write the SVD decomposition of X as:

X = U DV (18)

where U and V are N × D and D × D orthogonal matrices and D is a diagonal matrix with entries

d 1 ≥ d 2 ≥ • • • ≥ d D > 0.

We next simplify the function y(x; w) by setting u = s Q, so that y(x; w) = s Qx = u x. We see that the Jacobian ∂y /∂x ≡ u is constant, and Problem 13 simplifies to standard L 2 -regularized least-squares linear regression with the following loss function:

(x,t)∈D (u x -t) 2 + 1 |D| exp(λ) u 2 (19)

It is well-known (see Hastie et al. (2001) , Chapter 3) that the optimal solution u * (λ) minimizing Equation 19 is given by:

u * (λ) = (X X + exp(λ)I) -1 X t = V (D 2 + exp(λ)I) -1 DU t (20)

Furthermore, the optimal solution u * to the unregularized version of Problem 19 is given by:

u * = V D -1 U t (21)

Recall that we defined Q 0 = V , i.e., the change-of-basis matrix from the standard basis to the principal components of the data matrix, and we defined s 0 to solve the unregularized regression problem given Q 0 . Thus, we require that Q 0 s 0 = u * which implies s 0 = D -1 U t.

There are not unique solutions to Problem 13, so we take any functions Q(λ), s(λ) which satisfy Q(λ) s(λ) = v * (λ) as "best-response functions". We will show that our chosen functions Q * (λ) = σ(λv + c) row Q 0 and s * (λ) = s 0 , where v = -1 and c i = 2 log(d i ) for i = 1, . . . , D, meet this criteria. We start by noticing that for any d ∈ R + , we have:

EQUATION

It follows that:

Q * (λ) s * (λ) = [σ(λv + c) row Q 0 ] s 0 (23) =     diag     σ(-λ + 2 log(d 1 )) . . . σ(-λ + 2 log(d D ))     Q 0     s 0 (24) = Q 0      diag      d 2 1 d 2 1 +exp(λ)

. . .

d 2 D d 2 D +exp(λ)           s 0 (25) = V      diag      d 2 1 d 2 1 +exp(λ)

. . .

d 2 D d 2 D +exp(λ)           D -1 U t (26) = V      diag      d1 d 2 1 +exp(λ)

. . .

EQUATION

B.3 THEOREM 3:

By assumption f is quadratic, so there exist

EQUATION

One can easily compute that:

EQUATION

EQUATION

Since we assume ∂ 2 f /∂w 2 0, we must have C 0. Setting the derivative equal to 0 and using second-order sufficient conditions, we have:

EQUATION

Hence, we find:

EQUATION

We let ŵφ (λ) = U λ + b, and define f to be the function given by:

EQUATION

Substituting and simplifying:

f (λ 0 , U , b, σ) = E ∼p( |σ) 1 2 (λ 0 + ) A(λ 0 + ) + (λ 0 + ) B(U (λ 0 + ) + b) + 1 2 (U (λ 0 + ) + b) C(U (λ 0 + ) + b) + d (λ 0 + ) + e (U (λ 0 + ) + b) (36)

Expanding, we find that equation 36 is equal to:

E ∼p( |σ) 1 + 2 + 3 + 4 (37)

where we have:

EQUATION

We can simplify these expressions considerably by using linearity of expectation and that ∼ p( |σ) has mean 0:

EQUATION

EQUATION

We can use the cyclic property of the Trace operator, E ∼p( |σ) [ ] = σ 2 I, and commutability of expectation and a linear operator to simplify the expectations of 2 and 3 :

E ∼p( |σ) 2 = λ 0 BU λ 0 + λ 0 Bb + Tr σ 2 BU (46) E ∼p( |σ) 3 = 1 2 (λ 0 U CU λ 0 + λ 0 U Cb+ Tr σ 2 U CU + b CU λ 0 + b Cb) (47)

We can then differentiate f by making use of various matrix-derivative equalities (Duchi, 2007) to find:

EQUATION

Setting the derivative ∂ f/∂b(λ 0 , U , b, σ) equal to 0, we have:

EQUATION

Setting the derivative for ∂ f/∂U (λ 0 , U , b, σ) equal to 0, we have:

EQUATION

51) Substituting the expression for b given by equation 50 into equation 51 and simplifying gives:

EQUATION

EQUATION

54) This is exactly the best-response Jacobian ∂w * /∂λ(λ) as given by Equation 34. Substituting U = C -1 B into the equation 50 gives:

b = C -1 B λ 0 -C -1 B λ 0 -C -1 e

(55) This is w * (λ 0 ) -∂w * /∂λ(λ 0 ), thus the approximate best-response is exactly the first-order Taylor series of w * about λ 0 .

B.4 BEST-RESPONSE GRADIENT LEMMA:

Lemma 4. Under the same conditions as Lemma 1 and using the same notation, for all λ ∈ U , we have that:

∂w * ∂λ (λ) = - ∂ 2 f ∂w 2 (λ, w * (λ)) -1 ∂ 2 f ∂λ∂w (λ, w * (λ)) (56) Proof. Define ι * : U → R n × R m by ι * (λ) = (λ, w * (λ))

. By first-order optimality conditions, we know that:

EQUATION

Hence, for all λ ∈ U :

EQUATION

Rearranging gives Equation 56.

C BEST-RESPONSE APPROXIMATIONS FOR CONVOLUTIONAL FILTERS:

We let L denote the number of layers, C l the number of channels in layer l's feature map, and K l the size of the kernel in layer l. We let W l,c ∈ R C l-1 ×K l ×K l and b l,c ∈ R denote the weight and bias respectively of the c th convolution kernel in layer l (so c ∈ {1, . . . , C l }). For u l,c , a l,c ∈ R n , we define best-response approximations Ŵ l,c φ and bl,c φ by: Ŵ hyper , b l,c elem }. We can compute the number of parameters used as 2n + 2(|W l,c | + |b l,c |). Summing over channels c, we find the total number of parameters is 2nC l + 2p, where p is the total number of parameters in the normal CNN layer. Hence, we use twice the number of parameters in a normal CNN, plus an overhead that depends on the number of hyperparameters.

For an implementation in code, see Appendix G.

D LANGUAGE MODELING EXPERIMENT DETAILS:

Here we present additional details on the setup of our LSTM language modeling experiments on PTB, and on the role of each hyperparameter we tune.

We trained a 2-layer LSTM with 650 hidden units per layer and 650-dimensional word embeddings (similar to (Zaremba et al., 2014; Gal & Ghahramani, 2016) ) on sequences of length 70 in minibatches of size 40. To optimize the baseline LSTM, we used SGD with initial learning rate 30, which was decayed by a factor of 4 based on the non-monotonic criterion introduced by Merity et al. (2018) (i.e., whenever the validation perplexity fails to improve for 5 epochs). Following Merity et al. (2018) , we used gradient clipping 0.25.

To optimize the ST-LSTM, we used the same optimization setup as for the baseline LSTM. For the hyperparameters, we used Adam with learning rate 0.01. We used an alternating training schedule in which we updated the model parameters for 2 steps on the training set and then updated the hyperparameters for 1 step on the validation set. We used one epoch of warm-up, in which we updated the model parameters, but did not update hyperparameters. We terminated training when the learning rate dropped below 0.0003.

We tuned variational dropout (re-using the same dropout mask for each step in a sequence) on the input to the LSTM, the hidden state between the LSTM layers, and the output of the LSTM. We also tuned embedding dropout, which sets entire rows of the word embedding matrix to 0, effectively removing certain words from all sequences. We regularized the hidden-to-hidden weight matrix using DropConnect (zeroing out weights rather than activations) (Wan et al., 2013) . Because DropConnect operates directly on the weights and not individually on the mini-batch elements, we cannot use independent perturbations per example; instead, we sample a single DropConnect rate per mini-batch. Finally, we used activation regularization (AR) and temporal activation regularization (TAR). AR penalizes large activations, and is defined as:

EQUATION

)

where m is a dropout mask and h t is the output of the LSTM at time t. TAR is a slowness regularizer, defined as:

EQUATION

For AR and TAR, we tuned the scaling coefficients α and β. For the baselines, the hyperparameter ranges were: [0, 0.95] for the dropout rates, and [0, 4] for α and β. For the ST-LSTM, all the dropout rates and the coefficients α and β were initialized to 0.05 (except in Figure 3 , where we varied the output dropout rate).

E IMAGE CLASSIFICATION EXPERIMENT DETAILS:

Here, we present additional details on the CNN experiments. For all results, we held out 20% of the training data for validation.

We trained the baseline CNN using SGD with initial learning rate 0.01 and momentum 0.9, on mini-batches of size 128. We decay the learning rate by 10 each time the validation loss fails to decrease for 60 epochs, and end training if the learning rate falls below 10 -5 or validation loss has not decreased for 75 epochs. For the baselines-grid search, random search, and Bayesian optimization-the search spaces for the hyperparameters were as follows: dropout rates were in the range [0, 0.75]; contrast, saturation, and brightness each had range [0, 1]; hue had range [0, 0.5]; the number of cutout holes had range [0, 4], and the length of each cutout hole had range [0, 24].

We trained the ST-CNN's elementary parameters using SGD with initial learning rate 0.01 and momentum of 0.9, on mini-batches of size 128 (identical to the baselines). We use the same decay schedule as the baseline model. The hyperparameters are optimized using Adam with learning rate 0.003. We alternate between training the best-response approximation and hyperparameters with the same schedule as the ST-LSTM, i.e. T train = 2 steps on the training step and T valid = 1 steps on the validation set. Similarly to the LSTM experiments, we used five epochs of warm-up for the model parameters, during which the hyperparameters are fixed. We used an entropy weight of τ = 0.001 in the entropy regularized objective (Eq. 15). The cutout length was restricted to lie in {0, . . . , 24} while the number of cutout holes was restricted to lie in {0, . . . , 4}. All dropout rates, as well as the continuous data augmentation noise parameters, are initialized to 0.05. The cutout length is initialized to 4, and the number of cutout holes is initialized to 1. Overall, we found the ST-CNN to be relatively robust to the initialization of hyperparameters, but starting with low regularization aided optimization in the first few epochs.

F ADDITIONAL DETAILS ON HYPERPARAMETER SCHEDULES:

Here, we draw connections between hyperparameter schedules and curriculum learning. Curriculum learning (Bengio et al., 2009) is an instance of a family of continuation methods (Allgower & Georg, 2012), which optimize non-convex functions by solving a sequence of functions that are ordered by increasing difficulty. In a continuation method, one considers a family of training criteria C λ (w) with a parameter λ, where C 1 (w) is the final objective we wish to minimize, and C 0 (w) represents the training criterion for a simpler version of the problem. One starts by optimizing C 0 (w) and then gradually increases λ from 0 to 1, while keeping w at a local minimum of C λ (w) (Bengio et al., 2009) . This has been hypothesized to both aid optimization and improve generalization. In this section, we explore how hyperparameter schedules implement a form of curriculum learning; for example, a schedule that increases dropout over time increases stochasticity, making the learning problem more difficult. We use the results of grid searches to understand the effects of different hyperparameter settings throughout training, and show that greedy hyperparameter schedules can outperform fixed hyperparameter values.

First, we performed a grid search over 20 values each of input and output dropout, and measured the validation perplexity in each epoch. Figure 7 shows the validation perplexity achieved by different combinations of input and output dropout, at various epochs during training. We see that at the start of training, the best validation loss is achieved with small values of both input and output dropout.

As we train for more epochs, the best validation performance is achieved with larger dropout rates. Figure 9 shows the perturbed values for output dropout we used to investigate whether the improved performance yielded by STNs is due to the regularization effect, and not the schedule, in Section 4.1. 

G CODE LISTINGS:

In this section, we provide PyTorch code listings for the approximate best-response layers used to construct ST-LSTMs and ST-CNNs: the HyperLinear and HyperConv2D classes. We also provide a simplified version of the optimization steps used on the training set and validation set. 

H SENSITIVITY STUDIES:

In this section, we present experiments that show how sensitive our STN models are to different meta-parameters.

In particular, we investigate the effect of using alternative schedules (Figure 10 ) for the number of optimization steps performed on the training and validation sets.

Additionally, we investigate the effect of using different initial perturbation scales for the hyperparameters, which are either fixed or tuned (Figure 11 ). 

:

We illustrate how this is done for the PyTorch library(Paszke et al., 2017) in Appendix G.

We describe modifications for convolutional filters in Appendix C.

We illustrate how this is done for the PyTorch library(Paszke et al., 2017) in Appendix G.

This is an abuse of notation since there is not a unique solution to Problem 13 for each λ in general.

For grid search and random search we used the Ray Tune libraries (https://github.com/ ray-project/ray/tree/master/python/ray/tune).For Bayesian optimization, we used Spearmint (https://github.com/HIPS/Spearmint).

Figures and Tables:

Figure 1: Best-response architecture for an L2-Jacobian regularized two-layer linear network.

|σ) is Gaussian with mean 0 and variance σ 2 I, and ŵφ is affine. Fix λ 0 ∈ R n and let φ * = arg min φ E ∼p( |σ) [f (λ + , ŵφ (λ + ))]. Then we have ŵφ/∂λ(λ 0 ) = ∂w * /∂λ(λ 0 ). Proof. See Appendix B.3.

Figure2: The effect of the sampled neighborhood. Left: If the sampled neighborhood is too small (e.g., a point mass) the approximation learned will only match the exact best-response at the current hyperparameter, with no guarantee that its gradient matches that of the best-response. Middle: If the sampled neighborhood is not too small or too wide, the gradient of the approximation will match that of the best-response. Right: If the sampled neighborhood is too wide, the approximation will be insufficiently flexible to model the best-response, and again the gradients will not match.

Figure 4: (a) A comparison of the best validation perplexity achieved on PTB over time, by grid search, random search, Bayesian optimization, and STNs. STNs achieve better (lower) validation perplexity in less time than the other methods. (b) The hyperparameter schedule found by the STN for each type of dropout. (c) The hyperparameter schedule found by the STN for the coefficients of activation regularization and temporal activation regularization.

Hyperparameter and elementary parameter dimension f (λ, w), F (λ, w) Lower-level & upper-level objective r Function mapping unconstrained hyperparameters to the appropriate restricted space L T (λ, w), L V (λ, w) Training loss & validation loss -(L T (r(λ), w), L V (r(λ), w)) = (f (λ, w), F (λ, w)) w * (λ) Best-response of the parameters to the hyperparameters F * (λ) Single-level objective from best-response, equals F (λ, w * (λ)) λ * Optimal hyperparameters ŵφ (λ) Parametric approximation to the best-response function φ Approximate best-response parameters σ Scale of the hyperparameter noise distribution σ The sigmoid function Sampled perturbation noise, to be added to hyperparameters p( |σ), p(λ|σ) The noise distribution and induced hyperparameter distribution α A learning rate T train , T valid Number of training steps on the training and validation data x, t An input datapoint and its associated target D A data set consisting of tuples of inputs and targets D The dimensionality of input data y(x, w) Prediction function for input data x and elementary parameters w row Row-wise rescaling -not elementwise multiplication Q, s First and second layer weights of the linear network in Problem 13 Q 0 , s 0 The basis change matrix and solution to the unregularized Problem 13 Q * (λ), s * (λ) The best response weights of the linear network in Problem 13 a(x; w) Activations of hidden units in the linear network of Problem 13 W , b A layer's weight matrix and bias D out , D in A layer's input dimensionality and output dimensionality ∂L T (λ,w) ∂λ

l,c φ (λ) = (λ u l,c ) W l,c hyper + W l{u l,c , a l,c , W l,c hyper , W l,celem , b l,c

Figure 7: Validation performance of a baseline LSTM given different settings of input and output dropout, at various epochs during training. (a), (b), and (c) show the validation performance on PTB given different hyperparameter settings, at epochs 1, 10, and 25, respectively. Darker colors represent lower (better) validation perplexity.Next, we present a simple example to show the potential benefits of greedy hyperparameter schedules. For a single hyperparameter-output dropout-we performed a fine-grained grid search and constructed a dropout schedule by using the hyperparameter values that achieve the best validation perplexity at each epoch in training. As shown in Figure8, the schedule formed by taking the best output dropout value in each epoch yields better generalization than any of the fixed hyperparameter values from the initial grid search. In particular, by using small dropout values at the start of training, the schedule achieves a fast decrease in validation perplexity, and by using larger dropout later in training, it achieves better overall validation perplexity.

Figure 9: Comparison of output dropout schedules. (a) Gaussian-perturbed output dropout rates around the best value found by grid search, 0.68; (b) sinusoid-perturbed output dropout rates with amplitude 0.1 and a period of 1200 mini-batches; (c) the output dropout schedule found by the ST-LSTM.

HyperLinear, used as a drop-in replacement for Linear modules class HyperLinear(nn.Module): def __init__(self, input_dim, output_dim, n_hparams): super(HyperLinear, self).__init__() self.input_dim = input_dim self.output_dim = output_dim self.n_hparams = n_hparams self.n_scalars = output_dim self.elem_w = nn.Parameter(torch.Tensor(output_dim, input_dim)) self.elem_b = nn.Parameter(torch.Tensor(output_dim)) self.hnet_w = nn.Parameter(torch.Tensor(output_dim, input_dim)) self.hnet_b = nn.Parameter(torch.Tensor(output_dim)) self.htensor_to_scalars = nn.Linear(self.n_hparams, self.n_scalars * 2, bias=False) self.init_params() def forward(self, input, hnet_tensor): output = F.linear(input, self.elem_w, self.elem_b) if hnet_tensor is not None: hnet_scalars = self.htensor_to_scalars(hnet_tensor) hnet_wscalars = hnet_scalars[:, :self.n_scalars] hnet_bscalars = hnet_scalars[:, self.n_scalars:] hnet_out = hnet_wscalars * F.linear(input, self.hnet_w) hnet_out += hnet_bscalars * self.hnet_b output += hnet_out return output Listing 2: HyperConv2d, used as a drop-in replacement for Conv2d modules class HyperConv2d(nn.Module): self.htensor_to_scalars.weight.data.normal_(std=0.01) Listing 3: Stylized optimization step on the training set for updating elementary parameters # Perturb hyperparameters around current value in unconstrained # parametrization. batch_htensor = perturb(htensor, hscale) # Apply necessary reparametrization of hyperparameters. hparam_tensor = hparam_transform(batch_htensor) # Sets data augmentation hyperparameters in the data loader. dataset.set_hparams(hparam_tensor) # Get next batch of examples and apply any input transformation # (e.g. input dropout) as dictated by the hyperparameters. images, labels = next_batch(dataset) images = apply_input_transform(images, hparam_tensor) # Run everything through the model and do gradient descent. pred = hyper_cnn(images, batch_htensor, hparam_tensor) xentropy_loss = F.cross_entropy(pred, labels) xentropy_loss.backward() cnn_optimizer.step() Listing 4: Stylized optimization step on the validation set for updating hyperparameters/noise scale # Perturb hyperparameters around current value in unconstrained # parametrization, so we can assess sensitivity of validation # loss to the scale of the noise. batch_htensor = perturb(htensor, hscale) # Apply necessary reparametrization of hyperparameters. hparam_tensor = hparam_transform(batch_htensor) # Get next batch of examples and run through the model. images, labels = next_batch(valid_dataset) pred = hyper_cnn(images, batch_htensor, hparam_tensor) xentropy_loss = F.cross_entropy(pred, labels) # Add extra entropy weight term to loss. entropy = compute_entropy(hscale) loss = xentropy_loss -args.entropy_weight * entropy loss.backward() # Tune the hyperparameters. hyper_optimizer.step() # Tune the scale of the noise applied to hyperparameters. scale_optimizer.step()

Figure10: The effect of using a different number of train/val steps. For the CNN, we include Bayesian Optimization and the reported STN parameters for comparison. During these experiments we found schedules which achieve better final loss with CNNs.

Comparing an LSTM trained with fixed and perturbed output dropouts, an STN, and LSTM trained with the STN schedule.

Table content: <table><tr><td/><td>1.0</td><td/><td/><td/></tr><tr><td>Method p = 0.68, Fixed p = 0.68 w/ Gaussian Noise 85.87 82.29 Val Test 85.83 83.19 p = 0.68 w/ Sinusoid Noise 85.29 82.15 p = 0.78 (Final STN Value) 89.65 86.90 STN 82.58 79.02 LSTM w/ STN Schedule 82.87 79.93</td><td>0.0 0.2 0.4 0.6 0.8 Output Dropout Rate</td><td>0</td><td>5k</td><td>10k 15k 20k 25k Iteration Init=0.05 Init=0.3 Init=0.5 Init=0.7 Init=0.9</td></tr></table>

Final validation and test performance of each method on the PTB word-level language modeling task, and the CIFAR-10 image-classification task.

Table content: <table/>

The hyperparameter schedule prescribed by the STN while training for image classification. The dropouts are indexed by the convolutional layer they are applied to. FC dropout is for the fully-connected layers.ods, achieving lower validation perplexity more quickly. The final validation and test perplexities achieved by each method are shown in Table2. We show the schedules the STN finds for each hyperparameter in Figures4b and 4c; we observe that they are nontrivial, with some forms of dropout used to a greater extent at the start of training (including input and hidden dropout), some used throughout training (output dropout), and some that are increased over the course of training (embedding and weight dropout).

Table content: <table><tr><td>0.0 0.2 0.4 0.6 0.8 1.0 Dropout Rate</td><td>0</td><td>40k</td><td>80k Iteration</td><td>120k Layer 0 Layer 1 Layer 2 Layer 3 Layer 4 Input FC 0 FC 1</td><td>0.0 0.1 0.2 0.3 0.4 0.5 Amount of Noise</td><td>0</td><td>40k</td><td>80k Iteration</td><td colspan="2">120k Hue Contrast Saturation Brightness Inscale</td><td>0 2 4 6 Number/Length of Holes</td><td>0</td><td>40k</td><td>80k Iteration Cutout holes 120k Cutout length</td></tr><tr><td colspan="5">Figure 6: 4.3 IMAGE CLASSIFICATION</td><td/><td/><td/><td/><td/><td/><td/><td/><td/></tr><tr><td/><td/><td/><td/><td/><td/><td/><td/><td/><td>0.5 0.6 0.7 0.8 0.9 1.0 1.1 1.2 1.3 1.4 Best Val Loss</td><td>0</td><td colspan="2">20k</td><td colspan="2">40k Time (s)</td><td>60k Grid Random BayesOpt 80k STN</td></tr><tr><td/><td/><td/><td/><td/><td/><td/><td/><td/><td colspan="6">Figure 5: A comparison of the best validation</td></tr><tr><td/><td/><td/><td/><td/><td/><td/><td/><td/><td colspan="6">loss achieved on CIFAR-10 over time, by grid</td></tr><tr><td/><td/><td/><td/><td/><td/><td/><td/><td/><td colspan="6">search, random search, Bayesian optimization,</td></tr><tr><td/><td/><td/><td/><td/><td/><td/><td/><td/><td colspan="6">and STNs. STNs outperform other methods for</td></tr><tr><td/><td/><td/><td/><td/><td/><td/><td/><td/><td colspan="6">many computational budgets.</td></tr></table>

Table of Notation

Table content: <table/>

Peer review (iclr2019):
 ========. Summary. ========. The paper deals with hyper-parameter optimization of neural networks. The authors formulate the problem as a bilevel optimization problem: minimizing the validation loss over the hyperparameters, subject to the parameters being at the minimum of the training loss. The authors propose an approximation of the so-called best-response function, that maps the hyperparameters to the corresponding optimal parameters (w.r.t the minimization of the training loss), allowing a formulate as a single-level optimization problem and the use gradient descent algorithm. The proposed approximation is based on shifting and scaling the weights and biases of the network. There are no guarantee on its quality except in some very simple cases. The approach assumes a distribution on the hyperparameters, governed by a parameter, which is adapted during the course of the training to achieve a compromise between the flexibility of the best-response function and the quality of its local approximation around the current hyperparameters. The authors show that their approach beats grid-search, random search and Bayesian optimization on the CIFAR-10 and PTB datasets. They point out that the dynamic update of the hyperparameters during the training allows to reach a better performance than any fixed hyperparameter. . ======================. Comments and questions. ======================. Can cross-validation be adapted to this approach? . Can this be used to optimize the learning rate? Which is of course a crucial hyperparameter and that needs an update schedule during the training. . Section 3.2:. "If the entries are too large, then θ̂ φ will not be flexible enough to capture the best- response over the sampled neighborhood. However, its entries must remain sufficiently large so that θ̂ φ captures the local shape around the current hyperparameter values." Not clear why -- more explanations would be helpful. . "minimizing the first term eventually moves all probability mass towards an optimum λ∗ ,resulting in σ = 0". I can*t see how minimizing the first term w.r.t .phi (as in section "2.2.Local approximation") would alter .sigma. . "τ must be set carefully to ensure...". The authors still do not explain how to set .tau. . Section 3.3: . If the hyperparameter is discrete and falls in Case 2, then REINFORCE gradient estimator is used. What about the quality of this gradient? . Section 5, paragraph Gradient-Based HO: "differentiating gradient descent" needs reformulation -- an algorithm cannot be differentiated. . Pros . - The paper is pretty clear . - Generalizes a previous idea and makes it handle discrete hyperparameters and scale better. . - I like the idea of hyperparameters changing dynamically during the training which allows to explore a much larger space than one value . - Although limited, the experimental results are convincing . Cons . - The method itself depends on some parameters and it is not clear how to choose them. Therefore it might be tricky to make it work in practice. I feel like there is a lot of literature around HO but very often people still use the very simple grid/random search, because the alternative methods are often quite complex to implement and make really work. So the fact that the method depends on "crucial" parameters but that are not transparently managed may be a big drawback to its applicability. . - No theoretical guarantee on the quality of the used approximation for neural networks . - Does not handle the learning rate which is a crucial hyperparameter (but maybe it could) .