Table of Contents Author Guidelines Submit a Manuscript
The Scientific World Journal
Volume 2014, Article ID 314728, 20 pages
Research Article

Self-Adaptive MOEA Feature Selection for Classification of Bankruptcy Prediction Data

1Institute of Polymers and Composites-I3N, University of Minho, Guimarães, Portugal
2Department of Computer Science, Universidad Carlos III de Madrid, Leganes, Madrid, Spain
3Department of Production and Systems Engineering, University of Minho, Braga, Portugal

Received 1 October 2013; Accepted 25 December 2013; Published 23 February 2014

Academic Editors: Z. Cui and X. Yang

Copyright © 2014 A. Gaspar-Cunha et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


Bankruptcy prediction is a vast area of finance and accounting whose importance lies in the relevance for creditors and investors in evaluating the likelihood of getting into bankrupt. As companies become complex, they develop sophisticated schemes to hide their real situation. In turn, making an estimation of the credit risks associated with counterparts or predicting bankruptcy becomes harder. Evolutionary algorithms have shown to be an excellent tool to deal with complex problems in finances and economics where a large number of irrelevant features are involved. This paper provides a methodology for feature selection in classification of bankruptcy data sets using an evolutionary multiobjective approach that simultaneously minimise the number of features and maximise the classifier quality measure (e.g., accuracy). The proposed methodology makes use of self-adaptation by applying the feature selection algorithm while simultaneously optimising the parameters of the classifier used. The methodology was applied to four different sets of data. The obtained results showed the utility of using the self-adaptation of the classifier.

1. Introduction

Bankruptcy prediction has become an important economic phenomenon [1, 2]. The high individual, economical, and social costs arising from bankruptcies have motivated further effort in understanding the problem and finding better prediction methods. In finances, bankruptcy prediction is an important topic of research as it provides a way of identifying business failure, that is, situations in which a firm or particular cannot pay lenders, preferred stock shareholders, suppliers, and so forth. An organisation which is unable to meet its scheduled payments when estimations of future cash show that the current financial situation will not change in the near future is said to undergo into financial distress. Signs of financial distress are evident long before bankruptcy occurs. Research in bankruptcy prediction started in [3] where a univariate discriminant model was used. This was followed by studies using traditional statistical methods which include correlation, regression, logistic models, and factor analysis [4, 5]. More recently, an overview of the classic statistical divided them into four types: univariate analysis, risk index models, multivariate discriminant analysis, and conditional probability models [6].

Modern bankruptcy prediction models combine both statistical analysis and artificial intelligence techniques improving then the decision support tools and decision making [79]. In this manner, back propagation artificial neural networks have been applied to bankruptcy prediction [10] whose results revealed better accuracy than predictions made using some other techniques (recursive portioning, -nearest neighbours, C4.5, etc.). Consequently, research has focused on the combination of artificial neural networks with other soft computing tools such as fuzzy sets, genetic programming, ant colony optimisation, or particle swarm optimisation [1114].

Support vector machines (SVMs) have been largely used for classification and pattern recognition applications. SVMs are a family of generalised linear classifiers widely used for classification of financial data. In particular, several studies have been published on the application of SVMs to the problem of bankruptcy prediction [1518]. A survey on support vector machines applied to the problem of bankruptcy prediction can be found in [19]. It is worth mentioning that support vector machines require solving a quadratic programming problem which is time consuming when considering large dimensional problems and also that it requires the optimisation of algorithm parameters which may affect its performance. The aim behind this research is to overcome the above limitations which will be accomplished by using feature selection and self-adaptation of the classification algorithm parameters.

Feature selection can be described as one of the initial stages of a classification process by which the complexity of the problem is reduced by elimination of irrelevant features [20]. Feature selection must be approached with the minor lose of information of the original set after the noisy or irrelevant features are removed; that is, the elimination of irrelevant features should not reduce the overall classification accuracy. Being the original set of features for a given classification task, the continuous feature selection problem consists in assigning weights to each feature in such a way that the order corresponding to its theoretical relevance is preserved. In a similar way, the binary feature selection problem refers to the assignment of binary weights that leads to a reduced subset of features (with ). In the general case, all features take part in the learning process, each one with a particular contribution. In binary feature selection, only a subset of the features is considered in the learning process for which all of them contribute in the same manner. For the purpose of this work, binary feature selection will be used. In [21], the problem of binary feature selection was formally defined, which, for the general case, consists in finding a compromise between minimising the number of features in and maximising an evaluation measure over the subset . Notice that an optimal subset of features is not necessarily unique which has motivated further research into this field. Also, there are many potential benefits of feature selection [22], that is, facilitating data visualisation and understanding, reducing the measurement and storage requirements, reducing training and using times, and so forth. Traditional feature selection methods used in bankruptcy prediction consist on applying statistical methods, such as -test, correlation matrix, stepwise regression, principal component analysis, or factor analysis to examine their prediction performance [23]. The application of artificial intelligence techniques, such as evolutionary computation, to the problem of feature selection is now emerging in order to enhance the effectiveness of traditional methods [20].

The general case for feature selection fits into a multi-objective optimisation approach where the aim is to simultaneously optimise two or more conflicting objectives. In addition, identifying a set of solutions representing the best possible trade-offs among objectives of the problem instead of a single solution might be of interest in many cases. Within this context, evolutionary algorithms constitute a preferred choice as they simultaneously deal with a set of solutions, referred to as population, which allows several different solutions to be generated in a single run. Several evolutionary multi-objective approaches (MOEAs) have been applied to finances and economics. The most popular application of MOEAs in the literature deals with the portfolio optimisation problem [2426], although MOEAs have also been successfully applied to stock ranking [27], risk-return analysis [28, 29], and economic modelling [30, 31]. In a sense, this work constitutes a study on the consequences of simultaneously optimised two or three objective functions over real-world benchmark problems.

