Wasserstein Task Embedding for Measuring Task Similarities

Xinran Liu¹, Yikun Bai¹, Yuzhe Lu¹, Andrea Soltoggio², Soheil Kolouri¹

Abstract

Measuring similarities between different tasks is critical in a broad spectrum of machine learning problems, including transfer, multi-task, continual, and meta-learning. Most current approaches to measuring task similarities are architecture-dependent: 1) relying on pre-trained models, or 2) training networks on tasks and using forward transfer as a proxy for task similarity. In this paper, we leverage the optimal transport theory and define a novel task embedding for supervised classification that is model-agnostic, training-free, and capable of handling (partially) disjoint label sets. In short, given a dataset with ground-truth labels, we perform a label embedding through multi-dimensional scaling and concatenate dataset samples with their corresponding label embeddings. Then, we define the distance between two datasets as the 2-Wasserstein distance between their updated samples. Lastly, we leverage the 2-Wasserstein embedding framework to embed tasks into a vector space in which the Euclidean distance between the embedded points approximates the proposed 2-Wasserstein distance between tasks. We show that the proposed embedding leads to a significantly faster comparison of tasks compared to related approaches like the Optimal Transport Dataset Distance (OTDD). Furthermore, we demonstrate the effectiveness of our proposed embedding through various numerical experiments and show statistically significant correlations between our proposed distance and the forward and backward transfer between tasks on a wide variety of image recognition datasets.

\affiliations

¹Department of Computer Science, Vanderbilt University
²Department of Computer Science, Loughborough University

Introduction

Learning from a broad spectrum of tasks and transferring knowledge between them is a cornerstone of intelligence, and primates perfectly exemplify this characteristic. Modern Machine Learning (ML) is rapidly moving toward multi-task learning, and there is great interest in methods that can integrate, rapidly adapt, and seamlessly transfer knowledge between tasks. When learning from multiple possibly heterogeneous tasks, it is essential to understand the relationships between the tasks and their fundamental properties. It is, therefore, highly desirable to define (dis)similarity measures between tasks that will allow one to cluster tasks, have better control over the forward and backward transfer, and ultimately require less supervision for learning tasks.

There has been an increasing interest in assessing task similarities and their relationship with forward and backward knowledge transfer among tasks. For instance, various recent works look into the selection of good source tasks/models for a given target task to maximize the forward transfer to the target task achille2019task2vec; zamir2018taskonomy; bao2019information; bhattacharjee2020p2l; fifty2021efficiently. Others have demonstrated the relationship between negative backward transfer (i.e., catastrophic forgetting) and task similarities nguyen2019toward.

Many existing methods for measuring task similarities depend on the choice of model(s), architecture(s), and the training process leite2005predicting; zamir2018taskonomy; achille2019task2vec; khodak2019adaptive; nguyen2020leep; venkitaraman2020task; gao2021information. For example, zamir2018taskonomy; venkitaraman2020task; gao2021information use pre-trained task specified models to measure a notion of forward transfer and define it as task similarity. achille2019task2vec embed tasks into a vector space that relies on a partially-trained network. khodak2019adaptive use the optimal parameters as a proxy for each task and leite2005predicting use the learning curves of a pre-specified model to measure task similarities. Besides being model-dependent, these approaches are often computationally expensive as they involve training deep models (or require pre-trained models).

