Meta-Learning with Latent Embedding Optimization

Andrei A. Rusu, Dushyant Rao, Jakub Sygnowski, Oriol Vinyals, Razvan Pascanu, Simon Osindero, Raia Hadsell

Introduction

Humans have a remarkable ability to quickly grasp new concepts from a very small number of examples or a limited amount of experience, leveraging prior knowledge and context. In contrast, traditional deep learning approaches (LeCun et al., 2015; Schmidhuber, 2015) treat each task independently and hence are often data inefficient – despite providing significant performance improvements across the board, such as for image classification (Simonyan & Zisserman, 2014; He et al., 2016), reinforcement learning (Mnih et al., 2015; Silver et al., 2017), and machine translation (Cho et al., 2014; Sutskever et al., 2014). Just as humans can efficiently learn new tasks, it is desirable for learning algorithms to quickly adapt to and incorporate new and unseen information.

Few-shot learning tasks challenge models to learn a new concept or behaviour with very few examples or limited experience (Fei-Fei et al., 2006; Lake et al., 2011). One approach to address this class of problems is meta-learning, a broad family of techniques focused on learning how to learn or to quickly adapt to new information. More specifically, optimization-based meta-learning approaches (Ravi & Larochelle, 2017; Finn et al., 2017) aim to find a single set of model parameters that can be adapted with a few steps of gradient descent to individual tasks. However, using only a few samples (typically 1 or 5) to compute gradients in a high-dimensional parameter space could make generalization difficult, especially under the constraint of a shared starting point for task-specific adaptation.

In this work we propose a new approach, named Latent Embedding Optimization (LEO), which learns a low-dimensional latent embedding of model parameters and performs optimization-based meta-learning in this space. Intuitively, the approach provides two advantages. First, the initial parameters for a new task are conditioned on the training data, which enables a task-specific starting point for adaptation. By incorporating a relation network into the encoder, this initialization can better consider the joint relationship between all of the input data. Second, by optimizing in the lower-dimensional latent space, the approach can adapt the behaviour of the model more effectively. Further, by allowing this process to be stochastic, the ambiguities present in the few-shot data regime can be expressed.

We demonstrate that LEO achieves state-of-the-art results on both the miniImageNet and tieredImageNet datasets, and run an ablation study and further analysis to show that both conditional parameter generation and optimization in latent space are critical for the success of the method. Source code for our experiments is available at https://github.com/deepmind/leo.

Model

We define the NN-way KK-shot problem using the episodic formulation of Vinyals et al. (2016). Each task instance Ti\mathcal{T}_{i} is a classification problem sampled from a task distribution p(T)p(\mathcal{T}). The tasks are divided into a training meta-set Str\mathcal{S}^{{tr}}, validation meta-set Sval\mathcal{S}^{{val}}, and test meta-set Stest\mathcal{S}^{{test}}, each with a disjoint set of target classes (i.e., a class seen during testing is not seen during training). The validation meta-set is used for model selection, and the testing meta-set is used only for final evaluation.

Each task instance Ti∼p(T)\mathcal{T}_{i}\sim p\left(\mathcal{T}\right) is composed of a training set Dtr\mathcal{D}^{{tr}} and validation set Dval\mathcal{D}^{{val}}, and only contains NN classes randomly selected from the appropriate meta-set (e.g. for a task instance in the training meta-set, the classes are a subset of those available in Str\mathcal{S}^{{tr}}). In most setups, the training set \mathcal{D}^{{tr}}=\big{\{}({\mathbf{x}}_{n}^{k},{y}_{n}^{k})\mid k=1\dots K;n=1\dots N\big{\}} contains KK samples for each class. The validation set Dval\mathcal{D}^{{val}} can contain several other samples from the same classes, providing an estimate of generalization performance on the NN classes for this problem instance. We note that the validation set of a problem instance Dval\mathcal{D}^{{val}} (used to optimize a meta-learning objective) should not be confused with the held-out validation meta-set Sval\mathcal{S}^{{val}} (used for model selection).

2 Model-Agnostic Meta-Learning

Model-agnostic meta-learning (MAML) (Finn et al., 2017) is an approach to optimization-based meta-learning that is related to our work. For some parametric model fθf_{\theta}, MAML aims to find a single set of parameters θ\theta which, using a few optimization steps, can be successfully adapted to any novel task sampled from the same distribution. For a particular task instance Ti=(Dtr,Dval)\mathcal{T}_{i}=\left(\mathcal{D}^{{tr}},\mathcal{D}^{{val}}\right), the parameters are adapted to task-specific model parameters θi′\theta_{i}^{\prime} by applying some differentiable function, typically an update rule of the form:

