Multi-task Learning - Multi Task Learning

Multi-task learning (MTL) is a subfield of machine learning in which multiple learning tasks are solved at the same time, while exploiting commonalities and differences across tasks. This can result in improved learning efficiency and prediction accuracy for the task-specific models, when compared to training the models separately.. Early versions of MTL were called "hints"

In a widely cited 1997 paper, Rich Caruana gave the following characterization:

Multitask Learning is an approach to inductive transfer that improves generalization by using the domain information contained in the training signals of related tasks as an inductive bias. It does this by learning tasks in parallel while using a shared representation; what is learned for each task can help other tasks be learned better.

In the classification context, MTL aims to improve the performance of multiple classification tasks by learning them jointly. One example is a spam-filter, which can be treated as distinct but related classification tasks across different users. To make this more concrete, consider that different people have different distributions of features which distinguish spam emails from legitimate ones, for example an English speaker may find that all emails in Russian are spam, not so for Russian speakers. Yet there is a definite commonality in this classification task across users, for example one common feature might be text related to money transfer. Solving each user's spam classification problem jointly via MTL can let the solutions inform each other and improve performance. Further examples of settings for MTL include multiclass classification and multi-label classification.

Multi-task learning works because regularization induced by requiring an algorithm to perform well on a related task can be superior to regularization that prevents overfitting by penalizing all complexity uniformly. One situation where MTL may be particularly helpful is if the tasks share significant commonalities and are generally slightly under sampled. However, as discussed below, MTL has also been shown to be beneficial for learning unrelated tasks.

Methods

source : www.researchgate.net

Task grouping and overlap

Within the MTL paradigm, information can be shared across some or all of the tasks. Depending on the structure of task relatedness, one may want to share information selectively across the tasks. For example, tasks may be grouped or exist in a hierarchy, or be related according to some general metric. Suppose, as developed more formally below, that the parameter vector modeling each task is a linear combination of some underlying basis. Similarity in terms of this basis can indicate the relatedness of the tasks. For example with sparsity, overlap of nonzero coefficients across tasks indicates commonality. A task grouping then corresponds to those tasks lying in a subspace generated by some subset of basis elements, where tasks in different groups may be disjoint or overlap arbitrarily in terms of their bases. Task relatedness can be imposed a priori or learned from the data.

Exploiting unrelated tasks

One can attempt learning a group of principal tasks using a group of auxiliary tasks, unrelated to the principal ones. In many applications, joint learning of unrelated tasks which use the same input data can be beneficial. The reason is that prior knowledge about task relatedness can lead to sparser and more informative representations for each task grouping, essentially by screening out idiosyncrasies of the data distribution. Novel methods which builds on a prior multitask methodology by favoring a shared low-dimensional representation within each task grouping have been proposed. The programmer can impose a penalty on tasks from different groups which encourages the two representations to be orthogonal. Experiments on synthetic and real data have indicated that incorporating unrelated tasks can result in significant improvements over standard multi-task learning methods.

Transfer of knowledge

Related to multi-task learning is the concept of knowledge transfer. Whereas traditional multi-task learning implies that a shared representation is developed concurrently across tasks, transfer of knowledge implies a sequentially shared representation. Large scale machine learning projects such as the deep convolutional neural network GoogLeNet, an image-based object classifier, can develop robust representations which may be useful to further algorithms learning related tasks. For example, the pre-trained model can be used as a feature extractor to perform pre-processing for another learning algorithm. Or the pre-trained model can be used to initialize a model with similar architecture which is then fine-tuned to learn a different classification task.

Mathematics

source : mmlab.ie.cuhk.edu.hk

Reproducing Hilbert space of vector valued functions (RKHSvv)

The MTL problem can be cast within the context of RKHSvv (a complete inner product space of vector-valued functions equipped with a reproducing kernel). In particular, recent focus has been on cases where task structure can be identified via a separable kernel, described below. The presentation here derives from Ciliberto et al, 2015.

RKHSvv concepts

Suppose the training data set is ${\displaystyle {\mathcal {S}}_{t}=\{(x_{i}^{t},y_{i}^{t})\}_{i=1}^{n_{t}}}$, with ${\displaystyle x_{i}^{t}\in {\mathcal {X}}}$, ${\displaystyle y_{i}^{t}\in {\mathcal {Y}}}$, where ${\displaystyle t}$ indexes task, and ${\displaystyle t\in 1,...,T}$. Let ${\displaystyle n=\sum _{t=1}^{T}n_{t}}$. In this setting there is a consistent input and output space and the same loss function ${\displaystyle {\mathcal {L}}:\mathbb {R} \times \mathbb {R} \rightarrow \mathbb {R} _{+}}$ for each task: . This results in the regularized machine learning problem:

where ${\displaystyle {\mathcal {H}}}$ is a vector valued reproducing kernel Hilbert space with functions ${\displaystyle f:{\mathcal {X}}\rightarrow {\mathcal {Y}}^{T}}$ having components ${\displaystyle f_{t}:{\mathcal {X}}\rightarrow {\mathcal {Y}}}$.

The reproducing kernel for the space ${\displaystyle {\mathcal {H}}}$ of functions ${\displaystyle f:{\mathcal {X}}\rightarrow \mathbb {R} ^{T}}$ is a symmetric matrix-valued function ${\displaystyle \Gamma :{\mathcal {X}}\times {\mathcal {X}}\rightarrow \mathbb {R} ^{T\times T}}$ , such that ${\displaystyle \Gamma (\cdot ,x)c\in {\mathcal {H}}}$ and the following reproducing property holds:

The reproducing kernel gives rise to a representer theorem showing that any solution to equation 1 has the form:

Separable kernels

The form of the kernel ${\displaystyle \Gamma }$ induces both the representation of the feature space and structures the output across tasks. A natural simplification is to choose a separable kernel, which factors into separate kernels on the input space ${\displaystyle {\mathcal {X}}}$ and on the tasks ${\displaystyle \{1,...,T\}}$. In this case the kernel relating scalar components ${\displaystyle f_{t}}$ and ${\displaystyle f_{s}}$ is given by ${\textstyle \gamma ((x_{i},t),(x_{j},s))=k(x_{i},x_{j})k_{T}(s,t)=k(x_{i},x_{j})A_{s,t}}$. For vector valued functions ${\displaystyle f\in {\mathcal {H}}}$ we can write ${\displaystyle \Gamma (x_{i},x_{j})=k(x_{i},x_{j})A}$, where ${\displaystyle k}$ is a scalar reproducing kernel, and ${\displaystyle A}$ is a symmetric positive semi-definite ${\displaystyle T\times T}$ matrix. Henceforth denote ${\displaystyle S_{+}^{T}=\{{\text{PSD matrices}}\}\subset \mathbb {R} ^{T\times T}}$ .

This factorization property, separability, implies the input feature space representation does not vary by task. That is, there is no interaction between the input kernel and the task kernel. The structure on tasks is represented solely by ${\displaystyle A}$. Methods for non-separable kernels ${\displaystyle \Gamma }$ is an current field of research.

For the separable case, the representation theorem is reduced to ${\textstyle f(x)=\sum _{i=1}^{N}k(x,x_{i})Ac_{i}}$. The model output on the training data is then ${\displaystyle KCA}$ , where ${\displaystyle K}$ is the ${\displaystyle n\times n}$ empirical kernel matrix with entries ${\textstyle K_{i,j}=k(x_{i},x_{j})}$, and ${\displaystyle C}$ is the ${\displaystyle n\times T}$ matrix of rows ${\displaystyle c_{i}}$.

With the separable kernel, equation 1 can be rewritten as

where ${\displaystyle V}$ is a (weighted) average of ${\displaystyle {\mathcal {L}}}$ applied entry-wise to Y and KCA. (The weight is zero if ${\displaystyle Y_{i}^{t}}$ is a missing observation).

Note the second term in P can be derived as follows:

${\displaystyle ||f||_{\mathcal {H}}^{2}=\langle \sum _{i=1}^{n}k(\cdot ,x_{i})Ac_{i},\sum _{j=1}^{n}k(\cdot ,x_{j})Ac_{j}\rangle _{\mathcal {H}}}$

${\displaystyle =\sum _{i,j=1}^{n}\langle k(\cdot ,x_{i})Ac_{i},k(\cdot ,x_{j})Ac_{j}\rangle _{\mathcal {H}}}$ (bilinearity)

${\displaystyle =\sum _{i,j=1}^{n}\langle k(x_{i},x_{j})Ac_{i},c_{j}\rangle _{\mathbb {R} ^{T}}}$ (reproducing property)