Model-agnostic task similarity measures provide a fundamentally different approach to quantifying task relationships ben2006analysis; alvarez2020geometric; tran2019transferability; tan2021otce. These methods often measure the similarity between tasks as a function of the similarity between the joint or conditional input/output distributions, sometimes also taking the loss function into account. The classic theoretical results for such similarity measures ben2006analysis; batu2000testing focus on information theoretic divergences between the source and target distributions. More recently, Optimal Transport (OT) based approaches alvarez2020geometric; tan2021otce; xu2022selecting have shown promise in modeling task similarities. Notably, alvarez2020geometric approach measuring task similarities through the lens of a hierarchical OT yurochkin2019hierarchical where they solve an inner OT problem to calculate the label distance between the class-conditional distributions of two supervised learning tasks. The label distance is then incorporated into the transportation cost of an outer OT problem, resulting in a distance between two datasets that integrates both sample and label discrepancies. tan2021otce treats the optimal transport plan between the input distributions of two tasks as a joint probability distribution and use conditional entropy to measure the difference between the two tasks. One major shortcoming of these OT-based approaches is their computational complexity. These methods require the pairwise calculation of OT (or entropy regularized OT) between different tasks, which can be prohibitively expensive in applications requiring frequent evaluations of task similarities, e.g., in continual learning.

We propose a novel OT-based task embedding for supervised learning problems that is model-agnostic and computationally efficient. On the one hand, our proposed approach is similar to (achille2019task2vec) and (peng2020domain2vec), which embed datasets into a vector space in which one can easily measure the difference between tasks, e.g., via the Euclidean distance between embedded vectors. On the other hand, our approach is inspired by the Optimal Transport Dataset Distance (OTDD) alvarez2020geometric framework, and it essentially provides a Euclidean embedding for a hierarchical OT-based distance between tasks. To calculate such a task embedding, we use the Wasserstein embedding framework wang2013linear; kolouri2020wasserstein. Importantly, our approach alleviates the need for pairwise calculation of OT problems between tasks, turning it into a more desirable solution than previously proposed methods.

Contributions. We propose a computationally efficient and model-agnostic task embedding, denoted as Wasserstein Task Embedding (WTE), in which the Euclidean distance between embedded vectors approximates a hierarchical OT distance between the tasks. We provide extensive numerical experiments and demonstrate that: 1) our calculated distances between embedded tasks are highly correlated with the OTDD distance alvarez2020geometric, 2) our proposed embedding and similarity calculation is significantly faster than the OTDD distance, and 3) our proposed similarity measure provides strong and statistically significant correlation with both forward and backward transfer.

Related work

Model-based task similarity. Most existing approaches to measuring task similarity are model-dependent and use forward transferability as a proxy for similarity. For example, zamir2018taskonomy use pre-trained models on source tasks and measure their performance on a target task to obtain an asymmetric notion of similarity between source and target tasks. Following zamir2018taskonomy’s work, dwivedi2019representation measure the transferability in a more efficient manner by applying the Representation Similarity Analysis (RSA) between the trained models (e.g. DNNs) from different tasks. Similarly, nguyen2020leep assume the source and target tasks share the same set of inputs but have different sets of labels, and estimate the transferability by the empirical conditional distribution of target labels given the inputs computed by a pre-trained model on the source task.

Another class of approaches embed the tasks into a vector space and then define the (dis)similarity on the embedded vector representations. achille2019task2vec discuss processing data (images) through a partially trained “probe network” and obtain vector embedding by computing the Fisher information matrix (FIM). The (dis)similarity of two tasks is then computed from the difference between the the FIMs. Similarly, peng2020domain2vec propose a domain (labeled dataset) to vector technique. In particular, given a domain, they feed the data to a pre-trained CNN to compute the Gram matrices of the activations of the hidden convolutional layers, and apply feature disentanglement to extract the domain-specified features. Concatenation of the diagonal entries of Gram matrices and the domain-specified features gives the final domain embedding. These methods, however, highly rely on the pre-trained models and training process, and lack theoretical guarantees. On the opposite side of the spectrum is directly measuring the discrepancy between domains.

