Volume 10 Supplement 4
Proceedings of the IEEE International Conference on Bioinformatics and Biomedicine (BIBM) 2008
Mixture of experts models to exploit global sequence similarity on biomolecular sequence labeling
- Cornelia Caragea^{1, 2}Email author,
- Jivko Sinapov^{1, 2},
- Drena Dobbs^{3, 4} and
- Vasant Honavar^{1, 2, 4}
DOI: 10.1186/1471-2105-10-S4-S4
© Caragea et al; licensee BioMed Central Ltd. 2009
Published: 29 April 2009
Abstract
Background
Identification of functionally important sites in biomolecular sequences has broad applications ranging from rational drug design to the analysis of metabolic and signal transduction networks. Experimental determination of such sites lags far behind the number of known biomolecular sequences. Hence, there is a need to develop reliable computational methods for identifying functionally important sites from biomolecular sequences.
Results
We present a mixture of experts approach to biomolecular sequence labeling that takes into account the global similarity between biomolecular sequences. Our approach combines unsupervised and supervised learning techniques. Given a set of sequences and a similarity measure defined on pairs of sequences, we learn a mixture of experts model by using spectral clustering to learn the hierarchical structure of the model and by using bayesian techniques to combine the predictions of the experts. We evaluate our approach on two biomolecular sequence labeling problems: RNA-protein and DNA-protein interface prediction problems. The results of our experiments show that global sequence similarity can be exploited to improve the performance of classifiers trained to label biomolecular sequence data.
Conclusion
The mixture of experts model helps improve the performance of machine learning methods for identifying functionally important sites in biomolecular sequences.
Background
Advances in high throughput data acquisition technologies have resulted in rapid increase in the amount of data in biological sciences. For example, progress on sequencing technologies has resulted in the release of hundreds of complete genome sequences. With the exponentially growing number of biomolecular sequences from genome projects and high-throughput experimental studies, sequence annotations do not keep pace with sequencing.
The wet-lab experiments to determine the annotations (e.g., functional site annotations) are still difficult and time consuming. Hence, there is an urgent need for development of computational tools that can accurately annotate biomolecular data.
Machine learning methods currently offer one of the most cost-effective approaches to construction of predictive models in applications where representative training data are available. Biomolecular sequence labeling is an instance of a supervised learning problem. Given a data set (x_{ i }, y_{ i })_{i = 1,⋯, n}of pairs of sequences, x_{ i }= (x_{i,1}x_{i,2}⋯x_{ i, m }) and y_{ i }= (y_{i,1}y_{i,2}⋯y_{ i, m }), where y_{ i, j }in the output sequence is the label for x_{ i, j }in the input (or observation) sequence, j = 1,⋯, m, the task is to learn a classifier that can predict the labels for each element of a new input sequence, x_{ test }.
There is a large body of work on learning predictive models to label biomolecular sequence data. Terribilini et al. [1] trained Naïve Bayes classifiers to identify RNA-protein interface residues in a protein sequence. Yan et al. [2] developed a two-stage classifier to identify protein-protein interaction sites. Qian and Sejnowski [3] trained Neural Networks to predict protein secondary structure, i.e., classifying each residue in a protein sequence into one of the three classes: helix (H), strand (E) or coil (C). Caragea et al. [4] and Kim et al. [5] used Support Vector Machines to identify residues in a protein sequence that undergo post-translational modifications.
Typically, to solve the biomolecular sequence labeling problem using standard machine learning algorithms, each element in a sequence is encoded using a local, fixed-length window corresponding to the target element and its sequence context (an equal number of its sequence neighbors on each side) [6]. The classifier is trained to label the target element. This procedure can produce reliable results in settings where there exists a local sequence pattern that is predictive of the label for the target site. However, there are cases where the local amino acid distribution around functionally important sites in a given set of proteins is highly variable. For example, in identifying RNA-protein and DNA-protein interface residues from amino acid sequences, there is typically no consensus sequence around each site.
Classifiers trained using machine learning to distinguish "positive" examples from the "negative" ones, must "learn" to do so by learning the characteristics associated with known "positive" and "negative" examples. The greater the commonality among members of a subset, the more likely it is that a machine learning approach will be successful in identifying the predictive characteristics.
Against this background, we hypothesize that classifiers trained to label biomolecular sequence data can be improved by taking into account the global sequence similarity between the protein sequences in addition to the local features extracted around each site. The intuition behind this hypothesis is that the more similar two sequences are, the greater the likelihood that their functional sites have similar patterns. Therefore, we propose to improve the biomolecular sequence labeling problem by using a machine learning approach, that is, a mixture of experts model that considers the global similarity between protein sequences when building the model and making the predictions. We evaluate our approach to learning a mixture of experts model on two biomolecular sequence labeling tasks: RNA- and DNA-protein interface prediction tasks.
Results
The main result of our study is that taking into account global sequence similarity through the means of a mixture of experts model can improve the performance of the classifiers trained to label biomolecular sequence data.
The mixture of experts that exploits the global similarity between the protein sequences in a data set in addition to the local features extracted around each residue outperforms the baseline classifiers on the biomolecular sequence labeling task
We trained mixtures of Naïve Bayes (NB) and Logistic Regression (LR) classifiers on both RNA- and DNA-protein interface prediction tasks considered in this study to predict whether or not a residue in a protein sequence is an interface residue. We used various identity cutoffs to construct the data sets. The mixture of experts models have a hierarchical structure that is constructed using 2-way spectral clustering based on a global similarity functions, i.e., we computed the entries in the similarity matrix W by applying the Needleman-Wunsch global alignment algorithm on each pair of sequences. The Blosum62 substitution matrix was used for costs. The resulting entries in the matrix W are normalized and scaled so that each value is between 0 and 1. The mixture of experts models consist of NB and LR at the leaves, respectively (see Methods section for further details).
Experimental results on the RNA-protein sequence data set. Experimental results with Naive Bayes (NB) and Logistic Regression (LR) models, and Mixture of Experts (ME) models on the non-redundant RNA-protein sequence data set, where the identity cutoffs are 30% and 90%. The results are shown for default threshold θ = 0.5. ME-NB-global and ME-LR-global use NB and LR at the leaves and exploits the global sequence similarity to construct the hierarchical structure. ME-NB-local exploits the local sequence similarity to construct the hierarchical structure. ME-NB-random randomizes the global similarity matrix and constructs the hierarchical structure based on the randomized matrix.
RNA-protein 30% | RNA-protein 90% | |||||||||
---|---|---|---|---|---|---|---|---|---|---|
Classifier | Precision | Recall | CC | FM | AUC | Precision | Recall | CC | FM | AUC |
NB | 0.58 | 0.25 | 0.31 | 0.35 | 0.75 | 0.58 | 0.30 | 0.33 | 0.40 | 0.77 |
ME-NB-global | 0.61 | 0.27 | 0.34 | 0.38 | 0.77 | 0.61 | 0.32 | 0.36 | 0.42 | 0.78 |
ME-NB-local | 0.62 | 0.25 | 0.33 | 0.35 | 0.76 | 0.61 | 0.30 | 0.34 | 0.40 | 0.77 |
ME-NB-random | 0.59 | 0.24 | 0.31 | 0.35 | 0.75 | 0.59 | 0.30 | 0.33 | 0.40 | 0.77 |
LR | 0.62 | 0.18 | 0.28 | 0.29 | 0.76 | 0.63 | 0.23 | 0.31 | 0.34 | 0.77 |
ME-LR-global | 0.60 | 0.23 | 0.31 | 0.34 | 0.77 | 0.61 | 0.27 | 0.33 | 0.38 | 0.78 |
Experimental results on the DNA-protein sequence data set. Experimental results with Naive Bayes (NB) and Logistic Regression (LR) models, and Mixture of Experts (ME) models on the non-redundant DNA-protein sequence data set, where the identity cutoffs are 30% and 90%. The results are shown for default threshold θ = 0.5. ME-NB-global and ME-LR-global use NB and LR at the leaves and exploits the global sequence similarity to construct the hierarchical structure. ME-NB-local exploits the local sequence similarity to construct the hierarchical structure. ME-NB-random randomizes the global similarity matrix and constructs the hierarchical structure based on the randomized matrix.
Classifier | DNA-protein 30% | DNA-protein 90% | ||||||||
---|---|---|---|---|---|---|---|---|---|---|
Precision | Recall | CC | FM | AUC | Precision | Recall | CC | FM | AUC | |
NB | 0.59 | 0.05 | 0.16 | 0.10 | 0.75 | 0.56 | 0.07 | 0.18 | 0.13 | 0.75 |
ME-NB-global | 0.62 | 0.12 | 0.25 | 0.20 | 0.77 | 0.65 | 0.15 | 0.29 | 0.25 | 0.78 |
ME-NB-local | 0.65 | 0.06 | 0.18 | 0.12 | 0.76 | 0.64 | 0.08 | 0.21 | 0.15 | 0.76 |
ME-NB-random | 0.58 | 0.05 | 0.15 | 0.09 | 0.75 | 0.56 | 0.07 | 0.18 | 0.13 | 0.75 |
LR | 0.57 | 0.07 | 0.18 | 0.12 | 0.79 | 0.57 | 0.08 | 0.18 | 0.14 | 0.79 |
ME-LR-global | 0.57 | 0.14 | 0.26 | 0.23 | 0.80 | 0.63 | 0.17 | 0.29 | 0.26 | 0.81 |
The mixture of experts that exploits the global similarity between protein sequences outperforms a mixture of experts that exploits the local similarity between protein sequences
In order to verify that indeed global sequence similarity is helpful in improving the performance of classifiers, and that the improvement does not come from the more sophisticated structure of the model, we computed the entries in the similarity matrix W by applying Smith-Waterman local alignment algorithm with Blosum62, thus taking into account local sequence similarity (the matrix W is normalized and scaled as before). We also randomized the global similarity matrix computed previously and use this randomized matrix to construct the hierarchical structure of the mixture of experts models.
In Tables 1 and 2 we show the performance of NB and mixture of NB models using global (ME-NB-global) and local (ME-NB-local) sequence similarities, as well as a random (ME-NB-random) sequence similarity for the default threshold θ = 0.5. The results of our experiments show that the mixture of experts models that capture global sequence similarity outperform the other models in terms of the majority of standard measures for comparing the performance of classifiers used in this study (the results are similar for the mixture of LR models, data not shown). For example, for 30% identity cutoff, Correlation Coefficient increases from 0.33 (local similarity) to 0.34 (global similarity) on the RNA-protein data set (Table 1), and from 0.18 (local similarity) to 0.25 (global similarity) on the DNA-protein data set (Table 2). Hence, we conclude that global similarity is helpful in improving the performance of classifiers trained to label biomolecular sequence data.
The mixture of experts has consistently higher performance than the baseline classifier for all identity cutoffs
The mixture of experts that exploits the global sequence similarity offers a higher precision than the ensemble of classifiers for the same Recall
We also trained ensembles of NB and LR classifiers on both RNA- and DNA-protein interface prediction tasks to predict whether or not a residue in a protein sequence is an interface residue. An ensemble of classifiers [7, 8] is simply a collection of classifiers, each trained on a balanced subsample of the training data. The prediction of the ensemble is computed from the predictions of the individual classifiers (see Methods section for further details).
We compared the performance of the mixtures of NB and LR with that of ensembles of NB and LR, respectively. In Figures 1, 2, 3, and 4 we show the PR curves for the mixture and the ensemble models on both RNA- and DNA-protein sequence data sets for 30% identity cutoff. As can be seen from the figures, the mixtures of experts consistently offer a higher Precision than the ensembles of classifiers for the same Recall. Note that the PR curves of the ensembles are closer to those of the baseline classifiers.
Discussion
Reliable methods for identifying putative functional sites in protein sequences is an important problem with broad applications in computational biology, e.g., rational drug design. Computational tools for identifying functional sites from sequences are especially important because of the cost and efforts involved in structure determination.
In this work we sought to improve the performance of classifiers that make predictions on residues in protein sequences by taking into account the global similarity between the protein sequences in the data set in addition to the local features extracted around each residue. We evaluated mixture of experts models that consider the global similarity between protein sequences when building the model and making the predictions on the RNA-protein and DNA-protein interface prediction tasks. Two closely related models are the Hierarchical Mixture of Experts model [9] and the ensemble of classifiers model [7]
Hierarchical Mixture of Experts
The Hierarchical Mixture of Experts model (HME) was first proposed by Jordan and Jacobs (1994) [9] to solve nonlinear classification and regression problems by combining linear models: the input space is divided into a set of nested regions and simple (e.g., linear) models are fit to the data that fall in these regions. Hence, instead of using a "hard" partitioning of the data, the authors use a "soft" partitioning, i.e., the data is allowed to simultaneously lie in more than one region.
The HME has a tree-structured architecture that is known a priori. The internal nodes of the tree correspond to gating networks and the leaf nodes correspond to expert networks. The expert networks output class probabilities for each input x, while the gating networks learn how to combine the predictions of the experts up the tree with the final prediction output by the root. The parameters of the gating networks are learned using Expectation Maximization algorithm [10]. The gating and the expert networks are generalized linear models.
Ensemble of classifiers
An ensemble of classifiers is a collection of independent classifiers, each classifier being trained on a subsample of the training data [7]. The prediction of the ensemble of classifiers is computed from the predictions of the individual classifiers using majority voting. An example is misclassified by the ensemble if a majority of the classifiers misclassifies it. When the errors made by the individual classifiers are uncorrelated, the predictions of the ensemble of classifiers are often more reliable.
Mixture of experts – our approach
Our approach to learning a mixture of experts model takes into account the global similarity between biomolecular sequences in a data set. Unlike the HME model [9], we assume that the structure of our model is not known a priori. Hence, to learn the hierarchical structure of the model, we use hierarchical clustering of the sequences in the data set. The leaf nodes consist of expert classifiers, while the gating nodes combine the output of each classifier to the root of the tree which makes the final prediction. The gating nodes combine the predictions of the expert classifiers based on an estimate of the cluster membership of a test protein sequence. Following the approach taken by Jordan and Jacobs [9], we considered a "soft" partitioning of the data, i.e., each sequence in the training set simultaneously lies in all clusters of the hierarchical structure with a different weight in each cluster. The combination scheme of the predictions of the expert classifiers and the "soft" partitioning of the data that considers the global sequence similarity differentiate our model from an ensemble of classifiers model.
Conclusion
Identification of functionally important sites in biomolecular sequences has broad applications ranging from rational drug design to the analysis of metabolic and signal transduction networks. With the rapid increase in the amount of data (e.g., protein sequences) there is a growing need for reliable procedures to accurately identify such sites.
In this study, we have presented a mixture of experts approach to identification of functionally important sites from amino acid sequence of proteins that takes into account global similarity between the protein sequences. Specifically, we systematically evaluated Naive Bayes and Logistic Regression classifiers, as well as mixtures of Naive Bayes and Logistic Regression in a sequence-based 10-fold cross-validation setup. The results of our experiments show that global sequence similarity through the means of the mixture of experts approach can be exploited to improve the performance of classifiers trained to label biomolecular sequence data.
Methods
Data sets and parameter settings
We used two datasets to perform experiments: RNA-protein and DNA-protein interface data sets that are available online at http://www.cs.iastate.edu/~cornelia/rna_dna. RNA- and DNA-protein interactions play a pivotal role in protein function. Reliable identification of such interaction sites from protein sequences has broad applications ranging from rational drug design to the analysis of metabolic and signal transduction networks.
The RNA- and DNA-protein interface data sets consist of RNA- and DNA-binding protein sequences, respectively, extracted from structures in the Protein Data Bank (PDB) [11]. We downloaded all the protein structures of known RNA- and DNA-protein complexes from PDB solved by X-ray crystallography and having X-ray resolution between 0 and 3.5 Å. As of May 2008, the number of RNA-protein complexes was 435 and DNA-protein complexes was 1259. A residue was identified as interface residue using Entangle with the default parameters [12].
Furthermore, to remove redundancy in each data set, we used BlastClust, a toolkit that clusters sequences with statistically significant matches, available at http://toolkit.tuebingen.mpg.de/blastclust[13]. While constructing our non-redundant sequence data sets, we applied various identity cutoffs, starting from 30% and ending at 90% in steps of 10. For example, in the 30% identity cutoff sequence data set, two sequences were pairwise matched if they were 30% or more identical over an area covering 90% of the length of each sequence. We randomly selected a sequence from each cluster returned by BlastClust. Thus, the resulting non-redundant RNA-protein sequence data set for 30% identity cutoff has 180 protein sequences. The total number of amino acid residues is 33,235.
We represented residues identified as interface residues in a protein sequence as positive instances (+) and those not identified as interface residues as negative instances (-). Furthermore, we encoded each residue by a local window of fixed length, winLength = 21, corresponding to the target residue and ten neighboring residues on each side.
Number of sequences, as well as positive and negative instances used in our experiments for the RNA- and DNA-protein data sets. Number of sequences as well as number of positive (+) and negative (-) instances in the non-redundant RNA- and DNA-protein sequence data sets for 30%, 60%, and 90% identity cutoffs.
Data Sets | Number of Sequences | Number of + Instances | Number of - Instances |
---|---|---|---|
RNA-prot 30% | 180 | 5398 | 27837 |
RNA-prot 60% | 215 | 6689 | 32073 |
RNA-prot 90% | 246 | 7798 | 34675 |
DNA-prot 30% | 257 | 5326 | 53494 |
DNA-prot 60% | 289 | 5974 | 58031 |
DNA-prot 90% | 317 | 6551 | 60877 |
Learning mixture of experts models
Here we present our approach to learning a mixture of experts model that takes into account the global similarity between biomolecular sequences. Unlike the Hierarchical Mixture of Experts model [9], we assume that the structure of our model is not known a priori. Hence, to learn the hierarchical structure of the model, we use hierarchical clustering of the sequences in the data set. The leaf nodes consist of expert classifiers, while the gating nodes combine the output of each classifier to the root of the tree which makes the final prediction. The gating nodes combine the predictions of the expert classifiers based on an estimate of the cluster membership of a test protein sequence. Similar to Jordan and Jacobs [9], we considered a "soft" partitioning of the data, i.e., each sequence in the training set simultaneously lies in all clusters of the hierarchical structure with a different weight in each cluster.
Learning the structure of the mixture of experts model
We first compute the pairwise similarity matrix W_{n × n}for the protein sequences in the training set based on a common global sequence alignment method. Second, using this similarity matrix, we apply 2-way spectral clustering algorithm, described in the next subsection, to recursively bipartition the training set of protein sequences until a splitting criterion is met.
The output of the algorithm is a hierarchical clustering of the protein sequences, i.e., a tree $\mathcal{T}$ such that each node (cluster) consists of a subset of sequences. The root node is the largest cluster containing all the protein sequences in the training set. Once a cluster is partitioned into its two subclusters, it becomes their parent in the resulting tree structure. We store all the intermediate clusters computed by the algorithm. If the number of sequences at a given cluster falls below some percentage of the total sequences in the training set, then the node becomes a leaf and thus is not further partitioned (we used 10% in our experiments).
Figure 9 shows the tree structure produced by the 2-way spectral clustering algorithm when applied to a set of 147 RNA-protein sequences. The similarity matrix is computed based on the Needleman-Wunsch global alignment algorithm. In the figure, to keep the tree smaller, we stopped bipartitioning a node when the number of sequences at a given cluster falls below 30% of the total sequences in the training set.
2-Way spectral clustering
Spectral clustering has been successfully applied in many applications, including image segmentation [15], document clustering [16], grouping related proteins according to their structural SCOP classification [17].
Spectral clustering falls within the category of graph partitioning algorithms that partition the data into disjoint clusters by exploiting the eigenstructure of a similarity matrix. In general, to find an optimal graph partitioning is NP complete. Shi and Malik [15] proposed an approximate spectral clustering algorithm that optimizes the normalized cut (NCut) objective function. It is a divisive, hierarchical clustering algorithm that recursively bi-partitions the graph until some criterion is reached, producing a tree structure.
- 1.
Let W_{n × n}= [S(i, j)] be the symmetrical matrix containing the similarity score for each pair of sequences.
- 2.
Let D_{n × n}be the degree matrix of W, i.e., a diagonal matrix such that D_{ ii }= ∑_{ j }S(i, j).
- 3.
Solve the eigenvalue system (D - W)x = λ D x for the eigenvector corresponding to the second smallest eigenvalue and use it to bipartition the graph.
- 4.
Recursively bipartition each subgraph obtained at Step 3. if necessary.
Note that the quality of the clusters found by the 2-way spectral clustering algorithm depends on the choice of the similarity function S.
Estimating the parameters of the mixture of experts model
Following the approach taken by Jordan and Jacobs [9], we make use of the "soft" partitioning of the biomolecular sequence data. Thus, having the hierarchical clustering $\mathcal{T}$ stored, we devise a procedure that allows each sequence in the training set to simultaneously lie in all clusters, with a different weigth in each cluster.
- 1.
Find the K closest sequences to x_{ i }at the parent node based on the similarity function used to construct the hierarchical clustering $\mathcal{T}$ (in our experiments we used K equal to 20% of the sequences at the parent node).
- 2.
Let K_{0} out of K sequences go to the left child node, and K_{1} out of K go to the right child node.
- 3.
The estimated probability of x_{ i }for being in child node j is computed as p(x_{ i }∈ V_{ j }|x_{ i }∈ par(V_{ j })) = K_{ j }/K, where j = 0, 1.
We recursively place the sequence x_{ i }in all the nodes of $\mathcal{T}$ with different weights, starting from the root, based on its estimated cluster membership computed above. Thus, the sequence weight at the root is 1 (all the sequences in the training set lie at the root of the tree), and the weight at any other node in the tree is the product of the sequence weights on the path from the root to that node.
Let ${V}_{1}^{l},{V}_{2}^{l},\cdots ,{V}_{M}^{l}$ be the leaf nodes and ${V}_{1}^{g},{V}_{2}^{g},\cdots ,{V}_{N}^{g}$ be the internal or gating nodes in the hierarchical clustering $\mathcal{T}$. During learning, we train either a collection of M Naïve Bayes classifiers or a collection of M Logistic Regression classifiers, one classifier at each leaf node ${V}_{k}^{l}$, k = 1,⋯, M. Naïve Bayes and Logistic Regression are briefly described in the next section.
To solve the biomolecular sequence labeling problem, one approach is to predict each element x_{ i, j }in the sequence x_{ i }independently, i.e., to assume that the observation-label pairs (x_{ i, j }, y_{ i, j })_{j = 1, m}are independent of each other (the label independence assumption). However, x_{ i, j }may not contain all the information necessary to predict y_{ i, j }. Hence, it is fairly common to encode each element x_{ i, j }in the sequence x_{ i }based on a local, fixed-length window corresponding to the target element and its sequence context (an equal number of its sequence neighbors on each side) ${{x}^{\prime}}_{i,j}$ = x_{i, j-t},⋯, x_{ i, j },⋯, x_{i, j+t}. The classifier is trained to label the target element x_{ i, j }[6].
Machine learning classifiers
Naïve Bayes
Naïve Bayes (NB) [18] is a supervised learning algorithm that belongs to the class of generative models, in which the probabilities p(x|y) and p(y) of the input x and the class label y are estimated from the training data using maximum likelyhood estimates. Typically, the input x is high-dimensional, represented as a set of features (attributes), x = (x_{1}, x_{2}, ⋯, x_{ d }), making it impossible to estimate p(x|y) for large values of d.
Therefore, training a Naïve Bayes classifier reduces to estimating probabilities p(x_{ i }|y), i = 1,⋯, d, and p(y), from the training data, for all class labels y.
The class label with the highest posterior probability is assigned to the new input x_{ test }.
Logistic Regression
where [β, θ] are the parameters of the model that can be estimated either by maximizing the conditional likelihood on the training data or by minimizing the loss function.
During classification, Logistic Regression predicts a new input x_{ test }as 1 if and only if
β^{ T }x_{ test }+ θ > 0
Ensemble of classifiers
In our experiments, we used m = 300. Each individual classifier in the collection was trained on approximately $\frac{l}{10}$ instances, where l represents the total number of training instances available to the ensemble.
The implementation of all the models considered in this study is built on Weka, an open source machine learning software available at http://www.cs.waikato.ac.nz/ml/weka/[20].
Performance evaluation
To obtain the estimates for TP, FN, FP and TN, we performed 10-fold sequence-based cross-validation [21] wherein the set of sequences is partitioned into 10 disjoint subsets (folds). At each run of a cross-validation experiment, 9 subsets are used for training and the remaining one is used for testing the classifier. The values for TP, FN, FP and TN are obtained using the default threshold θ = 0.5, i.e., an instance is classified as positive if the probability of being in the positive class returned by the classifier is greater than or equal to 0.5, and as negative otherwise.
With any classifier, it is possible to tradeoff the Precision against Recall. Hence, it is more informative to compare the Precision-Recall curves which show the tradeoff over their entire range of possible values than to compare the performance of the classifiers for a particular choice of the tradeoff.
The Precision-Recall curve is a good indicator of the performance of classifiers when the data sets are highly unbalanced, as is the case with our both RNA- and DNA-protein data sets [22]. It has also been shown that if a curve dominates in PR space, it also dominates in ROC space [22].
To evaluate how good a classifier is at discriminating between the positive and negative examples, we also report the Area Under the Receiver Operating Characteristic Curve (AUC) on the test set, which represents the probability of correct classification [23].
Declarations
Acknowledgements
This research has been suppported in part by a grant from the National Institutes of Health (GM066387) to VH and DD.
This article has been published as part of BMC Bioinformatics Volume 10 Supplement 4, 2009: Proceedings of the IEEE International Conference on Bioinformatics and Biomedicine (BIBM) 2008. The full contents of the supplement are available online at http://www.biomedcentral.com/1471-2105/10?issue=S4.
Authors’ Affiliations
References
- Terribilini M, Lee JH, Yan C, Jernigan RL, Honavar V, Dobbs D: Predicting RNA-binding Sites from Amino Acid Sequence. RNA Journal. 2006, 12: 1450-1462.View ArticleGoogle Scholar
- Yan C, Dobbs D, Honavar V: A Two-Stage Classifier for Identification of Protein-Protein Interface Residues. Bioinformatics. 2004, 20 (Suppl 1): i371-i378.View ArticlePubMedGoogle Scholar
- Qian N, Sejnowski T: Predicting the secondary structure of globular proteins using neural networks models. J Mol Biol. 1988, 202: 865-884.View ArticlePubMedGoogle Scholar
- Caragea C, Sinapov J, Silvescu A, Dobbs D, Honavar V: Glycosylation site prediction using ensembles of Support Vector Machine classifiers. BMC Bioinformatics. 2007, 8: 438-PubMed CentralView ArticlePubMedGoogle Scholar
- Kim JH, Lee J, Oh B, Kimm K, Koh I: Prediction of phosphorylation sites using SVMs. Bioinformatics. 2004, 20 (17): 3179-3184.View ArticlePubMedGoogle Scholar
- Diettrich TG: Machine Learning for Sequential Data: A Review. Proceedings Joint IAPR International Workshop on Structural, Syntactic, and Statistical Pattern Recognition. 2002, 15-30.View ArticleGoogle Scholar
- Dietterich TG: Ensemble Methods in Machine Learning. Lecture Notes in Computer Science. 2000, 1857: 1-15.View ArticleGoogle Scholar
- Russell S, Norvig P: Artificial Intelligence: A Modern Approach. 2003, Prentice HallGoogle Scholar
- Jordan MI, Jacobs RA: Hierarchical mixtures of experts and the EM algorithm. Neural Computation. 1994, 6: 181-214.View ArticleGoogle Scholar
- Dempster AP, Laird NM, Rubin DB: Maximum likelihood from incomplete data via the EM algorithm. Journal of the Royal Statistical Society. 1977, 39: 1-38.Google Scholar
- Berman H, Westbrook J, Feng Z, Gilliland G, Bhat T, Weissig H, Shindyalov I, Bourne P: The Protein Data Bank. Nucleic Acid Res. 2000, 28: 235-242.PubMed CentralView ArticlePubMedGoogle Scholar
- Allers J, Shamoo Y: Structure-based analysis of protein-RNA interactions using the program ENTANGLE. J mol Biol. 2001, 311: 75-86.View ArticlePubMedGoogle Scholar
- Using BLASTClust to Make Non-redundant Sequence Sets. [http://toolkit.tuebingen.mpg.de/blastclust]
- Duda R, Hart E, Stork D: Pattern Classification. 2001, Wiley, SecondGoogle Scholar
- Shi J, Malik J: Normalized cuts and image segmentation. Pattern Analysis and Machine Intelligence. 2000, 22 (8): 888-905.View ArticleGoogle Scholar
- Dhillon IS: Co-clustering documents and words using bipartite spectral graph partitioning. Proceedings of SIGKDD International Conference on Knowledge Discovery and Data Mining. 2001, 269-274.Google Scholar
- Paccanaro A, Casbon JA, Saqi MAS: Spectral clustering of protein sequences. Nucleic Acids Research. 2006, 34 (5): 1571-1580.PubMed CentralView ArticlePubMedGoogle Scholar
- Mitchell TM: Machine Learning. 1997, McGraw HillGoogle Scholar
- Ng AY, Jordan MI: On discriminative vs. generative classifiers: A comparison of logistic regression and naive Bayes. Advances in Neural Information Processing Systems (NIPS), NIPS. 2002Google Scholar
- Weka 3: Data Mining Software in Java. [http://www.cs.waikato.ac.nz/ml/weka/]
- Caragea C, Sinapov J, Dobbs D, Honavar V: Assessing the Performance of Macromolecular Sequence Classifiers. IEEE 7th International Symposium on Bioinformatics and Bioengineering. 2007, 320-326.Google Scholar
- Davis J, Goadrich M: The Relationship Between Precision-Recall and ROC Curves. Proceedings of the 23rd International Conference on Machine Learning. 2006, 233-240.Google Scholar
- Baldi P, Brunak S, Chauvin Y, Andersen C, Nielsen H: Assessing the accuracy of prediction algorithms for classification: an overview. Bioinformatics. 2000, 16 (5): 412-424.View ArticlePubMedGoogle Scholar
Copyright
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.