${\displaystyle =\sum _{i,j=1}^{n}k(x_{i},x_{j})c_{i}^{\top }Ac_{j}=tr(KCAC^{\top })}$

Known task structure

Task structure representations

There are three largely equivalent ways to represent task structure: through a regularizer; through an output metric, and through an output mapping.

Regularizer - With the separable kernel, it can be shown (below) that ${\textstyle ||f||_{\mathcal {H}}^{2}=\sum _{s,t=1}^{T}A_{t,s}^{\dagger }\langle f_{s},f_{t}\rangle _{{\mathcal {H}}_{k}}}$, where ${\displaystyle A_{t,s}^{\dagger }}$ is the ${\displaystyle t,s}$ element of the pseudoinverse of ${\displaystyle A}$, and ${\displaystyle {\mathcal {H}}_{k}}$ is the RKHS based on the scalar kernel ${\displaystyle k}$, and ${\textstyle f_{t}(x)=\sum _{i=1}^{n}k(x,x_{i})A_{t}^{\top }c_{i}}$. This formulation shows that ${\displaystyle A_{t,s}^{\dagger }}$ controls the weight of the penalty associated with ${\textstyle \langle f_{s},f_{t}\rangle _{{\mathcal {H}}_{k}}}$. (Note that ${\textstyle \langle f_{s},f_{t}\rangle _{{\mathcal {H}}_{k}}}$ arises from ${\textstyle ||f_{t}||_{{\mathcal {H}}_{k}}=\langle f_{t},f_{t}\rangle _{{\mathcal {H}}_{k}}}$.)

Proof:

${\displaystyle ||f||_{\mathcal {H}}^{2}=\langle \sum _{i=1}^{n}\gamma ((x_{i},t_{i}),\cdot )c_{i}^{t_{i}},\sum _{j=1}^{n}\gamma ((x_{j},t_{j}),\cdot )c_{j}^{t_{j}}\rangle _{\mathcal {H}}}$

${\displaystyle =\sum _{i,j=1}^{n}c_{i}^{t_{i}}c_{j}^{t_{j}}\gamma ((x_{i},t_{i}),(x_{j},t_{j}))}$

${\displaystyle =\sum _{i,j=1}^{n}\sum _{s,t=1}^{T}c_{i}^{t}c_{j}^{s}k(x_{i},x_{j})A_{s,t}}$

${\displaystyle =\sum _{i,j=1}^{n}k(x_{i},x_{j})\langle c_{i},Ac_{j}\rangle _{\mathbb {R} ^{T}}}$

${\displaystyle =\sum _{i,j=1}^{n}k(x_{i},x_{j})\langle c_{i},AA^{\dagger }Ac_{j}\rangle _{\mathbb {R} ^{T}}}$

${\displaystyle =\sum _{i,j=1}^{n}k(x_{i},x_{j})\langle Ac_{i},A^{\dagger }Ac_{j}\rangle _{\mathbb {R} ^{T}}}$

${\displaystyle =\sum _{i,j=1}^{n}\sum _{s,t=1}^{T}(Ac_{i})^{t}(Ac_{j})^{s}k(x_{i},x_{j})A_{s,t}^{\dagger }}$

${\displaystyle =\sum _{s,t=1}^{T}A_{s,t}^{\dagger }\langle \sum _{i=1}^{n}k(x_{i},\cdot )(Ac_{i})^{t},\sum _{j=1}^{n}k(x_{j},\cdot )(Ac_{j})^{s}\rangle _{{\mathcal {H}}_{k}}}$

${\displaystyle =\sum _{s,t=1}^{T}A_{s,t}^{\dagger }\langle f_{t},f_{s}\rangle _{{\mathcal {H}}_{k}}}$

Output metric - an alternative output metric on ${\displaystyle {\mathcal {Y}}^{T}}$ can be induced by the inner product ${\displaystyle \langle y_{1},y_{2}\rangle _{\Theta }=\langle y_{1},\Theta y_{2}\rangle _{\mathbb {R} ^{T}}}$. With the squared loss there is an equivalence between the separable kernels ${\displaystyle k(\cdot ,\cdot )I_{T}}$ under the alternative metric, and ${\displaystyle k(\cdot ,\cdot )\Theta }$, under the canonical metric.

