Bonsai -- Diverse and Shallow Trees for Extreme Multi-label Classification

TL;DR

Bonsai introduces shallow, diverse trees for extreme multi-label classification, leveraging novel label representations to improve accuracy and training speed, especially on tail labels, on very large datasets.

Contribution

The paper proposes a new algorithm, Bonsai, that generalizes label representation and learns shallow trees for efficient and accurate extreme multi-label classification.

Findings

Outperforms state-of-the-art methods in accuracy on large datasets.

Achieves approximately 200 times faster training than one-vs-rest methods.

Effectively improves tail-label prediction accuracy.

Abstract

Extreme multi-label classification (XMC) refers to supervised multi-label learning involving hundreds of thousand or even millions of labels. In this paper, we develop a suite of algorithms, called Bonsai, which generalizes the notion of label representation in XMC, and partitions the labels in the representation space to learn shallow trees. We show three concrete realizations of this label representation space including : (i) the input space which is spanned by the input features, (ii) the output space spanned by label vectors based on their co-occurrence with other labels, and (iii) the joint space by combining the input and output representations. Furthermore, the constraint-free multi-way partitions learnt iteratively in these spaces lead to shallow trees. By combining the effect of shallow trees and generalized label representation, Bonsai achieves the best of both worlds - fast training which is comparable to state-of-the-art tree-based methods in XMC, and much better prediction accuracy, particularly on tail-labels. On a benchmark Amazon-3M dataset with 3 million labels, \bonsai outperforms a state-of-the-art one-vs-rest method in terms of prediction accuracy, while being approximately 200 times faster to train. The code for Bonsai is available at \url{https://github.com/xmc-aalto/bonsai}

Tables

Table 1. Table 1: Multi-label datasets used in the experiment. APpL and ALpP represent average points per label and average labels per point respectively.

Dataset	# Training	# Test	# Labels	# Features	APpL	ALpP
EURLex-4K	15,539	3,809	3993	5000	25.7	5.3
Wikipedia-31K	14,146	6,616	30,938	101,938	8.5	18.6
WikiLSHTC-325K	1,778,351	587,084	325,056	1,617,899	17.4	3.2
Wikipedia-500K	1,813,391	783,743	501,070	2,381,304	24.7	4.7
Amazon-670K	490,499	153,025	670,091	135,909	3.9	5.4
Amazon-3M	1,717,899	742,507	2,812,281	337,067	31.6	36.1

Table 2. Table 2: p ​ r ​ e ​ c ​ @ ​ k 𝑝 𝑟 𝑒 𝑐 @ 𝑘 prec@k (P@k) on benchmark datasets for k = 1 , 3 ​  and  ​ 5 𝑘 1 3  and  5 k=1,3\text{ and }5 . For each case of P@k and dataset, the best performed score is highlighted in bold. Entries marked ”-” imply the corresponding method could not scale to the particular dataset, thus the scores are unavailable.

Dataset	Our Approach (Bonsai)			Embedding based		Tree based		Linear one-vs-rest
	Bonsai-i	Bonsai-o	Bonsai-io	SLEEC	RobustXML	Fast-XML	Parabel	PD-Sparse	DiSMEC
EURLex-4K
P@1	83.0	82.5	82.9	79.3	78.7	71.4	82.2	76.4	82.4
P@3	69.7	69.4	69.4	64.3	63.5	59.9	68.7	60.4	68.5
P@5	58.4	58.1	58.0	52.3	51.4	50.4	57.5	49.7	57.7
Wikipedia-31K
P@1	84.7	84.70	84.8	85.5	85.5	82.5	84.2	73.8	84.1
P@3	73.6	73.57	73.6	73.6	74.0	66.6	72.5	60.9	74.6
P@5	64.7	64.81	64.8	63.1	63.8	56.7	63.4	50.4	65.9
WikiLSHTC-325K
P@1	66.6	63.4	65.8	55.5	53.5	49.3	65.0	58.2	64.4
P@3	44.5	42.8	44.1	33.8	31.8	32.7	43.2	36.3	42.5
P@5	33.0	32.0	32.7	24.0	29.9	24.0	32.0	28.7	31.5
Wikipedia-500K
P@1	69.2	68.7	69.1	48.2	41.3	54.1	68.7	-	70.2
P@3	49.8	48.8	49.7	29.4	30.1	35.5	49.6	-	50.6
P@5	38.8	37.6	38.8	21.2	19.8	26.2	38.6	-	39.7
Amazon-670K
P@1	45.5	44.5	45.7	35.0	31.0	33.3	44.9	-	44.7
P@3	40.3	39.8	40.6	31.2	28.0	29.3	39.8	-	39.7
P@5	36.5	36.4	36.9	28.5	24.0	26.1	36.0	-	36.1
Amazon-3M
P@1	48.4	47.5	48.5	-	-	44.2	47.5	-	47.8
P@3	45.6	44.7	45.5	-	-	40.8	44.6	-	44.9
P@5	43.4	42.6	43.5	-	-	38.6	42.5	-	42.8

Table 3. Table 3: Coverage@k (C@k) statistics comparing Parabel and Bonsai . Along each C@k and dataset configuration, the best performing score is highlighted in bold.

Dataset	Methods	C@1	C@3	C@5
EUR-Lex	Parabel	31.46	43.11	54.38
	Bonsai	31.38	44.09	55.61
Wiki10	Parabel	7.00	5.77	6.76
	Bonsai	7.52	6.82	8.01
WikiLSHTC	Parabel	22.73	35.94	43.18
	Bonsai	24.14	38.49	46.37
Amazon-670k	Parabel	32.73	33.77	38.82
	Bonsai	33.28	34.76	40.11
Amazon-3M	Parabel	21.16	20.49	21.81
	Bonsai	22.27	21.89	23.36

Equations

$$N_r=N_1{r}^{-\beta }$$

$${\mathcal{V}}_i={\text{Y}}^T\text{X}={\left[\begin{array}{c}{\text{v}}_1^T\\ {\text{v}}_2^T\\ ⋮\\ {\text{v}}_L^T\end{array}\right]}_{L\times D}$$

$$\mathrm{where}\text{X}={\left[\begin{array}{c}{\text{x}}_1^T\\ {\text{x}}_2^T\\ ⋮\\ {\text{x}}_N^T\end{array}\right]}_{N\times D},\text{Y}={\left[\begin{array}{c}{\text{y}}_1^T\\ {\text{y}}_2^T\\ ⋮\\ {\text{y}}_N^T\end{array}\right]}_{N\times L}.$$

$${\mathcal{V}}_o={\text{Y}}^T\text{Y}={\left[\begin{array}{c}{\text{v}}_1^T\\ {\text{v}}_2^T\\ ⋮\\ {\text{v}}_L^T\end{array}\right]}_{L\times L},\mathrm{where}\text{Y}={\left[\begin{array}{c}{\text{y}}_1^T\\ {\text{y}}_2^T\\ ⋮\\ {\text{y}}_N^T\end{array}\right]}_{N\times L}.$$

$${\mathcal{V}}_j={\text{Y}}^T\text{Z}={\left[\begin{array}{c}{\text{v}}_1^T\\ {\text{v}}_2^T\\ ⋮\\ {\text{v}}_L^T\end{array}\right]}_{L\times (D+L)}$$

$$\mathrm{where}\text{Z}={\left[\begin{array}{c}{\text{z}}_1^T\\ {\text{z}}_2^T\\ ⋮\\ {\text{z}}_N^T\end{array}\right]}_{N\times (D+L)}\text{Y}={\left[\begin{array}{c}{\text{y}}_1^T\\ {\text{y}}_2^T\\ ⋮\\ {\text{y}}_N^T\end{array}\right]}_{N\times L}$$

$$\underset{{\text{c}}_1,\dots ,{\text{c}}_K{\in }^{\eta }}{\min }\left[\sum _{k=1}^K\sum _{\ell \in {\text{c}}_i}d({\text{v}}_{\ell },{\text{c}}_k)\right]$$

$$\underset{{\text{w}}_k}{\min }\left[||{\text{w}}_k|{|}_2^2+C\sum _{i=1}^{|{\text{X}}_k|}\mathcal{L}(s_{k_i}{\text{w}}_k^T{\text{x}}_i)\right]$$

$$\text{Pr}({\text{y}}_{\ell }=1\mid \text{x})$$

$$\times \prod _{n\in \mathcal{A}(e)}\text{Pr}(z_n=1\mid z_{p(n)}=1,\text{x})$$

$$prec@k(\hat{\text{y}},\text{y})$$

$$nDCG@k(\hat{\text{y}},\text{y})$$

$$PSP@k(\hat{\text{y}},\text{y})$$

$$PSnDCG@k(\hat{\text{y}},\text{y})$$

$$C@k=|\mathbf{P}|/|\mathbf{L}|$$

Full text

By combining the effect of shallow trees and generalized label representation, Bonsai achieves the best of both worlds - fast training which is comparable to state-of-the-art tree-based methods in XMC, and much better prediction accuracy, particularly on tail-labels. On a benchmark Amazon-3M dataset with 3 million labels, Bonsai outperforms a state-of-the-art one-vs-rest method in terms of prediction accuracy, while being approximately 200 times faster to train. The code for Bonsai is available at https://github.com/xmc-aalto/bonsai.

1 Introduction

Extreme Multi-label Classification (XMC) refers to supervised learning of a classifier which can automatically label an instance with a small subset of relevant labels from an extremely large set of all possible target labels. Machine learning problems consisting of hundreds of thousand labels are common in various domains such as product categorization for e-commerce [27, 32, 7, 1], hash-tag suggestion in social media [11], annotating web-scale encyclopedia [29], and image-classification [21, 10]. It has been demonstrated that, the framework of XMC can also be leveraged to effectively address ranking problems arising in bid-phrase suggestion in web-advertising and suggestion of relevant items for recommendation systems [31].

From the machine learning perspective, building effective extreme classifiers is faced with the computational challenge arising due to large number of (i) output labels, (ii) input training instances, and (iii) input features. Another important statistical characteristic of the datasets in XMC is that a large fraction of labels are tail labels, i.e., those which have very few training instances that belong to them (also referred to as power-law, fat-tailed distribution and Zipf’s law). Formally, let $N_{r}$ denote the size of the $r$-th ranked label, when ranked in decreasing order of number of training instances that belong to that label, then :

[TABLE]

where $N_{1}$ represents the size of the 1-st ranked label and $\beta>0$ denotes the exponent of the power law distribution. This distribution is shown in Figure 1 for a benchmark dataset, WikiLSHTC-325K from the XMC repository [8]. In this dataset, only $\sim$150,000 out of 325,000 labels have more than 5 training instances in them. Tail labels exhibit diversity of the label space, and contain informative content not captured by the head or torso labels. Indeed, by predicting well the head labels, yet omitting most of the tail labels, an algorithm can achieve high accuracy [36]. However, such behavior is not desirable in many real world applications, where fit to power-law distribution has been observed [2, 6].

1.1 Related work

Various works in XMC can be broadly categorized into one of the four strands :

One-vs-rest : As the name suggests, these methods learn a classifier per label which distinguishes it from rest of the labels. In terms of prediction accuracy and label diversity, these methods have been shown to be among the best performaning ones for XMC [5, 41, 6]. However, due to their reliance on a distributed training framework, it remains challenging to employ them in resource constrained environments. 2. 2.

Tree-based : Tree-based methods implement a divide-and-conquer paradigm and scale to large label sets in XMC by partitioning the labels space. As a result, these scheme of methods have the computational advantage of enabling faster training and prediction [31, 15, 17, 26, 38]. Approaches based on decision trees have also been proposed for multi-label classification and those tailored to XMC regime [18, 33]. However, tree-based methods suffer from error propagation in the tree cascade as also observed in hierarchical classification [3, 4]. As a result, these methods tend to perform particularly worse on metrics which are sensitive for tail-labels [30]. 3. 3.

Label embedding : Label-embedding approaches assume that, despite large number of labels, the label matrix is effectively low rank and therefore project it to a low-dimensional sub-space. These approaches have been at the fore-front in multi-label classification for small scale problems with few tens or hundred labels [14, 35, 37, 23]. For power-law distributed labels in XMC settings, the crucial assumption made by the embedding-based approaches of a low rank label space breaks down [39, 9, 34]. Under this condition, embedding based approaches leads to high prediction error. 4. 4.

Deep learning : Deeper architectures on top of word-embeddings have also been explored in recent works [24, 19, 28]. However, their performance still remains sub-optimal compared to the methods discussed above which are based on bag-of-words feature representations. This is mainly due to the data scarcity in tail-labels which is substantially below the sample complexity required for deep learning methods to reach their peak performance.

Therefore, a central challenge in XMC is to build classifiers which retain the accuracy of one-vs-rest paradigm while being as efficiently trainable as the tree-based methods. Recently, there have been efforts for speeding up the training of existing classifiers by better initialization and exploiting the problem structure [13, 22, 16]. In a similar vein, a recently proposed tree-based method, Parabel [30], partitions the label space recursively into two child nodes using 2-means clustering. It also maintains a balance between these two label partitions in terms of number of labels. Each intermediate node in the resulting binary label-tree is like a meta-label which captures the generic properties of its constituent labels. The leaves of the tree consist of the actual labels from the training data. During training and prediction each of these labels is distinguished from other labels under the same parent node through the application of a binary classifier at internal nodes and one-vs-all classifier for the leaf nodes. By combination of tree-based partitioning and one-vs-rest classifier, it has been shown to give better performance than previous tree-based methods [31, 15, 17] while simultaneously allowing efficient training.

However, in terms of prediction performance, Parabel remains sub-optimal compared to one-vs-rest approaches. In addition to error propagation due to cascading effect of the deep trees, its performance is particularly worse on tail labels. This is the result of two strong constraints in its label partitioning process, (i) each parent node in the tree has only two child nodes, and (ii) at each node, the labels are partitioned into equal sized parts, such that the number of labels under the two child nodes differ by at most one. As a result of the coarseness imposed by the binary partitioning of labels, the tail labels get subsumed by the head labels.

1.2 Bonsai overview

In this paper, we develop a family of algorithms, called Bonsai. At a high level, Bonsai follows a similar paradigm which is common in most tree-based approaches, i.e., label partitioning followed by learning classifiers at the internal nodes. However, it has two main features, which distinguish it from state-of-the-art tree based approaches. These are summarized below :

• •

Generalized label representation - In this work, we argue that the notion of representing the labels is quite general, and their exist various meaningful manifestations of the label representation space. As three concrete examples, we show the applicability of the following representations of labels : (i) input space representation as a function of feature vectors (ii) output space representation based on their co-occurrence with other labels, and (iii) a combination of the output and input representations. In this regard, our work generalizes the view taken in most earlier works, which have represented labels only in the input space such as by representing them as sum of the training instances in which their are active  [30, 38]. We show that these representations, when combined with shallow trees (described next), surpass existing methods demonstrating the efficacy of the proposed generalization representation.

• •

Shallow trees - To avoid error propagation in the tree cascade, we propose to construct a shallow tree architecture. This is achieved by enabling (i) a flexible clustering via $K-$means for $K>2$, and (ii) relaxing balancedness constraints in the clustering step. Multi-way partitioning initializes diverse sub-groups of labels, and the unconstrained nature maintains the diversity during the entire process. These are in contrast to tree-based methods which impose such constraints for a balanced tree construction. As we demonstrate in our empirical findings, by relaxing the constraints, Bonsai leads to prediction diversity and significantly better tail-label coverage.

By synergizing the effect of a richer label representation and shallow trees, Bonsai achieves the best of both worlds - prediction diversity better than state-of-the-art tree-based methods with comparable training speed, and prediction accuracy at par with one-vs-rest methods. The code for Bonsai is available at https://github.com/xmc-aalto/bonsai.

2 Formal description of Bonsai

We assume to be given a set of $N$ training points $\left\{(\textbf{x}_{i},\textbf{y}_{i})\right\}_{i=1}^{N}$ with $D$ dimensional feature vectors $\textbf{x}_{i}\in^{D}$ and $L$ dimensional label vectors $\textbf{y}_{i}\in\{0,1\}^{L}$. Without loss of generality, let the set of labels is represented by $\{1,\ldots,\ell,\ldots,L\}$ Our goal is to learn a multi-label classifier in the form of a vector-valued output function $f:\mathbb{R}^{D}\mapsto\{0,1\}^{L}$. This is typically achieved by minimizing an empirical estimate of $\mathbb{E}_{(\textbf{x},\textbf{y})\sim\mathcal{D}}[\mathcal{L}(\textbf{W};(\textbf{x},\textbf{y}))]$ where $\mathcal{L}$ is a loss function, and samples $(\textbf{x},\textbf{y})$ are drawn from some underlying distribution $\mathcal{D}$. The desired parameters W can take one of the myriad of choices. In the simplest (and yet effective) of setups for XMC such as linear classification, W can be in the form of matrix. In other cases, it can be representative of a deeper architecture or a cascade of classifiers in a tree structured topology. Due to their scalability to extremely large datasets, Bonsai follows a tree-structured partitioning of labels.

In this section, we next present in detail the two main components of Bonsai: (i) generalized label representation and (ii) shallow trees.

2.1 Label representation

In the extreme classification setting, labels can be represented in various ways. To motivate this, as an analogy in terms of publications and their authors, one can think of labels as authors, the papers they write as their training instances, and multiple co-authors of a paper as the multiple labels. Now, one can represent authors (labels) either solely based on the content of the papers they authored (input space representation), or based only on their co-authors (output space representation) or as a combination of the two.

Formally, let each label $\ell$ be represented by $\eta$-dimensional vector $\textbf{v}_{\ell}\in\mathbb{R}^{\eta}$. Now, $\textbf{v}_{\ell}$ can be represented as a function only of input instances $\left\{\textbf{x}_{i}\right\}_{i=1}^{N}$, only of output instances $\left\{\textbf{y}_{i}\right\}_{i=1}^{N}$ or as a combination of both $\left\{(\textbf{x}_{i},\textbf{y}_{i})\right\}_{i=1}^{N}$. We now present three concrete realizations of the label representation $\textbf{v}_{\ell}$. We later show that these representations can be seamlessly combined with shallow tree cascade of classifiers, and yield state-of-the-art performance on XMC tasks.

• a.

Input space representation of $\textbf{v}_{\ell}$ - The label representation for label $\ell$ can be arrived at by summing all the training examples for which it is active. Let $\mathcal{V}_{i}$ be the label representation matrix given by

[TABLE]

[TABLE]

We follow the notation that each bold letter such as x is a vector in column format and $\textbf{x}^{T}$ represents the correponding row vector. Hence, each row $\textbf{v}_{\ell}$ of matrix $\mathcal{V}_{i}$ which represents the label $\ell$, is given by the sum of all the training instances for which label $\ell$ is active. This can also be represented as, $\textbf{v}_{\ell}=\sum_{i=1}^{N}\textbf{y}_{i{\ell}}\textbf{x}_{i}$. Note that even though $\textbf{v}_{\ell}$ also depends on the label vectors, it is still in the same space as the input instance and has dimensionality $D$. Furthermore, each $\textbf{v}_{\ell}$ can be normalized to unit length in euclidean norm as follows : $\textbf{v}_{\ell}:=\textbf{v}_{\ell}/\|\textbf{v}_{\ell}\|_{2}$.

• b.

Output space representation of $\textbf{v}_{\ell}$ - In the multi-label setting, another way to represent the labels is to represent them solely as a function of the degree of their co-occurence with other labels. That is, if two labels co-occur with similar set of labels, then these are bound to be related to each other, and hence should have similar representation. In this case, the label representation matrix $\mathcal{V}_{o}$ is given by

[TABLE]

Here $\mathcal{V}_{o}$ is an $L\times L$ symmetric matrix, where each row $\textbf{v}_{\ell}^{T}$, corresponds to the number of times the label $\ell$ co-occurs with all other labels. Hence these label co-occurrence vectors $\textbf{v}_{\ell}$ give us another way of representing the label $\ell$. It may be noted that in contrast to the previous case, being an output space representation, the dimensionality of the label vector is same as that of the output space having the same dimensionality, i.e. $\eta=L$.

• c.

Joint input-output representation of $\textbf{v}_{\ell}$ - Given the previous input and output space representations of labels, a natural way to extend it is by combining these representations via concatenation. This is achieved as follows, for a training instance $i$ with feature vector $\textbf{x}_{i}$ and corresponding label vector $\textbf{y}_{i}$, let $\textbf{z}_{i}$ be the concatenated vector given by, $\textbf{z}_{i}=[\textbf{x}_{i}\odot\textbf{y}_{i}]$. Then, the joint representation can be computed in the matrix $\mathcal{V}_{j}$ as follows

[TABLE]

[TABLE]

Here each row $\vec{v}_{\ell}$ of the label representation matrix $\mathcal{V}_{j}$ which is the label representation in the joint space, is therefore a concatenation of representations obtained from $\mathcal{V}_{i}$ and $\mathcal{V}_{o}$, hence being of length $(D+L)$. Since both the input vectors $\textbf{x}_{i}$ and output vectors $\textbf{y}_{i}$ are highly sparse, this does not lead to any major computational burden in training.

It may be noted that label representation based solely on the input as considered by recent works [30, 38], can be considered as a special case of our more general formulation of label representation. As also shown later in our empirical findings, in combination with shallow tree cascade of classifiers, partitioning of :

• •

output space representation ($\mathcal{V}_{o}$) yields competitive results compared to state-of-the-art classifiers in XMC such as Parabel.

• •

joint representation ($\mathcal{V}_{j}$) further surpasses the state-of-the-art methods in terms of prediction performance and label diversity.

2.2 Label partitioning via $K$-means clustering

Once we have obtained the representation $\vec{v}_{\ell}$ for each label $\ell$ in the set $S=\left\{1,\ldots,L\right\}$, the next step is to iteratively partition $S$ into disjoint subsets. This is achieved by $K$-means clustering, which also presents many choices such as number of clusters and degree of balancedness among the clusters. Our goal, in this work, is to avoid propagation error in a deep tree cascade. We, therefore, choose a relatively large value of $K$ (e.g. $\geq 100$) which leads to shallow trees.

The clustering step in Bonsai first partitions $S$ into $K$ disjoint sets $\left\{S_{1},\ldots,S_{K}\right\}$. Each of the elements, $S_{k}$, of the above set can be thought of as a meta-label which semantically groups actual labels together in one cluster. Then, $K$ child nodes of the root are created, each contains one of the partitions, $\left\{S_{k}\right\}_{k=1}^{K}$. The same process is repeated on each of the newly-created $K$ child nodes in an iterative manner. In each sub-tree, the process terminates either when the node’s depth exceeds pre-defined threshold $d_{\text{max}}$ or the number of associated labels is no larger than $K$, e.g, $|S_{k}|\leq K$.

Formally, without loss of generality, we assume a non-leaf node has labels $\left\{1,\ldots,L\right\}$. We aim at finding $K$ cluster centers $\textbf{c}_{1},\ldots,\textbf{c}_{K}\in\mathbb{R}^{\eta}$, i.e., in an appropriate space (input, output, or joint) by optimizing the following :

[TABLE]

where $d(.,.)$ represents a distance function and $\textbf{v}_{\ell}$ represents the vector representation of the label ${\ell}$. The distance function is defined in terms of the dot product as follows : $d(\textbf{v}_{\ell},\textbf{c}_{k})=1-\textbf{v}_{\ell}^{T}\cdot\textbf{c}_{k}$. The above problem is NP-hard and we use the standard $K$-means algorithm (also known as Lloyd’s algorithm) [25]111We also tried $K\text{-means++}$ and observed that faster convergence did not out-weigh extra computation time for seed initialization. for finding an approximate solution to equation (2.5).

The $K$-way unconstrained clustering in Bonsai has the following advantages over Parabel which enforces binary and balanced partitioning :

Initializing label diversity in partitioning : By setting $K>2$, Bonsai allows a varied partitioning of the labels space, rather than grouping all labels in two clusters. This facet of Bonsai is especially favorable for tail labels by allowing them to be part of separate clusters if they are indeed very different from the rest of the labels. Depending on the similarity to other labels, each label can choose to be part of one of the $K$ clusters. 2. 2.

Sustaining label diversity : Bonsai sustains the diversity in the label space by not enforcing the balanced-ness constraint of the form, $\lvert|S_{k}|-|S_{k^{\prime}}|\rvert\leq 1,\forall 1\leq k,k^{\prime}\leq K$ (where $|.|$ operator is overloaded to mean set cardinality for the inner one and absolute value for the outer ones) among the partitions. This makes the Bonsai partitions more data-dependent since smaller partitions with diverse tail-labels are very moderately penalized under this framework. 3. 3.

Shallow tree cascade : Furthermore, $K$-way unconstrained partitioning leads to shallower trees which are less prone propagation error in deeper trees constructed by Parabel. As we will show in Section 3, the diverse partitioning reinforced by shallower architecture leads to better prediction performance, and significant improvement is achieved on tail labels.

A pictorial description of the partitioning scheme of Bonsai and its difference compared to Parabel is also illustrated in Figure 3.

2.3 Learning node classifiers

Once the label space is partitioned into a diverse and shallow tree structure, we learn a $K$-way One-vs-All linear classifier at each node. These classifiers are trained independently using only the training examples that have at least one of the node labels. We distinguish the leaf nodes and non-leaf nodes in the following way : (i) for non-leaf nodes, the classifier learns $K$ linear classifiers separately, each maps to one of the $K$ children. During prediction, the output of each classifier determines whether the test point should traverse down the corresponding child. (ii) for leaf nodes, the classifier learns to predict the actual labels on the node.

Without loss of generality, given a node in the tree, denote by $\{c_{k}\}_{k=1}^{K}$ as the set of its children. For the special case of leaf nodes, the set of children represent the final labels. We learn $K$ linear classifiers parameterized by $\left\{\textbf{w}_{1},\ldots,\textbf{w}_{K}\right\}$, where $\textbf{w}_{k}\in^{D}$ for $\forall k=1,\ldots,K$. Each output label determines if the corresponding $K$ children should be traversed or not.

For each of the child node $c_{k}$, we define the training data as $T_{k}=(\textbf{X}_{k},\textbf{s}_{k})$, where $\textbf{X}_{k}=\left\{\textbf{x}_{i}\mid\textbf{y}_{ik}=1,i=1,\ldots,N\right\}$. Let $\textbf{s}_{k}\in\left\{+1,-1\right\}^{N}$ represent the vector of signs depending on whether $\textbf{y}_{ik}=1$ corresponds to +1 and $\textbf{y}_{ik}=0$ for -1. We consider the following optimization problem for learning linear SVM with squared hinge loss and $\ell_{2}$-regularization

[TABLE]

where $\mathcal{L}(z)=(\max(0,1-z))^{2}$. This is solved using the Newton method based primal implementation in LIBLINEAR [12]. To restrict the model size, and remove spurious parameters, thresholding of small weights is performed as in [5].

The tree-structured architecture of Bonsai is illustrated in Figure 2. The details of Bonsai’s training procedure in the form of an algorithm is shown in Algorithm 1. The partitioning process in Section 2.1 is described as the procedure GROW in the algorithm. The One-vs-All procedure is shown as one-vs-all in Algorithm 1.

2.4 Prediction error propagation in shallow versus deep trees

During prediction, a test point x traverses down the tree. At each non-leaf node, the classifier narrows down the search space by deciding which subset of child nodes x should further traverse. If the classifier decides not to traverse down some child node $c$, all descendants of $c$ will not be traversed. Later, as x reaches to one or more leaf nodes, One-vs-All classifiers are evaluated to assign probabilities to each label. Bonsai uses beam search to avoid the possibility of evaluating all nodes.

The above search space pruning strategy implies errors made at non-leaf nodes could propagate to their descendants. Bonsai sets relatively large values to the branching factor $K$ (typically $100$), resulting in much shallower trees compared to Parabel, and hence significantly reducing error propagation, particularly for tail-labels.

More formally, given a data point x and a label $\ell$ that is relevant to x, we denote $e$ as the leaf node $\ell$ belongs to and $\mathcal{A}(e)$ as the set of ancestor nodes of $e$ and $e$ itself. Note that $|\mathcal{A}(e)|$ is path length from root to $e$. Denote the parent of $n$ as $p(n)$. We define the binary indicator variable $z_{n}$ to take value 1 if node $n$ is visited during prediction and 0 otherwise. From the chain rule, the probability that $\ell$ is predicted as relevant for x is as follows:

[TABLE]

Consider the Amazon-3M dataset with $L\approx 3\times 10^{6}$, setting $K=2$ produces a tree of depth 16. Assuming $\text{Pr}(z_{n}=1\mid z_{p(n)}=1,\textbf{x})=0.95$, for $\forall n\in p(n)$ and $\text{Pr}(\textbf{y}_{\ell}=1\mid z_{e}=1,\textbf{x})=1$, it gives $\text{Pr}(\textbf{y}_{\ell}=1\mid\textbf{x})=(0.95)^{16}\approx 0.46$. This is to say, even if $\text{Pr}(z_{n}=1\mid z_{p(n)}=1,\textbf{x})$ is high (e.g, 0.95) at each $n\in\mathcal{A}(e)$, multiplying them together can result in small probability (e.g, 0.46) if the depth of the tree, i.e., $|\mathcal{A}(e)|$ is large. We choose to mitigate this issue by increasing $K$, and hence limiting the propagation error.

3 Experimental Evaluation

In this section, we detail the dataset description, and the set up for comparison of the proposed approach against state-of-the-art methods in XMC.

3.1 Dataset and evaluation metrics

We perform empirical evaluation on publicly available datasets from the XMC repository 222http://manikvarma.org/downloads/XC/XMLRepository.html curated from sources such as Amazon for item-to-item recommendation tasks and Wikipedia for tagging tasks. The datasets of various scales in terms of number of labels are used, EURLex-4K consisting of approximately 4,000 labels to Amazon-3M consisting of 3 million labels. The datasets also exhibit a wide range of properties in terms of number of training instances, features, and labels. The detailed statistics of the datasets are shown in Table 1.

With applications in recommendation systems, ranking and web-advertising, the objective of the machine learning system in XMC is to correctly recommend/rank/advertise among the top-k slots. We therefore use evaluation metrics which are standard and commonly used to compare various methods under the XMC setting - Precision@$k$ ($prec@k$) and normalised Discounted Cumulative Gain ($nDCG@k$). Given a label space of dimensionality $L$, a predicted label vector $\hat{\textbf{y}}\in\mathbb{R}^{L}$ and a ground truth label vector $\textbf{y}\in\{0,1\}^{L}$ :

[TABLE]

where $DCG@k=\frac{\textbf{y}_{\ell}}{\sum_{l=1}{\frac{1}{\log(\ell+1)}}}$, and $rank_{k}(\hat{\textbf{y}})$ returns the $k$ largest indices of $\hat{\textbf{y}}$.

For better readability, we report the percentage version of above metrics (multiplying the original scores by 100). In addition, we consider $k\in\left\{1,3,5\right\}$.

3.2 Methods for comparison

We consider three different variants of the proposed family of algorithms, Bonsai, which is based on the generalized label representations (discussed in Section 2.1) combined with the shallow tree cascades. We refer the algorithms learnt by partitioning the input space, output space and the joint space as Bonsai-i, Bonsai-o, and Bonsai-io respectively. These are compared against six state-of-the-art algorithms from each of the three main strands for XMC namely, label-embedding, tree-based and one-vs-all methods :

• •

Label-embedding methods: Due to the fat-tailed distritbution of instances among labels, SLEEC [9] makes a locally low-rank assumption on the label space, RobustXML [40] decomposes the label matrix into tail labels and non tail labels so as to enforce an embedding on the latter without the tail labels damaging the embedding. LEML [43] makes a global low-rank assumption on the label space and performs a linear embedding on the label space. As a result, it gives much worse results, and is not compared explicitly in the interest of space.

• •

Tree-based methods: FastXML [31] learns an ensemble of trees which partition the label space by directly optimizing an nDCG based ranking loss function, PFastXML [15] replaces the nDCG loss in FastXML by its propensity scored variant which is unbiased and assigns higher rewards for accurate tail label predictions, Parabel [30] which has been described earlier in the paper.

• •

One-vs-All methods: PD-Sparse [42] enforces sparsity by exploiting the structure of a margin-maximizing loss with L1-penalty, DiSMEC [5] learns one-vs-rest classifiers for every label with weight pruning to control model size.

Since we are considering only bag-of-words representation across all datasets, we do not compare against deep learning methods explicitly. However, it may be noted that despite using raw data and corresponding word-embeddings, deep learning methods in XMC are still sub-optimal in terms of prediction performance in XMC  [24, 19, 20]. More details on the performance of deep methods can be found in [38].

Bonsai is implemented in C++ on a 64-bit Linux system. For all the datasets, we set the branching factor $K=100$ at every tree depth. We will explore the effect of tree depth in details later. This results in depth-1 trees (excluding the leaves which represent the final labels) for smaller datasets such as EURLex-4K, Wikipedia-31K and depth-2 trees for larger datasets such as WikiLSHTC-325K and Wikipedia-500K. Bonsai learns an ensemble of three trees similar to Parabel. For other approaches, the results were reproduced as suggested in the respective papers.

4 Experimental results

In this section, we report the main findings of our empirical evaluation.

4.1 Precision@$k$

The comparison of Bonsai against other baselines is shown in Table 2. The results are averaged over five runs with different initializations of the clustering algorithm. The important findings from these results are the following :

• •

The competitive performance of the different variants of Bonsai shows the success and applicability of the notion of generalized label representation, and their concrete realization discussed in section 2.1. It further highlights that it is possible to enrich these representations further, and achieve better partitioning.

• •

The consistent improvement of Bonsai over Parabel on all datasets validates the choice of higher fanout and advantages of using shallow trees.

• •

Another important insight from the above results is that when the average number of labels per training point are higher such as in Wikipedia-31K, Amazon-670K and Amazon-3M, the joint space label representation, used in Bonsai-io, leads to better partitioning and further improves the strong performance of input only label representation in Bonsai-i.

• •

Even though DiSMEC performs slightly better on Wiki-500K and Wikipedia-31K, its computational complexity of training and prediction is orders of magnitude higher than Bonsai. As a result, while Bonsai can be run in environments with limited computational resources, DiSMEC requires a distributed infrastructure for training and prediction.

4.2 Performance on tail labels

We also evaluate prediction performance on tail labels using propensity scored variants of $prec@k$ and $nDCG@k$. For label $\ell$, its propensity $p_{\ell}$ is related to number of its positive training instances $N_{\ell}$ by $p_{\ell}\propto 1/\left(1+e^{-\log(N_{\ell})}\right)$. With this formulation, $p_{\ell}\approx 1$ for head labels and $p_{\ell}\ll 1$ for tail labels. Let $\textbf{y}\in\{0,1\}^{L}$ and $\hat{\textbf{y}}\in\mathbb{R}^{L}$ denote the true and predicted label vectors respectively. As detailed in [15], propensity scored variants of $P@k$ and $nDCG@k$ are given by

[TABLE]

where $PSDCG@k:=\sum_{\ell\in rank_{k}{(\hat{\textbf{y}})}}{[\frac{\textbf{y}_{\ell}}{p_{\ell}\log(\ell+1)}]}$ , and $rank_{k}(\textbf{y})$ returns the $k$ largest indices of y.

To match against the ground truth, as suggested in [15], we use $100\cdot\mathbb{G}(\{\hat{\textbf{y}}\})/\mathbb{G}(\{\textbf{y}\})$ as the performance metric. For $M$ test samples, $\mathbb{G}(\{\hat{\textbf{y}}\})=\frac{-1}{M}\sum_{i=1}^{M}\mathbb{L}(\hat{\textbf{y}}_{i},\textbf{y})$, where $\mathbb{G}(.)$ and $\mathbb{L}(.,.)$ signify gain and loss respectively. The loss $\mathbb{L}(.,.)$ can take two forms, (i)$\mathbb{L}(\hat{\textbf{y}}_{i},\textbf{y})=-PSnDCG@k$, and (ii) $\mathbb{L}(\hat{\textbf{y}}_{i},\textbf{y})=-PSP@k$. This leads to the two metrics which are sensitive to tail labels and are denoted by $prec\_wt@k$, and $nDCG\_wt@k$.

Figure 4 shows the result w.r.t $prec\_wt@k$, and $nDCG\_wt@k$ among the tree-based approaches. Again, Bonsai-i shows consistent improvement over Parabel. For instance, on WikiLSHTC-325K, the relative improvement over Parabel is approximately 6.7% on $prec\_wt@5$. This further validates the applicability of the shallow tree architecture resulting from the design choices of $K$-way partitioning along with flexibility to allow unbalanced partitioning in Bonsai, which allows tail labels to be assigned into different partitions w.r.t the head ones.

4.3 Unique label coverage

We also evaluate coverage@k, denoted $C@k$, which is the percentage of normalized unique labels present in an algorithm’s top-$k$ labels. Let $\mathbf{P}=P_{1}\cup P_{2}\cup...\cup P_{M}$ where $P_{i}=\{l_{i1},l_{i2},..,l_{ik}\}$ i.e the set of top-$k$ labels predicted by the algorithm for test point $i$ and $M$ is the number of test points. Also, let $\mathbf{L}=L_{1}\cup L_{2}\cup...\cup L_{M}$ where $L_{i}=\{g_{i1},g_{i2},..,g_{ik}\}$ i.e the top-k propensity scored ground truth labels for test point $i$, then, coverage@k is given by

[TABLE]

The comparison between Bonsai-i and Parabel of this metric on five different datasets is shown in Table 3. It shows that the proposed method is more effective in discovering correct unique labels. These results further reinforce the results in the previous section on the diversity preserving feature of Bonsai.

4.4 Impact of tree depth

We next evaluate prediction performance produced by Bonsai trees with different depth values. We set the fan-out parameter $K$ appropriately to achieve the desired tree depth. For example, to partition 4,000 labels into a hierarchy of depth two, we set $K=64$.

In Figure 5, we report the result on three datasets, averaged over ten runs under each setting. The trend is consistent - as the tree depth increases, prediction accuracy tends to drop, though it is not very stark for Wikipedia-31K.

Furthermore, in Figure 6, we show that the shallow architecture is an integral part of the success of the Bonsai family of algorithms. To demonstrate this, we plugged in the label representation used in Bonsai-o into Parabel, called Parabel-o in the figure. As can be seen, Bonsai-o outperforms Parabel-o by a large margin showing that shallow trees substantially alleviate the prediction error.

4.5 Training and prediction time

Growing shallower trees in Bonsai comes at a slight price in terms of training time. It was observed that Bonsai leads to approximately 2-3x increase in training time compared to Parabel. For instance, on three cores, Parabel take one hour for training on WikiLSHTC-325K dataset, while Bonsai takes approximately three hours for the same task. However, it may also be noted that the training process can be performed in an offline manner. Though, unlike Parabel, Bonsai does not come with logarithmic dependence on the number of labels for the computational complexity of prediction. However, its prediction time is typically in milli-seconds, and hence it remains quite practical in XMC applications with real-time constraints such as recommendation systems and advertising.

5 Conclusion

In this paper, we present Bonsai, which is a class of algorithms for learning shallow trees for label partitioning in extreme multi-label classification. Compared to the existing tree-based methods, it improves this process in two fundamental ways. Firstly, it generalizes the notion of label representation beyond the input space representation, and shows the efficacy of output space representation based on its co-occurrence with other labels, and by further combining these in a joint representation. Secondly, by learning shallow trees which prevent error propagation in the tree cascade and hence improving the prediction accuracy and tail-label coverage. The synergizing effects of these two ingredients enables Bonsai to retain the training speed comparable to tree-based methods, while achieving better prediction accuracy as well as significantly better tail-label coverage. As a future work, the generalized label representation can be further enriched by combining with embeddings from raw text. This can lead to the amalgamation of methods studied in this paper with those that are based on deep learning.