Another issue that will be considered in this work is the self-adaptation of the classifier algorithm parameters. Self-adaptation aims at finding suitable adjustment of the algorithm parameters efficiently [32]. In general, the definition of self-adaptation in evolutionary algorithms refers to the adjustment of control parameters that are related to evolutionary routines [33], that is, mutation or crossover rates, population size, and selection strategy. In this work, the scope of this definition will be modified and the aim will be the automatic adjustment of the classification process parameters, which in the present case include the training method, the training fraction, and the specific SVM parameters (e.g., kernel and regularisation parameters). Some other recent works that might be of interest for the reader are [3438].

The aim of this work is to further investigate into the feature selection problem in bankruptcy prediction using a multi-objective approach, including self-adaptation of the classification algorithm parameters. This work is expected to contribute by introducing a novel multi-objective methodology for feature selection which provides a solution to the problem of bankruptcy prediction compromising both the minimisation of the number of features selected and the maximisation/minimisation of a quality measure of the classifier, for example, accuracy or error. Also, this paper will help to create a better understanding of the application of SVMs to real-world data. The proposed methodology will be validated using bankruptcy prediction datasets found in the literature.

The remaining of this paper is organised as follows. The proposed methodology will be described in detail in Section 2, for which the corresponding expertise areas of classification, using SVMs, feature selection in classification and multi-objective evolutionary optimisation will be introduced. Section 3 describes the datasets used during the experimental part of this research. Discussion on the performance of the algorithm will follow in Section 4. The paper finalises in Section 5 by pointing out the main contributions, limitations and further extensions to this work.

2. Multiobjective Feature Selection

2.1. Feature Selection

As stated above, the feature selection problem consists in finding the minimum number of features that are necessary to evaluate correctly a set of data. Considering as the original set of features with a cardinality , the following definition applies [39].

Definition 1 (feature selection). Let be an evaluation measure to be optimised (considering here a maximisation problem, without loss of generality) defined as . The selection of a feature subset can be seen under three considerations.(i)Set . Find , such that is maximum.(ii)Set a value , that is, the minimum that is going to be tolerated. Find the with smaller , such that .(iii)Find a compromise among minimising and maximising (general case).

In the present work, a wrapper approach was used [40]. Usually, the existing data is divided in two sets, the training and the test data. For that purpose, the existence of (i) a representative set of data, capable of allowing the identification of the relations between the features and the classification of such data, (ii) an algorithm able to classify the data accurately (classification algorithm), (iii) and an optimisation algorithm able to find the best set (or the minimum) of features that classify the data with the best accuracy and/or the minimum error is necessary.

Figure 1 illustrates the well-known confusion matrix, for a situation with two classes. TP (true positives) are the positive instances correctly classified, TN (true negatives) are the negative instances correctly classified, FP (false positives) are the positive instances incorrectly classified, and FN (false negatives) are the negative instances incorrectly classified. Based on this taxonomy, different measures can be defined to quantify the accuracy and the error achieved by the classifier as follows: where is the accuracy, is the recall or sensitivity, is the precision, and are the classification errors of types I and II, respectively, and is the harmonic mean of the sensitivity () and precision (). After the above formalism, the problem consists in maximising , , , and and minimising the errors. There are other type of classification measures that can be also applied. However, the problem to be addressed is the simultaneous optimisation of some of these measures. For example, in bankruptcy prediction, the maximisation of the profits, but, simultaneously, the minimisation of losses is desired. In the present situation the profits can be quantified by recall (), since it is a direct measure of the positives correctly classified (TP), that is, the companies that test well and are healthy, and the losses can be quantified by the error of type I (), a measure of the positives incorrectly classified (FP), that is, the companies that test well but actually are in bankruptcy. The trade-off between and is known as the Receiver Operating Characteristics (ROC) curve [41, 42]. Figure 2 illustrates this concept. The ideal point is identified by “1’’ and means a perfect classification.

Figure 1: Confusion matrix.
Figure 2: ROC curve.

The above example illustrates the importance of optimising more than one objective simultaneously. In fact, in the case of feature selection, the first objective to be optimised (minimised) is the number of features that are necessary to get an accurate classification, which can be taken into account by maximising or , for example, but also by obtaining the best trade-off between and , as illustrated in Figure 2. In the first case, there are two objectives to be satisfied simultaneously, while, in the second, three objectives are considered. Therefore, the use of a multi-objective optimisation algorithm together with an accurate classifier is of primordial importance.

2.2. Support Vector Machines

There are available in the literature a large number of algorithms/methods for classification of data. For example, the WEKA software offers a great number of different methods ready to be used in an straightforward way [43]. A good survey about the best classification algorithms can be found in [44].

The method adopted here is the Support vector machines (SVMs). SVMs are a set of supervised learning methods based on the use of a kernel, which can be applied to classification and regression [45]. In the SVM, a hyperplane or set of hyperplanes are constructed in a high-dimensional space. The initial step consists in transforming the data points, through the use of nonlinear mapping, into the high-dimensional space. In this case, a good separation is achieved by the hyperplane that has the largest distance to the nearest training data points of any class. Thus, the larger this margin, the smaller the generalisation error of the classifier. SVMs can be seen as an extension to nonlinear models of the generalised portrait algorithm developed in [46]. For the purpose of this work, the SVM package from LIBSVM was used [47].