Discrepancy measures of domains. Over the years, numerous notions of discrepancy to measure the (dis)similarity of datasets (domains) were proposed, including $L_{1}$ -norm (batu2000testing), generalized Kolmogorov-Smirnov distance (devroye1996parametric; kifer2004detecting), and loss-oriented discrepancy distance (mansour2009domain). In the context of domain adaptation, generalized Kolmogorov-Smirnov distance (later known as the $A$ -distance) is a principled notation of discrepancy, which is a relaxation of total variation. ben2006analysis show that the target performance (generalization error) is controlled by the empirical estimate of the source domain error and the $A$ -distance between source and target domains. Another widely used distance is the Maximum Mean Discrepancy (MMD) gretton2006kernel, which captures the (dis)similarity of the embedding of distribution measures in a reproducing kernel Hilbert space. pan2010domain propose to learn transfer components across domains in reproducing kernel Hilbert space using MMD, and show that the subspace spanned by these transfer components preserves data properties. Such domain discrepancy methods, however, can not take labels into account, and thus may not be enough to reflect the similarity of tasks.

Optimal transport based task similarity. In recent years, metrics rooted in the optimal transport problem, e.g., the “Wasserstein distance” villani2009optimal; villani2021topics (or the “earth mover’s distance” rubner2000earth; solomon2014earth), have attracted growing interest in the machine learning community. Wasserstein distance is a rigorous metric of probability measures endowed with desired statistical convergence behavior, in contrast to other classical discrepancies (e.g. KL-divergence, total variation, JS-divergence, Hellinger distance, Maximum mean discrepancy, etc). OT based metrics are widely used in generative modeling (arjovsky2017wasserstein; liu2019wasserstein), domain adaptation (courty2017joint; courty2014domain; alvarez2018gromov), graph embedding kolouri2020wasserstein; xu2019gromov, and neural architecture search (kandasamy2018neural).

alvarez2020geometric propose a notion of distance between two datasets in a supervised learning setting. They introduce Optimal Transport Dataset Distance (OTDD) based on the OT theory, which can be thought as a hierarchical OT distance where the transportation cost measures the distance between samples as well as labels. With the assumption that the label-induced distributions can be approximated by Gaussians, the distance between labels is defined as the Bures-Wasserstein distance.

tan2021otce introduce another OT-based method to measure the transferability, named OTCE (Optimal Transport Conditional Entropy) score. In particular, they first use the entropic optimal transport to estimate domain differences and then use the optimal coupling between the source and target distributions to compute the conditional entropy of the target task given source task. The OTCE is defined by the linear combination of the OT distance and the conditional entropy. Both OTDD and OTCE were shown to be effectively aligned with forward transfer, however, the computation of the pairwise Wasserstein distances among increasing number of datasets remains expensive. This hinders the application of these methods to problems where one needs to perform nearest dataset retrieval frequently (e.g., memory replay approaches in continual learning).

Computation Cost of OT Distance. Calculating the Wasserstein distance involves solving an $n^{2}$ dimension linear programming and the computational cost is $O (n^{3} log (n))$ for a pair of $n$ -size empirical distributions (pele2009fast). To facilitate the computation, one common method is adding entropic regularization (cuturi2013sinkhorn; peyre2017computational), by which the original linear programming problem is converted into a strictly convex problem. By applying the Sinkhorn-Knopp algorithm (peyre2017computational; chizat2018scaling) to find an $ϵ$ - accurate solution, the computational complexity reduces to $O (n^{2} log (n) / ϵ^{3})$ (altschuler2017near). However, this technique suffers a stability-accuracy trade-off. When the regularity coefficient is high, the objective is biased toward the entropy term; when it is small, the Sinkhorn algorithm will not be numerically stable.

Figure 1: Wasserstein Task Embedding framework. Given tasks $τ_{1}$ and $τ_{2}$ with input space $R^{d}$ , WTE first map them into $R^{d + L}$ as probability distributions $μ_{1}$ and $μ_{2}$ by MDS, then apply WE to get vector $v_{1}$ and $v_{2}$ with respect to a fixed reference measure $μ_{0}$ . Here $M$ is the size of reference set.

Preliminaries

Multidimensional scaling (MDS)