where G\mathcal{G} is typically implemented as a step of gradient descent on the few-shot training set Dtr\mathcal{D}^{{tr}}, θi′=θ−α∇θLTitr(fθ)\theta^{\prime}_{i}=\theta-\alpha\nabla_{\theta}\mathcal{L}_{{\mathcal{T}_{i}}}^{{tr}}\left(f_{\theta}\right). Generally, multiple sequential adaptation steps can be applied. The learning rate α\alpha can also be meta-learned concurrently, in which case we refer to this algorithm as Meta-SGD (Li et al., 2017). During meta-training, the parameters θ\theta are updated by back-propagating through the adaptation procedure, in order to reduce errors on the validation set Dval\mathcal{D}^{{val}}:

The approach includes the main ingredients of optimization-based meta-learning with neural networks: initialization is done by maintaining an explicit set of model parameters θ\theta; the adaptation procedure, or “inner loop”, takes θ\theta as input and returns θi′\theta_{i}^{\prime} adapted specifically for task instance Ti{\mathcal{T}_{i}}, by iteratively using gradient descent (Eq. 1); and termination, which is handled simply by choosing a fixed number of optimization steps in the “inner loop”. MAML updates θ\theta by differentiating through the “inner loop” in order to minimize errors of instance-specific adapted models fθi′f_{\theta_{i}^{\prime}} on the corresponding validation set (Eq. 2). We refer to this process as the “outer loop” of meta-learning. In the next section we use the same stages to describe Latent Embedding Optimization (LEO).

3 Latent Embedding Optimization for Meta-Learning

The primary contribution of this paper is to show that it is possible, and indeed beneficial, to decouple optimization-based meta-learning techniques from the high-dimensional space of model parameters. We achieve this by learning a stochastic latent space with an information bottleneck, conditioned on the input data, from which the high-dimensional parameters are generated.

Instead of explicitly instantiating and maintaining a unique set of model parameters θ\theta, as in MAML, we learn a generative distribution of model parameters which serves the same purpose. This is a natural extension: we relax the requirement of finding a single optimal θ∗∈Θ\theta^{\ast}\in\Theta to that of approximating a data-dependent conditional probability distribution over Θ\Theta, which can be more expressive. The choice of architecture, composed of an encoding process, and decoding (or parameter generation) process, enables us to perform the MAML gradient-based adaptation steps (or “inner loop”) in the learned, low-dimensional embedding space of the parameter generative model (Figure 1).

Figure 2 shows the architecture of the resulting network. Intuitively, the decoder is akin to a generative model, mapping from a low-dimensional latent code to a distribution over model parameters. The encoding process ensures that the initial latent code and parameters before gradient-based adaptation are already data-dependent. This encoding process also exploits a relation network that allows the latent code to be context-dependent, considering the pairwise relationship between all classes in the problem instance. In the following sections, we explain the LEO procedure more formally.

3.2 Initialization: Generating Parameters Conditioned on a Few Examples

The first stage is to instantiate the model parameters that will be adapted to each task instance. Whereas MAML explicitly maintains a single set of model parameters, LEO utilises a data-dependent latent encoding which is then decoded to generate the actual initial parameters. In what follows, we describe an encoding scheme which leverages a relation network to map the few-shot examples into a single latent vector. This design choice allows the approach to consider context when producing a parameter initialization. Intuitively, decision boundaries required for fine-grained distinctions between similar classes might need to be different from those for broader classification.

The encoding process involves a simple feed-forward mapping of each data point, followed by a relation network that considers the pair-wise relationship between the data in the problem instance. The overall encoding process is defined in Eq. 3, and proceeds as follows. First, each example from a problem instance Ti=(Dtr,Dval)∼p(T)\mathcal{T}_{i}=\left(\mathcal{D}^{{tr}},\mathcal{D}^{{val}}\right)\sim p\left(\mathcal{T}\right) is processed by an encoder network gϕe ⁣:Rnx→Rnhg_{\phi_{{e}}}\colon\mathcal{R}^{n_{x}}\to\mathcal{R}^{n_{h}}, which maps from input space to a code in an intermediate hidden-layer code space H\mathcal{H}. Then, codes in H\mathcal{H} corresponding to different training examples are concatenated pair-wise (resulting in (NK)2(NK)^{2} pairs in the case of KK-shot classification) and processed by a relation network gϕrg_{\phi_{{r}}}, in a similar fashion to Oreshkin et al. (2018) and Sung et al. (2017). The (NK)2(NK)^{2} outputs are grouped by class and averaged within each group to obtain the (2×N2\times N) parameters of a probability distribution in a low-dimensional space Z=Rnz\mathcal{Z}=\mathcal{R}^{n_{z}}, where nz≪dim(θ)n_{z}\ll\textrm{dim}(\theta), for each of the NN classes.

