Supervised redundant feature detection for tumor classification
© Zeng and Li; licensee BioMed Central Ltd. 2014
Published: 22 October 2014
As a high dimensional problem, analysis of microarray data sets is a challenging task, where many weakly relevant or redundant features affect overall performance of classifiers.
The previous works used redundant feature detection methods to select discriminative compact gene set, which only considered the relationship among features, not the redundancy of classification ability among features. This study propose a novel algorithm named RESI (Redundant fEature Selection depending on Instance), which considers label information in the measure of feature subset redundancy.
Experimental results on benchmark data sets show that RESI performs better than the previous state-of-the-art algorithms on redundant feature selection methods like mRMR.
We propose an effective supervised redundant feature detection method for tumor classification.
Rapid advances in gene expression microarray technology enable simultaneous measurement of the expression levels for thousands or tens of thousands of genes in a single experiment. Analysis of microarray data presents unprecedented opportunities and challenges for data mining in areas such as gene clustering, class discovery, and tumor classification . A typical data set may contain thousands of genes but only a small number of samples (often less than a hundred). The number of samples is likely to remain small at least for the near future due to the expense of collecting microarray samples . The nature of relatively high dimensionality but small sample size in microarray data causes the known problem of "curse of dimensionality". Therefore, selecting a small number of discriminative ones from thousands of genes is essential for successful sample classification [3–5].
Feature selection, a process of choosing a subset of features from the original, is frequently used as a preprocessing technique in data mining. It has been proved effective in reducing dimensionality, improving mining efficiency, increasing mining accuracy, and enhancing result comprehensibility [6–9]. The goal of the feature selection algorithm is to select the minimum set of features that are strongly related to the desired decision variable and have the least redundancy among them [10, 11]. Existing efficient feature selection algorithms (feature ranking methods) usually assume feature independence, and assign a discriminative score to each feature. Because the interactions and correlations among features are not considered, these algorithms fail to remove redundant features completely.
Hall  pointed out that the prediction capability and the inter-correlation of feature subset are two important aspects in feature selection. However, the optimal feature selection requires an exponentially large search space (O(2 m ), where m is the number of features) . In practise, researchers often resort to various approximations to determine the optimal feature subset. The prediction capability is usually estimated by the interrelation of individual feature with the target. For the issue of redundancy, many researchers attempt to explicitly address feature interactions by finding some low-order interactions, i.e. 2-way (one feature and the label) and 3-way (two features and the label) interactions.
In recent years, researchers have proposed many techniques to eliminate redundant features according to the above principle. Battiti proposed the Mutual Information based Feature Selection (MIFS) method , and then improved versions like MIFS-U  and mMIFS-U  were proposed. Hall proposed the Correlation-based Feature Selection (CFS) method , Ding and Peng proposed the minimum Redundancy-Maximum Relevance (mRMR) method , Bontempi and Meyer proposed a Causal filter selection method, called min-Interaction Max-Relevance (mIMR) , Fleuret used the Conditional Mutual Information Maximization (CMIM) to select feature . These methods consider the problem of feature selection from different motivations and various solutions have been proposed. But the same point is that they are trying to find the optimal discriminative feature subset by considering to remove feature redundancy, where feature redundancy is computed by various pair-wise similarity measures, i.e. mutual information or conditional mutual information.
However, the used traditional pair-wise similarity metrics only consider the numerical values of given variables, but not the similarity of discriminative ability among them. Therefore, feature redundancy can not be measured correctly in terms of feature prediction ability by existing methods. How to measure feature (subset) redundancy is a fundamental problem, which should be reconsidered in the field of feature selection. In the application of classification, we consider a feature redundant only when its predictive power is redundant, not that its numerical value is similar with some selected features. Because two highly similar features are obviously not redundant to each other when the minor difference happen to be critical to the classification. In this paper, we proposed a novel method to measure feature (subset) redundancy by comparing features' predictive powers directly. Feature's prediction power is recorded by its instances' distribution explicitly, which includes clear-discerned instances and blur-discerned instances. Based on the predictive power distributions, a new feature redundancy metric is defined on the ground of comparing predictive powers. Furthermore, we proposed a novel algorithm named RESI (Redundant fEature Selection depending on Instance). Experiments on several benchmark microarray data sets demonstrate the outstanding performance of our proposed method.
Supervised redundant feature detection
Given a data set (X = [X1, X2, . . . , X m ], Y) with n instances represented by m features (genes), where X and Y are instances matrix and class label respectively. The task of classification is to tag instances with a label in low probability of error. Theoretically, having more features implies more discriminative power in classification. However, many features are relevant to each other and they have no contribution to classification, except to degrade performance of classifiers . Thus, it is necessary to remove irrelevant and redundant features as more as possible, without losing information greatly. Generally, the task of feature selection in classification issue is formally defined as the process of selecting the optimal feature subset S ⊆ X to have high classification performance.
Thus, qualities of relevant and non-redundant features of the selected subset are two key points considered in feature selection methods. In existing methods, predictive power of a feature subset has often been divided into computation of the interrelation score of individual features with the label. Due to computational limitation in high dimensional feature space, redundancy degree of a feature subset is usually estimated by pair-wise similarity measures, i.e. mutual information or conditional mutual information of two given variables.
However, to measure the feature subset redundancy is a fundamental problem, which should be reconsidered in the field of feature selection. Without including learning model (i.e. wrapper methods), measuring the redundancy degree of a feature subset as a whole is impractical. So, approximately measuring individual pair-wise feature similarity is a reasonable way which has been adopted by various existing methods. But, the previous pair-wise similarity measures only consider the numerical values of given features, but not the discriminative ability among them. Therefore, feature redundancy is not measured correctly. For instance, two highly similar features are usually considered as redundant to each other by the previous pair-wise similarity measures. But this is not correct when the minor difference of the two features happens to be critical to the classification.
In the context of classification, it is reasonable to address the feature set redundancy from the view of predictive power. So, we define the pair-wise feature redundancy as follow.
Theorem 1 In the context of classification, a feature Xp is redundant to a given feature Xq if its predictive power to label Y has already been expressed by feature Xq.
As argued in Theorem 1, the similarity should be measured by comparing the distribution of predictive power between two features. Thus, the predictive power needs to be measured in a comparable way. More concretely, we need to record those instances which are clear-discerned by the given feature and which are not.
Theorem 2 Given a feature Xp, an instance Xu is clear-discerned if the majority of its k nearest neighbour instances on feature Xp have the same class labels with that of Xu.
where δ(·) is the Kronecker delta function i.e. δ(·) = 1 when Y v == Y u is true, otherwise δ(·) = 0, Y u is the target label of instance X u , is the value of instance X u on feature X p , N k () is the neighborhood of X u defined by the k closest instances X v on feature X p ( is close to ), k is a predefined parameters. Equation 1 measures how many instance labels are the same in the neighborhood of X u via feature X p .
The value of varies from 0 to 1. High means the corresponding instance has more neighbors with the same class. An instance is defined as clear-discerned when its is higher than a given threshold µ, i.e. , otherwise it is blur-discerned.
The numerator of Equation 2 gives the count of instances which are blur-discerned on both feature X p and feature X q , the denominator is the blur-discerned instances which count on feature X p . The value of REMI(X p ; X q |Y) varies from 0 to 1. It is 0, when all blur-discerned instances on feature X p are clear-discerned on feature X q , which means the predictive power of X q is completely complementary to that of X p and we consider their corresponding redundancy is zero. The value is 1 when no blur-discerned instance on feature X p is clear-discerned on feature X q , which means X q has no predictive power contribution to X p , i.e. X q is redundant to X p .
where , n1 and var1 are the mean value of features, the number of examples and variance of one class, , n2 and var2 have the similar meaning for the other class.
The meaning of t-statistic value is intuitive, which measures the weighted distance between the centroid of one class and the other class on the feature p. It is usually believed that the value of t(p|Y) largely represents the discriminative power of the giving feature, higher absolute score means greater discriminative power. A feature is usually regarded as irrelevant when its t-statistic value is trivial. The t-statistic measure has demonstrated as an effective feature selection method and widely used in the field of bioinformatics .
which denotes the merit score of candidate feature X p given the selected features S (S = X1:p−1) and the label Y. Furthermore, we propose a novel algorithm named RESI (Redundant fEature Selection depending on Instance), which combines REMI with the sequential forward search strategy. The detail algorithm is described as in Algorithm 1.
Algorithm 1 The RESI Algorithm
Input: Feature set X = [X1, X2, . . . , X m ] Target label Y
Output: Selected feature subset S
1: F ⇐ X;
2: S ⇐ Ø;
3: for all feature X i in F do
4: if abs(t(X i |Y)) < ∈ then
5: F ⇐ F \ X i ;
6: end if
7: end for
9: S ⇐ S ∪ X p ;
10: F ⇐ F \ X p ;
11: while pre-defined stopping criteria is not satified do
13: S ⇐ S ∪ X p ;
14: F ⇐ F \ X p ;
15: end while
The RESI algorithm works in two stages, including irrelevant features removing and redundant feature elimination. Firstly, the t-statistic value is computed for each feature. The irrelevant features, whose absolute t-statistic values are trival i.e. abs(t(X i |Y)) < ∈, are removed from the candidate set. Then, the most discriminative features with the highest abs(t(X p |Y)) are selected. In the stage of redundant feature elimination, RESI starts with the current selected subset, and adds one important feature at a time. Given a selected feature subset of p−1 features S , a new feature X p is chosen from the rest of the feature subset X \ S by evaluating the merit function J REMI (X p |S , Y). It will be terminated when pre-defined stopping criteria is satisfied. For example, the number of selected features is larger than a threshold or the score of J REMI (X p |S, Y) is trivial.
In our experiments, the parameter ∈ = 0.1, the number of neighborhood instances k = 3, the threshold of clear-discerned instance µ = 0.66.
Related works on redundant feature detection
Related feature selection methods
Merit function of candidate feature
I(X i ; Y)
I(X p ; X i )
I(X p ; Y)
I(X p ; X i )
I(X p ; Y)
I(X p ; X i )
I(X p ; Y)
I(X p ; Y)
I(X i ; X p ; Y)
I(X p ; Y)
I(X i ; X p ; Y)
I(X p ; Y) − maxi ∈[1:p−1][I(X p ; X i ) − I(X p ; X i |Y)]
Hall proposed the Correlation-based Feature Selection (CFS) method , where some correlation measures are used to evaluate the goodness of a subset by considering the individual predictive ability of each feature and the degree of correlation between them. The symmetrical uncertainty (a normalized version of mutual information), for discrete data, and the standard linear correlation, for continuous data, are used by Hall to measure I(X i ; Y) and I(X i ; X j ).
Ding and Peng proposed the minimum Redundancy-Maximum Relevance (mRMR) method in 2005 , which requires that selected discriminative features are maximally dissimilar to each other. mRMR is almost as the same as the MIFS method , except that the parameter β is set as in mRMR. Both MIFS and mRMR use mutual information to measure I(X i ; Y) and I(X i ; X j ). Kwak and Choi proposed an improvement to MIFS, called MIFS-U , which uses a re-weighted mutual information to measure the feature redundancy. Without explicit claim, I(·;· ) is measured by the mutual information in the rest of paper.
Bontempi and Meyer proposed a causal filter selection method, called min-Interaction Max-Relevance (mIMR) . Bontempi and Meyer try to maximize the mutual information between X1:pand Y directly. Due to the number of the subset of X1:pis O(2 p ), it is impractical to compute I(X1:p; Y) in a precise way. Only low-order interactions are considered in the approximate solution given by mIMR. And in the final solution, the merit function of mIMR can also be divided into two parts: predictive power and feature set redundancy. Obviously, the difference between mRMR and mIMR is the two-way mutual information I(X i ; X p ) is replaced by the three-way mutual information I(X i ; X p ; Y), and I(X i ; X p ; Y) = I(X p ; X i ) − I(X p ; X i |Y).
Fleuret used the Conditional Mutual Information Maximization (CMIM) to select feature , which examines the information between a feature and the target, conditioned on each current feature. It is clear that CMIM is very similar with mIMR. The only difference is the sum function is replaced by the maximum function to measure the feature redundancy.
Yang and Moody proposed using Joint Mutual Information (JMI) to select feature , which tries to maximize joint mutual information . This is the information between the targets and a joint random variable, defined by pairing the candidate X p with each current selected feature. But after deduction, the merit function of JMI is exactly equivalent to that of mIMR.
Vidal-Naquet and Ullman proposed another criterion used for computer vision, which is refereed as Informative Fragments (IF) . The authors motivate the criterion mini∈[1:p−1][I(X p , X i ; Y) − I(X i ; Y)] by noting that it measures the predictive ability gain of combining a new feature X p with each existing feature X i , over simply using X i by itself. The X i with the least "gain" from being paired with X p is taken as the score for X p . Interestingly, using the chain rule I(X p , X i ; Y) = I(X i ; Y) + I(X p ; Y|X i ), therefore IF is equivalent to CMIM.
From the short review, we can easily find the common point of these methods. Although the motivations are various, the merit functions of all these methods are divided into two parts: predictive power and feature redundancy. And due to the practical limitation, one certain pair-wise similarity measure is adopted to compute the predictive power and feature redundancy. However, the traditional pair-wise similarity measures, i.e. mutual information or conditional mutual information, only consider the numerical values of given variables, but not the the similarity of discriminative ability between them. Therefore, feature redundancy can not be measured correctly in existing methods.
Experimental data sets
We use the stratified 10-fold cross-validation procedure, where each data set is split into ten subsets of equal size. Each subset is used as test set once, and the corresponding left subsets are combined together and used as training set. Within each cross-validation fold, the gene expression data are standardized. The expressions of the training set are transformed to zero mean and unit standard deviation across samples, and the test set are transformed according to the means and standard deviations of the corresponding training set. The Irani's MDL method is applied when discretization is required . The 10-fold cross-validation is repeated 10 times, which is also denoted as the 10 × 10 cross-validation measuring procedure.
We should note that the 10 × 10 cross-validation measuring procedure is more reliable than the randomized re-sampling testing strategy and the leave-one-out cross-validation due to the correlations between the test and training sets, some detailed discussions can be found at  Even in the small sample problem like gene expression data, 10 × 10 cross-validation is still one of the most reliable measuring way .
where TP, TN, FP, and FN, stand for the number of true positive, true negative, false positive, and false negative samples, respectively. Without explicit clarification, all the scores are averaged on 10 × 10 cross-validation.
To make conclusions sound, six widely used classifiers are used, including Support Vector Machine (SVM) with linear kernel and c = 1, non-linear support vector machine using Sequential Minimal Optimization (SMO) and polynomial kernel, k Nearest Neighbor (kNN) with k = 3, Logistic Regression (LR), Naïve Bayes classifier (NB) and decision tree with J48 algorithm (J48). All these classifiers are trained on the training set to predict the label of the test samples on the same cross-validation partition.
The algorithms are implemented in JAVA language based on WEKA , and carried out on a DELL PC workstation with 24 × X5680 3.33GHz CPU and 64G RAM.
Results and discussion
From Figure 1 and Figure 2, it's clear there is no method has overwhelming performance over all classifiers, although the plotted lines in Figure 2 fluctuate much more gently than the drawings on individual data set. However, our RESI is relatively the best one among all the feature selection methods. This is especially obvious by using the classifier of LR, and we will demonstrate the point by t-test later.
The performances of mIMR, mRMR and CMIM are often tied, and mRMR is the best one in most cases. We suggest it's because all these methods use mutual information to represent the predictive power, and the traditional mutual information of mRMR is the most robust one to estimate feature redundancy in our experiments.
It's also interesting to find that the performances of RESI and Ranking are obviously better than that of others, when the selected dimension is small, i.e. below 5. This is true with almost all classifiers. We believe the superiority is caused by the difference of two metrics, mutual information and the absolute t-statistic, used to represent the predictive power. The absolute t-statistic measures the weighted distance between the centroid of two classes, which is directly related to the discriminative power of the giving feature. On the contrary, mutual information represents the general information between two variables. So from our experimental results, the absolute t-statistic performs better than mutual information to represent the predictive power.
Investigating the difference between RESI and Ranking, we find that their performances are similar when the selected dimension is small. This is because the absolute t-statistic is used to represent the predictive power by both of them. And when only a few features are selected, the influence of feature redundancy is not obvious. But as the selected dimension grows, the impact of feature redundancy becomes more and more critical to the BACC performance. Eliminating redundant features is meaningful when the feature dimension is not small, which is coincided with Figure 2, where RESI is much better than Ranking when the dimension grows.
Comparative t-test results (wins/ties/losses) summarized on data sets and classifiers
Redundant feature selection is an important topic in the field of bioinformatics. Here, we propose a novel redundant feature subset measure REMI by comparing feature predictive powers directly, which is recorded by its instance distribution explicitly including clear-discerned instances and blur-discerned instances. Furthermore, a novel feature selection method RESI based on REMI was proposed. Experimental results on benchmark microarray data sets demonstrate that RESI performs better than the state-of-the-art algorithms like mRMR on fourteen benchmark data sets.
Future works include improving its efficiency and applying it to more scientific fields.
This work was supported by the Natural Science Foundation of China under grant nos. 61105053, 61273305 and 61463033, the Nature Science Foundation of Jiangxi Province under grant no. 20132BAB201043, and China Postdoctoral Science Foundation under grant no. 2013M541540.
Funding for open access charge: Natural Science Foundation of China.
This article has been published as part of BMC Medical Genomics Volume 7 Supplement 2, 2014: IEEE International Conference on Bioinformatics and Biomedicine (BIBM 2013): Bioinformatics in Medical Genomics. The full contents of the supplement are available online at http://www.biomedcentral.com/bmcmedgenomics/supplements/7/S2.
- Ge L, Du N, Zhang A: Finding informative genes from multiple microarray experiments: A graph-based consensus maximization model. Bioinformatics and Biomedicine (BIBM), 2011 IEEE International Conference on. 2011, 506-511. 10.1109/BIBM.2011.34.View ArticleGoogle Scholar
- Dougherty ER: Small sample issue for microarray-based classification. Comparative and Functional Genomics. 2001, 2: 28-34. 10.1002/cfg.62.PubMed CentralView ArticlePubMedGoogle Scholar
- Zheng S, Liu W: Selecting informative genes by lasso and dantzig selector for linear classifiers. Bioinformatics and Biomedicine (BIBM), 2010 IEEE International Conference on. 2010, 677-680. 10.1109/BIBM.2010.5706651.View ArticleGoogle Scholar
- Zeng XQ, Li GZ: Incremental partial least squares analysis of big streaming data. Pattern Recognition. 2014, 47 (11): 3726-3735. 10.1016/j.patcog.2014.05.022.View ArticleGoogle Scholar
- You M, Li GZ: Feature selection for multi-class problems by using pairwise-class and all-class techniques. International Journal of General Systems. 2011, 40 (4): 381-394. 10.1080/03081079.2010.530027.View ArticleGoogle Scholar
- Li GZ, Meng HH, Yang MQ, Yang JY: Combining support vector regression with feature selection for multivariate calibration. Neural Computing and Applications. 2009, 18 (7): 813-820. 10.1007/s00521-008-0202-6.View ArticleGoogle Scholar
- Shao H, Li G, Liu G, Wang Y: Symptom selection for multi-label data of inquiry diagnosis in traditional chinese medicine. Science China Information Sciences. 2013, 56 (5): 1-13.View ArticleGoogle Scholar
- Li GZ, Sun S, You M, Wang YL, Liu GP: Inquiry diagnosis of coronary heart disease in chinese medicine based on symptom-syndrome interactions. Chinese Medicine. 2012, 7 (1): 9-10.1186/1749-8546-7-9.PubMed CentralView ArticlePubMedGoogle Scholar
- You M, Liu J, Li GZ, Chen Y: Embedded feature selection for multi-label classification of music emotions. International Journal of Computational Intelligence Systems. 2012, 5 (4): 668-678. 10.1080/18756891.2012.718113.View ArticleGoogle Scholar
- Guyon I, Elisseeff A: An introduction to variable and feature selection. The Journal of Machine Learning Research. 2003, 3: 1157-1182.Google Scholar
- Zeng XQ, Li GZ, Yang JY, Yang MQ, Wu GF: Dimension reduction with redundant gene elimination for tumor classification. BMC Bioinformatics. 2008, 9 (Suppl 6): 8-10.1186/1471-2105-9-S6-S8.View ArticleGoogle Scholar
- Hall MA: Correlation-based feature selection for discrete and numeric class machine learning. Proceedings of the 2000 International Conference on Machine Learning (ICML'00). 2000, 359-366.Google Scholar
- Battiti R: Using mutual information for selecting features in supervised neural net learning. IEEE Transactions on Neural Networks. 1994, 5 (4): 537-550. 10.1109/72.298224.View ArticlePubMedGoogle Scholar
- Kwak N, Choi CH: Input feature selection for classification problems. IEEE Transactions on Neural Networks. 2002, 13 (1): 143-159. 10.1109/72.977291.View ArticlePubMedGoogle Scholar
- Novovičová J, Somol P, Haindl M, Pudil P: Conditional mutual information based feature selection for classification task. Progress in Pattern Recognition, Image Analysis and Applications. 2007, Springer, 4756: 417-426.View ArticleGoogle Scholar
- Peng H, Long F, Ding C: Feature selection based on mutual information criteria of max-dependency, max-relevance, and min-redundancy. IEEE Transactions on Pattern Analysis and Machine Intelligence. 2005, 27 (8): 1226-1238.View ArticlePubMedGoogle Scholar
- Bontempi G, Meyer P: Causal filter selection in microarray data. Proceedings of the 2010 International Conference on Machine Learning (ICML'10). 2010, 95-102.Google Scholar
- Fleuret F: Fast binary feature selection with conditional mutual information. The Journal of Machine Learning Research. 2004, 5: 1531-1555.Google Scholar
- Zeng XQ, Li GZ: Dimension reduction for p53 protein recognition by using incremental partial least squares. IEEE Transactions on NanoBioscience. 2014, 13 (2): 73-79.View ArticlePubMedGoogle Scholar
- Liu H, Li J, Wong L: A comparative study on feature selection and classification methods using gene expression profiles and proteomic patterns. Genome informatics. 2002, 13: 51-60.PubMedGoogle Scholar
- Brown G: A new perspective for information theoretic feature selection. Proceedings of the 2009 International Conference on Artificial Intelligence and Statistics (ICAIS'09). 2009, 5: 49-56.Google Scholar
- Liu H, Liu L, Zhang H: Feature selection using mutual information: An experimental study. Proceedings of the 2008 Pacific Rim International Conference on Artificial Intelligence (PRICAI'08). 2008, 235-246.View ArticleGoogle Scholar
- Bolón-Canedo V, Sánchez-Maroño N, Alonso-Betanzos A: A review of feature selection methods on synthetic data. Knowledge and Information Systems. 2013, 34 (3): 483-519. 10.1007/s10115-012-0487-8.View ArticleGoogle Scholar
- Yang HH, Moody J: Data visualization and feature selection: New algorithms for nongaussian data. Advances in Neural Information Processing Systems. 1999, 12: 687-693.Google Scholar
- Vidal-Naquet M, Ullman S: Object recognition with informative features and linear classification. Proceedings of the 2003 IEEE International Conference on Computer Vision (ICCV'03). 2003, 281-288.Google Scholar
- Guyon I, Li J, Mader T, Pletscher PA, Schneider G, Uhr M: Competitive baseline methods set new standards for the NIPS 2003 feature selection benchmark. Pattern Recognition Letters. 2007, 28 (12): 1438-1444. 10.1016/j.patrec.2007.02.014.View ArticleGoogle Scholar
- West M, Blanchette C, Dressman H, Huang E, Ishida S, Spang R, Zuzan H, Olson JA, Marks JR, Nevins JR: Predicting the clinical status of human breast cancer by using gene expression profiles. Proceedings of the 2001 National Academy of Sciences. 2001, 11462-11467.Google Scholar
- Li J, Liu H, Wong L: Mean-entropy discretized features are effective for classifying high-dimensional biomedical data. Proceedings of the 3rd ACM SIGKDD Workshop on Data Mining. 2003, 17-24.Google Scholar
- Fayyad U, Irani K: Multi-interval discretization of continuous-valued attributes for classification learning. Proceedings of the 1993 International Joint Conference on Articial Intelligence (IJCAI'93). 1993, 1022-1027.Google Scholar
- Dietterich TG: Approximate statistical tests for comparing supervised classification learning algorithms. Neural Computation. 1998, 10: 1895-1923. 10.1162/089976698300017197.View ArticlePubMedGoogle Scholar
- Molinaro AM, Simon R, Pfeiffer RM: Prediction error estimation: a comparison of resampling methods. Bioinformatics. 2005, 21 (15): 3301-3307. 10.1093/bioinformatics/bti499.View ArticlePubMedGoogle Scholar
- Witten IH, Frank E: Data Mining: Practical Machine Learning Tools and Techniques. Morgan Kaufmann Publishers Inc., San Francisco, 2Google Scholar
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. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.