GaKCo: a Fast Gapped k-mer string Kernel using Counting - supplementary information, code and data
2017-12-29T18:18:38Z (GMT) by
This data record consists of supplementary information files, all datasets and code used in the ECML PKDD 2017 paper GaKCo: a Fast Gapped k-mer string Kernel using Counting (README included).
GaKCo is a fast and naturally parallelizable algorithm for gapped k-mer based string kernel calculation. GaKCo uses associative arrays to calculate the co-occurrence of substrings using cumulative counting. This algorithm is scalable to larger dictionary size and more number of mismatches.
GaKCo_ECML17_Supplementary.pdf - presents schematics of the GakCo algorithm, formal proof regarding Hamming Distance Property, justification of GaKCo's Sort and Count Method, connections to previous studies, details of the datasets, an overview of empirical Performance of GaKCo versus Neural Networks and other related experiments.
data.zip - contains 38 training and testing nucloetide and peptide datasets in .fasta format: a text-based format in which base pairs or amino acids are represented using single-letter codes. A sequence in FASTA format begins with a single-line description, followed by lines of sequence data. The description line is distinguished from the sequence data by a greater-than (">") symbol in the first column. Protein, DNA and text dictionaries are provided in .txt format. Both formats are openly accessible via text edit software.
**Datasets for GaKCo:**
We perform 19 different classification tasks to evaluate the performance of GaKCo. These tasks belong to the discussed three categories: (1) TF binding site prediction (DNA dataset), (2) Remote Protein Homology prediction (protein dataset), and (3) Character-based English text classification (text dataset).
code.zip - contains code files in C++ format: .cpp, .h, and bash file: .sh to compile GaKCo using the openMP g++ compiler, and to obtain kernel output using data files and dictionaries above. See below for further detail.
Compiling GaKCo (with openMP) :
g++ -c GaKCo.cpp -o GaKCo -fopenmp
To get kernel output:
#User options: > 5, =(0,..,g-1), = 0 for single-thread/ 1 for multithread
Bash script to run end-to-end kernel calculation:
String Kernel (SK) techniques, especially those using gapped k-mers as features (gk), have obtained great success in classifying sequences like DNA, protein, and text. However, the state-of-the-art gk-SK runs extremely slow when we increase the dictionary size (Σ) or allow more mismatches (M). This is because current gk-SK uses a trie-based algorithm to calculate co-occurrence of mismatched substrings resulting in a time cost proportional to O(ΣM). We propose a fast algorithm for calculating Gapped k-mer Kernel using Counting (GaKCo). GaKCo uses associative arrays to calculate the co-occurrence of substrings using cumulative counting. This algorithm is fast, scalable to larger Σ and M, and naturally parallelizable. We provide a rigorous asymptotic analysis that compares GaKCo with the state-of-the-art gk-SK. Theoretically, the time cost of GaKCo is independent of the ΣM term that slows down the trie-based approach. Experimentally, we observe that GaKCo achieves the same accuracy as the state-of-the-art and outperforms its speed by factors of 2, 100, and 4, on classifying sequences of DNA (5 datasets), protein (12 datasets), and character-based English text (2 datasets).