Output mapping - Outputs can be mapped as ${\displaystyle L:{\mathcal {Y}}^{T}\rightarrow {\mathcal {\tilde {Y}}}}$ to a higher dimensional space to encode complex structures such as trees, graphs and strings. For linear maps ${\displaystyle L}$, with appropriate choice of separable kernel, it can be shown that ${\displaystyle A=L^{\top }L}$.

Task structure examples

Via the regularizer formulation, one can represent a variety of task structures easily.

• Letting ${\textstyle A^{\dagger }=\gamma I_{T}+(\gamma -\lambda ){\frac {1}{T}}{\mathbf {1}}{\mathbf {1}}^{\top }}$ (where ${\displaystyle I_{T}}$ is the TxT identity matrix, and ${\textstyle {\mathbf {1}}{\mathbf {1}}^{\top }}$ is the TxT matrix of ones) is equivalent to letting ${\displaystyle \gamma }$ control the variance ${\textstyle \sum _{t}||f_{t}-{\bar {f}}||_{{\mathcal {H}}_{k}}}$ of tasks from their mean ${\textstyle {\frac {1}{T}}\sum _{t}f_{t}}$. For example, blood levels of some biomarker may be taken on ${\displaystyle T}$ patients at ${\displaystyle n_{t}}$ time points during the course of a day and interest may lie in regularizing the variance of the predictions across patients.
• Letting ${\displaystyle A^{\dagger }=\alpha I_{T}+(\alpha -\lambda )M}$ , where ${\displaystyle M_{t,s}={\frac {1}{|G_{r}|}}\mathbb {I} (t,s\in G_{r})}$ is equivalent to letting ${\displaystyle \alpha }$ control the variance measured with respect to a group mean: ${\displaystyle \sum _{r}\sum _{t\in G_{r}}||f_{t}-{\frac {1}{|G_{r}|}}\sum _{s\in G_{r})}f_{s}||}$. (Here ${\displaystyle |G_{r}|}$ the cardinality of group r, and ${\displaystyle \mathbb {I} }$ is the indicator function). For example people in different political parties (groups) might be regularized together with respect to predicting the favorability rating of a politician. Note that this penalty reduces to the first when all tasks are in the same group.
• Letting ${\displaystyle A^{\dagger }=\delta I_{T}+(\delta -\lambda )L}$, where ${\displaystyle L=D-M}$ is the Laplacian for the graph with adjacency matrix M giving pairwise similarities of tasks. This is equivalent to giving a larger penalty to the distance separating tasks t and s when they are more similar (according to the weight ${\displaystyle M_{t,s}}$,) i.e. ${\displaystyle \delta }$ regularizes ${\displaystyle \sum _{t,s}||f_{t}-f_{s}||_{{\mathcal {H}}_{k}}^{2}M_{t,s}}$.
• All of the above choices of A also induce the additional regularization term ${\textstyle \lambda \sum _{t}||f||_{{\mathcal {H}}_{k}}^{2}}$ which penalizes complexity in f more broadly.

Learning tasks together with their structure

Learning problem P can be generalized to admit learning task matrix A as follows:

Choice of ${\displaystyle F:S_{+}^{T}\rightarrow \mathbb {R} _{+}}$ must be designed to learn matrices A of a given type. See "Special cases" below.

Optimization of Q

Restricting to the case of convex losses and coercive penalties Ciliberto et al have shown that although Q is not convex jointly in C and A, a related problem is jointly convex.

Specifically on the convex set ${\displaystyle {\mathcal {C}}=\{(C,A)\in \mathbb {R} ^{n\times T}\times S_{+}^{T}|Range(C^{\top }KC)\subseteq Range(A)\}}$, the equivalent problem

is convex with the same minimum value. And if ${\displaystyle (C_{R},A_{R})}$ is a minimizer for R then ${\displaystyle (C_{R}A_{R}^{\dagger },A_{R})}$ is a minimizer for Q.

R may be solved by a barrier method on a closed set by introducing the following perturbation:

The perturbation via the barrier ${\displaystyle \delta ^{2}tr(A^{\dagger })}$ forces the objective functions to be equal to ${\displaystyle +\infty }$ on the boundary of ${\displaystyle R^{n\times T}\times S_{+}^{T}}$ .

S can be solved with a block coordinate descent method, alternating in C and A. This results in a sequence of minimizers ${\displaystyle (C_{m},A_{m})}$ in S that converges to the solution in R as ${\displaystyle \delta _{m}\rightarrow 0}$, and hence gives the solution to Q.