Thus, given the KK-shot training samples corresponding to a class nn: \mathcal{D}_{n}^{{tr}}=\big{\{}\big{(}{\mathbf{x}}_{n}^{k},{y}_{n}^{k}\big{)}\mid k=1\dots K\big{\}} the encoder gϕeg_{\phi_{{e}}} and relation network gϕrg_{\phi_{{r}}} together parameterize a class-conditional multivariate Gaussian distribution with a diagonal covariance, which we can sample from in order to output a class-dependent latent code zn∈Z{\mathbf{z}}_{n}\in\mathcal{Z} as follows:

Intuitively, the encoder and relation network define a stochastic mapping from one or more class examples to a single code in the latent embedding space Z\mathcal{Z} corresponding to that class. The final latent code can be obtained as the concatenation of class-dependent codes: z=[z1,z2,…,zN]{\mathbf{z}}=[{\mathbf{z}}_{1},{\mathbf{z}}_{2},\dots,{\mathbf{z}}_{N}].

In other words, codes zn{\mathbf{z}}_{n} are mapped independently to the top-layer parameters θi\theta_{i} of a softmax classifier using the decoder gϕdg_{\phi_{{d}}}, which is essentially a stochastic generator of model parameters.

3.3 Adaptation by Latent Embedding Optimization (LEO) (The “Inner Loop”)

Given the decoded parameters, we can then define the “inner loop” classification loss using the cross-entropy function, as follows:

It is important to note that the decoder gϕdg_{\phi_{{d}}} is a differentiable mapping between the latent space Z\mathcal{Z} and the higher-dimensional model parameter space Θ\Theta. Primarily, this allows gradient-based optimization of the latent codes with respect to the training loss, with zn′=zn−α∇znLTitr{\mathbf{z}}^{\prime}_{n}={\mathbf{z}}_{n}-\alpha\nabla_{{\mathbf{z}}_{n}}\mathcal{L}^{{tr}}_{{\mathcal{T}_{i}}}. The decoder gϕdg_{\phi_{{d}}} will convert adapted latent codes zn′{\mathbf{z}}^{\prime}_{n} to effective model parameters θi′\theta^{\prime}_{i} for each adaptation step, which can be repeated several times, as in Algorithm 1. In addition, by backpropagating errors through the decoder, the encoder and relation net can learn to provide a data-conditioned latent encoding z{\mathbf{z}} that produces an appropriate initialization point θi\theta_{i} for the classifier model.

3.4 Meta-Training Strategy (The “Outer Loop”)

For each task instance Ti\mathcal{T}_{i}, the initialization and adaptation procedure produce a new classifier fθi′f_{\theta_{i}^{\prime}} tailored to the training set Dtr\mathcal{D}^{{tr}} of the instance, which we can then evaluate on the validation set of that instance Dval\mathcal{D}^{{val}}. During meta-training we use that evaluation to differentiate through the “inner loop” and update the encoder, relation, and decoder network parameters: ϕe\phi_{{e}}, ϕr\phi_{{r}}, and ϕd\phi_{{d}}. Meta-training is performed by minimizing the following objective:

where p(zn)=N(0,I)p({\mathbf{z}}_{n})=\mathcal{N}(0,\mathcal{I}). Similar to the loss defined in (Higgins et al., 2017) we use a weighted KL-divergence term to regularize the latent space and encourage the generative model to learn a disentangled embedding, which should also simplify the LEO “inner loop” by removing correlations between latent space gradient dimensions. The third term in Eq. (6) encourages the encoder and relation net to output a parameter initialization that is close to the adapted code, thereby reducing the load of the adaptation procedure if possible.

L2L_{2} regularization was used with all weights of the model, as well as a soft, layer-wise orthogonality constraint on decoder network weights, which encourages the dimensions of the latent code as well as the decoder network to be maximally expressive. In the case of linear encoder, relation, and decoder networks, and assuming that Cd\mathcal{C}_{{d}} is the correlation matrix between rows of ϕd\phi_{{d}}, then the regularization term takes the following form:

3.5 Beyond Classification and Linear Output Layers

Thus far we have used few-shot classification as a working example to highlight our proposed method, and in this domain we generate only a single linear output layer. However, our approach can be applied to any model fθif_{\theta_{i}} which maps observations to outputs, e.g. a nonlinear MLP or LSTM, by using a single latent code z{\mathbf{z}} to generate the entire parameter vector θi\theta_{i} with an appropriate decoder. In the general case, z{\mathbf{z}} is conditioned on Dtr\mathcal{D}^{{tr}} by passing both inputs and labels to the encoder. Furthermore, the loss LTi\mathcal{L}_{{\mathcal{T}_{i}}} is not restricted to be a classification loss, and can be replaced by any differentiable loss function which can be computed on Dtr\mathcal{D}^{{tr}} and Dval\mathcal{D}^{{val}} sets of a task instance Ti{\mathcal{T}_{i}}.