The SVMs depend on the definition of some important parameters. First, it is necessary to select the the type of kernel. In the present work, the Radial Basis Function (RBF) kernel was adopted due to its efficiency. Then, it becomes necessary to select the SVM type which depends on its usage, that is, if it is used for classification or regression. Since this work deals with classification, the -SVC and -SVC methods were selected. Both, the kernel and the type of SVM, depend on the value defined for some parameters that must be carefully set, the kernel parameter and the regularisation parameter that depends on the type of SVM chosen ( and ). Finally, some other parameters were studied including the training method and the training fraction. Two different training methods were tested, the holdout method, where a fraction (training fraction) of the instances are used to train the SVM and the remaining are used for testing and the -fold method, that consists in dividing the set of instances in subsets. Then subsets are used to train the SVM and the remaining set is used for validation. The process is repeated times, accounting for all subsets used for validation, and the accuracy is obtained as the average of the training/testing steps [47].

Due to the large number of parameters that must be set before applying the optimisation algorithm, it makes sense to apply the feature selection algorithm and the optimisation of these parameters simultaneously. This is what is done in the present work. Therefore, the following parameters were optimised simultaneously with the process of feature selection: training method (holdout, H; or 10-fold, , validation), training fraction (TF), kernel , and regularisation parameters ( or ). More details about the implementation of this strategy are given in the next subsection.

2.3. Multiobjective Evolutionary Algorithm

In order to deal with multiple objectives multiobjective optimisation algorithms (MOOA) must accomplish two basic functions simultaneously: (i) they need to guide the population towards the optimal Pareto set. This can be done by using a fitness assignment operator that takes into account the non-dominance concept. (ii) The nondominated set must be maintained as diverse as possible; that is, the solutions must be well distributed along the entire optimal Pareto front. Additionally, it is also necessary to maintain an archive of the best solutions found during the various generations in order to prevent some nondominated solutions from being lost. Therefore, generally in MOEAs, it is only necessary to replace the selection phase of a traditional EA by a routine able to deal with multiple objectives [48, 49].

In this work, the MOEA adopted is the reduced Pareto set genetic algorithm (RPSGA) [50]. However, any other multi-objective algorithm can be used for the same purpose. The main steps of this algorithm are described below (Algorithm 1). The algorithm starts by the random creation of an internal population of size and an empty external population of size . Then, at each generation (i.e., while the stopping criteria are not met), the following operations are performed: (i) the internal population is evaluated using the SVM routine; (ii) a clustering technique is applied to reduce the number of solutions on the efficient frontier and to calculate the ranking of the individuals of the internal population; (iii) the fitness of the individuals is calculated using a ranking function; (iv) a fixed number of the best individuals is copied to the external population; (v) if the external population is not totally full, the genetic operators of selection, crossover, and mutation are applied to the internal population to generate a new population; (vi) when the external population becomes full, the clustering technique is applied to sort the individuals of the external population, and a predefined number of the best individuals is incorporated in the internal population by replacing lowest fitness individuals.

Algorithm 1: Reduced Pareto set genetic algorithm (RPSGA).

Detailed information about this algorithm can be found in [50, 51]. The influence of some important parameters of the algorithm, such as size of internal and external populations, number of individuals copied to the external population in each generation and from the external to the internal population, and the limits of the indifference of the clustering technique, had already been studied and the best values have been suggested [50].

2.4. Methodology for Feature Selection

The linkage between the problem to solve (the selection of features), the SVM, and the MOEA is done as follows. During the generation of the initial population, the chromosome (generated randomly) is constituted by a binary string identifying if the corresponding feature is present (value equal to 1) or not (value equal to 0) and the values of the classification algorithm/process (TF, H or , and , or ), which are used for self-adaptation. These chromosome values are then passed to the SVM during the evaluation of the population. The SVM returns the achieved values of accuracy and errors (1) obtained with the selected features and parameter values that are present in the chromosome of each individual.

The RPSGA algorithm was adapted to deal with the above feature selection problem. With respect to the classifier parameters, two approaches were considered. Initially, a pure feature selection problem was analysed where these parameters were not allowed to vary after being set up at the beginning of the algorithm. In a second approach, these parameters were included in the chromosome as variables to be optimised. The latter approach has the advantage of obtaining in a single run the best features and, simultaneously, fine tuning the classifier parameters (self-adaptation). Each candidate solution generated by the RPSGA will be externally evaluated by the SVM whose result will be returned to the RPSGA to be used as fitness in the genetic routine. New solutions will be generated based on the performance of the previous generation. As usual, the fittest solutions have more possibilities of survival.

3. Datasets

In the present study, the four datasets presented below will be used to validate the proposed methodology. Note that the DIANE data consists of two datasets from different years.

3.1. Industrial French Companies Data

In the present work, two samples, from the years 2005 and 2006, respectively, obtained from the DIANE database were selected. The original database comprised financial ratios of about 60 000 industrial French companies with at least 10 employees. The dataset includes information about 30 financial ratios of companies covering a wide range of industrial sectors (see Table 1). Since the original database contained many instances with missing values, especially, concerning defaults companies, the default cases were sorted by the number of missing values and only samples with less than 10 missing values were selected. A final set of 600 default examples was obtained. In order to create a balanced dataset, 600 random nondefault examples were selected and added to the dataset, thus resulting in a set of 1200 examples. Similar preprocessing of this dataset can be found in [31, 52, 53].

Table 1: Set of features considered for the industrial French companies.
3.2. German Credit Data

The German Credit database was created at the University of Hamburg and is publicly accessible at the UCI Machine Learning Repository [54]. It consists of 1000 instances of credit applications which are described by the 20 attributes shown in Table 2. Examples of previous usage of the German Credit dataset can be found in [55, 56]. There are two versions of the German dataset available, the original German Credit dataset which consists of numerical and nominal attributes and its numeric version produced at the Strathclyde University. As the method proposed in this paper only accepts numerical attributes, the numeric version of the data will be used.

Table 2: Set of features considered for the German Credit.
3.3. Australian Credit Data

