# Omikuji [![Build Status](https://dev.azure.com/yubingdong/omikuji/_apis/build/status/tomtung.omikuji?branchName=master)](https://dev.azure.com/yubingdong/omikuji/_build/latest?definitionId=1&branchName=master) [![Crate version](https://img.shields.io/crates/v/omikuji)](https://crates.io/crates/omikuji) [![PyPI version](https://img.shields.io/pypi/v/omikuji)](https://pypi.org/project/omikuji/) An efficient implementation of Partitioned Label Trees (Prabhu et al., 2018) and its variations for extreme multi-label classification, written in Rust🦀 with love💖. ## Features & Performance Omikuji has has been tested on datasets from the [Extreme Classification Repository](http://manikvarma.org/downloads/XC/XMLRepository.html). All tests below are run on a quad-core Intel® Core™ i7-6700 CPU, and we allowed as many cores to be utilized as possible. We measured training time, and calculated precisions at 1, 3, and 5. (Note that, due to randomness, results might vary from run to run, especially for smaller datasets.) ### Parabel, better parallelized Omikuji provides a more parallelized implementation of Parabel (Prabhu et al., 2018) that trains faster when more CPU cores are available. Compared to the [original implementation](http://manikvarma.org/code/Parabel/download.html) written in C++, which can only utilize the same number of CPU cores as the number of trees (3 by default), Omikuji maintains the same level of precision but trains 1.3x to 1.7x faster on our quad-core machine. **Further speed-up is possible if more CPU cores are available**. | Dataset | Metric | Parabel | Omikuji
(balanced,
cluster.k=2) | |----------------- |------------ |--------- |------------------------------------------ | | EURLex-4K | P@1 | 82.2 | 82.1 | | | P@3 | 68.8 | 68.8 | | | P@5 | 57.6 | 57.7 | | | Train Time | 18s | 14s | | Amazon-670K | P@1 | 44.9 | 44.8 | | | P@3 | 39.8 | 39.8 | | | P@5 | 36.0 | 36.0 | | | Train Time | 404s | 234s | | WikiLSHTC-325K | P@1 | 65.0 | 64.8 | | | P@3 | 43.2 | 43.1 | | | P@5 | 32.0 | 32.1 | | | Train Time | 959s | 659s | ### Regular k-means for shallow trees Following Bonsai (Khandagale et al., 2019), Omikuji supports using regular k-means instead of balanced 2-means clustering for tree construction, which results in wider, shallower and unbalanced trees that train slower but have better precision. Comparing to the [original Bonsai implementation](https://github.com/xmc-aalto/bonsai), Omikuji also achieves the same precisions while training 2.6x to 4.6x faster on our quad-core machine. (Similarly, further speed-up is possible if more CPU cores are available.) | Dataset | Metric | Bonsai | Omikuji
(unbalanced,
cluster.k=100,
max\_depth=3) | |----------------- |------------ |--------- |-------------------------------------------------------------- | | EURLex-4K | P@1 | 82.8 | 83.0 | | | P@3 | 69.4 | 69.5 | | | P@5 | 58.1 | 58.3 | | | Train Time | 87s | 19s | | Amazon-670K | P@1 | 45.5* | 45.6 | | | P@3 | 40.3* | 40.4 | | | P@5 | 36.5* | 36.6 | | | Train Time | 5,759s | 1,753s | | WikiLSHTC-325K | P@1 | 66.6* | 66.6 | | | P@3 | 44.5* | 44.4 | | | P@5 | 33.0* | 33.0 | | | Train Time | 11,156s | 4,259s | *\*Precision numbers as reported in the paper; our machine doesn't have enough memory to run the full prediction with their implementation.* ### Balanced k-means for balanced shallow trees Sometimes it's desirable to have shallow and wide trees that are also balanced, in which case Omikuji supports the balanced k-means algorithm used by HOMER (Tsoumakas et al., 2008) for clustering as well. | Dataset | Metric | Omikuji
(balanced,
cluster.k=100) | |----------------- |------------ |------------------------------------------ | | EURLex-4K | P@1 | 82.1 | | | P@3 | 69.4 | | | P@5 | 58.1 | | | Train Time | 19s | | Amazon-670K | P@1 | 45.4 | | | P@3 | 40.3 | | | P@5 | 36.5 | | | Train Time | 1,153s | | WikiLSHTC-325K | P@1 | 65.6 | | | P@3 | 43.6 | | | P@5 | 32.5 | | | Train Time | 3,028s | ### Layer collapsing for balanced shallow trees An alternative way for building balanced, shallow and wide trees is to collapse adjacent layers, similar to the tree compression step used in AttentionXML (You et al., 2019): intermediate layers are removed, and their children replace them as the children of their parents. For example, with balanced 2-means clustering, if we collapse 5 layers after each layer, we can increase the tree arity from 2 to 2⁵⁺¹ = 64. | Dataset | Metric | Omikuji
(balanced,
cluster.k=2,
collapse 5 layers) | |----------------- |------------ |--------------------------------------------------------------- | | EURLex-4K | P@1 | 82.4 | | | P@3 | 69.3 | | | P@5 | 58.0 | | | Train Time | 16s | | Amazon-670K | P@1 | 45.3 | | | P@3 | 40.2 | | | P@5 | 36.4 | | | Train Time | 460s | | WikiLSHTC-325K | P@1 | 64.9 | | | P@3 | 43.3 | | | P@5 | 32.3 | | | Train Time | 1,649s | ## Build & Install Omikuji can be easily built & installed with [Cargo](https://doc.rust-lang.org/cargo/getting-started/installation.html) as a CLI app: ``` cargo install omikuji --features cli --locked ``` Or install from the latest source: ``` cargo install --git https://github.com/tomtung/omikuji.git --features cli --locked ``` The CLI app will be available as `omikuji`. For example, to reproduce the results on the EURLex-4K dataset: ``` omikuji train eurlex_train.txt --model_path ./model omikuji test ./model eurlex_test.txt --out_path predictions.txt ``` ### Python Binding A simple Python binding is also available for training and prediction. It can be install via `pip`: ``` pip install omikuji ``` Note that you might still need to install Cargo should compilation become necessary. You can also install from the latest source: ``` pip install git+https://github.com/tomtung/omikuji.git -v ``` The following script demonstrates how to use the Python binding to train a model and make predictions: ```python import omikuji # Train hyper_param = omikuji.Model.default_hyper_param() # Adjust hyper-parameters as needed hyper_param.n_trees = 5 model = omikuji.Model.train_on_data("./eurlex_train.txt", hyper_param) # Serialize & de-serialize model.save("./model") model = omikuji.Model.load("./model") # Optionally densify model weights to trade off between prediction speed and memory usage model.densify_weights(0.05) # Predict feature_value_pairs = [ (0, 0.101468), (1, 0.554374), (2, 0.235760), (3, 0.065255), (8, 0.152305), (10, 0.155051), # ... ] label_score_pairs = model.predict(feature_value_pairs) ``` ## Usage ``` $ omikuji train --help Train a new omikuji model USAGE: omikuji train [OPTIONS] ARGS: Path to training dataset file The dataset file is expected to be in the format of the Extreme Classification Repository. OPTIONS: --centroid_threshold Threshold for pruning label centroid vectors [default: 0] --cluster.eps Epsilon value for determining linear classifier convergence [default: 0.0001] --cluster.k Number of clusters [default: 2] --cluster.min_size Labels in clusters with sizes smaller than this threshold are reassigned to other clusters instead [default: 2] --cluster.unbalanced Perform regular k-means clustering instead of balanced k-means clustering --collapse_every_n_layers Number of adjacent layers to collapse This increases tree arity and decreases tree depth. [default: 0] -h, --help Print help information --linear.c Cost coefficient for regularizing linear classifiers [default: 1] --linear.eps Epsilon value for determining linear classifier convergence [default: 0.1] --linear.loss Loss function used by linear classifiers [default: hinge] [possible values: hinge, log] --linear.max_iter Max number of iterations for training each linear classifier [default: 20] --linear.weight_threshold Threshold for pruning weight vectors of linear classifiers [default: 0.1] --max_depth Maximum tree depth [default: 20] --min_branch_size Number of labels below which no further clustering & branching is done [default: 100] --model_path Optional path of the directory where the trained model will be saved if provided If an model with compatible settings is already saved in the given directory, the newly trained trees will be added to the existing model") --n_threads Number of worker threads If 0, the number is selected automatically. [default: 0] --n_trees Number of trees [default: 3] --train_trees_1_by_1 Finish training each tree before start training the next This limits initial parallelization but saves memory. --tree_structure_only Build the trees without training classifiers Might be useful when a downstream user needs the tree structures only. ``` ``` $ omikuji test --help Test an existing omikuji model USAGE: omikuji test [OPTIONS] ARGS: Path of the directory where the trained model is saved Path to test dataset file The dataset file is expected to be in the format of the Extreme Classification Repository. OPTIONS: --beam_size Beam size for beam search [default: 10] -h, --help Print help information --k_top Number of top predictions to write out for each test example [default: 5] --max_sparse_density Density threshold above which sparse weight vectors are converted to dense format Lower values speed up prediction at the cost of more memory usage. [default: 0.1] --n_threads Number of worker threads If 0, the number is selected automatically. [default: 0] --out_path Path to the which predictions will be written, if provided ``` ### Data format Our implementation takes dataset files formatted as those provided in the [Extreme Classification Repository](http://manikvarma.org/downloads/XC/XMLRepository.html). A data file starts with a header line with three space-separated integers: total number of examples, number of features, and number of labels. Following the header line, there is one line per each example, starting with comma-separated labels, followed by space-separated feature:value pairs: ``` label1,label2,...labelk ft1:ft1_val ft2:ft2_val ft3:ft3_val .. ftd:ftd_val ``` ## Trivia The project name comes from [o-mikuji](https://en.wikipedia.org/wiki/O-mikuji) (御神籤), which are predictions about one's future written on strips of paper (labels?) at jinjas and temples in Japan, often tied to branches of pine trees after they are read. ## References - Y. Prabhu, A. Kag, S. Harsola, R. Agrawal, and M. Varma, “Parabel: Partitioned Label Trees for Extreme Classification with Application to Dynamic Search Advertising,” in Proceedings of the 2018 World Wide Web Conference, 2018, pp. 993–1002. - S. Khandagale, H. Xiao, and R. Babbar, “Bonsai - Diverse and Shallow Trees for Extreme Multi-label Classification,” Apr. 2019. - G. Tsoumakas, I. Katakis, and I. Vlahavas, “Effective and efficient multilabel classification in domains with large number of labels,” ECML, 2008. - R. You, S. Dai, Z. Zhang, H. Mamitsuka, and S. Zhu, “AttentionXML: Extreme Multi-Label Text Classification with Multi-Label Attention Based Recurrent Neural Networks,” Jun. 2019. ## License Omikuji is licensed under the MIT License.