Multidimensional scaling (MDS) cox2008multidimensional is a non-linear dimensionality reduction approach that embeds $N$ samples into an $l$ -dimensional Euclidean space while preserving their pairwise distances. Given a set of high-dimensional data $X = {x_{n}}_{n = 1}^{N}$ and the proximity matrix $D \in R^{N \times N}$ , where $D_{i, j} = d (x_{i}, x_{j})$ , and $d (\cdot, \cdot)$ denotes the metric in $X$ , the goal of MDS is to construct a distance-preserving map from $X$ to a lower-dimensional Euclidean space $R^{l}$ . Depending on the objective and inputs, MDS can be classified into metric MDS and non-metric MDS. Specifically, metric MDS aims to find a map $ψ : X \to R^{l}$ such that

min ψ \sqrt{\frac{\sum_{i, j} {(d (x_{i}, x_{j}) - ∥ ψ (x_{i}), ψ (x_{j}) ∥)}_{i}^{2}}{\sum_{i, j} d (x_{i}, x_{j})^{2}}},

(1)

which can be solved by Algo. 1.

procedure MDS(

X = {x_{n}}_{n = 1}^{N}, D = [d (x_{i}, x_{j})]_{i, j}, l

)

B = - \frac{1}{2} (i d - \frac{1}{N} 1_{N}) D (i d - \frac{1}{N} 1_{N})

Eigen-decomposition

B = V Λ V^{T}

Rearrange

Λ

into

^Λ

with descending order of variances

Rearrange

V

into

^V

in correspondence with

^Λ

{^Λ}_{(l)} =^Λ [: l, : l]

;

{^V}_{(l)} =^V [: l]

return

ψ (X) = {^Λ}_{(l)}^{\frac{1}{2}} {^V}_{(l)}

end procedure

Algorithm 1 Multidimensional Scaling

Note that MDS not only works for Euclidean distances, but also for other dissimilarities such as Wasserstein distances 5609205; hamm2022wassmap.

Wasserstein Distances

Let $μ, ν$ be Borel probability measures on $X \subseteq R^{d}$ with finite $p^{th}$ moment, and the corresponding probability density functions are $p_{μ}$ and $p_{ν}$ , i.e. $d μ = p_{μ} d x$ , $d ν = p_{ν} d x$ . The 2-Wasserstein distance between $μ$ and $ν$ is defined as villani2009optimal:

W_{2} (μ, ν) = {(inf γ \in Γ (μ, ν) \int_{X \times X} ∥ x - x^{'} ∥^{2} d γ (x, x^{'}))}_{X \times X}^{\frac{1}{2}},

(2)

where $Γ (μ, ν)$ is the set of all transport plans between $μ$ and $ν$ , i.e. probability measures on $X \times X$ with marginals $μ$ and $ν$ . We also note that by Brenier theorem Brenier1991PolarFA, given two absolutely continuous probability measures $μ, ν$ on $R^{d}$ with densities $p_{μ}, p_{ν}$ , there exists a convex function $ϕ$ such that $T = \nabla ϕ$ is a transport map sending $μ$ to $ν$ . Moreover, it is the optimal map in the Monge-Kantorovitch optimal transport problem with quadratic cost:

(3)