The Australian Credit database originates from [57] and concerns data form 690 credit card applications. The data are publicly available in the UCI Machine Learning Repository [54]. Each instance consists of 14 attributes and one of two possible classes (all attribute names and values were changed to meaningless symbols to protect the confidentiality of the data). The class distribution is similar for both, 44.5% versus 55.5%. Examples of previous usage of this dataset can be found in [58].

3.4. Data Normalisation

In general, a large amount of data is available and often these data are inconsistent and redundant being necessary considerable manipulation to make it useful for problems like credit risk analysis. It becomes important to identify the ratios or ranges of data that are relevant to the problem. Restricting the data to the relevant ranges represents an advantage to reduce the complexity of the problem.

Due to the large diversity of data concerning the type of data (e.g., real or integer values, numeric or categorical) and the range of variation of the values for each feature, some normalisation of the data becomes necessary. Therefore, the data was transformed as follows:(1)logarithmic transformation: (2)centering and standardizing the data: (3)normalisation of the data in the interval : where represents the instance, stands for feature, is the original data in a matrix form (which is transformed successively in and ), AVG and STD are the average and the standard deviation of all instances for feature , respectively, and is the final value used by the classifier. The data used by the classifier is restricted to the interval as recommended in [44].

4. Results and Discussion

4.1. Computational Experiments

Table 3 presents the set of experiments carried out to test the proposed approach. Due to the stochastic nature of the evolutionary algorithm, 10 runs for each experiment were performed, using different seed values (as required by the random number generator). In the case of Experiment 1, the C-SVC method was used with the following fixed parameters: holdout (H) validation as training method, , training fraction, , equal to 0.7, kernel parameter, , equal to 0.1, and the regularisation parameter, , equal to 10. In the remaining experiments, these parameters were allowed to range in the following intervals: , , , , and . In Table 3, represents the maximum number of features allowed in the initial generation, that is, if is equal to 5 means that in the initial generation the individuals of the population have at the most 5 features. In consecutive generations, the number of selected features was allowed to grow until the maximum of features for each database is reached: for French industrial companies, subscript F, ; for German Credit data, subscript G, ; and for Australian Credit data, subscript A, . Besides, Figures 1 and 2 should not become a problem (with respect to the dataset dimension) for standard SVMs experimentation; this work tries to demonstrate that feature selection is useful for the application of SVMs over datasets of high dimension.

Table 3: Set of computational experiments (N1 is the maximum number of features in the initial population; X means experiment not done).

The aim of Experiment 1 is to compare the performance of the feature selection method proposed when the classifier parameters are fixed to that of the same method when the parameters are allowed to vary. This will be done by comparing Experiments 1, 2, and 12. Experiments 2 to 7 are thought to illustrate the influence of the method when different classification measures are applied. In the case of Experiments 8 to 11, the aim is to study the influence of the maximum number of features of the initial population () in the evolution of ROC curves (i.e., versus ). Finally, Experiments 12, 13, and 14 were intended to show the influence of the SVM method used. In all runs, the following RPSGA parameters were used (see [50] for more details): the main and elitist population sizes were 100 and 200 individuals, respectively; fitness proportional selection was adopted; crossover rate of 0.8 and mutation probability of 0.05 were used; the number of ranks was set to 30 and the limit of indifference of the clustering technique was set to 0.01, whereas the number of generations was set to 100 for all runs.

4.2. Analysis of a Standard Experiment

This section is aimed at showing the type of results that can be obtained using the proposed methodology. For that purpose, Figure 3 shows the entire initial population and the nondominated solutions corresponding to generations 25, 50, 75, and 100 for Run-a of Experiment 2. This graph presents the trade-off between Acc (to be maximised) and NF (to be minimised). It can be easily observed that the algorithm is able to evolve the population significantly, from the initial population (randomly generated), located predominantly at the bottom right corner, towards the top left corner. It is also noticeable that only 50 generations are needed to reach a reasonable approximation of the Pareto front. The use of 100 generations was only used to guarantee the convergence of the algorithm.

Figure 3: Run-a of Experiment (initial population and nondominated solutions of the final population).

Table 4 shows the obtained results corresponding to the decision variable domain for the above run after 100 generations. The accuracy is ranged between 76.3% and 83.5%, when considering a minimum number of 3 features and a maximum of 17, respectively. In all cases, the holdout (H) cross validation training method was selected and the training fraction lies around 52% and is ranged between 0.13 and 0.55, whereas fluctuate between 10 and 211. This indicates that decision variables (, , , and ) converge for a small interval when compared to the initial range where they are allowed to vary.

Table 4: Optimal solutions for Run-a of Experiment 2 (Diane 2005 database).

However, the target consists in finding better solutions than those obtained over a single run. Figure 4 shows the optimal Pareto curves of the 10 runs that were performed for Experiment 2. It can be seen that there is one of these runs that dominates the others, Run-f, except when , where the best solution is obtained for Run-e. Table 5 shows the decision variable values of the corresponding Pareto front, for which is ranged between 75.6% and 85.8%, the obtained is hold out for all cases, and the lies around 75%. On the other hand, the SVM parameters have a large variation which indicates that and play an important role in acquiring best accuracies. Similar conclusions can be drawn when analysing the results obtained using the remaining datasets.

Table 5: Optimal solutions for Experiment 2 (Diane 2005 database).
Figure 4: All runs of Experiment 2 (nondominated solutions of the final population).
4.3. Analysis and Comparison of Results

Figures 5 and 6 represent the nondominated solutions of the 10 different runs carried out in Experiments 2 to 7 using to the French industrial companies in 2005 dataset. These plots allow to assess the efficiency of the proposed optimisation methodology when dealing with all the objective function measures presented in Section 2. As expected, and since the common objective used in these experiments is the minimisation of , the solutions evolve nicely towards the region where the true Pareto front is supposed to be; that is, when simultaneously maximising a second objective (e.g., , , , and ) the solutions evolve towards the top left corner, while when simultaneously minimising a second objective (e.g., and ), the solutions evolve towards the bottom left corner.