Special cases

Spectral penalties - Dinnuzo et al suggested setting F as the Frobenius norm ${\displaystyle {\sqrt {tr(A^{\top }A)}}}$. They optimized Q directly using block coordinate descent, not accounting for difficulties at the boundary of ${\displaystyle \mathbb {R} ^{n\times T}\times S_{+}^{T}}$.

Clustered tasks learning - Jacob et al suggested to learn A in the setting where T tasks are organized in R disjoint clusters. In this case let ${\displaystyle E\in \{0,1\}^{T\times R}}$ be the matrix with ${\displaystyle E_{t,r}=\mathbb {I} ({\text{task }}t\in {\text{group }}r)}$. Setting ${\displaystyle M=I-E^{\dagger }E^{T}}$, and ${\displaystyle U={\frac {1}{T}}{\mathbf {11}}^{\top }}$, the task matrix ${\displaystyle A^{\dagger }}$ can be parameterized as a function of ${\displaystyle M}$: ${\displaystyle A^{\dagger }(M)=\epsilon _{M}U+\epsilon _{B}(M-U)+\epsilon (I-M)}$ , with terms that penalize the average, between clusters variance and within clusters variance respectively of the task predictions. M is not convex, but there is a convex relaxation ${\displaystyle {\mathcal {S}}_{c}=\{M\in S_{+}^{T}:I-M\in S_{+}^{T}\land tr(M)=r\}}$. In this formulation, ${\displaystyle F(A)=\mathbb {I} (A(M)\in \{A:M\in {\mathcal {S}}_{C}\})}$.

Generalizations

Non-convex penalties - Penalties can be constructed such that A is constrained to be a graph Laplacian, or that A has low rank factorization. However these penalties are not convex, and the analysis of the barrier method proposed by Ciliberto et al does not go through in these cases.

Non-separable kernels - Separable kernels are limited, in particular they do not account for structures in the interaction space between the input and output domains jointly. Future work is needed to develop models for these kernels.

Applications

source : manufacturingscience.asmedigitalcollection.asme.org

Spam filtering

Using the principles of MTL, techniques for collaborative spam filtering that facilitates personalization have been proposed. In large scale open membership email systems, most users do not label enough messages for an individual local classifier to be effective, while the data is too noisy to be used for a global filter across all users. A hybrid global/individual classifier can be effective at absorbing the influence of users who label emails very diligently from the general public. This can be accomplished while still providing sufficient quality to users with few labeled instances.

Web search

Using boosted decision trees, one can enable implicit data sharing and regularization. This learning method can be used on web-search ranking data sets. One example is to use ranking data sets from several countries. Here, multitask learning is particularly helpful as data sets from different countries vary largely in size because of the cost of editorial judgments. It has been demonstrated that learning various tasks jointly can lead to significant improvements in performance with surprising reliability.

RoboEarth

In order to facilitate transfer of knowledge, IT infrastructure is being developed. One such project, RoboEarth, aims to set up an open source internet database that can be accessed and continually updated from around the world. The goal is to facilitate a cloud-based interactive knowledge base, accessible to technology companies and academic institutions, which can enhance the sensing, acting and learning capabilities of robots and other artificial intelligence agents.

Software package

source : 0xmachine.com

The Multi-Task Learning via StructurAl Regularization (MALSAR) Matlab package implements the following multi-task learning algorithms:

• Mean-Regularized Multi-Task Learning
• Multi-Task Learning with Joint Feature Selection
• Robust Multi-Task Feature Learning
• Trace-Norm Regularized Multi-Task Learning
• Alternating Structural Optimization
• Incoherent Low-Rank and Sparse Learning
• Robust Low-Rank Multi-Task Learning
• Clustered Multi-Task Learning
• Multi-Task Learning with Graph Structures

See also

source : www.computer.org

• Artificial Intelligence
• Artificial neural network
• Evolutionary computation
• Human-based genetic algorithm
• Kernel methods for vector output
• Machine Learning
• Robotics

References

External links

• The Biosignals Intelligence Group at UIUC
• Washington University at St. Louis Depart. of Computer Science

Software

• The Multi-Task Learning via Structural Regularization Package
• Online Multi-Task Learning Toolkit (OMT) A general-purpose online multi-task learning toolkit based on conditional random field models and stochastic gradient descent training (C#, .NET)