where $T = \nabla ϕ$ pushes $μ$ to $ν$ , denoted by $T_{#} μ = ν$ .

Wasserstein Embedding (WE)

Wasserstein Embedding (wang2013linear; kolouri2016continuous; courty2017learning; kolouri2021wasserstein) provides a Hilbertian embedding for probability distributions such that the Euclidean distance between the embedded vectors approximates the 2-Wasserstein distance between the two distributions. Let ${μ_{i}}_{i = 0}^{I}$ be a set of a probability distributions over $X \subseteq R^{d}$ with densities ${p_{i}}_{i = 0}^{I}$ . We fix $μ_{0}$ as the reference measure. Assume $T_{i}$ is the optimal transport map that pushes $μ_{0}$ to $μ_{i}$ , the Wasserstein embedding of $μ_{i}$ is through a function $Φ$ defined as

Φ (μ_{i}) = (T_{i} - i d) \sqrt{p_{0}}

(4)

where the $i d$ is the identity function, i.e., $i d (x) = x$ . $Φ$ admits nice properties including but not limited to kolouri2021wasserstein:

$d (μ_{i}, μ_{j}) ≜ ∥ Φ (μ_{i}) - Φ (μ_{j}) ∥_{2}$ is a true metric between $μ_{i}$ and $μ_{j}$ , moreover, it approximates the 2-Wasserstein distance: $d (μ_{i}, μ_{j}) \approx W_{2} (μ_{i}, μ_{j})$ .
In particular, $∥ Φ (μ_{i}) ∥_{2} = ∥ Φ (μ_{i}) - Φ (μ_{0}) ∥_{2} \approx W_{2} (μ_{i}, μ_{0})$ . Here we leveraged the fact $Φ (μ_{0}) = 0$ .

Although these hold true for both continuous and discrete measures ${μ_{i}}_{i = 0}^{I}$ , we focus on the (uniformly distributed) discrete setting in this paper and provide the following numerical computation details. Let $p_{i} = \frac{1}{N_{i}} \sum_{n = 1}^{N_{i}} δ_{x_{n}^{i}}$ , where $δ_{x}$ is the Dirac delta function centered at $x \in X$ and $X_{i} = {x_{n}^{i}}_{n = 1}^{N_{i}}$ is the set of locations of non-negative mass for $μ_{i}$ . Then the Kantorovich problem with quadratic cost between $μ_{i}$ and $μ_{0}$ can be formulated as

min π \in Π_{i} N_{0} \sum n = 1 N_{i} \sum k = 1 π_{n k} ∥ x_{n}^{i} - x_{k}^{0} ∥_{2}^{2}

(5)

where the feasible set is

Π_{i} ≜ {π \in R^{N_{0} \times N_{i}} | N_{0} N_{i} \sum n = 1 π_{n k} = N_{i} N_{0} \sum k = 1 π_{n k} = 1, \forall n, k} .

(6)

The optimal transport plan $π_{i}^{*}$ is the minimizer of the above optimization problem, which is solved by linear program at cost $O (N^{3} log (N))$ , $N$ being the number of input samples. To avoid mass splitting, the barycentric projection wang2013linear assigns each $x_{j}^{0}$ in the reference distribution to the center of mass it is sent to and thus outputs an approximated Monge map $T_{i}$ . Then the Wasserstein Embedding for input $X_{i}$ is calculated by

Φ (X_{i}) = (T_{i} - X_{0}) / \sqrt{N_{0}} \in R^{N_{0} \times d} .

(7)

One of the motivations behind Wasserstein embedding is to reduce the need for computing pairwise Wasserstein distances. Given $M$ datasets, computation of $\frac{M (M - 1)}{2}$ Wasserstein distances across all distinct pairs is impractically expensive especially when $M$ is large, while leveraging Wasserstein embedding, it suffices to calculate only $M$ Wasserstein distances and the pairwise Euclidean distances between the embedded distributions.

Method

In this section, we specify the problem setting, review the OTDD framework, and then propose our Wasserstein task embedding (WTE).

Problem Setting

In supervised classification problems, tasks are represented by input-label pairs and can be denoted as $τ = X \times Y = {(x_{n}, y_{n})}_{n = 1}^{N}$ , where $X \subseteq R^{d}$ is the data/inputs and $Y$ is the labels. We aim to define a similarity/dissimilarity measure for tasks that enables task clustering and allows for better control over the forward and backward transfer.

Figure 2: Label-to-label Bures-Wasserstein distance (left) and label MDS embedding Euclidean distances (middle) between MNIST and USPS datasets, squared error is provided on the right.

Optimal Transport Dataset Distance (OTDD)

Let $X = {x_{n} \in R^{d}}_{n = 1}^{N}$ be the input set with labels (classes) $Y = {y_{j}}_{j = 1}^{J}$ . For each $j$ , let $C_{y_{j}} \subseteq X$ denote the class with label $y_{j}$ . Following the OTDD framework, let $τ = {(x_{i}, y_{j}) \in X \times Y | x_{i} \in C_{y_{j}}}_{i, j}$ denote the set of data-label pairs. OTDD encodes the label $y_{j}$ as distribution $ν_{y_{j}}$ , where $ν_{y_{j}} = \frac{1}{| C_{y_{j}} |} \sum_{x_{i} \in C_{y_{j}}} δ_{x_{i}}$ . The ground distance in $τ$ is then defined by combining the Euclidean distance between the data points and the 2-Wasserstein distance between label distributions:

d_{τ} ((x, y), (x^{'}, y^{'})) := {(∥ x - x^{'} ∥^{2} + W_{2}^{2} (ν_{y}, ν_{y^{'}}))}_{y}^{\frac{1}{2}} .