Related Work

The problem of few-shot adaptation has been approached in the context of fast weights (Hinton & Plaut, 1987; Ba et al., 2016), learning-to-learn (Schmidhuber, 1987; Thrun & Pratt, 1998; Hochreiter et al., 2001; Andrychowicz et al., 2016), and through meta-learning. Many recent approaches to meta-learning can be broadly categorized as metric-based methods, which focus on learning similarity metrics for members of the same class (e.g. Koch et al., 2015; Vinyals et al., 2016; Snell et al., 2017); memory-based methods, which exploit memory architectures to store key training examples or directly encode fast adaptation algorithms (e.g. Santoro et al., 2016; Ravi & Larochelle, 2017); and optimization-based methods, which search for parameters that are conducive to fast gradient-based adaptation to new tasks (e.g. Finn et al., 2017; 2018).

Related work has also explored the use of one neural network to produce (some fraction of) the parameters of another (Ha et al., 2016; Krueger et al., 2017), with some approaches focusing on the goal of fast adaptation. Munkhdalai et al. (2017) meta-learn an algorithm to change additive biases across deep networks conditioned on the few-shot training samples. In contrast, Gidaris & Komodakis (2018) use an attention kernel to output class conditional mixing of linear output weights for novel categories, starting from a pre-trained deep model. Qiao et al. (2017) learn to output top linear layer parameters from the activations provided by a pre-trained feature embedding, but they do not make use of gradient-based adaptation. None of the aforementioned approaches to fast adaptation explicitly learn a probability distribution over model parameters, or make use of latent variable generative models to characterize it.

Approaches which use optimization-based meta-learning include MAML (Finn et al., 2017) and REPTILE (Nichol & Schulman, 2018). While MAML backpropagates the meta-loss through the “inner loop”, REPTILE simplifies the computation by incorporating an L2L_{2} loss which updates the meta-model parameters towards the instance-specific adapted models. These approaches use the full, high-dimensional set of model parameters within the “inner loop”, while Lee & Choi (2018) learn a layer-wise subspace in which to use gradient-based adaptation. However, it is not clear how these methods scale to large expressive models such as residual networks (especially given the uncertainty in the few-shot data regime), since MAML is prone to overfitting (Mishra et al., 2018). Recognizing this issue, Zhou et al. (2018) train a deep input representation, or “concept space”, and use it as input to an MLP meta-learner, but perform gradient-based adaptation directly in its parameter space, which is still comparatively high-dimensional. As we will show, performing adaptation in latent space to generate a simple linear layer can lead to superior generalization.

Probabilistic meta-learning approaches such as those of Bauer et al. (2017) and Grant et al. (2018) have shown the advantages of learning Gaussian posteriors over model parameters. Concurrently with our work, Kim et al. (2018) and Finn et al. (2018) propose probabilistic extensions to MAML that are trained using a variational approximation, using simple posteriors. However, it is not immediately clear how to extend them to more complex distributions with a more diverse set of tasks. Other concurrent works have introduced deep parameter generators (Lacoste et al., 2018; Wu et al., 2018) that can better capture a wider distribution of model parameters, but do not employ gradient-based adaptation. In contrast, our approach employs both a generative model of parameters, and adaptation in a low-dimensional latent space, aided by a data-dependent initialization.

Finally, recently proposed Neural Processes (Garnelo et al., 2018a; b) bear similarity to our work: they also learn a mapping to and from a latent space that can be used for few-shot function estimation. However, coming from a Gaussian processes perspective, their work does not perform “inner loop” adaptation and is trained by optimizing a variational objective.

Evaluation

We evaluate the proposed approach on few-shot regression and classification tasks. This evaluation aims to answer the following key questions: (1) Is LEO capable of modeling a distribution over model parameters when faced with uncertainty? (2) Can LEO learn from multimodal task distributions and is this reflected in ambiguous problem instances, where multiple distinct solutions are possible? (3) Is LEO competitive on large-scale few-shot learning benchmarks?