Figure 5: Optimal Pareto fronts for Diane 2005 data (10 runs).
Figure 6: Optimal Pareto fronts for Diane 2005 data (10 runs).

Further analysis of Figures 5 and 6 helps to identify the ranges that can be accomplished when using the different objective functions (for the French datasets): and , , , , and . However, when considering the best values in a particular run, the following values were found: , , , , , and , corresponding to NF equal to 10, 11, 5, 13, 3, and 13, respectively. Considering a given number of features, for example, , the following best values are found: , , , , , and . On the other hand, when considering all ten runs of each experiment, the variation range for each objective function can be graphically observed. Such a variation enforces the use of several runs with different seed values in order to select the best set of features as well as the best classifier parameters. Since the final accuracy will depend certainly on the combination of the right features, the methodology adopted cannot be based on selecting the features that appear more frequently in the 10 runs performed for each experiment [59].

The above reasoning was used to select the best solution of the front when comparing the results from Experiments 1, 2, and 12 over all datasets studied. Note that Experiments 1, 2, and 12 consist on simultaneously optimising NF and (see Figures 7, 8, 9, and 10). Furthermore, the above analysis allowed to create Table 6 which summarises the solutions found for three different cases: solutions with best accuracy (Best) and best solutions using only 5 (NF ≤ 5) and 10 (NF ≤ 10) features, respectively.

Table 6: Results summary for Experiments 1, 2, and 12 and for all databases.
Figure 7: Optimal Pareto fronts for Diane 2005 data (10 runs).
Figure 8: Optimal Pareto fronts for Diane 2006 data (10 runs).
Figure 9: Optimal Pareto fronts for German Credit data (10 runs).
Figure 10: Optimal Pareto fronts for Australian Credit data (10 runs).

As expected and in general, the results of Table 6 show that the best accuracy is obtained when the classifier parameters are also optimised (Experiments 2 and 12). Concerning the use of the C-SVC or the -SVC kernels, no definitive conclusion can be drawn, since the C-SVC kernel yielded the best result for Diane05, whereas the -SVC kernel yielded the best result the the Australian data and for some other cases the best result depends on the number of features (Diane06 and German data). With respect to the runs where the “best” results were obtained for each of the three conditions that were analysed, again there is some variability; in some cases, the results were taken from the same run but in most of the cases they were not. Again, this fact was expected after the analysis made in the previous section. In all cases, the holdout validation method is selected, TF ranges between 70% and 80% in most cases (except in the case of the German database), and the kernel and regularisation parameters have a high variability to maximise the accuracy. This was also expected after the analysis of the previous section.

The analysis or results show that the desired accuracy can be achieved using several combinations of features. Results coming from the same run tend to select the same features (this fact was also observed in the results presented in Tables 4 and 5). An interesting finding came from Experiment 2 over Diane05 database; it was observed that when the number of features was reduced to 5 at the most (), four out of five of the features selected were identical to one of the features that were selected for the best solution condition (features 7, 13, 21, and 23), but the last feature selected when using this constraint was not included in the best solution (feature 29). Many valuable information can be obtained from Table 6. As an example, if the problem consists on obtaining the best accuracy using five features at the most (), the features identified in bold should be selected to be used in future classifications together with their corresponding parameter for each dataset considered.

Figures 11 and 12 show the best results achieved in Experiments 8, 9, 10, 11, and 14. Note that these experiments consist in optimising three different objectives (, , and ) and were aimed to obtain the results that best fit in a ROC curve; that is, . Besides the optimisation that was carried out considering all three objectives, only nondominated solution with respect to objectives and are presented (best of 10 runs for each experiment). Table 7 shows a summary of results from the above experiments for all databases using two different conditions ( and ). The area under ROC was computed at first for all cases and then best results were presented for each condition. Identical conclusions, to that of the beginning of this section, Section 4.3, can be made here concerning the algorithm parameters, that is, best kernel (which depend on the database), best validation method, training fraction, and kernel and regularisation parameters. Similarly, there exist various combinations of features that allow the obtention of the best and values. As before, the best solutions using five features at the most, , can be selected for each database. Such features are identified in bold in Table 7 and can be used in future classification together with their corresponding classifier parameters.

Table 7: Results summary for Experiments 8, 9, 10, 11, and 14 and for all databases.
Figure 11: ROC curves for Diane 05 data (best of 10 runs).
Figure 12: ROC curves for Diane 05 data (best of 10 runs).

In [60], clustering feature selection methods were used to identify the most relevant features on several datasets. The Australian Credit dataset was used to test three versions of a clustering based algorithm with different optimisation strategies. The structure of clusters, found by the optimisation version of the algorithm proposed in the above paper, indicates the subset of three relevant features for classification, features 8, 9, and 14. Almost all solutions obtained in the experiments carried out in this paper (see Tables 6 and 7) using the Australian data include features 8 and 9. Feature 14 is also present in most of the solutions obtained in this work.

Prediction of financial distress of companies using Diane dataset, was previously analysed using several machine learning approaches [61]. Support vector machines achieved the highest accuracy considering five features selected by SVM attribute evaluation method. The five features selected for predicting failures during 2007 using historical data from 2006, 2005, and 2004 were 1, 4, 11, 16, and 28 (which differs from the solution obtained here). However, it should be noted that these results were obtained using historical data and, therefore, they are not comparable to the results obtained in this work.

