A Corporate Credit Rating Model Using Support Vector Domain Combined with Fuzzy Clustering Algorithm
Corporate credit-rating prediction using statistical and artificial intelligence techniques has received considerable attentions in the literature. Different from the thoughts of various techniques for adopting support vector machines as binary classifiers originally, a new method, based on support vector domain combined with fuzzy clustering algorithm for multiclassification, is proposed in the paper to accomplish corporate credit rating. By data preprocessing using fuzzy clustering algorithm, only the boundary data points are selected as training samples to accomplish support vector domain specification to reduce computational cost and also achieve better performance. To validate the proposed methodology, real-world cases are used for experiments, with results compared with conventional multiclassification support vector machine approaches and other artificial intelligence techniques. The results show that the proposed model improves the performance of corporate credit-rating with less computational consumption.
Techniques of credit ratings have been applied by bond investors, debt issuers, and governmental officials as one of the most efficient measures of risk management. However, company credit ratings are too costly to obtain, because agencies including Standard and Poor’s (S&P), and Moody’s are required to invest lots of time and human resources to accomplish critical analysis based on various aspects ranging from strategic competitiveness to operational level in detail [1–3]. Moreover, from a technical perspective, credit rating constitutes a typical multiclassification problem, because the agencies generally have much more than two categories of ratings. For example, ratings from S&P range from AAA for the highest-quality bonds to D for the lowest-quality ones.
The final objective of credit rating prediction is to develop the models, by which knowledge of credit risk evaluation can be extracted from experiences of experts and to be applied in much broader scope. Besides prediction, the studies can also help users capture fundamental characteristics of different financial markets by analyzing the information applied by experts.
Although rating agencies take emphasis on experts’ subjective judgment in obtaining ratings, many promising results on credit rating prediction based on different statistical and Artificial Intelligence (AI) methods have been proposed, with a grand assumption that financial variables extracted from general statements, such as financial ratios, contain lots of information about company’s credit risk, embedded in their valuable experiences [4, 5].
Among the technologies based on AI applied in credit rating prediction, the Artificial Neural Networks (ANNs) have been applied in the domain of finance because of the ability to learn from training samples. Moreover, in terms of defects of ANN such as overfitting, Support Vector Machine (SVM) has been regarded as one of the popular alternative solutions to the problems, because of its much better performance than traditional approaches such as ANN [6–11]. That is, an SVM’s solution can be globally optimal because the models seek to minimize the structural risk . Conversely, the solutions found by ANN tend to fall into local optimum because of seeking to minimize the empirical risk.
However, SVM, which was originally developed for binary classification, is not naturally modified for multiclassification of many problems including credit ratings. Thus, researchers have tried to extend original SVM to multiclassification problems , with some techniques of multiclassification SVM (MSVM) proposed, which include approaches that construct and combine several binary classifiers as well as the ones that directly consider all the data in a single optimization formulation.
In terms of multiclassification in the domain of credit rating containing lots of data, current approaches applied in MSVM still have some drawbacks in integration of multiple binary classifiers as follows.(1)Some unclassifiable regions may exist if a data point belongs to more than one class or to none.(2)Training binary classifiers based on two-class SVM multiple times for the same data set often result in a highly intensive time complexity for large-scale problems including credit ratings prediction to improve computational consumption.
To overcome the drawbacks associated with current MSVM in credit rating prediction, a novel model based on support vector domain combined with kernel-based fuzzy clustering is proposed in the paper to accomplish multiclassification involved in credit ratings prediction.
2. Literature Review
2.1. Credit Rating Using Data Mining Techniques
Major researches applying data mining techniques for bond rating prediction can be found in the literature.
Early investigations of credit rating techniques mainly focused on the applicability of statistical techniques including multiple discriminant analysis (MDA) [14, 15] and logistic regression analysis (LRA) , and so forth, while typical techniques of AI including ANN [17, 18] and case-based reasoning (CBR) , and so forth are applied in the second phase of research.
The important researches applying AI techniques in bond-rating prediction are listed in Table 1. In summary, the most prior ones accomplish prediction using ANN with comparison to other statistical methods, with general conclusions that neural networks outperformed conventional statistical methods in the domain of bond rating prediction.
On the other hand, to overcome the limitations such as overfitting of ANN, techniques based on MSVM are applied in credit rating in recent years. Among the models based on MSVM in credit rating, method of Grammar and Singer was early proposed by Huang et al., with experiments based on different parameters so as to find the optimal model . Moreover, methodologies based on One-Against-All, One-Against-One, and DAGSVM are also proposed to accomplish S&P’s bond ratings prediction, with kernel function of Gaussian RBF applied and the optimal parameters derived form a grid-search strategy . Another automatic-classification model for credit rating prediction based on One-Against-One approach was also applied . And Lee applied MSVM in corporate credit rating prediction , with experiments showing that model based on MSVM outperformed other AI techniques such as ANN, MDA, and CBR.
2.2. Multiclassification by Support Vector Domain Description
Support Vector Domain Description (SVDD), proposed by Tax and Duin in 1999  and extended in 2004 , is a method for classification with the aim to accomplish accurate estimation of a set of data points originally. The methods based on SVDD differ from two or multiclass classification in that a single object type is interested rather than to be separated from other classes. The SVDD is a nonparametric method in the sense that it does not assume any particular form of distribution of the data points. The support of unknown distribution of data points is modeled by a boundary function. And the boundary is “soft” in the sense that atypical points are allowed outside it.
The boundary function of SVDD is modeled by a hypersphere rather than a hyperplane applied in standard SVM, which can be made with less constrains by mapping the data points to a high-dimensional space using methodology known as kernel trick, where the classification is performed.
SVDD has been applied in a wide range as a basis for new methodologies in statistical and machine learning, whose application in anomaly detection showed that the model based on it can improve accuracy and reduce computational complexity . Moreover, ideas of improving the original SVDD through weighting each data point by an estimate of its corresponding density were also proposed  and applied in area of breast cancer, leukemia, and hepatitis, and so forth. Other applications including pump failure detection , face recognition , speaker recognition , and image retrieval  are argued by researchers.
The capability of SVDD in modeling makes it one of the alternative to large-margin classifiers such as SVM. And some novel methods applied in multiclass classification were proposed based on SVDD  combined with other algorithms such as fuzzy theories [41, 42] and Bayesian decision .
3. The Proposed Methodology
In terms of SVDD, which is a boundary-based method for data description, it needs more boundary samples to construct a closely fit boundary. Unfortunately, more boundary ones usually imply that more target objects have to be rejected with the overfitting problem arising and computational consumption increased. To accomplish multiclassification in corporate credit rating, a method using Fuzzy SVDD combined with fuzzy clustering algorithm is proposed in the paper. By mapping data points to a high-dimensional space by Kernel Trick, the hypersphere applied to every category is specified by training samples selected as boundary ones, which are more likely to be candidates of support vectors. After preprocessing using fuzzy clustering algorithm, rather than by original ones directly in standard SVDD [32, 33], one can improve accuracy and reduce computational consumption. Thus, testing samples are classified by the classification rules based on hyperspheres specified for every class. And the thoughts and framework of the proposed methodology can be illustrated in Figures 1 and 2, respectively.
3.1. Fuzzy SVDD
3.1.1. Introduction to Hypersphere Specification Algorithm
The hypersphere, by which SVDD models data points, is specified by its center and radius . Let denote the data matrix with data points and variables, which implies that is p-dimensional while is scalar. The geometry of one solution to SVDD in two dimensions is illustrated in Figure 3, where represents the perpendicular distance from the boundary to an exterior points . In terms of interior points, and the ones positioned on the boundary, is to be assigned as 0. Hence, can be calculated using the following equation:
In the following, another closely related measure can be obtianed in (3.2) in terms of exterior points
To obtain an exact and compact representation of the data points, the minimization of both the hypersphere radius and to any exterior point is required. Moreover, inspired by fuzzy set theory, matrix can be extended to with coefficients representing fuzzy membership associated with introduced. So, the data domain description can be formulated as (3.3), where nonnegative slack variables are a measure of error in SVDD, and the term is the one with different weights based on fuzzy set theory
To solve the problem, the Lagrange Function is introduced, where are Lagrange Multipliers shown as follows:
Setting (3.4) to 0, the partial derivates of leads to the following equations: That is,
The Karush-Kuhn-Tucker complementarities conditions result in the following equations:
Therefore, the dual form of the objective function can be obtained as follows:
And the problem can be formulated as follows:
The center of the hypersphere is a linear combination of data points with weighting factors obtained by optimizing (3.9). And the coefficients , which are nonzero, are thus selected as support vectors, only by which the hypersphere is specified and described. Hence, to judge whether a data point is within a hypersphere, the distance to the center should be calculated with (3.10) in order to judge whether it is smaller than the radius . And the decision function shown as (3.12) can be concluded from
3.1.2. Introduction to Fuzzy SVDD Based on Kernel Trick
Similarly to the methodology based on kernel function proposed by Vapnik , the Fuzzy SVDD can also be generalized to high-dimensional space by replacing its inner products by kernel functions .
For example, Kernel function of RBF can be introduced to SVDD algorithm, just as shown as follows:
3.2. Kernel-Based Fuzzy Clustering Algorithm
3.2.1. Introduction to Fuzzy Attribute C-Means Clustering
Suppose denote any finite sample set, where , and each sample is defined as . The category of attribute space is , where is the cluster number. For , let denote the attribute measure of , with .
Let denote the kth prototype of cluster , where .
Let denote the attribute measure of the th sample belonging to the kth cluster. That is, , , . The task of fuzzy clustering is to calculate the attribute measure , and determine the cluster which belongs to according to the maximum cluster index.
Fuzzy C-means (FCM) is an inner-product-induced distance based on the least-squared error criterion. A brief review of FCM can be found in Appendix based on coefficients definitions mentioned above.
Attribute Means Clustering (AMC) is an iterative algorithm by introducing the stable function . Suppose is a positive differential function in . Let , if , called as weight function, is a positive nonincreasing function, is called as stable function. And can be adopted as follows:
Hence, the relationship of objective function and its weight function is described by sable function, which was introduced to propose AMC.
According to current researches, some alternative functions including squared stable function, Cauchy stable function, and Exponential stable function are recommended.
Based on previous researches, AMC and FCM are extended to FAMC, which is also an iterative algorithm to minimize the following objective function shown as (3.16), where , which is a coefficient of FCM introduced in Appendix
3.2.2. Introduction to Kernel-Based Fuzzy Clustering
To gain a high-dimensional discriminant, FAMC can be extended to Kernel-based Fuzzy Attribute C-means Clustering (KFAMC). That is, the training samples can be first mapped into high-dimensional space by the mapping Φ using kernel function methods addressed in Section 3.1.2.
Since when Kernel function of RBF is introduced, (3.19) can be given as follows
And parameters in KFAMC can be estimated by where .
3.2.3. Algorithms of Kernel-Based Fuzzy Attribute C-Means Clustering
Based on theorem proved in , the updating procedure of KFAMC can be summarized in the following iterative scheme.
Step 1. Set , , and , and initialize , .
Step 2. For , calculate fuzzy cluster centers , . and .
Step 3. If or , stop, else go to Step 4.
Step 4. For step , update , , and , turn to Step 3,
where denotes iterate step, represents the maximum iteration times, and denotes the weighting matrix, respectively, which can be seen in  in detail.
3.3. The Proposed Algorithm
3.3.1. Classifier Establishment
In terms of SVDD, only support vectors are necessary to specify hyperspheres. But in the original algorithms [32, 33, 41], all the training samples are analyzed and thus computational cost is high consumption. Hence, if the data points, which are more likely to be candidates of support vectors, can be selected as training samples, the hypersphere will be specified with much less computational consumption.
Just as illustrated in Figure 4, only the data points, such as M, N positioned in fuzzy areas, which are more likely to be candidates of support vectors, are necessary to be classified with SVDD, while the ones in deterministic areas can be regarded as data points belonging to certain class.
So, the new methodology applied in SVDD is proposed as follows.(1)Preprocess data points using FAMC to reduce amount of training samples. That is, if fuzzy membership of a data point to a class is great enough, the data point can be ranked to the class directly. Just as shown in Figure 5, the data points positioned in deterministic area (shadow area A) are to be regarded as samples belonging to the class, while the other ones are selected as training samples.(2) Accomplish SVDD specification with training samples positioned in fuzzy areas, which has been selected using KFAMC. That is, among the whole data points, only the ones in fuzzy area, rather than all the data points, are treated as candidates of support vectors. And the classifier applied in multiclassification can be developed based on Fuzzy SVDD by specifying hypersphere according to every class. Hence, the main thoughts of Fuzzy SVDD establishment combined with KFAMC can be illustrated in Figure 6.
(a) Training data points obtained by preprocessing
(b) Hypersphere specification after data points preprocessing
The process of methods proposed in the paper can be depicted as follows.
In high-dimensional space, the training samples are selected according to their fuzzy memberships to clustering centers. Based on preprocessing with KFAMC, a set of training samples is given, which is represented by , where , and denote the number of training data, input pattern, and membership to class , respectively.
Hence, the process of Fuzzy SVDD specification can be summarized as follows.
Step 1. Set a threshold , and apply KFAMC to calculate the membership of each , to each class. If , is to be set as 1 and , is to be set as 0.
Step 2. Survey the membership of each . If , is to be ranked to class directly and removed from the training set. And an updated training set can be obtained.
Step 3. With hypersphere specified for each class using the updated training set obtained in Step 2, classifier for credit rating can be established using the algorithm of Fuzzy SVDD, just as illustrated in Figure 6.
3.3.2. Classification Rules for Testing Data Points
To accomplish multiclassification for testing data points using hyperspheres specified in Section 3.3.1, the following two factors should be taken into consideration, just as illustrated in Figure 7:(1) distances from the data point to centers of the hyperspheres;(2) density of the data points belonging to the class implied with values of radius of each hypersphere.
Just as shown in Figure 7, , denote the distances from data point to center of class A and class B, respectively. Even if , data point is expected more likely to belong to class A rather than class B because of difference in distributions of data points. That is, data points circled by hypersphere of class A are sparser than the ones circled by hypersphere of class B since is greater than .
So, classification rules can be concluded as follows.
Let denote the numbers of hyperspheres containing the data point.
Case I (). Data point belongs to the class represented by the hypersphere.
Case II ( or ). Calculate the index of membership of the data point to each hypersphere using (3.24), where denotes the radius of hypersphere denotes the distance from data point to the center of hypersphere
And the testing data points can be classified according to the following rules represented with
4.1. Data Sets
For the purpose of this study, two bond-rating data sets from Korea and China market, which have been used in [46, 47], are applied, in order to validate the proposed methodology. The data are divided into the following four classes: A1, A2, A3, and A4.
4.2. Variables Selection
Methods including independent-samples -test and F-value are applied in variable selection.
In terms of Korea data set, 14 variables, which are listed in Table 2, are selected from original ones, which were known to affect bond rating. For better comparison, similar methods were also used in China data set, with 12 variables among them being selected.
4.3. Experiment Results and Discussions
Based on the two data sets, some models based on AI are introduced for experiments. To evaluate the prediction performance, 10-fold cross validation, which has shown good performance in model selection , is followed. In the research, all features, which are represented with variables listed in Table 2, of data points range from 0 to 1 after Min-max transformation. To validate the methodology oriented multiclassification problem in credit rating, ten percent of the data points for each class are selected as testing samples. And the results of experiments on proposed method, with 0.9 being chosen as the value of threshold intuitively, are shown in Table 3.
To compare with other methods, the proposed model is compared with some other MSVM techniques, namely, ANN, One-Against-All, One-Against-One, DAGSVM, Grammer & Singer, OMSVM , and standard SVDD. The results concluded in the paper are all shown as average values obtained following 10-fold cross validation based on platform of Matlab 7.0.
To compare the performance of each algorithm, hit-ratio, which is defined according to the samples classified correctly, is applied. And the experiment results are listed in Table 4.
As shown in Table 4, the proposed method based on thoughts of hypersphere achieves better performance than conventional SVM models based on thoughts of hyperplane. Moreover, as one of modified models, some results obtained imply that the proposed method has better generalization ability and less computational complexity, which can be partially measured with training time labeled with “Time,” than standard SVDD.
Furthermore, as one of modified models based on standard SVDD, the proposed method accomplishes data preprocessing using KFAMC. Since the fuzzy area is determined by threshold , greater value of will lead to bigger fuzzy area. Especially, when , the algorithm proposed will be transformed to standard SVDD because almost all data points are positioned in fuzzy area. Hence, a model with too large threshold may be little different from standard SVDD, while a too small value will have poor ability of sphere-based classifier establishment due to lack of essential training samples. Thus, issues on choosing the appropriate threshold are discussed by empirical trials in the paper.
In the following experiment, the proposed method with various threshold values is tested based on different data sets, just as shown in Figure 8.
(a) Experiments based on korea data set
(b) Experiments based on china data set
The results illustrated in Figure 8 showed that the proposed method achieved best performance with threshold of 0.9 based on Korea data set. But in terms of China market, it achieved best performance with the threshold of 0.8 rather than a larger one due to effects of more outliers existing in data set.
Moreover, training time of proposed method can be also compared with standard SVDD, just as illustrated in Figure 9.
(a) Experiments based on korea data set
(b) Experiments based on china data set
Just as shown in Figure 9, results of experiments based on different data sets are similar. That is, with decline of threshold, more samples were eliminated from training set through preprocessing based on KFAMC to reduce training time. Hence, smaller values of threshold will lead to less computational consumption partly indicated as training time, while classification accuracy may be decreased due to lack of necessary training samples. Overall, threshold selection, which involves complex tradeoffs between computational consumption and classification accuracy, is essential to the proposed method.
5. Conclusions and Directions for Future Research
In the study, a novel algorithm based on Fuzzy SVDD combined with Fuzzy Clustering for credit rating is proposed. The underlying assumption of the proposed method is that sufficient boundary points could support a close boundary around the target data but too many ones might cause overfitting and poor generalization ability. In contrast to prior researches, which just applied conventional MSVM algorithms in credit ratings, the algorithm based on sphere-based classifier is introduced with samples preprocessed using fuzzy clustering algorithm.
As a result, through appropriate threshold setting, generalization performance measured by hit-ratio of the proposed method is better than that of standard SVDD, which outperformed many kinds of conventional MSVM algorithms argued in prior literatures. Moreover, as a modified sphere-based classifier, proposed method has much less computational consumption than standard SVDD.
One of the future directions is to accomplish survey studies comparing different bond-rating processes, with deeper market structure analysis also achieved. Moreover, as one of the MSVM algorithms, the proposed method can be applied in other areas besides credit ratings. And some more experiments on data sets such as UCI repository  are to be accomplished in the future.
Brief Review of FCM
Bezdek-type FCM is an inner-product-induced distance-based least-squared error criterion nonlinear optimization algorithm with constrains, where is the measure of the th sample belonging to the kth cluster and is the weighting exponent. The distance between and the prototype of kth cluster is as follows:
The above formula is also called as Mahalanobis distance, where is a positive matrix. When is a unit matrix, is Euclidean distance. We denote it as and adopt Euclidean distance in the rest of the paper. So, the parameters of FCM are estimated by updating according to the formulas:
The paper was sponsored by 985-3 project of Xi’an Jiaotong University.
W. F. Treacy and M. Carey, “Credit risk rating systems at large US banks,” Journal of Banking & Finance, vol. 24, no. 1-2, pp. 167–201, 2000.View at: Google Scholar
H. Ahn, K. J. Kim, and I. Han, “Purchase prediction model using the support vector machine,” Journal of Intelligence and Information Systems, vol. 11, pp. 69–81, 2005.View at: Google Scholar
V. Vapnik, The Nature of Statistical Learning Theory, Springer, New York, NY, USA, 1995.
G. E. Pinches and K. A. A. Mingo, “Multivariate analysis of industrial bond ratings,” The Journal of Finance, vol. 28, no. 1, pp. 1–18, 1973.View at: Google Scholar
A. Belkaoui, “Industrial bond ratings: a new look,” Financial Management, vol. 9, no. 3, pp. 44–51, 1980.View at: Google Scholar
L. H. Ederington, “Classification models and bond ratings,” The Financial Review, vol. 20, no. 4, pp. 237–262, 1985.View at: Google Scholar
J. W. Kim, H. R. Weistroffer, and R. T. Redmond, “Expert systems for bond rating: a comparative analysis of statistical, rule-based and neural network systems,” Expert Systems, vol. 10, no. 3, pp. 167–172, 1993.View at: Google Scholar
R. Chaveesuk, C. Srivaree-Ratana, and A. E. Smith, “Alternative neural network approaches to corporate bond rating,” Journal of Engineering Valuation and Cost Analysis, vol. 2, no. 2, pp. 117–131, 1999.View at: Google Scholar
H. J. Kim and K. S. Shin, “A hybrid approach using case-based reasoning and fuzzy logic for corporate bond rating,” Journal of Intelligence and Information system, vol. 16, pp. 67–84, 2010.View at: Google Scholar
S. Dutta and S. Shekhar, “Bond rating: a non-conservative application of neural networks,” in Proceedings of the IEEE International Conference on Neural Networks, vol. 2, pp. 443–450, San Diego, Calif, USA, July 1988.View at: Google Scholar
J. C. Singleton and A. J. Surkan, “Neural networks for bond rating improved by multiple hidden layers,” in Proceedings of the the IEEE International Conference on Neural Networks, vol. 2, pp. 163–168, June 1990.View at: Google Scholar
S. Garavaglia, “An application of a counter-propagation neural networks: simulating the standard & poor's corporate bond rating system,” in Proceedings of the 1st International Conference on Artificial Intelligence on Wall Street, pp. 278–287, 1991.View at: Google Scholar
J. Moody and J. Utans, “Architecture selection strategies for neural networks application to corporate bond rating,” in Neural Networks in the Capital Markets, pp. 277–300, John Wiley & Sons, 1995.View at: Google Scholar
J. J. Maher and T. K. Sen, “Predicting bond ratings using neural networks: a comparison with logistic regression,” Intelligent Systems in Accounting, Finance and Management, vol. 6, no. 1, pp. 59–72, 1997.View at: Google Scholar
Y. S. Kwon, I. G. Han, and K. C. Lee, “Ordinal pairwise partitioning (OPP) approach to neural networks training in bond rating,” Intelligent Systems in Accounting, Finance and Management, vol. 6, no. 1, pp. 23–40, 1997.View at: Google Scholar
R. Chaveesuk, C. S. Ratana, and A. E. Smith, “Alternative neural network approaches to corporate bond rating,” Journal of Engineering Valuation and Cost Analysis, vol. 2, no. 2, pp. 117–131, 1999.View at: Google Scholar
L. Cao, L. K. Guan, and Z. Jingqing, “Bond rating using support vector machine,” Intelligent Data Analysis, vol. 10, no. 3, pp. 285–296, 2006.View at: Google Scholar
D. Tax, A. Ypma, and R. Duin, “Pump failure detection using support vector data descriptions,” in Proceedings of the 3rd International Symposium on Advances in Intelligent Data (IDA '99), pp. 415–425, 1999.View at: Google Scholar
X. Dong, W. Zhaohui, and Z. Wanfeng, “Support vector domain description for speaker recognition,” in Proceedings of the IEEE Signal Processing Society Workshop, Neural Networks for Signal Processing XI, pp. 481–488, Falmouth, Mass, USA, 2001.View at: Google Scholar
L. L. Wei, W. J. Long, and W. X. Zhang, “Fuzzy data domain description using support vector machines,” in Proceedings of the 2nd International Conference on Machine Learning and Cybernetics, pp. 3082–3085, November 2003.View at: Google Scholar
L. A. Zadeh, “Fuzzy sets,” Information and Control, vol. 8, no. 3, pp. 338–353, 1965.View at: Google Scholar
J. Liu and M. Xu, “Bezdek type fuzzy attribute C-means clustering algorithm,” Journal of Beijing University of Aeronautics and Astronautics, vol. 33, no. 9, pp. 1121–1126, 2007.View at: Google Scholar
Q. S. Cheng, “Attribute means clustering,” Systems Engineering—Theory & Practice, vol. 18, no. 9, pp. 124–126, 1998.View at: Google Scholar
S. Pang, Credit Rating and Stock Market Prediction Model, Science Press, Beijing, China, 2005.
S. M. Weiss and C. A. Kulikowski, Computer Systems That Learn: Classification and Prediction Methods from Statistics, Neural Networks, Machine Learning and Expert Systems, Machine Learning, Morgan Kaufmann, San Mateo, Calif , USA, 1991.