To answer the first two questions we adopt the simple regression task of Finn et al. (2018). 1D regression problems are generated in equal proportions using either a sine wave with random amplitude and phase, or a line with random slope and intercept. Inputs are sampled randomly, creating a multimodal task distribution. Crucially, random Gaussian noise with standard deviation 0.30.3 is added to regression targets. Coupled with the small number of training samples (5-shot), the task is challenging for 2 main reasons: (1) learning a distribution over models becomes necessary, in order to account for the uncertainty introduced by noisy labels; (2) problem instances may be likely under both modes: in some cases a sine wave may fit the data as well as a line. Faced with such ambiguity, learning a generative distribution of model parameters should allow several different likely models to be sampled, in a similar way to how generative models such as VAEs can capture different modes of a multimodal data distribution.

We used a 3-layer MLP as the underlying model architecture of fθf_{\theta}, and we produced the entire parameter tensor θ\theta with the LEO generator, conditionally on DtrD^{{tr}}, the few-shot training inputs concatenated with noisy labels. For further details, see Appendix A.

In Figure 3 we show samples from a single model trained on noisy sines and lines, with true regression targets in black and training samples marked with red circles and vertical dashed lines. Plots (a) and (b) illustrate how LEO captures some of the uncertainty in ambiguous problem instances within each mode, especially in parts of the input space far from any training samples. Conversely, in parts which contain data, models fit the regression target well. Interestingly, when both sines and lines could explain the data, as shown in panels (c) and (d), we see that LEO can sample very different models, from both families, reflecting its ability to represent parametric uncertainty appropriately.

2 Few-shot Classification

In order to answer the final question we scale up our approach to 1-shot and 5-shot classification problems defined using two commonly used ImageNet subsets.

The miniImageNet dataset (Vinyals et al., 2016) is a subset of 100100 classes selected randomly from the ILSVRC-12 dataset (Russakovsky et al., 2014) with 600600 images sampled from each class. Following the split proposed by Ravi & Larochelle (2017), the dataset is divided into training, validation, and test meta-sets, with 6464, 1616, and 2020 classes respectively.

The tieredImageNet dataset (Ren et al., 2018) is a larger subset of ILSVRC-12 with 608 classes (779,165779{,}165 images) grouped into 3434 higher-level nodes in the ImageNet human-curated hierarchy (Deng et al., 2009a). This set of nodes is partitioned into 2020, 66, and 88 disjoint sets of training, validation, and testing nodes, and the corresponding classes form the respective meta-sets. As argued in Ren et al. (2018), this split near the root of the ImageNet hierarchy results in a more challenging, yet realistic regime with test classes that are less similar to training classes.

2.2 Pre-trained Features

Two potential difficulties of using LEO to instantiate parameters with a generator network are: (1) modeling distributions over very high-dimensional parameter spaces; and (2) requiring meta-learning (and hence, gradient computation in the inner loop) to be performed with respect to a high-dimensional input space. We address these issues by pre-training a visual representation of the data and then using the generator to instantiate the parameters for the final layer - a linear softmax classifier operating on this representation. We train a 2828-layer Wide Residual Network (WRN-28-10) (Zagoruyko & Komodakis, 2016a) with supervised classification using only data and classes from the training meta-set. Recent state-of-the-art approaches use the penultimate layer representation (Zhou et al., 2018; Qiao et al., 2017; Bauer et al., 2017; Gidaris & Komodakis, 2018); however, we choose the intermediate feature representation in layer 2121, given that higher layers tend to specialize to the training distribution (Yosinski et al., 2014). For details regarding the training, evaluation, and network architectures, see Appendix B.

2.3 Fine-tuning

Following the LEO adaptation procedure (Algorithm 1) we also use fine-tuningIn this context, “fine-tuning” refers to final adaptation in parameter space, rather than fine-tuning the pre-trained feature extractor. by performing a few steps of gradient-based adaptation directly in parameter space using the few-shot set Dtr\mathcal{D}^{{tr}}. This is similar to the adaptation procedure of MAML, or Meta-SGD (Li et al., 2017) when the learning rates are learned, with the important difference that starting points of fine-tuning are custom generated by LEO for every task instance Ti{\mathcal{T}_{i}}. Empirically, we find that fine-tuning applies a very small change to the parameters with only a slight improvement in performance on supervised classification tasks.

3 Results

The classification accuracies for LEO and other baselines are shown in Table 1. LEO sets the new state-of-the-art performance on the 1-shot and 5-shot tasks for both miniImageNet and tieredImageNet datasets. We also evaluated LEO on the “multi-view” feature representation used by Qiao et al. (2017) with miniImageNet, which involves significant data augmentation compared to the approaches in Table 1. LEO is state-of-the-art using these features as well, with 63.97±0.20%63.97\pm 0.20\% and 79.49±0.70%79.49\pm 0.70\% test accuracies on the 1-shot and 5-shot tasks respectively.

4 Ablation Study

