<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">AJOR</journal-id><journal-title-group><journal-title>American Journal of Operations Research</journal-title></journal-title-group><issn pub-type="epub">2160-8830</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/ajor.2017.71003</article-id><article-id pub-id-type="publisher-id">AJOR-73375</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Physics&amp;Mathematics</subject></subj-group></article-categories><title-group><article-title>
 
 
  The Development of an Alternative Method for the Sovereign Credit Rating System Based on Adaptive Neuro-Fuzzy Inference System
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Hakan</surname><given-names>Pabuçcu</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Tuba</surname><given-names>Yakıcı Ayan</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib></contrib-group><aff id="aff1"><addr-line>Department of Business Administration, Bayburt University, Bayburt, Turkey</addr-line></aff><aff id="aff2"><addr-line>Department of Econometrics, Karadeniz Technical University, Trabzon, Turkey</addr-line></aff><author-notes><corresp id="cor1">* E-mail:<email>hpabuccu@bayburt.edu.tr(HP)</email>;</corresp></author-notes><pub-date pub-type="epub"><day>27</day><month>12</month><year>2016</year></pub-date><volume>07</volume><issue>01</issue><fpage>41</fpage><lpage>55</lpage><history><date date-type="received"><day>4,</day>	<month>November</month>	<year>2016</year></date><date date-type="rev-recd"><day>8,</day>	<month>January</month>	<year>2017</year>	</date><date date-type="accepted"><day>11,</day>	<month>January</month>	<year>2017</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  The main purpose of this article is to determine the factors affecting credit rating and to develop the credit rating system based on statistical methods, fuzzy logic and artificial neural network. Variables used in this study were determined by the literature review and then the number of them was reduced by using stepwise regression analysis. Resulting variables were used as independent variables in the logistic model and as input variables for ANN and ANFIS model. After evaluating the models and comparing with each other, the ANFIS model was chosen as the best model to forecast credit rating. Rating determination was made for the countries that haven’t had a credit rating. Consequently, the ANFIS model made consistent, reliable and successful rating forecasts for the countries.
 