An approach to solving classification problems by combining feature selection and neural networks was proposed in [62]. A feature selection algorithm based on the notion of entropy from the information theory was applied to the German Credit dataset yielding the selection of the following seven features: 1, 2, 3, 5, 6, 15, and 20. Authors found that the predictive accuracy was marginally larger with the exclusion of the 13 redundant features. Most of the solutions obtained by the approach presented in this work include some of these features (in particular, features 1, 2, 3, and 5). However, it should be noted that the experimental set-up of the two studies is rather different and, therefore, conclusions must be drawn carefully.

5. Conclusion

With the current global economic situation where several countries are getting through economic recession, bankruptcy prediction is acquiring importance as a financial topic of research. When the financial data to be analysed becomes large, the need for feature selection arises as a tool used to reduce both computational times and number of computations by getting rid of irrelevant features. Feature selection also gives a method to evaluate the importance of each feature within the studied dataset.

This work aimed at investigating the feature selection problem in bankruptcy prediction using a multi-objective approach which includes self-adaptation of the classification algorithm parameters. For that purpose, a new methodology has been proposed and its performance has been evaluated using real-world benchmark problem datasets for bankruptcy prediction. A large set of experiments using different objective functions, such as accuracies, error, and sensitivity measures, have been performed which provides a better understanding on the application of SVMs to real-world data. The performance of the proposed method was also studied using two- and three-objective approaches.

Results have shown that the method performs well using the benchmark datasets studied. Large accuracies have been obtained using a significantly reduced subset of features. Consequently, the more the considered features, the lager the accuracies. Also, being a multi-objective technique, instead of a single solution, a set of nondominated solutions is provided which may help the decision maker to evaluate the trade-off in making a sacrifice in one of the objective functions towards obtaining a gain in some others. The inability for the classifier to handle nominal features within the data turned out to be the main limitation of the proposed method. This limitation was inherent to the classifier used by the method; it was overcome by converting nominal attributes of the data to numerical.

A possible extension to this work could be made by taking advantage of the multi-objective nature of the set of solutions and analysing in detail the trade-off between them, thus helping decision makers to choose the preferred solution for their needs. The proposed method could also be extended to work with many objectives as real-world situations actually do.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.