To assess the effects of different components, we also performed an ablation study, with detailed results in Table 2. To ensure a fair comparison, all approaches begin with the same pre-trained features (Section 4.2.2). The Meta-SGD case performs gradient-based adaption directly in the parameter space in the same way as MAML, but also meta-learns the inner loop learning rate (as we do for LEO). The main approach, labeled as LEO in the table, uses a stochastic parameter generator for several steps of latent embedding optimization, followed by fine-tuning steps in parameter space (see subsection 4.2.3). All versions of LEO are at or above the previous state-of-the-art on all tasks.

The largest difference in performance is between Meta-SGD and the other cases (all of which exploit a latent representation of model parameters), indicating that the low-dimensional bottleneck is critical for this application. The “conditional generator only” case (without adaptation in latent space) yields a poorer result than LEO, and even adding fine-tuning in parameter space does not recover performance; this illustrates the efficacy of the latent adaptation procedure. The importance of the data-dependent encoding is highlighted by the “random prior” case, in which the encoding process is replaced by the prior p(zn)p({\mathbf{z}}_{n}), and performance decreases. We also find that incorporating stochasticity can be important for miniImageNet, but not for tieredImageNet, which we hypothesize is because the latter is much larger. Finally, the fine-tuning steps only yield a statistically significant improvement on the 5-shot tieredImageNet task. Thus, both the data-conditional encoding and latent space adaptation are critical to the performance of LEO.

5 Latent Embedding Visualization

To qualitatively characterize the learnt embedding space, we plot codes produced by the relational encoder before and after the LEO procedure, using a 5-way 1-shot model and 10001000 task instances from the validation meta-set of miniImageNet. Figure 4 shows a t-SNE projection of class conditional encoder outputs zn{\mathbf{z}}_{n} as well as their respective final adapted versions zn′{\mathbf{z}}^{\prime}_{n}. If the effect of LEO were minimal, we would expect latent codes to have roughly the same structure before and after adaptation. In contrast, Figure 4(a) clearly shows that latent codes change substantially during LEO, since encoder output codes form a large cluster (blue) to which adapted codes (orange) do not belong. Figure 4(b) shows the same t-SNE embedding as (a) colored by class label. Note that encoder outputs, on the left side of plot (b), have a lower degree of class conditional separation compared to zn′{\mathbf{z}}^{\prime}_{n} clusters on the right, suggesting that qualitatively different structure is introduced by the LEO procedure. We further illustrate this point by highlighting latent codes for the “Jellyfish” validation class in Figure 4(c), which are substantially different before and after adaptation.

The additional structure of adapted codes zn′{\mathbf{z}}^{\prime}_{n} may explain LEO’s superior performance over approaches predicting parameters directly from inputs, since the decoder may not be able to produce sufficiently different weights for different classes given very similar latent codes, especially when the decoder is linear. Conversely, LEO can reduce the uncertainty of the encoder mapping, which is inherent in the few-shot regime, by adapting latent codes with a generic, gradient-based procedure.

6 Curvature and coverage analysis

We hypothesize that by performing the inner-loop optimization in a lower-dimensional latent space, the adapted solutions do not need to be close together in parameter space, as each latent step can cover a larger region of parameter space and effect a greater change on the underlying function. To support this intuition, we compute a number of curvature and coverage measures, shown in Figure 5.

The curvature provides a measure of the sensitivity of a function with respect to some space. If adapting in latent space allows as much control over the function as in parameter space, one would expect similar curvatures. However, as demonstrated in Figure 5(a), the curvature for LEO in z{\mathbf{z}} space (the absolute eigenvalues of the Hessian of the loss) is 2 orders of magnitude higher than in θ\theta, indicating that a fixed step in z{\mathbf{z}} will change the function more drastically than taking the same step directly in θ\theta. This is also observed in the “gen+ft” case, where the latent embedding is still used, but adaptation is performed directly in θ\theta space. This suggests that the latent bottleneck is responsible for this effect. Figure 5(b) shows that this is due to the expansion of space caused by the decoder. In this case the decoder is linear, and the singular values describe how much a vector projected through this decoder grows along different directions, with a value of one preserving volume. We observe that the decoder is expanding the space by at least one order of magnitude. Finally, Figure 5(c) demonstrates this effect along the specific gradient directions used in the inner loop adaptation: the small gradient steps in z{\mathbf{z}} taken by LEO induce much larger steps in θ\theta space, larger than the gradient steps taken by Meta-SGD in θ\theta space directly. Thus, the results support the intuition that LEO is able to ‘transport’ models further during adaptation by performing meta-learning in the latent space.

Conclusions and Future Work