</p></abstract><kwd-group><kwd>Credit Rating</kwd><kwd> Logistic Regression (LR)</kwd><kwd> Neural Networks (ANN)</kwd><kwd> Adaptive Neuro-Fuzzy Inference System (ANFIS)</kwd><kwd> Comparative Studies</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>Problems which occur during a global recession or decline periods affect all countries first and foremost developing countries. Countries lose their debt discharging competence and economic woes have seen almost each point of economic life. Independent credit rating agencies evaluated for government bonds in the point of predicting these kinds of problems determine credit ratings for countries and evaluate possibilities debtors failed to pay.</p><p>For credit rating studies some statistical methods such as regression analysis [<xref ref-type="bibr" rid="scirp.73375-ref1">1</xref>] [<xref ref-type="bibr" rid="scirp.73375-ref2">2</xref>] , logit and probit regression analysis [<xref ref-type="bibr" rid="scirp.73375-ref3">3</xref>] and discriminant analysis [<xref ref-type="bibr" rid="scirp.73375-ref4">4</xref>] are often used. But nowadays, determining non-linear relationships [<xref ref-type="bibr" rid="scirp.73375-ref5">5</xref>] instead of determining the existence of only linear relationships between variables gained importance. For that reason, advanced techniques such as neural networks, support vector machine [<xref ref-type="bibr" rid="scirp.73375-ref6">6</xref>] , have been applied. When several studies in literature were examined, credit rating was handled on the basis of enterprises and classification problem not on the basis of countries. Besides, any credit rating estimation was made for the countries which haven’t got credit rating except [<xref ref-type="bibr" rid="scirp.73375-ref2">2</xref>] . In this study, rating application was carried out on the basis of countries in order to fill the gap in the literature. While performing this, assessments of three large credit rating agencies (Moody’s, S&amp;P and Fitch) were considered. Also statistic, math and econometric models were used together with both comparative and supporting each other for credit rating application. Finally, credit prediction of 21 countries of which credit rating was not determined until today was carried out. Consequently, it is thought that this study was original and important either in respect to the methods used or in respect to determining the supremacy and deficiency of a method.</p></sec><sec id="s2"><title>2. Sovereign Credit Rating and Related Works</title><p>There is a growing interest to sovereign credit ratings in recent years. The risk assessment performed by the rating agencies that represent the obligation of governments. A rating is a prospective forecast of the default risk. Sovereign ratings are not “country ratings”. There is an important differentiation between them. Sovereign rating is the credit risk of the national governments not the specific default risk of other issuers [<xref ref-type="bibr" rid="scirp.73375-ref7">7</xref>] . Governments usually look for credit ratings to ease their own reach to international capital markets. Sovereign credit ratings are determined by using some macroeconomic-qualitative factors. Although sovereign credit ratings are assigned by the credit rating agencies, lots of questions are in governments minds related with the ratings’ rationale and consistency. The real question is “how clear the factors affect sovereign credit ratings are” [<xref ref-type="bibr" rid="scirp.73375-ref1">1</xref>] . In this study, it is tried to determine the factors affecting the sovereign credit ratings and an alternative model to assign ratings for the countries is developed. Related works with the sovereign-credit ratings are presented as follows.</p><p>Basic variables such as income per capita, gross domestic product (GDP), growth rate, inflation, fiscal balance, external balance, external debt, default back- ground and development sorting which are used for credit valuableness in the study carried out by [<xref ref-type="bibr" rid="scirp.73375-ref1">1</xref>] which was shown as a reference point for almost all studies regarding the credit rating. Direction and volume of Country ratings and relationship between the factors determining the ratings were tried to determine by regression analysis. [<xref ref-type="bibr" rid="scirp.73375-ref8">8</xref>] used the data of S&amp;P and Moody’s in order to analyze the numeric data of credit ratings by dividing into two parts as countries have higher ratings and lower ratings.</p><p>As the most important result of the study, it has been expressed that ratings can’t be explained only with economic and financial indicators. In [<xref ref-type="bibr" rid="scirp.73375-ref3">3</xref>] classification methodology with neural networks and aligned probit analysis comparing was carried out for credit rating. Credit ratings were used as dependent variable, foreign debt/export, financial balance, external balance, inflation rate, gross national product (GNP) per capita, growth rate, development situation of the country was used as an independent variable. [<xref ref-type="bibr" rid="scirp.73375-ref9">9</xref>] used variables of political stability, government efficiency, superiority of law and illegality in order to determine how official bodies affect the political organization’s credit rating. [<xref ref-type="bibr" rid="scirp.73375-ref10">10</xref>] included corruption perception index, the default history, commercial gap, position of a country in the world, democracy index, source of commercial laws, population and petrol production situation along with macro-economic variables into their models and determined the existence of a negative relationship between corruption perception. In the study, where political violence, diplomatic pressure, illegality, military domination, religious trends, effectiveness of laws, structure of democracy and foreign direct investment data were used [<xref ref-type="bibr" rid="scirp.73375-ref11">11</xref>] not only the macro-economic variables but also the qualitative variables were provided. In this study, methodologies of credit rating agencies were handled and even GDP per capita variable could affect credit rating 80%. [<xref ref-type="bibr" rid="scirp.73375-ref12">12</xref>] discussed credit rating agencies and their effects on developing countries. Handled mentioned agencies within the frame of Basel 1 - 2 criteria, examined their methodologies and tried to explain qualitative and quantitative methods in details. [<xref ref-type="bibr" rid="scirp.73375-ref2">2</xref>] determined ratings by using regression analysis for the countries which were not included in to Credit rating by Moody’s Fitch and S&amp;P. Credit ratings and GNP per capita, GDP per capita, reserve ratios, GDP per capita volatility, inflation and superiority of law were used as explanatory variables. Most of the predictions made for the countries which have no rating were predicted as “B” and over. [<xref ref-type="bibr" rid="scirp.73375-ref13">13</xref>] examined the relationships between economic freedom, credit ratings and situation of default of country. Probit analysis and Tobit analysis were used for this study. [<xref ref-type="bibr" rid="scirp.73375-ref14">14</xref>] used the method to compare economic growth which was considered to affect credibility, human development and political stability. As a result of the multiple regression analysis, GDP per capita, internal debt, current account balance and human development which was adapted to inequity, affected debt discharging liability negatively and unemployment and political stability affected debt discharging liability positively. [<xref ref-type="bibr" rid="scirp.73375-ref15">15</xref>] searched the reasons of change of credit ratings in developing markets and how changes of credit ratings in countries affected the bank credibility. Credit rating changes were used as dependent variable and economic freedom index, corruption perception index, property rights, income per capita, inflation, current account balance, financial balance and external debt were used as independent variable.</p></sec><sec id="s3"><title>3. Description of Selected Models</title><sec id="s3_1"><title>3.1. Variable Selection Methodology</title><p>According to [<xref ref-type="bibr" rid="scirp.73375-ref16">16</xref>] , the most significant point in credit rating studies existing variables which symbolize the problems as possible as fewer and has strong representational ability. Due to the variable selection, the dimensions of variable space are decreased to provide effective working of algorithm [<xref ref-type="bibr" rid="scirp.73375-ref17">17</xref>] . When the studies carried out for variable selection were assessed, it is observed that many methods were used. E.g. [<xref ref-type="bibr" rid="scirp.73375-ref18">18</xref>] [<xref ref-type="bibr" rid="scirp.73375-ref19">19</xref>] used Kernel principal components analysis, [<xref ref-type="bibr" rid="scirp.73375-ref20">20</xref>] , [<xref ref-type="bibr" rid="scirp.73375-ref4">4</xref>] used linear discriminant analysis, [<xref ref-type="bibr" rid="scirp.73375-ref17">17</xref>] used genetic algorithm, [<xref ref-type="bibr" rid="scirp.73375-ref21">21</xref>] used one way variance analysis (ANOVA), [<xref ref-type="bibr" rid="scirp.73375-ref22">22</xref>] used one way ANOVA, factor analysis and stepwise regression methods. In this study four different methods were implemented separately for variable selection and best considered stepwise regression analysis were preferred.</p></sec><sec id="s3_2"><title>3.2. Logistic Regression Analysis</title><p>Simple and multiple linear regression models could be predicted by ordinary least squares (OLS). But due to the dependent variable is qualitative (discontinuous, categorical) OLS predicts are not reliable. For that reason, alternative models were developed such as logit and probit [<xref ref-type="bibr" rid="scirp.73375-ref23">23</xref>] . Logistic regression is a special condition of linear regression. But while a dependent variable could take any numeric value in logistic regression, this value should be dual or categorical [<xref ref-type="bibr" rid="scirp.73375-ref24">24</xref>] . In [<xref ref-type="bibr" rid="scirp.73375-ref16">16</xref>] and [<xref ref-type="bibr" rid="scirp.73375-ref25">25</xref>] , logistic regression analysis is mentioned which successfully implemented statistical technique on credit rating studies and many fields.</p></sec><sec id="s3_3"><title>3.3. Artificial Neural Network</title><p>Artificial neural network (ANN) is an advanced mathematical technique which uses intelligent learning paradigms and having several implementation fields such as social, science and engineering fields. The architectural structure of model consists of three layers such as input, hidden and output [<xref ref-type="bibr" rid="scirp.73375-ref26">26</xref>] . In [<xref ref-type="bibr" rid="scirp.73375-ref27">27</xref>] each node in input layer transfers the value belonging to independent variable to the intermediate node and data coming to intermediate layer are combined in determined rules and transformed then mapped to target value in output layer. There is only one node in artificial neural network output layer which has been founded for credit rating. Artificial neural networks not to require to provide independent variables distributive characteristics or assumptions and they could model all non-linear relationships between input-output variables [<xref ref-type="bibr" rid="scirp.73375-ref28">28</xref>] .</p></sec><sec id="s3_4"><title>3.4. Adaptive Neuro Fuzzy Inference System</title><p>Ordinary mathematical method and tools are insufficient to model the systems which haven’t been defined well or uncertain systems. In contrary, fuzzy inference systems (FIS) can obtain better results by including human knowledge and reasoning process by fuzzy “if-then” rules. Fuzzy modeling was developed firstly by [<xref ref-type="bibr" rid="scirp.73375-ref29">29</xref>] and implemented several supervision, predict mechanisms. Adaptive Neuro Fuzzy Inference System (ANFIS) structure is the condition where FIS adapted to adaptive networks. Due to ANFIS hybrid learning algorithm fuzzy “if, then” rules and human knowledge reveal input output structure. ANFIS model which was firstly developed by [<xref ref-type="bibr" rid="scirp.73375-ref30">30</xref>] was used then for modeling the non-linear problems, control systems and solution of chaotic time series problems. Double input (x and y) and double fuzzy “if, then” rule could be mentioned by Takagi Sugeno type ANFIS Equations (1) &amp; (2).</p><p>R 1 : If x is A 1 and y is B 1 then f 1 = p 1 x + q 1 y + r 1 (1)</p><p>R 2 : If x is A 2 and y is B 2 then f 2 = p 2 x + q 2 y + r 2 (2)</p><p>Functions of node in each layer of ANFIS architecture and so the functions of layers are following respectively [<xref ref-type="bibr" rid="scirp.73375-ref31">31</xref>] [<xref ref-type="bibr" rid="scirp.73375-ref32">32</xref>] [<xref ref-type="bibr" rid="scirp.73375-ref33">33</xref>] [<xref ref-type="bibr" rid="scirp.73375-ref34">34</xref>] [<xref ref-type="bibr" rid="scirp.73375-ref35">35</xref>] .</p><p>1<sup>st</sup> layer: Each node in this layer transfers the input signals to another layer without implementing any collecting or activation process.</p><p>2<sup>nd</sup> layer: Node shown by square in this layer represent A i and B i fuzzy clusters. The output values of these nodes are the membership levels bounded to input values and used membership functions (Equation (3)).</p><p>{ O i 2 = μ A i ( x ) O i 2 = μ B i ( y ) i = 1 , 2 (3)</p><p>There are total four nodes for both inputs in the second layer. In this layer generally continuous or partial triangle, trapezoidal or bell shaped curve membership functions is used as membership function.</p><p>Equation (4) which was formed by using a bell shape curve (Gaussian) membership function density function could be used in each node for A i and B i expressions in order to calculate 0-1 membership levels.</p><p>μ A i ( x ) = 1 1 + | x − m i σ | 2 (4)</p><p>Here m i and σ i show the mean and standard deviation of the bell shaped curve membership function respectively. Parameters which are used in the meaning of “premise parameter” are adjusted while the network was being trained in this layer.</p><p>3<sup>rd</sup> layer: Each node was labeled by π in this layer and represents the multiply of all input signals. Here each node outputs which represents the firing strength of each rule is calculated by Equation (5).</p><p>O i 3 = μ i = μ A i ( x ) μ B i ( y ) , i = 1 , 2 (5)</p><p>4<sup>th</sup> layer: Outputs of node shown by N labeled circle in this layer means the normalized threshold of the rules. Mentioned threshold could be calculated by Equation (6).</p><p>O i 4 = μ ′ i = μ i μ 1 + μ 2 , i = 1 , 2 (6)</p><p>5<sup>th</sup> layer: Each output of i node which was shown by square in this layer could be calculated Equation (7).</p><p>O i 5 = μ ′ i f i = μ ′ i ( p i x + q i y + r i ) (7)</p><p>In equity, μ ′ i is the output of fourth layer and shows normalized firing strength. p i , q i , r i Parameter sets are used in the meaning of consequent parameter.</p><p>6<sup>th</sup> layer: Node in this layer was expressed by circle and labeled by Σ . Total output consisted in this layer (f) is calculated by Equation (8) as the sum of all coming signals.</p><p>o 6 = f = ∑ i μ ′ i f i = ∑ i μ i f i ∑ i μ i (8)</p></sec></sec><sec id="s4"><title>4. Empirical Study</title><sec id="s4_1"><title>4.1. Data</title><p>Data of credit ratings of countries were compiled from the reports of three largest credit rating agencies as dependent variables. For that reason, the letter points of Moody’s, Standard &amp; Poor’s and Fitch were transformed to numeric ratings by using a scale transformation and their averages were taken.</p><p>A variable pool was organized for the factors affecting credit rating, so country risk after literature research result and data belonging to those variables and 2011-2013 years were collected for 180 countries over the world. Data belonging to three years is used as cross section data as separate units by not considering the time factor. The purpose of this increasing safe estimate ratio by providing ANFIS models better learning with more samples and neural networks. Due to deficiency of some data belonging to some countries, some countries were removed from the analysis. Data belonging to totally 230 units could be collected completely. Used variables and information where they obtained were submitted in Appendix A.</p></sec><sec id="s4_2"><title>4.2. Variable Selection Methodology</title><p>Factor analysis, linear discriminant analysis, square discriminant analysis, stepwise regression analysis implementations were carried out in order to decrease the numbers of variable by selecting from the variable pool which was formed after a literature search. As a result of factor analysis, determined four factors. But we realized that if the number of samples increases, factor loads will change. So we decided not to use factor analysis in variable selection. When common variables which were used in similar studies were considered study was continued by the variables obtained from the regression analysis. The data set should provide some assumption in parametric analysis. For that reason, variables were examined whether they are suitable for normal distribution and normality assumption was tried to provide by logarithmic transformations. Transformed variables were renamed by adding “LG” code. After the normality test, multi-collinearity problem examination was carried out and some variables which have higher Variance Inflation Factor (VIF) values were removed from the analysis.</p></sec><sec id="s4_3"><title>4.3. Credit Rating Prediction Models</title><p>Credit rating prediction models are generally used as classification or clustering except in this study. We investigate this problem as directly credit rating prediction as points, by using some determinants of ratings. There have been lots of determinants for the credit rating problem so we decided to decrease the number of determinants. After completing the variable selection process by using stepwise regression analysis; logistic regression analysis, artificial neural network and ANFIS model were implemented in order to determine the relationship of selected variables to credit rating. Finally, model performances were evaluated and predicts were made. In logistic regression analysis, which is implemented for credit rating estimate, the dependent variable was coded as (0 - 1) and possibility values which were calculated as analysis result were evaluated as default risks. 76% of 230 pcs data were shared for training, 13% were shared for test and 11% were shared for validation for ANN and ANFIS models. For assessing model performances, total correct classification percentage and mean absolute error (MAE) scales were used. After the analysis, credit rating estimates were made for the countries which haven’t got credit ratings and models were compared.</p></sec></sec><sec id="s5"><title>5. Empirical Results</title><sec id="s5_1"><title>5.1. Variable Selection Model Results</title><p>Variables determined by stepwise regression analysis in order to use for estimation models were submitted in <xref ref-type="table" rid="table1">Table 1</xref>.</p><p>Here; LGGCI is global competitiveness index, LGCAB% is current account balance % of export, LGRL is the rule of law, LGDebtS is debt service and LGDEF is GDP Deflator. All variables are logarithmic forms. The “grade” is the dependent variable in all models that represents credit rating of countries.</p></sec><sec id="s5_2"><title>5.2. Logistic Regression Results</title><p>At this stage, variables obtained from stepwise regression analysis were used as explanatory variables. The categorical dependent variable has been formed with the classification which was made by credit rating agencies such as “investible” and “non-investible”. For that reason, “1” has been designated for the units which has credit note over 60 and “0” has been designated for the other units. As a result of logistic regression analysis, group memberships of countries were estimated and non-default possibility (being investible country) has been calculated.</p><p>In <xref ref-type="table" rid="table2">Table 2</xref>, several explanatory coefficients used for model selection are seen and in <xref ref-type="table" rid="table3">Table 3</xref>, it is seen that the model ( χ 2 = 13.061 v e p &gt; 0.05 ) estimated</p><table-wrap id="table1" ><label><xref ref-type="table" rid="table1">Table 1</xref></label><caption><title> Stepwise regression coefficient</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Model (e)</th><th align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-1040522x21.png" xlink:type="simple"/></inline-formula></th><th align="center" valign="middle" >Std. Dev.</th><th align="center" valign="middle" >t</th><th align="center" valign="middle" >Sig.</th></tr></thead><tr><td align="center" valign="middle" >Constant</td><td align="center" valign="middle" >122,862</td><td align="center" valign="middle" >13,874</td><td align="center" valign="middle" >8856</td><td align="center" valign="middle" >0.000</td></tr><tr><td align="center" valign="middle" >LGGCI</td><td align="center" valign="middle" >−328,972</td><td align="center" valign="middle" >39,799</td><td align="center" valign="middle" >−8266</td><td align="center" valign="middle" >0.000</td></tr><tr><td align="center" valign="middle" >LGCAB%</td><td align="center" valign="middle" >11,193</td><td align="center" valign="middle" >2222</td><td align="center" valign="middle" >5038</td><td align="center" valign="middle" >0.000</td></tr><tr><td align="center" valign="middle" >LGRL</td><td align="center" valign="middle" >16,437</td><td align="center" valign="middle" >3460</td><td align="center" valign="middle" >4750</td><td align="center" valign="middle" >0.000</td></tr><tr><td align="center" valign="middle" >LGDebtS</td><td align="center" valign="middle" >−5931</td><td align="center" valign="middle" >1602</td><td align="center" valign="middle" >−3703</td><td align="center" valign="middle" >0.000</td></tr><tr><td align="center" valign="middle" >LGDEF</td><td align="center" valign="middle" >−7069</td><td align="center" valign="middle" >2914</td><td align="center" valign="middle" >−2426</td><td align="center" valign="middle" >0.016</td></tr></tbody></table></table-wrap><table-wrap id="table2" ><label><xref ref-type="table" rid="table2">Table 2</xref></label><caption><title> Model summary</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Step</th><th align="center" valign="middle" >−2 Log Likelihood</th><th align="center" valign="middle" >Cox &amp; Snell R<sup>2</sup></th><th align="center" valign="middle" >Nagelkerke R<sup>2</sup></th></tr></thead><tr><td align="center" valign="middle" >1</td><td align="center" valign="middle" >145,688</td><td align="center" valign="middle" >0.466</td><td align="center" valign="middle" >0.632</td></tr><tr><td align="center" valign="middle" >2</td><td align="center" valign="middle" >138,238</td><td align="center" valign="middle" >0.485</td><td align="center" valign="middle" >0.657</td></tr><tr><td align="center" valign="middle" >3</td><td align="center" valign="middle" >126,257</td><td align="center" valign="middle" >0.514</td><td align="center" valign="middle" >0.697</td></tr><tr><td align="center" valign="middle" >4</td><td align="center" valign="middle" >122,870</td><td align="center" valign="middle" >0.522</td><td align="center" valign="middle" >0.708</td></tr></tbody></table></table-wrap><table-wrap id="table3" ><label><xref ref-type="table" rid="table3">Table 3</xref></label><caption><title> Hosmer-Lemeshow test</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Step</th><th align="center" valign="middle" >Chi Square</th><th align="center" valign="middle" >D.f.</th><th align="center" valign="middle" >Sig.</th></tr></thead><tr><td align="center" valign="middle" >1</td><td align="center" valign="middle" >13,456</td><td align="center" valign="middle" >7</td><td align="center" valign="middle" >0.062</td></tr><tr><td align="center" valign="middle" >2</td><td align="center" valign="middle" >28,010</td><td align="center" valign="middle" >8</td><td align="center" valign="middle" >0.000</td></tr><tr><td align="center" valign="middle" >3</td><td align="center" valign="middle" >24,293</td><td align="center" valign="middle" >8</td><td align="center" valign="middle" >0.002</td></tr><tr><td align="center" valign="middle" >4</td><td align="center" valign="middle" >13,061</td><td align="center" valign="middle" >8</td><td align="center" valign="middle" >0.110</td></tr></tbody></table></table-wrap><p>in fourth step according to the results of the Hosmer Lomeshow test which was implemented in order to assess model data compatibleness.</p><p>In <xref ref-type="table" rid="table4">Table 4</xref>, below logistic regression model obtained by forward selection method is as Equations (9)-(10).</p><p>A correct classification percentage which was realized by using that model as 0.88. This value is significantly higher. (The correct classification percentage was realized 72% in [<xref ref-type="bibr" rid="scirp.73375-ref36">36</xref>] ANN model, 75% in [<xref ref-type="bibr" rid="scirp.73375-ref2">2</xref>] and 92.4% in [<xref ref-type="bibr" rid="scirp.73375-ref37">37</xref>] .) But the possibilities obtained from logistic regression model and credit notes comply with each other. In another saying, it is not proper to transform possibilities to credit notes. That reduces the advantage of the model.</p><p>L = ln p 1 − p = 11 , 122 − 1 , 982 X 1 + 3 , 328 X 2 + 2 , 908 X 3 − 68 , 59 X 4 (9)</p><p>p 1 − p = e 11 , 122 − 1 , 982 X 1 + 3 , 328 X 2 + 2 , 908 X 3 − 68 , 59 X 4 = e 11 , 122 e − 1 , 982 X 1 e 3 , 328 X 2 e 2 , 908 X 3 e − 68 , 59 X 4 (10)</p></sec><sec id="s5_3"><title>5.3. ANN Results</title><p>Here the relationship between dependent and independent variables was determined by the multi-layer perception (MLP). Graphic examinations revealed that credit note was in the same direction with RL, CAB and GCI variables and opposite direction with DebtS and DEF variables. Parameters about ANN architectural structure and obtained error values are provided in <xref ref-type="table" rid="table5">Table 5</xref>.</p></sec><sec id="s5_4"><title>5.4. ANFIS Results</title><p>The variables selected by stepwise regression analysis were used in ANFIS model. But ANFIS model doesn’t generate good results if there are more than 4 inputs due to the rule base enlarging as [<xref ref-type="bibr" rid="scirp.73375-ref30">30</xref>] . For that reason, in order to determine the variable which causes the highest error, tests were carried out with quart variable groups (five group) and it has been decided to use the GCI, CAB, RL and DebtS</p><table-wrap id="table4" ><label><xref ref-type="table" rid="table4">Table 4</xref></label><caption><title> Forward selection method results</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Step</th><th align="center" valign="middle" >Variable</th><th align="center" valign="middle" ><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/3-1040522x24.png" xlink:type="simple"/></inline-formula></th><th align="center" valign="middle" >Std. Dev.</th><th align="center" valign="middle" >Wald</th><th align="center" valign="middle" >D.f.</th><th align="center" valign="middle" >Sig</th><th align="center" valign="middle" >Exp( )</th></tr></thead><tr><td align="center" valign="middle" >4</td><td align="center" valign="middle" >LGDEF (X<sub>1</sub>) LGCAB (X<sub>2</sub>) LGRL (X<sub>3</sub>) LGGCI (X<sub>4</sub>) Constant</td><td align="center" valign="middle" >−1982 3328 2908 −68,590 11,122</td><td align="center" valign="middle" >1134 1028 1143 16,583 4966</td><td align="center" valign="middle" >3056 10,487 6477 17,113 5017</td><td align="center" valign="middle" >1 1 1 1 1</td><td align="center" valign="middle" >0.040 0.001 0.011 0.000 0.025</td><td align="center" valign="middle" >0.138 27,892 18,322 0.000 67,670,180</td></tr></tbody></table></table-wrap><table-wrap id="table5" ><label><xref ref-type="table" rid="table5">Table 5</xref></label><caption><title> ANN model parameters and error results</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Network type</th><th align="center" valign="middle" >MLP</th></tr></thead><tr><td align="center" valign="middle" >Number of layers</td><td align="center" valign="middle" >4</td></tr><tr><td align="center" valign="middle" >Number of hidden layers</td><td align="center" valign="middle" >2</td></tr><tr><td align="center" valign="middle" >Input activation function</td><td align="center" valign="middle" >Tangent sigmoid</td></tr><tr><td align="center" valign="middle" >Output activation function</td><td align="center" valign="middle" >Linear</td></tr><tr><td align="center" valign="middle" >Iteration</td><td align="center" valign="middle" >2000</td></tr><tr><td align="center" valign="middle" >Learning rate</td><td align="center" valign="middle" >0.1</td></tr><tr><td align="center" valign="middle" >Momentum constant</td><td align="center" valign="middle" >0.7</td></tr><tr><td align="center" valign="middle" >MAE% for training</td><td align="center" valign="middle" >12.80</td></tr><tr><td align="center" valign="middle" >MAE% for test</td><td align="center" valign="middle" >10.01</td></tr><tr><td align="center" valign="middle" >MAE% for validation</td><td align="center" valign="middle" >12.87</td></tr></tbody></table></table-wrap><p>variables. Also model with five variables were trained by two membership function, but good results were not taken.</p><p>When mean absolute percentage errors (%MAE) related to ANFIS model in <xref ref-type="table" rid="table6">Table 6</xref> were examined and compared with in <xref ref-type="table" rid="table5">Table 5</xref>, it is understood that it was trained well compared to ANN and could make better estimations with lower error percentage.</p><p><xref ref-type="fig" rid="fig1">Figure 1</xref> and <xref ref-type="fig" rid="fig2">Figure 2</xref> shows the training errors trend of ANFIS model addition to MAE statistics in <xref ref-type="table" rid="table6">Table 6</xref> which is an indicator of whether the network was trained or not. As seen both errors in percentage (E%) and absolute percentage errors (AE%) trends approach to zero by decrease rapidly. The horizontal axis represents the sample number.</p><p>Also, it is seen in the graphics that errors belong to 60 and 62 numbered units are extremely large. When these units were excluded to obtain a more realistic value, MAE% value for training data was reduced 0.0771 to 0.0684.</p><p>While prediction success was 87.13% in ANN model, the higher accuracy percentage was obtained as 90.66% in ANFIS model. Here it is seen that ANFIS model can comprehend nonlinear relationships between variables successfully.</p></sec></sec><sec id="s6"><title>6. Conclusions</title><p>The sovereign credit rating is very important for country and other issuers, macro and micro level. Investors want to make their investment in the country have</p><table-wrap id="table6" ><label><xref ref-type="table" rid="table6">Table 6</xref></label><caption><title> ANFIS model parameters and error results</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Network type</th><th align="center" valign="middle" >ANFIS (Sugeno type)</th></tr></thead><tr><td align="center" valign="middle" >Number of layers</td><td align="center" valign="middle" >6</td></tr><tr><td align="center" valign="middle" >Iteration</td><td align="center" valign="middle" >250</td></tr><tr><td align="center" valign="middle" >Input membership function</td><td align="center" valign="middle" >type</td></tr><tr><td align="center" valign="middle" >Output function</td><td align="center" valign="middle" >Constant</td></tr><tr><td align="center" valign="middle" >Number of membership functions</td><td align="center" valign="middle" >3-3-3-3</td></tr><tr><td align="center" valign="middle" >Number of fuzzy rules</td><td align="center" valign="middle" >81</td></tr><tr><td align="center" valign="middle" >Optimization algorithm</td><td align="center" valign="middle" >Hybrid (Back prob. and LSE)</td></tr><tr><td align="center" valign="middle" >“and” method</td><td align="center" valign="middle" >prod</td></tr><tr><td align="center" valign="middle" >“or” method</td><td align="center" valign="middle" >probor</td></tr><tr><td align="center" valign="middle" >Clarification method</td><td align="center" valign="middle" >wtaver</td></tr><tr><td align="center" valign="middle" >%MAE (training)</td><td align="center" valign="middle" >7.71</td></tr><tr><td align="center" valign="middle" >%MAE (test)</td><td align="center" valign="middle" >13.70</td></tr><tr><td align="center" valign="middle" >%MAE (validation)</td><td align="center" valign="middle" >9.34</td></tr></tbody></table></table-wrap><p>a good rating. Good rating means, lower stock cost, lower interest rate, lower investment cost and higher profit for any country or firm. In this research LR, ANN and ANFIS models were compared and credit rating predictions were made for the countries which have not. When several studies were examined about credit rating, these are seen as micro rating applications such as company or bond rating. Most of these studies handled the problem as a classification problem. Designed models here handled the credit rating problem as country risk graded and configured directly for predicting credit rating estimation. This study is considered to contribute to literature in respect to trying different statistical methods during variable selection, determining ANFIS model as the best among several prediction methods, success of credit rating prediction and risks of 21 countries firstly measured.</p><p>Besides the main targets of this study, another examined issue is whether group membership possibilities could have transformed into credit ratings. Even group membership shows significantly good predictions for countries which have higher credit ratings; it was not same for the countries which have lower credit ratings. It is clear that LR credit rating prediction which has 88% correct accuracy classification success is unsuccessful in respect to credit rating estimation.</p><p>ANN and ANFIS are the models which learn the relationships over case studies and generate predictions. For that reason, it should be noted that as the number of samples increases, the correct estimation rates increase. Also, it is important that obtained data should be reliable. Relationships learned from incorrect data cause incorrect predictions. These models are considered useful to determine more complicated relationships instead of simple linear relationships. But due to these complicated relationships, it cannot be expressed by a simple equation; prediction equation could not be revealed after analysis.</p><p>ANFIS method which was the best among three prediction models with 91% accuracy percentage was used for first credit rating predictions of 21 countries (Appendix B). Thereby the developed ANFIS model proved its prediction success on the verification data set. Consequently, it was revealed that ANFIS model is the most proper model in order to measure country risks and assign credit ratings and it could be used trustfully.</p><p>ANFIS model is selected as the best prediction model because of its mathematical hybrid structure and learning algorithm. All possible structure and parameters are tried for LR and ANN. For the future research, it's possible to apply another artificial intelligence technic or heuristic-metaheuristic search technic for variable selection and credit rating prediction.</p></sec><sec id="s7"><title>Cite this paper</title><p>Pabu&#231;cu, H. and Ayan, T.Y. (2017) The Development of an Alternative Method for the Sovereign Credit Rating System Based on Adaptive Neuro-Fuzzy Inference System. American Journal of Operations Research, 7, 41-55. http://dx.doi.org/10.4236/ajor.2017.71003</p></sec><sec id="s8"><title>Appendix A</title></sec><sec id="s9"><title>Appendix B</title></sec></body><back><ref-list><title>References</title><ref id="scirp.73375-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Cantor, R. and Packer, F. (1996) Determinants and Impact of Sovereign Credit Ratings. Journal of Fixed Income, 6, 76-91. https://doi.org/10.3905/jfi.1996.408185</mixed-citation></ref><ref id="scirp.73375-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Ratha, D., De, P.K. and Mohapatra, S. (2011) Shadow Sovereign Ratings for Unrated Developing Countries. World Development, 39, 295-307.  
https://doi.org/10.1016/j.worlddev.2010.08.006</mixed-citation></ref><ref id="scirp.73375-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Bennell, J.A., Crabbe, D., Thomas, S. and Gwilym, O.A. (2005) Modelling Sovereign Credit Ratings: Neural Networks versus Ordered Probit. Expert System with Applications, 30, 415-425. https://doi.org/10.1016/j.eswa.2005.10.002</mixed-citation></ref><ref id="scirp.73375-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Subasi, A. and Gursoy, I.M. (2010) EEG Signal Classification Using PCA, ICA, LDA and Support Vector Machines. Expert System with Applications, 37, 8659-8666.  
https://doi.org/10.1016/j.eswa.2010.06.065</mixed-citation></ref><ref id="scirp.73375-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Harris, T. (2015) Credit Scoring Using the Clustered Support Vector Machine. Expert System with Applications, 42, 741-750.  
https://doi.org/10.1016/j.eswa.2014.08.029</mixed-citation></ref><ref id="scirp.73375-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Jones, S., Johnstone, D. and Wilson, R. (2015) An Empirical Evaluation of the Performance of Binary Classifiers in the Prediction of Credit Ratings Changes. Journal of Banking &amp; Finance, 56, 72-85. https://doi.org/10.1016/j.jbankfin.2015.02.006</mixed-citation></ref><ref id="scirp.73375-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Beers, T.D. and Cavanaugh, M. (2004) Sovereign Credit Ratings: A Primer. Standard &amp; Poor’s, New York.</mixed-citation></ref><ref id="scirp.73375-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Bissoondoyal-Bheenick, E. (2005) An Analysis of the Determinants of Sovereign Ratings. Global Finance Journal, 15, 251-280.  
https://doi.org/10.1016/j.gfj.2004.03.004</mixed-citation></ref><ref id="scirp.73375-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Butler, A.W. and Fauver, L. (2006) Institutional Environment and Sovereign Credit Ratings. Financial Management, 35, 53-79.  
https://doi.org/10.1111/j.1755-053X.2006.tb00147.x</mixed-citation></ref><ref id="scirp.73375-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">Depken, C., LaFountain, C. and Butters, R. (2007) Corruption and Creditworthiness: Evidence from Sovereign Credit Ratings. Working Papers 0601, Department of Economics, University of Texas at Arlington.</mixed-citation></ref><ref id="scirp.73375-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">Busse, M. and Hefeker, C. (2007) Political Risk, Institutions and Foreign Direct Investment. European Journal of Political Economy, 23, 397-415.  
https://doi.org/10.1016/j.ejpoleco.2006.02.003</mixed-citation></ref><ref id="scirp.73375-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">Elkhoury, M. (2009) Credit Rating Agencies and Their Potential Impact on Developing Countries. Compendium on Debt Sustainability, 165-180.</mixed-citation></ref><ref id="scirp.73375-ref13"><label>13</label><mixed-citation publication-type="other" xlink:type="simple">Roychoudhury, S. and Lawson, R.A. (2010) Economic Freedom, the Credit Ratings and the Sovereign default spread. Journal of Financial Economic Policy, 2, 149-162.  
https://doi.org/10.1108/17576381011070201</mixed-citation></ref><ref id="scirp.73375-ref14"><label>14</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Bundala</surname><given-names> N.N. </given-names></name>,<etal>et al</etal>. (<year>2013</year>)<article-title>Do Economic Growth, Human Development and Political Stability Favour Sovereign Creditworthiness of a Country? A Cross Country Survey on Developed and Developing Countries</article-title><source> International Journal of Advances in Management and Economics</source><volume> 1</volume>,<fpage> 32</fpage>-<lpage>46</lpage>.<pub-id pub-id-type="doi"></pub-id></mixed-citation></ref><ref id="scirp.73375-ref15"><label>15</label><mixed-citation publication-type="other" xlink:type="simple">Williams, G., Alsakka, R. and Gwilym, O. (2013) The Impact of Sovereign Rating Actions on Bank Ratings in Emerging Markets. Journal of Banking &amp; Finance, 37, 563-577. https://doi.org/10.1016/j.jbankfin.2012.09.021</mixed-citation></ref><ref id="scirp.73375-ref16"><label>16</label><mixed-citation publication-type="other" xlink:type="simple">Leonard, K.J. (1995) The Development of Credit Scoring Quality Measures for Consumer Credit Applications. International Journal of Quality &amp; Reliability Management, 12, 79-85. https://doi.org/10.1108/02656719510087346</mixed-citation></ref><ref id="scirp.73375-ref17"><label>17</label><mixed-citation publication-type="other" xlink:type="simple">Hajek, P. and Michalak, K. (2013) Feature Selection in Corporate Credit Rating Prediction. Knowledge-Based Systems, 51, 72-84.  
https://doi.org/10.1016/j.knosys.2013.07.008</mixed-citation></ref><ref id="scirp.73375-ref18"><label>18</label><mixed-citation publication-type="other" xlink:type="simple">Li, W., He, Q., Fan, X. and Fei, Z. (2012) Evaluation of Driver Fatigue on Two Channels of EEG Data. Neuroscience Letters, 506, 235-239.  
https://doi.org/10.1016/j.neulet.2011.11.014</mixed-citation></ref><ref id="scirp.73375-ref19"><label>19</label><mixed-citation publication-type="other" xlink:type="simple">Zhao, C., Zheng, C., Zhao, M., Tu, Y. and Liu, J. (2011) Multivariate Autoregressive Models and Kernel Learning Algorithms for Classifying Driving Mental Fatigue Based on Electroencephalographic. Expert System with Applications, 38, 1859-1865.  
https://doi.org/10.1016/j.eswa.2010.07.115</mixed-citation></ref><ref id="scirp.73375-ref20"><label>20</label><mixed-citation publication-type="other" xlink:type="simple">Boostani, R., Sadatnezhad, K. and Sabeti, M. (2009) An Efficient Classifier to Diagnose of Schizophrenia Based on the EEG Signals. Expert System with Applications, 36, 6492-6499. https://doi.org/10.1016/j.eswa.2008.07.037</mixed-citation></ref><ref id="scirp.73375-ref21"><label>21</label><mixed-citation publication-type="other" xlink:type="simple">Huang, Z., Chen, H., Hsu, C.J., Chen, W.H. and Wu, S. (2004) Credit Rating Analysis with Support Vector Machines and Neural Networks: A Market Comparative Study. Decision Support Systems, 37, 543-558.  
https://doi.org/10.1016/S0167-9236(03)00086-1</mixed-citation></ref><ref id="scirp.73375-ref22"><label>22</label><mixed-citation publication-type="other" xlink:type="simple">Kim, K.S. and Han, I. (2001) The Cluster-Indexing Method for Case-Based Reasoning Using Self-Organizing Maps and Learning Vector Quantization for Bond Rating Cases. Expert System with Applications, 21, 147-156.  
https://doi.org/10.1016/S0957-4174(01)00036-7</mixed-citation></ref><ref id="scirp.73375-ref23"><label>23</label><mixed-citation publication-type="other" xlink:type="simple">Gujarati, D.N. (2003) Student Solutions Manual for Use with Basic Econometrics. 4th Edition, McGraw Hill, New York.</mixed-citation></ref><ref id="scirp.73375-ref24"><label>24</label><mixed-citation publication-type="other" xlink:type="simple">Hosmer, D.W. and Lemeshow, S. (2004) Applied Logistic Regression. 2th Edition, Wiley, New York.</mixed-citation></ref><ref id="scirp.73375-ref25"><label>25</label><mixed-citation publication-type="other" xlink:type="simple">Hand, D.J. and Henley, W.E. (1997) Statistical Classification Methods in Consumer Credit Scoring: A Review. Journal of Royal Statistical Society, 160, 523-541.  
https://doi.org/10.1111/j.1467-985X.1997.00078.x</mixed-citation></ref><ref id="scirp.73375-ref26"><label>26</label><mixed-citation publication-type="other" xlink:type="simple">Yu, L., Wang, S. and Lai, K.K. (2008) Forecasting Crude Oil Price with an EMD-Based Neural Network Ensemble Learning Paradigm. Energy Economics, 30, 2623-2635. https://doi.org/10.1016/j.eneco.2008.05.003</mixed-citation></ref><ref id="scirp.73375-ref27"><label>27</label><mixed-citation publication-type="other" xlink:type="simple">Smith, K.A. and Gupta, J.N.D. (2002) Neural Networks in Business: Techniques and Applications. IGI Global, New York.</mixed-citation></ref><ref id="scirp.73375-ref28"><label>28</label><mixed-citation publication-type="other" xlink:type="simple">Jardin, P. and Séverin, E. (2012) Forecasting Financial Failure Using a Kohonen Map: A Comparative Study to Improve Model Stability over Time. European Journal of Operational Research, 221, 378-396.  
https://doi.org/10.1016/j.ejor.2012.04.006</mixed-citation></ref><ref id="scirp.73375-ref29"><label>29</label><mixed-citation publication-type="other" xlink:type="simple">Takagi, T. and Sugeno, M. (1985) Fuzzy Identification of Systems and Its Applications to Modeling and Control. IEEE Transactions on Systems Man &amp; Cybernetics, 15, 116-132. https://doi.org/10.1109/TSMC.1985.6313399</mixed-citation></ref><ref id="scirp.73375-ref30"><label>30</label><mixed-citation publication-type="other" xlink:type="simple">Jang, J.S.R. (1993) ANFIS: Adaptive-Network-Based Fuzzy Inference system. IEEE Transactions on Systems Man &amp; Cybernetics, 23, 665-685.  
https://doi.org/10.1109/21.256541</mixed-citation></ref><ref id="scirp.73375-ref31"><label>31</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Jang</surname><given-names> J. </given-names></name>,<etal>et al</etal>. (<year>1991</year>)<article-title>Fuzzy Modeling Using Generalized Neural Networks and Kalman Filter Algorithm</article-title><source> Proceeding of 9th National Conference Artificial Intelligence</source><volume> 91</volume>,<fpage> 762</fpage>-<lpage>767</lpage>.<pub-id pub-id-type="doi"></pub-id></mixed-citation></ref><ref id="scirp.73375-ref32"><label>32</label><mixed-citation publication-type="other" xlink:type="simple">Jang, J.S.R. (1993) ANFIS: Adaptive-Network-Based Fuzzy Inference System. IEEE Transactions on Systems Man &amp; Cybernetics, 23, 665-685.  
https://doi.org/10.1109/21.256541</mixed-citation></ref><ref id="scirp.73375-ref33"><label>33</label><mixed-citation publication-type="other" xlink:type="simple">Jang, J.S.R. (1996) Input Selection for ANFIS learning. Proceeding IEEE 5th International Fuzzy Systems, 2, 1493-1499. https://doi.org/10.1109/FUZZY.1996.552396</mixed-citation></ref><ref id="scirp.73375-ref34"><label>34</label><mixed-citation publication-type="other" xlink:type="simple">Fullér, R. (1995) Neural Fuzzy System.</mixed-citation></ref><ref id="scirp.73375-ref35"><label>35</label><mixed-citation publication-type="other" xlink:type="simple">Lin, C.T. and Lee, C.S.G. (1991) Neural-Network-Based Fuzzy Logic Control and Decision System. IEEE Transaction on Computers, 40, 1320-1336.  
https://doi.org/10.1109/12.106218</mixed-citation></ref><ref id="scirp.73375-ref36"><label>36</label><mixed-citation publication-type="other" xlink:type="simple">Leshno, M. and Spector, Y. (1996) Neural Network Prediction Analysis: The Bankruptcy Case. Neurocomputing, 10, 125-147.  
https://doi.org/10.1016/0925-2312(94)00060-3</mixed-citation></ref><ref id="scirp.73375-ref37"><label>37</label><mixed-citation publication-type="other" xlink:type="simple">Blanco, A., Pino-Mejías, R., Lara, J. and Rayo, S. (2013) Credit Scoring Models for the Microfinance Industry Using Neural Networks: Evidence from Peru. Expert System with Applications, 40, 356-364. https://doi.org/10.1016/j.eswa.2012.07.051</mixed-citation></ref></ref-list></back></article>