(8)

Based on this metric, the OT distance between two distributions $μ_{i}$ and $μ_{j}$ on $τ$ is

d_{O T} (μ_{i}, μ_{j}) = inf π \in Π (μ_{i}, μ_{j}) \int_{τ \times τ} d_{τ} (z, z^{'})^{2} d π (z, z^{'}),

(9)

where $Π (μ_{i}, μ_{j})$ denotes the set of transport plans between $μ_{i}$ and $μ_{j}$ . Note that Eq. 9 is a hierarchical transport problem, as the transportation cost itself depends on calculation of the Wasserstein distance. To avoid the computational cost of a hierarchical optimal transport problem, alvarez2020geometric replace the Wasserstein distance in Eq. 8 with the Bures-Wasserstein distance (malago2018wasserstein; bhatia2019bures), which assumes that $ν_{y}$ s are Gaussian distributions. Throughout the paper, we consider only the exact-OTDD, as opposed to the entropy-regularized and other variants.

Wasserstein Task Embedding

We define a task-2-vec framework (Fig. 1) using Wasserstein embedding (WE) such that the (squared) Euclidean distance between two vectors approximates the OTDD between the original tasks, and denote this embedding by WTE. We later show in the experiment section that the Euclidean distance between the embedded task vectors is not only highly predictive of forward transferability, but also significantly correlates with the backward transferability (catastrophic forgetting).

Label Embedding via MDS. The combination of optimal transport metric with MDS technique was first introduced as an approach to characterize and contrast the distribution of nuclear structure in different tissue classes (normal, benign, cancer, etc.) 5609205, and further studied in image manifold learning hamm2022wassmap. In short, it seeks to isometrically map probability distributions to vectors in relatively low-dimensional space. We leverage the prior work and define an approximated isometry on the label distributions by 1) calculating the pairwise Wasserstein distances and 2) applying MDS to obtain embedded vectors. We adopt the same simplification as in OTDD, that is, assuming the label distributions are Gaussians to replace Wasserstein distances with the closed form Bures-Wasserstein distance:

W_{2}^{2} (ν_{y}, ν_{y^{'}}) = ∥ u_{y} - u_{y^{'}} ∥_{2}^{2} + Tr (Σ_{y} + Σ_{y^{'}} - 2 (Σ_{y}^{\frac{1}{2}} Σ_{y^{'}} Σ_{y}^{\frac{1}{2}})^{\frac{1}{2}})

(10)

where $u$ and $Σ$ denote the mean and covariance matrix of Gaussian distributions. In consistence with previous notations, let us denote the label (MDS) embedding operator by $ψ$ , then

W_{2}^{2} (ν_{y}, ν_{y^{'}}) \approx ∥ ψ (ν_{y}) - ψ (ν_{y^{'}}) ∥_{2}^{2},

(11)