We have introduced Latent Embedding Optimization (LEO), a meta-learning technique which uses a parameter generative model to capture the diverse range of parameters useful for a distribution over tasks, and demonstrated a new state-of-the-art result on the challenging 5-way 1- and 5-shot miniImageNet and tieredImageNet classification problems. LEO achieves this by learning a low-dimensional data-dependent latent embedding, and performing gradient-based adaptation in this space, which means that it allows for a task-specific parameter initialization and can perform adaptation more effectively.

Future work could focus on replacing the pre-trained feature extractor with one learned jointly through meta-learning, or using LEO for tasks in reinforcement learning or with sequential data.

References

Appendix A Experimental setup - regression

We used the experimental setup of Finn et al. (2018) for 1D 5-shot noisy regression tasks. Inputs were sampled uniformly from $.Amultimodaltaskdistributionwasused.Halfoftheprobleminstancesweresinusoidswithamplitudeandphasesampleduniformlyfrom. A multimodal task distribution was used. Half of the problem instances were sinusoids with amplitude and phase sampled uniformly from[0.1,5]andand[0,\pi]respectively.Theotherhalfwerelineswithslopeandinterceptsampleduniformlyfromtheintervalrespectively. The other half were lines with slope and intercept sampled uniformly from the interval.Gaussiannoisewithstandarddeviation. Gaussian noise with standard deviation0.3$ was added to regression targets.

A.2 LEO Network Architecture

As Table 3 shows, the underlying model fθf_{\theta} (for which parameters θ\theta were generated) was a 3-layer MLP with 4040 units in all hidden layers and rectifier nonlinearities. A single code z{\mathbf{z}} was used to generate θ\theta with the decoder, conditioned on concatenated inputs and regression targets from Dtr\mathcal{D}^{{tr}} which were passed as inputs to the encoder. Sampling of latent codes and parameters was used both during training and evaluation.

The encoder was a 3-layer MLP with 3232 units per layer and rectifier nonlinearities; the bottleneck embedding space size was: nz=16n_{z}=16. The relation network and decoder were both 33-layer MLPs with 3232 units per layer. For simplicity we did not use biases in any layer of the encoder, decoder nor the relation network. Note that the last dimension of the relation network and decoder outputs are two times larger than nzn_{z} and dim(θ)\textrm{dim}(\theta) respectively, as they are used to parameterize both the means and variances of the corresponding Gaussian distributions.

Appendix B Experimental setup - classification

We used the standard 5-way 1-shot and 5-shot classification setups, where each task instance involves classifying images from 55 different categories sampled randomly from one of the meta-sets, and Dtr\mathcal{D}^{{tr}} contains 11 or 55 training examples respectively. Dval\mathcal{D}^{{val}} contains 15 samples during meta-training, as decribed in Finn et al. (2017), and all the remaining examples during validation and testing, following Qiao et al. (2017).

We did not employ any data augmentation or feature averaging during meta-learning, or any other data apart from the corresponding training and validation meta-sets. The only exception is the special case of “multi-view” embedding results, where features were averaged over representations of 4 corner and central crops and their horizontal mirrored versions, which we provide for full comparison with Qiao et al. (2017). Apart from the differences described here, the feature training pipeline closely followed that of Qiao et al. (2017).

B.2 Feature Pre-Training

As described in Section 4.2.2, we trained dataset specific feature embeddings before meta-learning, in a similar fashion to Qiao et al. (2017) and Bauer et al. (2017). A Wide Residual Network WRN-28-10 (Zagoruyko & Komodakis, 2016b) with 33 steps of dimensionality reduction was used to classify images of 80×8080\times 80 pixels from only the meta-training set into the corresponding training classes (6464 in case of miniImageNet and 351351 for tieredImageNet). We used dropout (pkeep=0.5)(p_{keep}=0.5) inside residual blocks, as described in (Zagoruyko & Komodakis, 2016b), which is turned off during evaluation and for dataset export. An L2 regularization term of 5e−45e^{-4} was used, 0.90.9 Nesterov momentum, and SGD with a learning rate schedule. The initial learning rate was 0.10.1 and it was multiplied with 0.2 at the steps given in Table 4. Mini-batches were of size of 10241024. Data augmentation for pre-training was similar to the inception pipeline (Szegedy et al., ), with color distortions and image deformations and scaling in training mode. For 6464-way evaluation accuracy and dataset export we used only the center crop (with a ratio of 8092\frac{80}{92}: about 85.95%85.95\% of the image) which was then resized to 80×8080\times 80 and passed to the network.

Activations in layer 2121, with average pooling over spatial dimensions, were precomputed and saved as feature embeddings with nx=dim(x)=640n_{x}=\textrm{dim}({\mathbf{x}})=640, which substantially simplified the meta-learning process.