This work was partially supported by the Portuguese Foundation for Science and Technology under Grant PEst-C/CTM/LA0025/2011 (Strategic Project-LA 25-2011-2012) and by the Spanish Ministerio de Ciencia e Innovación, under the project “Gestión de movilidad efficiente y sostenible, MOVES’’ with Grant Reference TIN2011-28336.


  1. E. L. Altman, Corporate Financial Distress and Bankruptcy, John Willey & Son, 1993.
  2. C. Pate, Bankrupties: Recent Level and Implications for 2002, chapter 11, Pricewaterhouse Coopers, 2002.
  3. W. Beaver, “Alternative accounting measures as predictors of failure,” The Accounting Review, vol. 1, no. 1, pp. 113–122, 1968. View at Google Scholar
  4. D. Martin, “Early warning of bank failure, a logit regression approach,” Journal of Banking and Finance, vol. 1, no. 3, pp. 249–276, 1977. View at Google Scholar · View at Scopus
  5. J. A. Ohlson, “Financial ratios and the probabilistic prediction of bankruptcy,” Journal of Accounting Research, vol. 18, pp. 109–131, 1980. View at Publisher · View at Google Scholar
  6. S. Balcaen and H. Ooghe, “35 years of studies on business failure: an overview of the classic statistical methodologies and their related problems,” The British Accounting Review, vol. 38, no. 1, pp. 63–93, 2006. View at Publisher · View at Google Scholar · View at Scopus
  7. P. P. M. Pompe and A. J. Feelders, “Using machine learning, neural networks, and statistics to predict corporate bankruptcy,” Microcomputers in Civil Engineering, vol. 12, no. 4, pp. 267–276, 1997. View at Google Scholar · View at Scopus
  8. G. Zhang, M. Y. Hu, B. E. Patuwo, and D. C. Indro, “Artificial neural networks in bankruptcy prediction: general framework and cross-validation analysis,” European Journal of Operational Research, vol. 116, no. 1, pp. 16–32, 1999. View at Publisher · View at Google Scholar · View at Scopus
  9. K.-S. Shin, T. S. Lee, and H.-J. Kim, “An application of support vector machines in bankruptcy prediction model,” Expert Systems with Applications, vol. 28, no. 1, pp. 127–135, 2005. View at Publisher · View at Google Scholar · View at Scopus
  10. K. Y. Tam, “Neural network models and the prediction of bank bankruptcy,” Omega, vol. 19, no. 5, pp. 429–445, 1991. View at Publisher · View at Google Scholar · View at Scopus
  11. W. L. Tung, C. Quek, and P. Cheng, “GenSo-EWS: a novel neural-fuzzy based early warning system for predicting bank failures,” Neural Networks, vol. 17, no. 4, pp. 567–587, 2004. View at Publisher · View at Google Scholar · View at Scopus
  12. A. Tsakonas, G. Dounias, M. Doumpos, and C. Zopounidis, “Bankruptcy prediction with neural logic networks by means of grammar-guided genetic programming,” Expert Systems with Applications, vol. 30, no. 3, pp. 449–461, 2006. View at Publisher · View at Google Scholar · View at Scopus
  13. S. Wang, L. Wu, Y. Zhang, and Z. Zhou, “Ant colony algorithm used for bankruptcy prediction,” in Proceedings of the 2nd International Symposium on Information Science and Engineering (ISISE '09), pp. 137–139, December 2009. View at Publisher · View at Google Scholar · View at Scopus
  14. L. Rui, “A particle swarm optimized Fuzzy Neural Network for bankruptcy prediction,” in Proceedings of the International Conference on Future Information Technology and Management Engineering (FITME '10), pp. 557–560, October 2010. View at Publisher · View at Google Scholar · View at Scopus
  15. K.-S. Shin, T. S. Lee, and H.-J. Kim, “An application of support vector machines in bankruptcy prediction model,” Expert Systems with Applications, vol. 28, no. 1, pp. 127–135, 2005. View at Publisher · View at Google Scholar · View at Scopus
  16. J. H. Min and Y.-C. Lee, “Bankruptcy prediction using support vector machine with optimal choice of kernel function parameters,” Expert Systems with Applications, vol. 28, no. 4, pp. 603–614, 2005. View at Publisher · View at Google Scholar · View at Scopus
  17. A. Fan and M. Palaniswami, “Selecting bankruptcy predictors using a support vector machine approach,” in Proceedings of the International Joint Conference on Neural Networks, vol. 6, pp. 354–359, July 2000. View at Scopus
  18. B. Baesens, S. Viaene, T. van Gestel et al., “Bankruptcy prediction with least squares support vector machine classifiers,” in Proceedings of the IEEE International Conference on Computational Intelligence for Financial Engineering, pp. 1–8, September 2000. View at Scopus
  19. J. Jayanthi, J. K. Suresh, and J. Vaishnavi, “Bankruptcy prediction using svm and hybrid svm survey,” International Journal of Computer Applications, vol. 34, no. 7, pp. 39–45, 2011. View at Google Scholar
  20. M. J. Martín-Bautista and M. A. Vila, “A survey of genetic feature selection in mining issues,” in Proceedings of the Congress on Evolutionary Computation (CEC '99), vol. 2, pp. 306–313, 1999.
  21. M. Kudo and J. Sklansky, “A comparative evaluation of medium and largescale feature selectors for pttern classifiers,” in Proceedings of the 1st International Workshop on Statistical Techniques in Pattern Recognition, pp. 91–96, 1997.
  22. I. Guyon and A. Elisseeff, “An introduction to variable and feature selection,” Journal of Machine Learning Research, vol. 3, pp. 1157–1182, 2003. View at Google Scholar
  23. C.-F. Tsai, “Feature selection in bankruptcy prediction,” Knowledge-Based Systems, vol. 22, no. 2, pp. 120–127, 2009. View at Publisher · View at Google Scholar · View at Scopus
  24. F. Duran, C. Cotta, and A. Fernandez, “On the use of sharpe’s index in evolutionary portfolio optimization under Markowitz’s model,” in Proceedings of the Concention on Adaptive and Emergent Behaviour and Complex Systems, 2009.
  25. P. Skolpadungket, K. Dahal, and N. Harnpornchai, “Portfolio optimization using multi-objective genetic algorithms,” in Proceedings of the IEEE Congress on Evolutionary Computation (CEC '07), pp. 516–523, September 2007. View at Publisher · View at Google Scholar · View at Scopus
  26. L. Dioşan, “A multi-objective evolutionary approach to the portfolio optimization problem,” in Proceeding of the International Conference on Computational Intelligence for Modelling, Control ad Automation, vol. 2, pp. 183–187, November 2005. View at Scopus
  27. S. Mullei and P. Beling, “Hybrid evolutionary algorithms for a multiobjective financial problem,” in Proceedings of the IEEE International Conference on Systems, Man. and Cybernetics, vol. 4, pp. 3925–3930, October 1998. View at Scopus
  28. F. Schlottmann and D. Seese, “Hybrid multi-objective evolutionary computation of constrained downside risk-return efficeient sets for credit portfolio,” in Proceedings of the 8th International Conference of the Society for Computational Economics, Computing in Economics and Finances, 2002.
  29. A. Mukerjee, R. Biswas, K. Deb, and A. Mathur, “Multi-objective evolutionary algorithms for the risk-return trade-off in bank-loan management,” International Transactions in Operational Research, vol. 9, no. 5, pp. 583–597, 2002. View at Google Scholar
  30. S. Mardle, S. Pascoe, and M. Tamiz, “An investigation of genetic algorithms forthe optimization of multiobjective fisheries bioeconomic models,” International Transactions in Operational Research, vol. 7, no. 1, pp. 33–49, 2000. View at Google Scholar
  31. A. Gaspar-Cunha, F. Mendes, J. A. Duarte, A. Vieira, B. Ribeiro, and J. A. Neves, “Multi-objective evolutionary algorithms for feature selection: application in bankruptcy prediction,” in Proceedings of the 8th International Conference on Simulated Evolution and Learning, pp. 319–328, 2010.
  32. S. Meyer-Nieberg and H.-G. Beyer, “Self-adaptation in evolutionary algorithms,” History, vol. 75, pp. 47–75, 2007. View at Publisher · View at Google Scholar · View at Scopus
  33. T. Bäck, “Self-adaptation,” in Handbook of Evolutionary Computation, pp. C7.1:1–C7.1:15, Oxford University Press, 1997. View at Google Scholar
  34. S. Zeng, Z. Liu, C. Li, Q. Zhang, and W. Wang, “An evolutionary algorithm and its application in antenna design,” Journal of Bioinformatics and Intelligent Control, vol. 2, no. 1, pp. 129–137, 2012. View at Google Scholar
  35. M. P. Poland, C. D. Nugent, H. Wang, and L. M. Chen, “Genetic algorithm and pure random search for exosensor distribution optimisation,” International Journal of Bio-Inspired Computation, vol. 6, no. 4, pp. 359–372, 2012. View at Google Scholar
  36. A. F. Sheta, P. Rausch, and A. S. Al-Afeef, “A monitoring and control framework for lost foam casting manufacturing processes using genetic programming,” International Journal of Bio-Inspired Computation, vol. 2, no. 4, pp. 111–118, 2012. View at Google Scholar
  37. J. Muñuzuri, P. C. Achedad, M. Rodríguez, and R. Grosso, “Use of a genetic algorithm for building efficient choice designs,” International Journal of Bio-Inspired Computation, vol. 1, no. 4, pp. 27–32, 2012. View at Google Scholar
  38. B. B. Pal, D. Chakraborti, P. Biswas, and A. Mukhopadhyay, “An applicationof genetic algorithm method for solving patrol manpower deployment problems through fuzzy goal programming in traffic management system: a case study,” International Journal of Bio-Inspired Computation, vol. 1, no. 4, pp. 47–60, 2012. View at Google Scholar
  39. L. Carlos Molina, L. Belanche, and À. Nebot, “Feature selection algorithms: a survey and experimental evaluation,” in Proceedings of the 2nd IEEE International Conference on Data Mining (ICDM '02), pp. 306–313, December 2002. View at Scopus
  40. R. Kohavi and G. H. John, “Wrappers for feature subset selection,” Artificial Intelligence, vol. 97, no. 1-2, pp. 273–324, 1997. View at Google Scholar · View at Scopus
  41. F. Provost and T. Fawcet, “Analysis and verification of classifier performance: classification under imprecise class and cost distributions,” in Proceedings of the International Conference on Knowledge Discovery and Data Mining (KDD '97), pp. 43–48, 1997.
  42. T. Fawcet, “An introduction to roc analysis,” Pattern Recognition Letters, vol. 27, pp. 861–874, 2006. View at Google Scholar
  43. M. Hall, E. Frank, G. Holmes, B. Pfahringer, P. Reutemann, and I. H. Witten, “The weka data mining software: an update,” SIGKDD Explorations, vol. 1, no. 1, 2009. View at Google Scholar
  44. X. Wu, “Top 10 algorithms in data mining,” Knowledge and Information Systems, vol. 14, no. 1, pp. 1–37, 2007. View at Google Scholar
  45. J. Shawe-Taylor and N. Cristianini, Support Vector Machines and Other Kernel-Based Learning Methods, Cambridge University Press, 2000.
  46. C. Cortes and V. Vapnik, “Support-vector networks,” Machine Learning, vol. 20, no. 3, pp. 273–297, 1995. View at Publisher · View at Google Scholar · View at Scopus
  47. C. C. Chang and C. J. Lin, “Libsvm a library for support vector machines,” 2000,
  48. D. Kalyanmoy, Multi-Objective Optimization Using Evolutionary Algorithms, Wiley, 2001.
  49. A. Carlos Coello, B. Gary Lamont, and D. A. van Veldhuizen, Evolutionary Algorithms for Solving Multi-Objective Problems (Genetic and Evolutionary Computation), Springer, New York, NY, USA, 2006.
  50. A. Gaspar-Cunha, “RPSGAe-reduced pareto set genetic algorithm: application to polymer extrusion,” in Metaheuristics for Multiobjective Optimisation, pp. 221–249, Springer, 2004. View at Google Scholar
  51. A. Gaspar-Cunha and J. A. Covas, “Robustness in multi-objective optimization using evolutionary algorithms,” Computational Optimization and Applications, vol. 39, no. 1, pp. 75–96, 2008. View at Publisher · View at Google Scholar · View at Scopus
  52. B. Ribeiro, C. Silva, A. Vieira, A. Gaspar-Cunha, and J. C. das Neves, “Financial distress model prediction using SVM+,” in Proceedings of the International Joint Conference on Neural Networks (IJCNN '10), pp. 1–7, July 2010. View at Publisher · View at Google Scholar · View at Scopus
  53. B. Ribeiro, N. Lopes, and C. Silva, “High-performance bankruptcy prediction model using Graphics Processing Units,” in Proceedings of the International Joint Conference on Neural Networks (IJCNN '10), pp. 1–7, July 2010. View at Publisher · View at Google Scholar · View at Scopus
  54. A. Frank and A. Asuncion, “UCI machine learning repository,” 2010,
  55. K. Wang, S. Zhou, A. W. Fu, and J. X. Yu, “Mining changes of classification by correspondence tracing,” in SDM, 2003. View at Google Scholar
  56. J. A. Gonzalez, L. B. Holder, and D. J. Cook, “Graph based concept learning,” in AAAI/IAAI, 2000. View at Google Scholar
  57. J. R. Quinlan, “Simplifying decision trees,” International Journal of Man-Machine Studies, vol. 27, no. 3, pp. 221–234, 1987. View at Google Scholar · View at Scopus
  58. M. A. Hal, Correlation-based feature selection for machine learning [Ph.D. thesis], Department of Computer Science, University of Waikato, Waikato, New Zealand, 1999.
  59. K. Deb and A. Raji Reddy, “Reliable classification of two-class cancer data using evolutionary algorithms,” BioSystems, vol. 72, no. 1-2, pp. 111–129, 2003. View at Publisher · View at Google Scholar · View at Scopus
  60. A. Bagirov, A. Rubinov, N. Soukhoroukova, and J. Yearwood, “Unsupervised and supervised data classification via nonsmooth and global optimization,” TOP, vol. 11, pp. 1–75, 2003. View at Publisher · View at Google Scholar
  61. A. Vieira, J. Duarte, B. Ribeiro, and J. Neves, “Accurate prediction of financial distress of companies with machine learning algorithms,” in Adaptive and Natural Computing Algorithms, M. Kolehmainen, P. Toivanen, and B. Beliczynski, Eds., vol. 5495 of Lecture Notes in Computer Science, pp. 569–576, Springer, Berlin, Germany, 2009. View at Google Scholar
  62. N. Kwak and C.-H. Choi, “Input feature selection for classification problems,” IEEE Transactions on Neural Networks, vol. 13, no. 1, pp. 143–159, 2002. View at Publisher · View at Google Scholar · View at Scopus