where $ψ (ν_{y}), ψ (ν_{y^{'}}) \in R^{l}$ are vectors whose dimension $l$ is selected to balance the trade-off between accuracy and computation cost tenenbaum2000global. Having both inputs and labels represented as vectors, we concatenate these two components and map the data-label pairs $τ$ to $τ^{'} \subseteq R^{d + l}$ such that

$d_{τ} ((x, y), (x^{'}, y^{'}))$	$\approx (∥ x - x^{'} ∥_{2}^{2} + ∥ ψ (ν_{y}) - ψ (ν_{y^{'}}) ∥_{2}^{2})^{\frac{1}{2}}$
	$= (∥ [x, ψ (ν_{y})] - [x^{'}, ψ (ν_{y^{'}})] ∥_{2}^{2})^{\frac{1}{2}}$
	$= ∥ [x, ψ (ν_{y})] - [x^{'}, ψ (ν_{y^{'}})] ∥_{2}$	(12)

where $[\cdot, \cdot]$ denotes the concatenation operator and the domain $τ^{'} = {[x, ψ (ν_{y})]}_{x \in X}$ is equipped with $l_{2}$ norm. Fig. 2 shows this approximation performance among labels in MNIST LeCun2005TheMD and USPS 291440 datasets. MDS embeddings can capture the pairwise relationships with a maximum of $7.26 %$ error by $10$ -dimensional vectors.

Figure 3: (Top row) forward transfer error drop and (bottom row) catastrophic forgetting against WTE distance on *NIST (left), Split-CIFAR100 (middle) and Split-Tiny ImageNet (right) over five runs. Pearson’s $r$ and the corresponding $p$ -value are reported on top of each experiment setting.

Wasserstein Embedding. By Eq. Wasserstein Task Embedding and Eq. 9, OTDD can be approximated by the squared 2-Wasserstein distance between the distributions over input-(label MDS embedding) pairs, $τ^{'}$ . Then we leverage the Wasserstein embedding framework to embed the updated task distributions into a Hilbert space, with the goal of reducing the cost of computing pairwise Wasserstein distances. Again, we emphasize that this can bring down the cost from quadratic to linear with the number of task distributions.

The WTE algorithm is summarized in Algo. 2. The outputs are the vector representations of input tasks with respect to a pre-determined MDS dimension and WE reference distribution.

procedure WTE(

{X_{i} = {(x_{n}^{i}, y_{n}^{i})}_{n = 1}^{N_{i}}}_{i = 1}^{I}

)

Calculate label-to-label distance matrix W (Eq. 10)

Calculate

ψ (ν_{y})

for all distinct labels

y

(Algo. 1)

Stack each input with its label vector:

x \to [x, ψ (ν_{y})]

Calculate the WE

v_{i} = Φ ({[x_{n}^{i}, ψ (ν_{y_{n}^{i}})]}_{n = 1}^{N_{i}})

(Eq. 7) return

{v_{i}}_{i = 1}^{I}

end procedure

Algorithm 2 Wasserstein Task Embedding

Experiments

To assess the effectiveness of our WTE framework, we empirically validate the correlation between WTE distance and forward/backward transferability on several datasets. Moreover, we provide both qualitative and quantitative comparison results with OTDD, and show WTE distance is well aligned with OTDD, and meanwhile is notably faster to compute. We use the MDS toolkit in scikit-learn and the exact linear programming solver in Python Optimal Transport (POT) flamary2021pot library for implementing WE. We carry out the distance calculations on CPU and all the model training experiments on a 24GB NVIDIA RTX A5000 GPU.

Figure 4: Pairwise OTDD (left) and WTE distances (middle) on the *NIST task group, and their correlation diagram (right). Notice that OTDD (Eq. 9) is the squared $W_{2}$ , we report the squared WTE distances and adjust to the same scale according to the cost function. Adjusted WTE distance is strongly correlated with OTDD, with correlation coefficient $r = 0.91$ .

Datasets

We conduct experiments on the following three task groups:

*NIST task group consists of the handwritten digits dataset MNIST LeCun2005TheMD and its extensions EMNIST cohen2017emnist, FashionMNIST xiao2017fashion, KMNIST Clanuwat2018DeepLF along with USPS 291440. We choose the mnist split for EMNIST dataset and thus all tasks contain 10 classes of gray-scale images. All datasets have a training set of 60,000 samples and a test set of 10,000 samples, except USPS, with a total of 9,298 samples. We resize the images from USPS into $28 \times 28$ pixel level to match with the others.

Split-CIFAR100 task group is generated by randomly splitting the CIFAR-100 Krizhevsky2009LearningML dataset with 100 image categories into 10 smaller tasks, each of which is a classification with 10 classes. There are 600 $32 \times 32$ color images in the training set and 100 in the test set per class.

Split-Tiny ImageNet task group follows the same splitting scheme as in Split-CIFAR100. We randomly divide the Tiny ImageNet Le2015TinyIV into 10 disjoint tasks with 20 classes. Each class contains 500 training images, 50 validating images and 50 test images. For better model performance, we first rescale each sample to $256 \times 256$ and then perform a center crop to get $224 \times 224$ pixel images.

Results

To study the transfer behaviors against the WTE distances, we fix a model architecture for each task group. Specifically, we use ResNet18 he2016deep on both *NIST and Split-Tiny ImageNet, and ResNet34 he2016deep on Split-CIFAR100. In the forward transfer setting, for each source-target task pair, we first train the head (i.e., the classifier) of a randomly initialized backbone on the target task, and use the test performance as the baseline. Next, we adapt from a model pre-trained on the source task and finetune the head on the target task. We define the forward transferability of the source-target pair as the performance gain, i.e. error drop when adapting from the source task. To analyse backward transfer, all source tasks are trained jointly during the first phase to avoid task bias, then in the second phase the model learns only the target task and suffers from “forgetting” the previous tasks. We use the catastrophic forgetting, i.e., negative backward transfer as a measure of backward (in)transferability. In implementations of WE, the reference distribution is fixed for each task group, and is randomly generated by upsampling random images at a lower spatial resolution to entail some smooth structure.

Figure 5: Wall-clock computation time comparison on the *NIST (left) and Split-CIFAR100 (right) task groups.

Fig. 3 summarizes the correlation diagrams between our proposed WTE distance and the forward/backward transferability on the aforementioned three task groups. WTE distance is negatively correlated with forward transferability, and positively correlated with catastrophic forgetting. In all scenarios, the correlation is strong and statistically significant, which confirms the efficacy of WTE distance as a measure of task similarities. We also visualize the comparison between WTE distance and OTDD on the *NIST task group in Fig. 4, showing strong correlation between the two distances.

Computation Complexity

As we mentioned before, OTDD suffers from a prohibitive computational cost as the number of tasks grows large. The pairwise OTDD calculation for a set of $M$ tasks requires $O (M^{2} N^{3} log (N))$ time in the worst case, where $N$ is the largest number of samples among the tasks. WTE distance requires solving only $M$ optimal transport problem, leading to $O (M N^{3} log (N))$ complexity. To better demonstrate the efficiency of WTE distance, we report the wall-clock time comparison on the *NIST and Split-CIFAR100 in Fig. 5.

Conclusion

In this paper, we propose Wasserstein task embedding (WTE), a model-agnostic task embedding framework for measuring task (dis)similarities in supervised classification problems. We perform a label embedding through multi-dimensional scaling and leverage the 2-Wasserstein embedding framework to embed tasks into a vector space, in which the Euclidean distance between the embedded points approximates the 2-Wasserstein distance between tasks. We demonstrate that our proposed task embedding distance is correlated with forward and backward transfer on *NIST, Split-CIFAR100 and Split-Tiny ImageNet task groups while being significantly faster than existing methods. In particular, we show statistically significant negative correlation between the WTE distances and the forward transfer, and positive correlation with the catastrophic forgetting (i.e. negative backward transfer). Lastly, we show the alignment of WTE distance with OTDD, but with a significant computational advantage as the number of tasks grows.