B.3 LEO Network Architecture

We used the same network architecture of parameter generator for all datasets and tasks. The encoder and decoder networks were linear with the bottleneck embedding space of size nz=64n_{z}=64. The relation network was a 33-layer fully connected network with 128128 units per layer and rectifier nonlinearities. For simplicity we did not use biases in any layer of the encoder, decoder nor the relation network. Table 5 summarizes this information. Note that the last dimension of the relation network and decoder outputs are two times larger than nzn_{z} and dim(x)\textrm{dim}({\mathbf{x}}) respectively, as they are used to parameterize both the means and variances of the corresponding Gaussian distributions.

The “Meta-SGD (our features)” baseline used the same one-layer softmax classifier as base model.

B.4 Optimization

We used a parallel implementation similar to that of Finn et al. (2017), where the “inner loop” is performed in parallel on a batch 1212 problem instances for every meta-update. Using a relation network in the encoder has negligible computational cost given that k2k^{2} is small in typical kk-shot learning domains, and the relation network is only used once per problem instance, to get the initial model parameters before adaptation. Within the LEO “inner loop” we perform 55 steps of adaptation in latent space, followed by 55 steps of fine-tuning in parameter space. The learning rates for these spaces were meta-learned in a similar fashion to Meta-SGD (Li et al., 2017), after being initialized to 11 and 0.0010.001 for the latent and parameter spaces respectively. We applied dropout independently on the feature embedding in every step, with the probability of not being dropped out pkeepp_{keep} chosen (together with other hyperparameters) using random search based on the validation meta-set accuracy.

Parameters of the encoder, relation, and decoder networks as well as per-parameter learning rates in latent and parameter spaces were optimized jointly using Adam (Kingma & Ba, 2014) to minimize the meta-learning objective (Eq. 6) over problem instances from the training meta-set, iterating for up to 100 000100\,000 steps, with early stopping using validation accuracy.

Meta-learning objectives can lead to difficult optimization processes in practice, specifically when coupled with stochastic sampling in latent and parameters spaces. For ease of experimentation we clip the meta-gradient, as well as its norm, at an absolute value of 0.10.1. Please note this was only done for the encoder, relation, decoder networks and learning rates, not the inner loop latent space adaptation gradients.

B.5 Hyper-parameters

To find the best values of hyperparameters, we performed a random grid search and we choose the set which lead to highest validation meta-set accuracy. The reported performance of our models is an average (±\pm a standard deviation) over 55 independent runs (using different random seeds) with the best hyperparameters kept fixed. The result of a single run is an average accuracy over 50000 task instances. After choosing hyperparameters (given in Table 6) we used both meta-training and meta-validation sets for training, in line with recent state-of-the-art approaches, e.g. Qiao et al. (2017).

The evaluation of each of the LEO baselines follow the same procedure; in particular, we perform a separate random search for each of them.

B.6 Training time

Training of LEO took 1-2 hours for miniImageNet and around 5 hours for tieredImageNet on a multi-core CPU (for each of the 55 independent runs). Our approach allows for caching the feature embeddings before training LEO, which leads to a very efficient meta-learning process.

Training of the image extractor was more compute-intensive, taking 5 hours for miniImageNet and around a day for tieredImageNet using 3232 GPUs.

B.7 Overview of the training procedure

In summary, there are three stages in our approach to meta-training:

In the first stage we use 64-way classification to pre-train the feature embedding only on the meta-training set, hence without the meta-validation classes.

In the second stage we train LEO on the meta-training set with early stopping on meta-validation, and we choose the best hyperparameters using random grid search.

In the third stage we train LEO again from scratch 55 times using the embedding trained in stage 11 and the chosen set of hyperparameters from stage 22. However, in this stage we meta-learn on embeddings from both meta-train and meta-validation sets, with early-stopping on meta-validation.

While it may not be intuitive to use early stopping on meta-validation in stage 33, it is still a proxy for good generalization since it favors models with high performance on classes excluded during feature embedding pre-training.

B.8 Overview of the evaluation procedure

The procedure for evaluation is similar to meta-training, except that we disable stochasticity and dropout. Naturally, instead of computing the meta-training loss, the parameters (adapted based on LTitr\mathcal{L}^{{tr}}_{{\mathcal{T}_{i}}}) are only used for inference on that particular task. That is:

A problem instance is drawn from the evaluation meta-set.

The few-shot samples are encoded to latent space, then decoded; the means are used to initialize the parameters of the inference model.

A few steps of adaptation are performed in latent space, followed (optionally) by a few steps of adaptation in parameter space.

The resulting parameters are used as the final adapted model for that particular problem instance.