1 Introduction
While in standard supervised learning problems we seek the best hypothesis in a given space and with a given learning algorithm, in hyperparameter optimization (HO) and metalearning (ML) we seek a configuration so that the optimized learning algorithm will produce a model that generalizes well to new data. The search space in ML often incorporates choices associated with the hypothesis space and the features of the learning algorithm itself (e.g., how optimization of the training loss is performed). Under this common perspective, both HO and ML essentially boil down to nesting two search problems: at the inner level we seek a good hypothesis (as in standard supervised learning) while at the outer level we seek a good configuration (including a good hypothesis space) where the inner search takes place. Surprisingly, the literature on ML has little overlap with the literature on HO and in this paper we present a unified framework encompassing both of them.
Classic approaches to HO (see e.g. Hutter et al., 2015, for a survey) have been only able to manage a relatively small number of hyperparameters, from a few dozens using random search (Bergstra and Bengio, 2012) to a few hundreds using Bayesian or modelbased approaches (Bergstra et al., 2013; Snoek et al., 2012). Recent gradientbased techniques for HO, however, have significantly increased the number of hyperparameters that can be optimized (Domke, 2012; Maclaurin et al., 2015; Pedregosa, 2016; Franceschi et al., 2017)
and it is now possible to tune as hyperparameters entire weight vectors associated with a neural network layer. In this way, it becomes feasible to design models that possibly have more hyperparameters than parameters. Such an approach is well suited for ML, since parameters are learned from a small dataset, whereas hyperparameters leverage multiple available datasets.
HO and ML only differ substantially in terms of the experimental settings in which they are evaluated. While in HO the available data is associated with a single task and split into a training set (used to tune the parameters) and a validation set (used to tune the hyperparameters), in ML we are often interested in the socalled fewshot
learning setting where data comes in the form of short episodes (small datasets with few examples per class) sampled from a common probability distribution over supervised tasks.
Early work on ML dates back at least to the 1990’s (Schmidhuber, 1992; Baxter, 1995; Thrun and Pratt, 1998) but this research area has received considerable attention in the last few years, mainly driven by the need in reallife and industrial scenarios for learning quickly a vast multitude of tasks. These tasks, or episodes, may appear and evolve continuously over time and may only contain few examples (Lake et al., 2017). Different strategies have emerged to tackle ML. Although they do overlap in some aspects, it is possible to identify at least four of them. The metric strategy attempts to use training episodes to construct embeddings such that examples of the same class are mapped into similar representations. It has been instantiated in several variants that involve nonparametric (or instancebased) predictors (Koch et al., 2015; Vinyals et al., 2016; Snell et al., 2017). In the related memorization strategy, the metalearner learns to store and retrieve data points representations in memory. It can be implemented either using recurrent networks (Santoro et al., 2016) or temporal convolutions (Mishra et al., 2018). The use of an attention mechanism (Vaswani et al., 2017) is crucial both in (Vinyals et al., 2016) and in (Mishra et al., 2018). The initialization strategy (Ravi and Larochelle, 2017; Finn et al., 2017) uses training episodes to infer a good initial value for the model’s parameters so that new tasks can be learned quickly by fine tuning. The optimization strategy (Andrychowicz et al., 2016; Ravi and Larochelle, 2017; Wichrowska et al., 2017) forges an optimization algorithm that will find it easier to learn on novel related tasks.
A main contribution of this paper is a unified view of HO and ML within the natural mathematical framework of bilevel programming, where an outer optimization problem is solved subject to the optimality of an inner optimization problem. In HO the outer problem involves hyperparameters while the inner problem is usually the minimization of an empirical loss. In ML the outer problem could involve a shared representation among tasks while the inner problem could concern classifiers for individual tasks. Bilevel programming
(Bard, 2013)has been suggested before in machine learning in the context of kernel methods and support vector machines
(Keerthi et al., 2007; Kunapuli et al., 2008), multitask learning (Flamary et al., 2014), and more recently HO (Pedregosa, 2016), but never in the context of ML. The resulting framework outlined in Sec. 2 encompasses some existing approaches to ML, in particular those based on the initialization and the optimization strategies.
Bilevel
programming 
Hyperparameter
optimization 
Metalearning 
Inner variables  Parameters 
Parameters of
Ground models 
Outer variables  Hyperparameters 
Parameters of
Metalearner 
Inner objective  Training error 
Training errors
on tasks (Eq. 3) 
Outer objective  Validation error 
Metatraining
error (Eq. 4) 
A technical difficulty arises when the solution to the inner problem cannot be written analytically (for example this happens when using the logloss for training neural networks) and one needs to resort to iterative optimization approaches. As a second contribution, we provide in Sec. 3 sufficient conditions that guarantee good approximation properties. We observe that these conditions are reasonable and apply to concrete problems relevant to applications.
In Sec. 4, by taking inspiration on early work on representation learning in the context of multitask and metalearning (Baxter, 1995; Caruana, 1998), we instantiate the framework for ML in a simple way treating the weights of the last layer of a neural network as the inner variables and the remaining weights, which parametrize the representation mapping, as the outer variables. As shown in Sec. 5, the resulting ML algorithm performs well in practice, outperforming most of the existing strategies on MiniImagenet.
2 A bilevel optimization framework
In this paper, we consider bilevel optimization problems (see e.g. Colson et al., 2007) of the form
(1) 
where function is defined at as
(2) 
We call the outer objective and, for every , we call the inner objective. Note that is a class of objective functions parameterized by . Specific instances of this problem include HO and ML, which we discuss next. Table 1 outlines the links among bilevel programming, HO and ML.
2.1 Hyperparameter Optimization
In the context of hyperparameter optimization, we are interested in minimizing the validation error of a model parameterized by a vector , with respect to a vector of hyperparameters . For example, we may consider representation or regularization hyperparameters that control the hypothesis space or penalties, respectively. In this setting, a prototypical choice for the inner objective is the regularized empirical error
where is a set of input/output points,
is a prescribed loss function, and
a regularizer parameterized by . The outer objective represents a proxy for the generalization error of , and it may be given by the average loss on a validation setor, in more generality, by a crossvalidation error, as detailed in Appendix B. Note that in this setting, the outer objective does not depend explicitly on the hyperparameters , since in HO is instrumental in finding a good model , which is our final goal. As a more specific example, consider linear models, , let be the square loss and let
, in which case the inner objective is ridge regression (Tikhonov regularization) and the bilevel problem optimizes over the regularization parameter the validation error of ridge regression.
2.2 MetaLearning
In metalearning (ML) the inner and outer objectives are computed by averaging a training and a validation error over multiple tasks, respectively. The goal is to produce a learning algorithm that will work well on novel tasks^{1}^{1}1The ML problem is also related to multitask learning, however in ML the goal is to extrapolate from the given tasks.. For this purpose, we have available a metatraining set , which is a collection of datasets, sampled from a metadistribution . Each dataset with is linked to a specific task. Note that the output space is task dependent (e.g. a multiclass classification problem with variable number of classes). The model for each task is a function , identified by a parameter vectors and hyperparameters . A key point here is that is shared between the tasks. With this notation the inner and outer objectives are
(3) 
(4) 
respectively. The loss represents the empirical error of the pair on a set of examples . Note that the inner and outer losses for task use different train/validation splits of the corresponding dataset . Furthermore, unlike in HO, in ML the final goal is to find a good and the are now instrumental.
The cartoon in Figure 1 illustrates ML as a bilevel problem. The parameter indexes an hypothesis space within which the inner objective is minimized. A particular example, detailed in Sec. 4, is to choose the model , in which case parameterizes a feature mapping. Yet another choice would be to consider , in which case represents a common model around which task specific models are to be found (see e.g. Evgeniou et al., 2005; Finn et al., 2017; Khosla et al., 2012; Kuzborskij et al., 2013, and reference therein).
2.3 GradientBased Approach
We now discuss a general approach to solve Problem (1)(2) when the hyperparameter vector is realvalued. To simplify our discussion let us assume that the inner objective has a unique minimizer . Even in this simplified scenario, Problem (1)(2) remains challenging to solve. Indeed, in general there is no closed form expression , so it is not possible to directly optimize the outer objective function. While a possible strategy (implicit differentiation) is to apply the implicit function theorem to (Pedregosa, 2016; Koh and Liang, 2017; Beirami et al., 2017), another compelling approach is to replace the inner problem with a dynamical system. This point, discussed in (Domke, 2012; Maclaurin et al., 2015; Franceschi et al., 2017), is developed further in this paper.
Specifically, we let where is a prescribed positive integer and consider the following approximation of Problem (1)(2)
(5) 
where is a smooth scalar function, and^{2}^{2}2In general, the algorithm used to minimize the inner objective may involve auxiliary variables, e.g., velocities when using gradient descent with momentum, so should be intended as a larger vector containing both model parameters and auxiliary variables.
(6) 
with a smooth initialization mapping and, for every , a smooth mapping that represents the operation performed by the th step of an optimization algorithm. For example, the optimization dynamics could be gradient descent: where is a sequence of steps sizes.
The approximation of the bilevel problem (1)(2) by the procedure (5)(6) raises the issue of the quality of this approximation and we return to this issue in the next section. However, it also suggests to consider the inner dynamics as a form of approximate empirical error minimization (e.g. early stopping) which is valid in its own right. From this perspective – conversely to the implicit differentiation strategy – it is possible to include among the components of variables which are associated with the optimization algorithm itself. For example, may include the step sizes or momentum factors if the dynamics in Eq. (6) is gradient descent with momentum; in (Andrychowicz et al., 2016; Wichrowska et al., 2017) the mapping
is implemented as a recurrent neural network, while
(Finn et al., 2017) focus on the initialization mapping by letting .A major advantage of this reformulation is that it makes it possible to compute efficiently the gradient of , which we call hypergradient, either in time or in memory (Maclaurin et al., 2015; Franceschi et al., 2017), by making use of reverse or forward mode algorithmic differentiation (Griewank and Walther, 2008; Baydin et al., 2017). This allows us to optimize a number of hyperparameters of the same order of that of parameters, a situation which arise in ML.
3 Exact and Approximate Bilevel Programming
In this section, we provide results about the existence of solutions of Problem (1)(2) and the approximation properties of Procedure (5)(6) with respect to the original bilevel problem. Proofs of these results are provided in the supplementary material.
Procedure (5)(6), though related to the bilevel problem (1)(2), may not be, in general, a good approximation of it. Indeed, making the assumptions (which sound perfectly reasonable) that, for every , for some , and that is continuous, one can only assert that . This is because the optimization dynamics converge to some minimizer of the inner objective , but not necessarily to the one that also minimizes the function . This is illustrated in Figure 2. The situation is, however, different if the inner problem admits a unique minimizer for every . Indeed in this case, it is possible to show that the set of minimizers of the approximate problems converge, as and in an appropriate sense, to the set of minimizers of the bilevel problem. More precisely, we make the following assumptions:

is a compact subset of ;

is jointly continuous;

the map is jointly continuous and such that is a singleton for every ;

remains bounded as varies in .
(7) 
Under the above assumptions, in the following we give results about the existence of solutions of problem (7) and the (variational) convergence of the approximate problems (5)(6) towards problem (7) — relating the minima as well as the set of minimizers. In this respect we note that, since both and are nonconvex, and are in general nonsingleton, so an appropriate definition of set convergence is required.
Theorem 3.1 (Existence).
Under Assumptions (i)(iv) problem (7) admits solutions.
Proof See Appendix A.
The result below follows from general facts on the stability of minimizers in optimization problems (Dontchev and Zolezzi, 1993).
Theorem 3.2 (Convergence).
In addition to Assumptions (i)(iv), suppose that:

is uniformly Lipschitz continuous;

The iterates converge uniformly to on as .
Then

,

, meaning that, for every such that , we have that:

admits a convergent subsequence;

for every subsequence such that , we have .

Proof See Appendix A.
We stress that assumptions (i)(vi) are very natural and satisfied by many problems of practical interests. Thus, the above results provide full theoretical justification to the proposed approximate procedure (5)(6). The following remark discusses assumption (vi), while the subsequent example will be relevant to the experiments in Sec. 5.
Remark 3.3.
If is strongly convex, then many gradientbased algorithms (e.g., standard and accelerated gradient descent) yield linear convergence of the iterates ’s. Moreover, in such cases, the rate of linear convergence is of type , where and are the Lipschitz constant of the gradient and the modulus of strong convexity of respectively. So, this rate can be uniformly bounded from above by , provided that and . Thus, in these cases converges uniformly to on (at a linear rate).
Example 3.4.
Let us consider the following form of the inner objective:
(8) 
where is a fixed regularization parameter and is the hyperparameter, representing a linear feature map. is strongly convex, with modulus (independent on the hyperparameter ), and Lipschitz smooth with constant , which is bounded from above, if ranges in a bounded set of square matrices. In this case assumptions (i)(vi) are satisfied.
4 Learning HyperRepresentations
In this section, we instantiate the bilevel programming approach for ML outlined in Sec. 2.2 in the case of deep learning where representation layers are shared across episodes. Finding good data representations is a centerpiece in machine learning. Classical approaches (Baxter, 1995; Caruana, 1998) learn both the weights of the representation mapping and those of the ground classifiers jointly on the same data. Here we follow the bilevel approach and split each dataset/episode in training and validation sets.
Our method involves the learning of a crosstask intermediate representation (parametrized by a vector ) on top of which task specific models (parametrized by vectors ) are trained. The final ground model for task is thus given by . To find , we solve Problem (1)(2) with inner and outer objectives as in Eqs. (3) and (4), respectively. Since, in general, this problem cannot be solved exactly, we instantiate the approximation scheme in Eqs. (5)(6) as follows:
(9) 
(10) 
Starting from an initial value, the weights of the taskspecific models are learned by iterations of gradient descent. The gradient of can be computed efficiently in time by making use of an extended reversehypergradient procedure (Franceschi et al., 2017) which we present in Algorithm 1. Since, in general, the number of episodes in a metatraining set is large, we compute a stochastic approximation of the gradient of by sampling a minibatch of episodes. At test time, given a new episode , the representation is kept fixed, and all the examples in are used to tune the weights of the episodespecific model .
Like other initialization and optimization strategies for ML, our method does not require lookups in a support set as the memorization and metric strategies do (Santoro et al., 2016; Vinyals et al., 2016; Mishra et al., 2018). Unlike (Andrychowicz et al., 2016; Ravi and Larochelle, 2017) we do not tune the optimization algorithm, which in our case is plain empirical loss minimization by gradient descent, and rather focus on the hypothesis space. Unlike (Finn et al., 2017), that aims at maximizing sensitivity of new task losses to the model parameters, we aim at maximizing the generalization to novel examples during training episodes, with respect to . Our assumptions about the structure of the model are slightly stronger than in (Finn et al., 2017) but still mild, namely that some (hyper)parameters define the representation and the remaining parameters define the classification function. In (Munkhdalai and Yu, 2017) the metaknowledge is distributed among fast and slow weights and an external memory; our approach is more direct, since the metaknowledge is solely distilled by . A further advantage of our method is that, if the episodespecific models are linear (e.g. logistic regressors) and each loss is strongly convex in , the theoretical guarantees of Theorem 3.2 apply (see Remark 3.3). These assumptions are satisfied in the experiments reported in the next section.
5 Experiments
The aim of the following experiments is threefold. First, we investigate the impact of the number of iterations of the optimization dynamics on the quality of the solution on a simple multiclass classification problem. Second, we test our hyperrepresentation method in the context of fewshot learning on two benchmark datasets. Finally, we constrast the bilevel ML approach against classical approaches to learn shared representations ^{3}^{3}3The code for reproducing the experiments, based on the package FarHO (https://bit.ly/farho), is available at https://bit.ly/hyperrepr.
5.1 The Effect of
Motivated by the theoretical findings of Sec. 3, we empirically investigate how solving the inner problem approximately (i.e. using small ) affects convergence, generalization performances, and running time. We focus in particular on the linear feature map described in Example 3.4, which allows us to compare the approximated solution against the closedform analytical solution given by
In this setting, the bilevel problem reduces to a (nonconvex) optimization problem in .
We use a subset of 100 classes extracted from Omniglot dataset (Lake et al., 2017) to construct a HO problem aimed at tuning . A training set and a validation set , each consisting of three randomly drawn examples per class, were sampled to form the HO problem. A third set , consisting of fifteen examples per class, was used for testing. Instead of using raw images as input, we employ feature vectors computed by the convolutional network trained on oneshot fiveways ML setting as described in Sec. 5.2.
For the approximate problems we compute the hypergradient using Algorithm 1, where it is intended that . Figure 3 shows the values of functions and (see Eqs. (1) and (5), respectively) during the optimization of . As increases, the solution of the approximate problem approaches the true bilevel solution. However, performing a small number of gradient descent steps for solving the inner problem acts as implicit regularizer. As it is evident from Figure 4, the generalization error is better when is smaller than the value yielding the best approximation of the inner solution.
1  4  16  64  256  Exact  
Time (sec)  60  119  356  1344  5532  320 
This is to be expected since, in this setting, the dimensions of parameters and hyperparameters are of the same order, leading to a concrete possibility of overfitting the outer objective (validation error). An appropriate, problem dependent, choice of may help avoiding this issue (see also Appendix C). As increases, the number of hyperiterations required to reach the maximum test accuracy decreases, further suggesting that there is an interplay between the number of iterations used to solve the inner and the outer objective. Finally, the running time of Algorithm 1, is linear in and the size of and independent of the size of (see also Table 2), making it even more appealing to reduce the number of iterations.
5.2 Fewshot Learning
We new turn our attention to learningtolearn, precisely to fewshot supervised learning, implementing the ML strategy outlined in Sec. 4 on two different benchmark datasets:
Omniglot (Lake et al., 2015), a dataset that contains examples of 1623 different handwritten characters from 50 alphabets. We downsample the images to .
MiniImagenet (Vinyals et al., 2016)
, a subset of ImageNet
(Deng et al., 2009), that contains 60000 downsampled images from 100 different classes.Following the experimental protocol used in a number of recent works, we build a metatraining set , from which we sample datasets to solve Problem (9)(10), a metavalidation set for tuning ML hyperparameters, and finally a metatest set
which is used to estimate accuracy. Operationally, each metadataset consists of a pool of samples belonging to different (nonoverlapping between separate metadataset) classes, which can be combined to form ground classification datasets
with 5 or 20 classes (for Omniglot). The ’s contain 1 or 5 examples per class which are used to fit (see Eq. 10). The ’s, containing 15 examples per class, is used either to compute (see Eq. (9)) and its (stochastic) gradient if or to provide a generalization score if comes from either or . For MiniImagenet we use the same split and images proposed in (Ravi and Larochelle, 2017), while for Omniglot we use the protocol defined by (Santoro et al., 2016).As ground classifiers we use multinomial logistic regressors and as task losses we employ crossentropy. The inner problems, being strongly convex, admit unique minimizers, yet require numerical computation of the solutions. We initialize ground models parameters to and, according to the observation in Sec. 5.1, we perform gradient descent steps, where is treated as a ML hyperparameter that has to be validated. Figure 6 shows an example of metavalidation of for oneshot learning on MiniImagenet. We compute a stochastic approximation of with Algorithm 1 and use Adam with decaying learning rate to optimize .
Regarding the specific implementation of the representation mapping
, we employ for Omniglot a fourlayers convolutional neural network with strided convolutions and 64 filters per layer as in
(Vinyals et al., 2016) and other successive works. For MiniImagenet we tried two different architectures:C4L
, a fourlayers convolutional neural network with maxpooling and 32 filters per layer;
RN: a residual network (He et al., 2016) built of four residual blocks followed by two convolutional layers.
The first network architecture has been proposed in (Ravi and Larochelle, 2017) and then used in (Finn et al., 2017), while a similar residual network architecture has been employed in a more recent work (Mishra et al., 2018). Further details on the architectures of , as well as other ML hyperparameters, are specified in the supplementary material. We report our results, using RN for MiniImagenet, in Table 3, alongside scores from various recently proposed methods for comparison.
Omniglot 5 classes  Omniglot 20 classes  MiniImagenet 5 classes  
Method  1shot  5shot  1shot  5shot  1shot  5shot 
Siamese nets (Koch et al., 2015)  
Matching nets (Vinyals et al., 2016)  
Neural stat. (Edwards and Storkey, 2016)  
Memory mod. (Kaiser et al., 2017)  
MetaLSTM (Ravi and Larochelle, 2017)  
MAML (Finn et al., 2017)  
Metanetworks (Munkhdalai and Yu, 2017)  
Prototypical Net. (Snell et al., 2017)  
SNAIL (Mishra et al., 2018)  
Hyperrepresentation  98.6  99.5  95.5  98.4 
, of various methods on 1shot and 5shot classification problems on Omniglot and MiniImagenet. For MiniImagenet 95% confidence intervals are reported. For Hyperrepresentation the scores are computed over 600 randomly drawn episodes. For other methods we show results as reported by their respective authors.
The proposed method achieves competitive results highlighting the relative importance of learning a task independent representation, on the top of which logistic classifiers trained with very few samples generalize well. Moreover, utilizing more expressive models such as residual network as representation mappings, is beneficial for our proposed strategy and, unlike other methods, does not result in overfitting of the outer objective, as reported in (Mishra et al., 2018). Indeed, compared to C4L, RN achieves a relative improvement of 6.5% on oneshot and 4.2% on fiveshot. Figure 5 provides a visual example of the goodness of the learned representation, showing that MiniImagenet examples (the first from metatraining, the second from the metatesting sets) from similar classes (different dog breeds) are mapped near each other by and, conversely, samples from dissimilar classes are mapped afar.
5.3 On Variants of Representation Learning Methods
In this section, we show the benefits of learning a representation within the proposed bilevel framework compared to other possible approaches that involve an explicit factorization of a classifier as . The representation mapping is either pretrained or learned with different metalearning algorithms. We focus on the problem of oneshot learning on MiniImagenet and we use C4L as architecture for the representation mapping. In all the experiments the ground models are multinomial logistic regressor as in Sec. 5.2, tuned with 5 steps of gradient descent. We ran the following experiments:
Multiclass: the mapping is given by the linear outputs before the softmax operation of a network^{4}^{4}4The network is similar to C4L but has 64 filters per layer. pretrained on the totality of examples contained in the training metadataset (600 examples for each of the 64 classes). In this setting, we found that using the second last layer or the output after the softmax yields worst results;
Bileveltrain: we use a bilevel approach but, unlike in Sec. 4, we optimize the parameter vector of the representation mapping by minimizing the loss on the training sets of each episode. The hypergradient is still computed with Algorithm 1, albeit we set for each training episodes;
Approx and Approxtrain: we consider an approximation of the hypergradient by disregarding the optimization dynamics of the inner objectives (i.e. we set ). In Approxtrain we just use the training sets;
Classic: as in (Baxter, 1995), we learn by jointly optimize and treat the problem as standard multitask learning, with the exception that we evaluate on minibatches of 4 episodes, randomly sampled every 5 gradient descent iterations.
In settings where we do not use the validation sets, we let the training sets of each episode contain 16 examples per class. Using training episodes with just one example per class resulted in performances just above random chance. While the first experiment constitutes a standard baseline, the others have the specific aim of assessing () the importance of splitting episodes of metatraining set into training and validation and () the importance of computing the hypergradient of the approximate bilevel problem with Algorithm 1. The results reported in Table 4 suggest that both the training/validation splitting and the full computation of the hypergradient constitute key factors for learning a good representation in a metalearning context.
Method  # filters  Accuracy 1shot 
Multiclass  64  
Bileveltrain  32  
Approx  32  
Approxtrain  32  
Classictrain  32  
HyperrepresentationC4L  32 
On the other side, using pretrained representations, especially in a lowdimensional space, turns out to be a rather effective baseline. One possible explanation is that, in this context, some classes in the training and testing metadatasets are rather similar (e.g. various dog breeds) and thus ground classifiers can leverage on very specific representations.
6 Conclusions
We have shown that both HO and ML can be formulated in terms of bilevel programming and solved with an iterative approach. When the inner problem has a unique solution (e.g. is strongly convex), our theoretical results show that the iterative approach has convergence guarantees, a result that is interesting in its own right. In the case of ML, by adapting classical strategies (Baxter, 1995) to the bilevel framework with training/validation splitting, we present a method for learning hyperrepresentations which is experimentally effective and supported by our theoretical guarantees.
Our framework encompasses recently proposed methods for metalearning, such as learning to optimize, but also suggests different design patterns for the inner learning algorithm which could be interesting to explore in future work. The resulting inner problems may not satisfy the assumptions of our convergence analysis, raising the need for further theoretical investigations. An additional future direction of research is the study of the statistical properties of bilevel strategies where outer objectives are based on the generalization ability of the inner model to new (validation) data. Ideas from (Maurer et al., 2016; Denevi et al., 2018) may be useful in this direction.
References
 Andrychowicz et al. (2016) Andrychowicz, M., Denil, M., Gomez, S., Hoffman, M. W., Pfau, D., Schaul, T., and de Freitas, N. (2016). Learning to learn by gradient descent by gradient descent. In Advances in Neural Information Processing Systems (NIPS), pages 3981–3989.
 Bard (2013) Bard, J. F. (2013). Practical bilevel optimization: algorithms and applications, volume 30. Springer Science & Business Media. 01251.

Baxter (1995)
Baxter, J. (1995).
Learning internal representations.
In
Proceedings of the 8th Annual Conference on Computational Learning Theory (COLT)
, pages 311–320. ACM.  Baydin et al. (2017) Baydin, A. G., Pearlmutter, B. A., Radul, A. A., and Siskind, J. M. (2017). Automatic differentiation in machine learning: a survey. Journal of Machine Learning Research, 18:153:1–153:43.
 Beirami et al. (2017) Beirami, A., Razaviyayn, M., Shahrampour, S., and Tarokh, V. (2017). On optimal generalizability in parametric learning. In Advances in Neural Information Processing Systems (NIPS), pages 3458–3468.
 Bergstra and Bengio (2012) Bergstra, J. and Bengio, Y. (2012). Random search for hyperparameter optimization. Journal of Machine Learning Research, 13(Feb):281–305.
 Bergstra et al. (2013) Bergstra, J., Yamins, D., and Cox, D. D. (2013). Making a science of model search: Hyperparameter optimization in hundreds of dimensions for vision architectures. In Proceedings of the 30th International Conference on Machine Learning, (ICML), pages 115–123.
 Caruana (1998) Caruana, R. (1998). Multitask learning. In Learning to learn, pages 95–133. Springer. 02683.
 Colson et al. (2007) Colson, B., Marcotte, P., and Savard, G. (2007). An overview of bilevel optimization. Annals of operations research, 153(1):235–256.
 Denevi et al. (2018) Denevi, G., Ciliberto, C., Stamos, D., and Pontil, M. (2018). Incremental learningtolearn with statistical guarantees. arXiv preprint arXiv:1803.08089. To appear in UAI 2018.
 Deng et al. (2009) Deng, J., Dong, W., Socher, R., Li, L.J., Li, K., and FeiFei, L. (2009). Imagenet: A largescale hierarchical image database. In Computer Vision and Pattern Recognition (CVPR), pages 248–255.
 Domke (2012) Domke, J. (2012). Generic Methods for OptimizationBased Modeling. In AISTATS, volume 22, pages 318–326.
 Dontchev and Zolezzi (1993) Dontchev, A. L. and Zolezzi, T. (1993). Wellposed optimization problems, volume 1543 of Lecture Notes in Mathematics. SpringerVerlag, Berlin.
 Edwards and Storkey (2016) Edwards, H. and Storkey, A. (2016). Towards a Neural Statistician. arXiv:1606.02185 [cs, stat]. 00027 arXiv: 1606.02185.
 Evgeniou et al. (2005) Evgeniou, T., Micchelli, C. A., and Pontil, M. (2005). Learning multiple tasks with kernel methods. Journal of Machine Learning Research, 6(Apr):615–637.
 Finn et al. (2017) Finn, C., Abbeel, P., and Levine, S. (2017). Modelagnostic metalearning for fast adaptation of deep networks. In Proceedings of the 34th International Conference on Machine Learning, (ICML), pages 1126–1135.
 Flamary et al. (2014) Flamary, R., Rakotomamonjy, A., and Gasso, G. (2014). Learning constrained task similarities in graphregularized multitask learning. Regularization, Optimization, Kernels, and Support Vector Machines, page 103.
 Franceschi et al. (2017) Franceschi, L., Donini, M., Frasconi, P., and Pontil, M. (2017). Forward and reverse gradientbased hyperparameter optimization. In Proceedings of the 34th International Conference on Machine Learning, (ICML), pages 1165–1173.

Glorot and Bengio (2010)
Glorot, X. and Bengio, Y. (2010).
Understanding the difficulty of training deep feedforward neural
networks.
In
Proceedings of the 13th International Conference on Artificial Intelligence and Statistics (AIStat)
, pages 249–256.  Griewank and Walther (2008) Griewank, A. and Walther, A. (2008). Evaluating derivatives: principles and techniques of algorithmic differentiation. SIAM.
 He et al. (2016) He, K., Zhang, X., Ren, S., and Sun, J. (2016). Deep residual learning for image recognition. In Computer Vision and Pattern Recognition (CVPR), pages 770–778.
 Hutter et al. (2015) Hutter, F., Lücke, J., and SchmidtThieme, L. (2015). Beyond Manual Tuning of Hyperparameters. KI  Künstliche Intelligenz, 29(4):329–337.
 Kaiser et al. (2017) Kaiser, L., Nachum, O., Roy, A., and Bengio, S. (2017). Learning to remember rare events.
 Keerthi et al. (2007) Keerthi, S. S., Sindhwani, V., and Chapelle, O. (2007). An efficient method for gradientbased adaptation of hyperparameters in svm models. In Advances in Neural Information Processing Systems (NIPS), pages 673–680.
 Khosla et al. (2012) Khosla, A., Zhou, T., Malisiewicz, T., Efros, A. A., and Torralba, A. (2012). Undoing the damage of dataset bias. In European Conference on Computer Vision, pages 158–171. Springer.
 Koch et al. (2015) Koch, G., Zemel, R., and Salakhutdinov, R. (2015). Siamese neural networks for oneshot image recognition. In ICML Deep Learning Workshop, volume 2. 00127.
 Koh and Liang (2017) Koh, P. W. and Liang, P. (2017). Understanding blackbox predictions via influence functions. In Proceedings of the 34th International Conference on Machine Learning (ICML), pages 1885–1894.
 Kunapuli et al. (2008) Kunapuli, G., Bennett, K., Hu, J., and Pang, J.S. (2008). Classification model selection via bilevel programming. Optimization Methods and Software, 23(4):475–489.
 Kuzborskij et al. (2013) Kuzborskij, I., Orabona, F., and Caputo, B. (2013). From n to n+ 1: Multiclass transfer incremental learning. In Computer Vision and Pattern Recognition (CVPR), 2013 IEEE Conference on, pages 3358–3365.
 Lake et al. (2015) Lake, B. M., Salakhutdinov, R., and Tenenbaum, J. B. (2015). Humanlevel concept learning through probabilistic program induction. Science, 350(6266):1332–1338.
 Lake et al. (2017) Lake, B. M., Ullman, T. D., Tenenbaum, J. B., and Gershman, S. J. (2017). Building machines that learn and think like people. Behavioral and Brain Sciences, 40. 00152.
 Maclaurin et al. (2015) Maclaurin, D., Duvenaud, D. K., and Adams, R. P. (2015). Gradientbased hyperparameter optimization through reversible learning. In Proceedings of the 32nd International Conference on Machine Learning, (ICML, pages 2113–2122.
 Maurer et al. (2016) Maurer, A., Pontil, M., and RomeraParedes, B. (2016). The benefit of multitask representation learning. The Journal of Machine Learning Research, 17(1):2853–2884.
 Mishra et al. (2018) Mishra, N., Rohaninejad, M., Chen, X., and Abbeel, P. (2018). A simple neural attentive metalearner.
 Munkhdalai and Yu (2017) Munkhdalai, T. and Yu, H. (2017). Meta networks. In Proceedings of the 34th International Conference on Machine Learning, (ICML), pages 2554–2563.
 Pedregosa (2016) Pedregosa, F. (2016). Hyperparameter optimization with approximate gradient. In Proceedings of The 33rd International Conference on Machine Learning (ICML), pages 737–746.
 Ravi and Larochelle (2017) Ravi, S. and Larochelle, H. (2017). Optimization as a model for fewshot learning. In In International Conference on Learning Representations (ICLR).
 Santoro et al. (2016) Santoro, A., Bartunov, S., Botvinick, M., Wierstra, D., and Lillicrap, T. (2016). Metalearning with memoryaugmented neural networks. In Proceedings ofthe 33rd International Conference on Machine Learning, pages 1842–1850.
 Schmidhuber (1992) Schmidhuber, J. (1992). Learning to Control FastWeight Memories: An Alternative to Dynamic Recurrent Networks. Neural Computation, 4(1):131–139. 00082.
 Snell et al. (2017) Snell, J., Swersky, K., and Zemel, R. S. (2017). Prototypical networks for fewshot learning. In Advances in neural information processing systems, pages 4080–4090.
 Snoek et al. (2012) Snoek, J., Larochelle, H., and Adams, R. P. (2012). Practical bayesian optimization of machine learning algorithms. In Advances in Neural Information Processing Systems (NIPS), pages 2951–2959.
 Thrun and Pratt (1998) Thrun, S. and Pratt, L. (1998). Learning to learn. Springer.
 Vaswani et al. (2017) Vaswani, A., Shazeer, N., Parmar, N., Uszkoreit, J., Jones, L., Gomez, A. N., Kaiser, L., and Polosukhin, I. (2017). Attention is all you need. In Advances in Neural Information Processing Systems (NIPS), pages 6000–6010.
 Vinyals et al. (2016) Vinyals, O., Blundell, C., Lillicrap, T., Kavukcuoglu, K., and Wierstra, D. (2016). Matching networks for one shot learning. In Advances in Neural Information Processing Systems (NIPS), pages 3630–3638.
 Wichrowska et al. (2017) Wichrowska, O., Maheswaranathan, N., Hoffman, M. W., Colmenarejo, S. G., Denil, M., Freitas, N., and SohlDickstein, J. (2017). Learned optimizers that scale and generalize. In Proceedings of the 34th International Conference on Machine Learning (ICML), pages 3751–3760.
Appendix A Proofs of the Results in Sec. 3
Proof of Theorem 3.1. Since is compact, it follows from Weierstrass theorem that a sufficient condition for the existence of minimizers is that is continuous. Thus, let and let be a sequence in such that . We prove that . Since is bounded, there exists a subsequence such that for some . Now, since and the map is jointly continuous, we have
Therefore, is a minimizer of and hence . This prove that is a bounded sequence having a unique cluster point. Hence is convergence to its unique cluster point, which is . Finally, since and is jointly continuous, we have and the statement follows. ∎
We recal a fundamental fact concerning the stability of minima and minimizers in optimization problems (Dontchev and Zolezzi, 1993). We provide the proof for completeness.
Theorem A.1 (Convergence).
Let and be lower semicontinuous functions defined on a compact set . Suppose that converges uniformly to on as . Then

,

, meaning that, for every such that , we have that:

admits a convergent subsequence;

for every subsequence such that , we have .

Proof.
Let be a sequence in such that, for every , . We prove that

admits a convergent subsequence.

for every subsequence such that , we have and .

.
The first point follows from the fact that is compact.
Concerning the second point, let be a subsequence
such that . Since converge uniformly
to , we have
Therefore, using also the continuity of , we have
So, and , that is, .
Finally, as regards the last point, we proceed by contradiction.
If does not convergce to ,
then there exists an and a subsequence
such that
(11) 
Now, let be a convergent subsequence of . Suppose that . Clearly is also a subsequence of . Then, it follows from point 2) above that . This latter finding together with equation (11) gives a contradiction. ∎
Appendix B Crossvalidation and Bilevel Programming
We note that the (approximate) bilevel programming framework easily accommodates also estimations of the generalization error generated by a crossvalidation procedures. We describe here the case of fold crossvalidation, which includes also leaveoneout cross validation.
Let be the set of available data; fold cross validation, with consists in partitioning in subsets and fit as many models on training data . The models are then evaluated on . Denoting by the vector of stacked weights, the fold cross validation error is given by
where . can be treated as the outer objective in the bilevel framework, while the inner objective may be given by the sum of regularized empirical errors over each for the models. Under this perspective, a fold crossvalidation procedure closely resemble the bilevel problem for ML formulated in Sec. 2.2, where, in this case, the metadistribution collapses on the data (ground) distribution and the episodes are sampled from the same dataset of points.
By following the procedure outlined in Sec. 2.3 we can approximate the minimization of with steps of an optimization dynamics and compute the hypergradient of by training the models and proceed with either forward or reverse differentiation. The models may be fitted sequentially, in parallel or stochastically. Specifically, in this last case, one can repeatedly sample one fold at a time (or possibly a minibatch of folds) and compute a stochastic hypergradient that can be used in a SGD procedure in order to minimize . At last, we note that similar ideas for leaveone out crossvalidation error are developed in (Beirami et al., 2017), where the hypergradient of an approximation of the outer objective is computed by means of the implicit function theorem.
Appendix C The Effect of : Ridge Regression
In Sec. 5.1 we showed that increasing the number of iterations leads to a better optimization of the outer objective through the approximations , which converge uniformly to by Proposition 3.2. This, however, does not necessary results in better test scores due to possible overfitting of the training and validation errors, as we noted for the linear hyperrepresentation multiclass classification experiment in Sec. 5.1. The optimal choice of appears to be, indeed, problem dependent: if, in the aforementioned experiment, a quite small
led to the best test accuracy (after hyperparameter optimization), in this section we present a small scale linear regression experiment that benefits from an increasing number of inner iterations.
We generated 90 noisy synthetic data points with 30 features, of which only 5 were informative, and divided them equally to form training, validation and test sets. As outer objective we set the mean squared error , where and are the validation design matrix and targets respectively and is the weight vector of the model. We set as inner inner objective
and optimized the vector of regularization coefficients (equivalent to a diagonal Tikhonov regularization matrix). The results, reported in Table 5, show that in this scenario overfitting is not an issue and 250 inner iterations yield the best test result. Note that, as in Sec. 5.1, also this problem admits an analytical solution, from which it is possible to compute the hypergradient analytically and perform exact hyperparameter optimization, as reported in the last row of the Table 5.
Validation MAPE  Test MAPE  
Exact 
Appendix D Further Details on Fewshot Experiments
This appendix contains implementation details of the representation mapping and the metalearning hyperparameters used in the fewshot learning experiments of Sec. 5.2.
To optimize the representation mapping, in all the experiments, we use Adam with learning rate set to and a decayrate of . We used the initialization strategy proposed in (Glorot and Bengio, 2010) for all the weights of the representation.
For Omniglot experiments we used a metabatch size of 32 episodes for fiveway and of 16 episodes for 20way. To train the episodespecific classifiers we set the learning rate to 0.1.
For one set of experiments with Miniimagenet we used an hyperrepresentation (C4L
) consisting of 4 convolutional layers where each layer is composed by a convolution with 32 filters, a batch normalization followed by a ReLU activation and a 2x2 maxpooling. The classifiers were trained using minibatches of 4 episodes for oneshot and 2 episodes for fiveshot with learning rate set to 0.01.
The other set of experiments with Miniimagenet employed a Residual Network (RN) as the representation mapping, built of 4 residual blocks (with 64, 96, 128, 256 filters) and then the block that follows conv (2048 filters), avg pooling, conv (512 filters) . Each residual block repeats the following block 3 times conv, batch normalization, leaky ReLU (leak 0.1)
before the residual connection. In this case the classifiers were optimized using minibatches of 2 episodes for both one and fiveshot with learning rate set to 0.04.
Comments
There are no comments yet.