Estimating corporate bankruptcy forecasting models by maximizing discriminatory power

In this paper, we estimate coefficients of bankruptcy forecasting models, such as logistic and neural network models, by maximizing their discriminatory power as measured by the Area Under Receiver Operating Characteristics (AUROC) curve. A method is introduced and compared with traditional logistic and neural network models, using out-of-sample analysis, in terms of discriminatory power, information content and economic impact while we forecast bankruptcy one year ahead, two years ahead but also financial distress, which is a situation that precedes firm bankruptcy. Using US public firms over the period 1990–2015, in all, we find that training models to maximize AUROC, provides more accurate out-of-sample forecasts relative to training them with traditional methods, such as maximizing the log-likelihood function, highlighting the benefits arising by using models with maximized AUROC. Among all models, however, a neural network trained with our method is the best performing one, even when we compare it with other methods proposed in the literature to maximize AUROC. Finally, our results are more pronounced when we increase the forecasting difficulty, such as forecasting financial distress. The implementation of our method to train bankruptcy models is robust in various settings and therefore well-justified.


Background and motivation
Increased attention has been paid in recent years for the development of powerful bankruptcy forecasting models, mainly for two reasons. First, the recent global financial crisis in [2007][2008][2009] has left banks to experience huge losses from their credit portfolios and consequently their lending policies and decision-making processes have been seriously criticized from regulators, investors and other stakeholders. Second, since the reform of Basel Accord in 2006, banks can develop their own internal models to assess credit risks and protect themselves through the capital reserves that should withhold to face potential losses. Thus, for a matter of bank viability, 1 financial stability and investor protection, it would be of great interest to develop powerful bankruptcy forecasting models, which is the aim of this paper.
One of the most significant measures to evaluate the performance of bankruptcy forecasting models is their ability to discriminate bankrupt from healthy firms. It has been shown that models with higher discriminatory power are associated with higher economic benefits for a bank (Bloechlinger et al. 2006;Agarwal et al. 2008). Furthermore, Bauer et al. (2014) show that even small differences in the discriminatory power among bankruptcy forecasting models yield superior bank economic performance. In addition, commercial vendors and industry experts, such as Moody's KMV, extensively use discriminatory power as an integral part of their validation processes, especially when comparing their newly developed models with existing ones (see for instance the RiskCalc 3.1 model in Dwyer et.al. 2004). As it is stated in their paper: "The greatest contribution to profitability, efficiency and reduced losses comes from the models' powerful ability to rank-order firms by riskiness so that the bank can eliminate high risk prospects." Beyond that, Moody's KMV provides ample explanatory documentation on how to use various discriminatory power measures in practice (see for instance Keenan et al. 1999 andSobehart et al. 2000) and it is also extensively used in academic research to compare various bankruptcy forecasting models. This extensive use from practitioners and academics alike, in fact, highlights the importance of using discriminatory power as a leading measure to evaluate the performance of bankruptcy forecasting models.
Despite the empirical evidence on the economic benefits arising by using models with higher discriminatory power, it is somewhat surprising that a common practice in bankruptcy forecasting is to use discriminatory power only ex-post as an indication of model performance, rather than obtaining model coefficients directly by maximizing discriminatory power. Exceptions include Miura et al. (2010) and Kraus et al. (2014) in the related area of credit scoring which we discuss and compare with our method. We contribute to this limited literature by introducing a method that we use to train bankruptcy forecasting models, such as logistic and neural network models and comparing these models with traditional logistic and neural network models, such as models which maximize the log-likelihood function. Ultimately, our goal in this study is to highlight the importance of using models which are trained to maximize discriminatory power.
To measure discriminatory power, we use the Area Under Receiver Operating Characteristics curve (AUROC or AUC). This is a widely used statistic that has been employed by many studies recently to compare discriminatory power of various bankruptcy forecasting models (including Chava et al. 2004;Campbell et al. 2008;Tinoco et al. 2013;Filipe et al. 2016 and many others). Furthermore, it has been used in related areas, such as mortgage default prediction (Fitzpatrick et al. 2016) and generally when assessing the performance of credit scoring models (see for instance Lessmann et al. 2015). Moreover, the AUROC is an appealing measure because it is easy to interpret and compute empirically and importantly, it does not depend on cut-off values, such as those needed when constructing the standard confusion matrices. Instead, the AUROC simply summarizes discriminatory power in a single number, thus it is easy to compare across various models, without using such cut-off values, which is the main reason it is has received considerable attention in bankruptcy studies. Due to these reasons, we select AUROC as the optimization criterion and we develop a method which seeks to maximize AUROC.
For our main analysis we collect annual financial data and daily equity prices for a large sample of U.S. public bankrupt and healthy firms and construct variables to make one-year forecasts, two-year forecasts and finally, we forecast financial distress which is a situation prior to the formal bankruptcy filing, over the period 1990-2015. We keep approximately 70% of the whole sample as a training set and evaluate the performance of the models in the testing set using three distinct type of tests, following Bauer et al. (2014); 1) AUROC analysis 2) Information content tests 3) Economic performance, when banks use various bankruptcy forecasting models in a competitive loan market.

Main findings
First, we employ standard statistical analysis to select few predictive variables, from a pool of variables, which individually exhibit high discriminatory power, have low correlation from each other and are statistically significant. This is to eliminate insignificant variables that may add noise and helps us constructing parsimonious models. When we consider only financial variables, we find that several financial variables related to firm leverage, profitability, liquidity and coverage, are significant predictors of bankruptcy. When we also consider market-based variables in the analysis, however, the model with both financial and market variables outperforms the model with only financial variables, consistent with prior research (Shumway 2001;Chava et al. 2004;Campbell et al. 2008;Wu et al. 2010 andTinoco et al. 2013). These two selected sets of variables (financial variables and financial with market variables) are the inputs to all models (logistic and neural networks trained to maximize AUROC and the log-likelihood function).
We begin our analysis by evaluating and comparing the out-of-sample performance of logistic and neural network models trained with our method, to those trained to maximize the log-likelihood function, one year ahead, two years ahead and finally, when we forecast financial distress. Overall, we find that our proposed method yields logistic and neural network models which outperform, out-of-sample, traditional logistic and neural network models. The results with respect to the three testing approaches suggest that models with maximized AUROC 1) Significantly outperform the traditional models in terms of their ability to discriminate bankrupt from healthy firms, 2) They provide significantly more information about future bankruptcy-financial distress relative to traditional methods 3) Banks using models with maximized AUROC earn superior returns on a risk-adjusted basis relative to banks that use traditional models to forecast bankruptcy-financial distress. From all models, however, our neural network is the best performing one. In addition, the results are more pronounced in the case of financial distress.
Next, we compare our method with other methods proposed in the literature to maximize AUROC. Using our proposed neural network model as a representative since it is the best performing model in all tests, we find that it outperforms the alternative AUROC maximization methods proposed by Miura et al. (2010) and Kraus et al. (2014). This result is more pronounced when forecasts are performed two years ahead and to the case of financial distress.
Finally, we compare, out-of-sample, the discriminating ability of logistic and neural networks trained by maximizing the AUROC, to the models trained with the traditional approach but this time, the input variables are constructed using quarterly data. In this way, we update the models as new information becomes available with higher frequency. In all, findings advocate the implementation of our estimation method since it provides better prediction performance, in terms of discriminatory power, relative to traditional estimation methods.
Our paper has implications in the way bankruptcy analysis is conducted and aims towards better decision-making through more accurate bankruptcy forecasts. First, our study can be viewed as a way towards improving the general practice of bankruptcy forecasting by providing an alternative estimation technique to obtain model coefficients relative to obtaining them with traditional methods. To this end, our proposed estimation method significantly improves performance, out-of-sample, especially when we increase the forecasting difficulty, such as forecasting bankruptcy two years in the future but also forecasting financial distress, which is a situation before the formal firm bankruptcy. In addition, our paper provides an extended methodological framework to commonly used traditional bankruptcy models such as Altman (1968); Ohlson (1980) but also more recent ones, such as Campbell et al. (2008) and many other similar models, by introducing a new optimization method to obtain their coefficients and increase forecasting accuracy. Finally, the advantage of our method is that it works well using any modelling approach where the output is a probability, thus it retains the same interpretability with the outputs of traditional estimation methods. This is also in contrast with the methodologies proposed by Miura et al. (2010) and Kraus et al. (2014) as these methods cannot be used by logistic and neural network models which are two of the most popular bankruptcy forecasting models (se for instance Kumar et al. 2007 and references therein). To the best of our knowledge, this is the first time such extensive work is performed to compare maximizing AUROC with the traditional maximizing the log-likelihood function for bankruptcy forecasting models and highlighting the benefits arising by AUROC-maximized models.
The remainder of the paper proceeds as follows: In Sect. 2 we discuss data collection, in Sect. 3 we present the methodology to maximize AUROC as well as the three distinct type of tests we use to evaluate performance, in Sect. 4 we discuss the results and Sect. 5 concludes.

Sample
Our sample consists of 11,096 non-financial U.S. firms from which 422 filed for bankruptcy under Chapter 7 or Chapter 11 between 1990 and 2015. We have a total of 97,133 2 firm-year observations with non-missing data to forecast bankruptcies using the corresponding data which are lagged by one or two years for our one or two year ahead forecasts respectively but also, to forecast financial distress. 3 Bankrupt firms and the date of their bankruptcy filing were identified from BankruptcyData, which is a comprehensive database containing corporate bankruptcy and distressed information for firms in the US. Table 1 reports the frequency of bankrupt and healthy firms (i.e. non-bankrupt firms) collected each year over the sample period spanning the years 1990-2015.
Since our main "distress" event is bankruptcy, we treat exits unrelated to bankruptcy as non-bankrupt observations 4 (i.e. healthy firms) and we report them in Table 1. In particular, Table 1 breaks down the healthy firms, each year, into three categories; 1) active firms are firms that survived during the year, 2) firms that stopped filing information due being merged or acquired (M&As) and 3) firms that stopped filing information for other reasons including conversion to private company, engaged in a levarage buyout etc. The delisting reasons was found in COMPUSTAT using the DLRSN variable which provides codes for each delisting reason.
Next, Fig. 1 presents graphically the yearly number of bankruptcies, to visualize the variation of the number of bankruptcies over our sample period. Figure 1 shows that bankruptcies peak in three major time-periods: 1) during the 1990-1991 US crisis, 2) during the dot-com bubble occurred around 2000 and 3) during the financial crisis period when bankruptcies peaked in 2009. Overall, the plot shows that the sample period we use captures the prevailing market conditions with higher (lower) number of bankruptcies during crisis (normal) periods.

Variables construction
We collect annual financial data and market (equity) data from Compustat and CRSP respectively and we construct several variables based on related studies in the literature. For example, in our analysis we consider variables used in traditional corporate bankruptcy studies, such as Altman (1968), Ohlson (1980), Zmijewski (1984) but also in more recent studies, such as Shumway (2001), Chava et al. (2004), Campbell et al. (2008) etc.
First, we construct financial ratios capturing aspects of a firm's financial performance, such as leverage, profitability, liquidity, coverage, activity, cash flows, as presented in panel A of Table 2. A limitation of financial variables is that by their nature look backwards and the quality of information they carry depends on accounting practices (Hillegeist et al. 2004;Agarwal et al. 2008). Market variables, instead, constructed from equity prices, are forward-looking since they carry market perceptions about the prospects of the firm. For publicly traded firms it would be more appropriate to incorporate market variables in the 2 Possibly, there are two main reasons that make our sample to differ from other studies, such as Campbell et al. (2008). First, we exclude financial firms with SIC codes in the range of 6000 which is something commonly done in the literature (corresponds to more than 8,000 unique firms over our sample period). Second, we delete observations with missing data whereas Campbell et al. (2008) replace missing data with sample averages. 3 We present more details about the financial distress case in Sect. 4.5. 4 In other words, unhealthy firms that did not file for bankruptcy but eventually exited the sample, are considered as healthy observation. Failure to account for this, it will overestimate the predictive performance. models. To this end, we collect daily equity prices from CRSP for the entire fiscal year and several market-based variables are constructed, as reported in panel B of Table 2. Annualized volatility of daily equity returns (VOLE) refers to the fluctuations of firm's equity value returns, expecting to be higher for bankrupt firms. Next, excess return (EXRET) refers to the difference between firm's annualized equity return and the annualized valueweighted return of a portfolio with NYSE, AMEX, NASDAQ stocks, expecting to be lower for bankrupt firms. Further, we consider the relative size of the firm (RSIZE), the logarithm of stock price (LOGPRICE) and the Market-to-Book ratio (MB), expecting a negative association with bankruptcy risk. Finally, we include three financial variables scaled by firm's market value. More precisely, Campbell et al. (2008) show that scaling financial variables with a market-based measure of firm's value i.e. market equity + liabilities (MTA), compared to total assets as reported in the balance sheet, increases the predictive accuracy of bankruptcy forecasting models. These variables are cash over MTA (CASHMTA), net income over MTA (NIMTA), expecting a negative association with bankruptcy risk and lastly, total liabilities over MTA (TLMTA). Following common practice, we winsorize the variables between 1st and 99th percentile to avoid problems induced by outliers. Table 2 presents an extensive list of variables that previous studies found to be significant predictors of bankruptcy risk. Out of these variables, a smaller set should be selected in order to construct parsimonious models with few variables but with high forecasting power. We establish a three-step approach to select the most powerful variables (see for instance Altman et al. 2007 andFilipe et al. 2016) and summarized in the following three steps:

Variables selection
Step 1: Removing variables with low discriminating ability (as a cut-off, we use AUROC equal to 0.60). The idea of this step is to qualify the variables that individually exhibit a satisfactory ability to discriminate bankrupt from healthy firms.
Step 2: Removing highly correlated variables using the Variance Inflation Factor (VIF) criterion. The idea of this step is to remove the variables that are highly correlated with others, since multicollinearity may yield misleading results regarding the significance of the variables in the final model. Beyond that, we end up with variables that provide  different information and explain bankruptcy uniquely. We use 5 as cut-off (variables with VIF ≥ 5 are removed).
Step 3: Performing a stepwise multivariate logistic regression to the remaining variables in order to obtain the most significant variables from a statistical point of view (we use a significance level of α = 5%). The logistic regression program estimates coefficients assuming independent observations, which is an invalid assumption, since the data contains information for firms over multiple periods. In such case, an appropriate correction measure which we adopt in our study, is to use clustered robust standard errors (also used by Filipe et al. 2016).
Using the three-step approach, we develop two types of models. The first one is a "private firm" type of model, including only financial variables. We further develop a "public firm" type of model, including both financial and market variables. For example, the private firm model includes five financial variables (TLTA, STDTA, NITA, CASHTA, EBITCL), while the public firm model includes six variables (TLTA, STDTA, LOGPRICE, CASHMTA, NIMTA, EXRET). Notice that two financial-based variables (CASHTA and NITA) are replaced with CASHMTA and NIMTA. Generally, the majority of variables that are found to be significant for the public firm model are market variables, which is consistent with the perception that market-based variables are better bankruptcy risk measures, due to their forward-looking nature. These two sets of variables are the inputs to all models (i.e. used in the models which we train to maximize AUROC and in the models trained to maximize the log-likelihood function). For simplicity, these two sets of variables are also used in the models when forecasting bankruptcy two years ahead and when forecasting financial distress. Table 3 reports descriptive statistics for the accounting and market variables that we find to be significant predictors of bankruptcy. As expected, bankrupt firms are more levered on average relative to healthy firms (TLTA and STDTA for bankrupt firms are higher), they are also less profitable (NITA and NIMTA are lower for bankrupt firms). Furthermore, bankrupt firms are more constrained in terms of cash available (CASHTA and CASHMTA are lower) as opposed to healthy firms. Going to the market variables, it is evident that the stock price of bankrupt firms (LOGPRICE) on average is lower than healthy firms, possibly due to their deteriorating financial position that is priced by investors, leading to a depreciation of their stock prices at the year prior to bankruptcy. Finally, bankrupt firms exhibit lower and negative market performance relative to the market (EXRET is lower one year prior to bankruptcy), as opposed to healthy firms.

Measuring discriminatory power
Discriminatory power refers to the ability of a model to discriminate bankrupt from healthy firms. According to a cut-off score, firms whose bankruptcy score exceeds the cut-off are classified as bankrupt and healthy otherwise. Therefore, a way to measure the discriminating ability of a model is, for a given cut-off score, to count the true forecasts (percentage of bankrupt firms correctly classified as bankrupt) and the false forecasts (percentage of healthy firms incorrectly classified as bankrupt). Doing this process with multiple cut-offs, we get a set of true and false forecasts. A graph made from this set is the ROC curve with false forecasts on the x-axis and true forecasts on the y-axis. A perfect model would always (never) make true (false) forecasts and thus its ROC curve would pass through the point (0,1). Generally, the closer the ROC curve to the top-left corner, the better the discriminatory power of the model. The ROC curve provides a graphical way to visualize discriminatory power. A quantitative assessment of the discriminatory power is given by the Area under ROC curve (AUROC) which is calculated as follows 5 : where I(x) is an indicator function, defined to be 1 if x is true and 0 otherwise, s i B and s j H denote the bankruptcy scores of a model for the i-th bankrupt firm and for the j-th healthy firm observation respectively. Finally, n is the number of bankrupt firms and m is the number of healthy firm observations. Note that Eq. (1) is discontinuous and non-differentiable.

Maximizing discriminatory power
In this section, we present a methodology to maximize the discriminatory power (AUROC) when the bankruptcy score, s, is a probability, meaning that the model has a probabilistic response function 6 which is the case of popular bankruptcy forecasting models such as logistic and neural network models. Ideally, we should have used Eq. (1) directly as the objective function in the optimization. However, traditional gradient-based optimization methods cannot be used to maximize Eq. (1) directly because it is discontinuous and non-differentiable. For this reason, we introduce a surrogate function that seeks to maximize the discriminatory power. We define: as the difference between the probability of bankruptcy for the i-th bankrupt firm, p i B X i , and the probability of bankruptcy for the j-th healthy firm observation, p j H X j , , conditional on the predictor variables in X which could be a set of financial and market variables. From Eq. (1), to obtain the coefficients, β, that maximize the discriminatory power of a model we would like as many as possible d i,j 's to be positive because AUROC increases in this way. A way to achieve this is through the minimization of the following surrogate merit function: where 0 ≤ γ ≤ 1. The above merit function ignores the terms where d i,j ( )> (meaning that the difference in bankruptcy probabilities between the i-th bankrupt firm and j-th healthy firm observation is relatively high, as specified by the parameter γ) and penalizes the terms where d i,j ( )≤ . In other words, the parameter can be considered as a parameter which controls the magnitude of the d i,j 's that are to be penalized. For instance, if γ = 0, we penalize only the negative d i,j 's (i.e. only the cases where the model assigned a higher probability of bankruptcy for a healthy firm than a bankrupt firm) while if γ = 1, we penalize all d i,j 's.
Based on the optimality conditions of minimizing F(β), at the optimal solution, a number of d i,j 's must satisfy the condition d i,j = γ. 7 Hence, by selecting γ (close) to zero, we force a number of d i,j 's to be close to zero in absolute terms. In that case, a small change of the input data can easily induce d i,j 's to change signs which in turn will cause a change in the AUROC. This may be particularly evident in the case of out-of-sample data. That is, by training a model to produce d i,j 's close to zero, may yield a model with poor generalization ability and consequently the out-of-sample AUROC will be very sensitive. On the other hand, selecting γ (close) to one, coefficient estimates can blow up and provide unreasonable results. Thus, theoretically, the parameter value should be in between 0 and 1 (we explain later in this section how we compute the parameter empirically).
However, the surrogate function in Eq. (3) is non-differentiable when z = − d i,j ( ) = 0 . To overcome this problem and thus being able to use traditional gradient-based optimization algorithms, we replace the term max (0, z) with a differentiable function. Note that, we can minimize F(β) given by Eq. (3) using linear programming provided that the response function is linear with respect to the coefficients, β. Here, the probability is a non-linear function and as such we should use non-linear optimization algorithms to obtain the coefficients. We replace the term max (0, z) by the following ε-smoothed differentiable approximation, h (z): where ε is a small positive number close to zero. Here we set ε = 0.001. The ε-smoothed function h (z) , which we graphically present in Fig. 2, is a shifted version of the smoothed function used previously by Charalambous et al. (2007) to value call options. As can be The function max(0,z) is a surrogate function aiming to maximize the AUROC. However, this function is non-differentiable when z = 0. Thus, we replace it with the differentiable ε-smoothed function, h (z)  Fig. 3 This figure summarizes the work in our study and specifically how we train the models. The input vector, x, which can be financial and market data, enter the bankruptcy model (logistic or neural network model). The output of the model is the probability of bankruptcy, p(β), which depends on the coefficients imposed by the model and enters the merit function along with the target, t. The merit function can be the log-likelihood function (LL) or the AUROC, which is the one we propose in this study to use as the optimization criterion to obtain model coefficients. At each iteration, the optimization algorithm updates the coefficients until the merit function is optimized. For training we use data from the period 1990-2006 seen from the graph, the h (z) function has similar properties with the max (0, z) except that h (z) is differentiable when z = 0. Hence, the merit function to be minimized is replaced by: The next step is to estimate the coefficients, β, by training the model to minimize F(β) given by Eq. (5). Figure 3 summarizes the work in our study.
Consider that we have N training input samples (i.e. observations). Each input sample, x n = x 1n , x 2n , … , x kn , is associated with a known target,t n , where n = 1,2,…, N and k is the number of variables. In the context of bankruptcy forecasting, the input sample x n can be information characterizing the n-th firm, such as financial and market information, whereas t n is an indicator variable which equals 1 if the corresponding firm-observation goes bankrupt and 0 otherwise. The inputs enter the bankruptcy model (logistic or neural network model) to produce the probability of bankruptcy which is a function of the coefficients imposed by the model. The output of the bankruptcy model, p(β), with the associated target, t, are used in the merit function. Traditionally, the log-likelihood function is used to obtain the coefficients. In this study, we propose another way to obtain coefficients and specifically we use the merit function given by Eq. (5) which is optimized in order to obtain the coefficients of the bankruptcy models and consequently the probability of bankruptcy. Note that the target, t, is indirectly used in the merit function given by Eq. (5) in order to identify the bankrupt and healthy firms and to estimate their probability of bankruptcy. In this paper, the training sample spans the period 1990-2006. To solve the problem, we formulate a nonlinear unconstrained optimization process using MATLAB. Specifically, we use the fminunc command and the trust-region optimization algorithm to obtain the coefficients of the logistic and neural network models. At each iteration, the optimization algorithm updates the coefficients and the probability of bankruptcy (as shown in Fig. 3) until the merit function we propose is optimized.
As far as the parameter γ is concerned, we compute it empirically based on validation-a straightforward and easy to implement approach, which makes use only the training data to determine the parameter γ while the testing data remain intact. Also, validation is a frequently used method implemented by many studies to determine parameters underlying the models. We further divide our training sample into training (70%) and validation (30%) sets. 8 We train the models by choosing from the set of parameter values γ = {0, 0.1, 0.2, …, 1} and keep the value that gives the highest AUROC on the validation set. For instance, using our private and public firm models we find that γ equals 0.3 and 0.1 respectively, consistent with our conjecture that the γ parameter should be between 0 and 1. Then we merge the training and validation sets, to train the models as explained before and test their performance on the testing set 2007-2015.
We further illustrate the role of γ by providing an example using our data to provide an idea of how our method works and why it increases AUROC. First, we estimate the coefficients of a logistic model by maximizing the log-likelihood function and we calculate the d i,j 's . Second, we estimate the coefficients of a logistic model by minimizing given by Eq. (5) and we calculate the d i,j 's . Figure 4 shows a sample of those d i,j 's , produced by logistic regression i.e. by the model trained to maximize the log-likelihood function (top plot) and by maximizing AUROC with the ε-smoothed function, setting γ = 0 (middle plot) and γ = 0.3 (bottom plot). Recall that we would like as many as possible of d i,j 's to be greater than zero. Hence, they should lie above the solid straight line. For the logistic regression, some lie above and some below. Using the ε-smoothed function, we want to make as many as possible negative d i,j 's to move above the straight line. Setting γ = 0, we observe that all d i,j 's are close to zero. Some cases, 21 in particular, that were negative according to the logistic regression became positive (denoted with green crosses) and one case that was positive became negative (denoted with a red star), highlighting the limitation of producing d i,j 's that are close to zero. Setting γ = 0.3, not only more d i,j 's that were negative according to logistic regression became positive (59 in particular), but now the majority lie well above the solid straight line, several also passing the γ parameter which are the points that lie above the dashed line. Notice now that none of the d i,j 's that were positive according to logistic regression became negative because the higher value of γ, causes d i,j 's to be well above zero and as a consequence, AUROC will not be sensitive. Finally, as a benchmark, we obtain the coefficients of the logistic and neural networks models by maximizing the log-likelihood function, LL. Assuming that we have N training samples, LL is defined as follows: where p n x n , is the bankruptcy probability of the n-th observation, given the input vector of variables,x n and the coefficients, .

Information content tests
We further consider information content tests, also employed by related studies (see for instance Hillegeist et al. 2004;Agarwal et al. 2008;Charitou et al. 2013;Bauer et al. 2014).
In such tests the out-of-sample bankruptcy probabilities produced by various models, such as by models with maximized AUROC, enter as inputs to logistic regression models and we are interested to assess their explanatory power. In particular, we estimate the following panel logit specification: where p i,t is the probability of bankruptcy at time t, that the i-th firm will go bankrupt the next year and Y i, t+1 is the status of the i-th firm the next year (1 if it goes bankrupt and 0 if it is solvent). The variable of interest is prob i,t , which is the out-of-sample bankruptcy probability of the i-th firm at time t, produced by a model, for instance with maximized AUROC. Finally, is the coefficient estimate and a t is the baseline hazard rate that is only time-dependent and it is common to all firms at time t. Similar with prior studies, we proxy the baseline hazard rate with the actual bankruptcy rate at time t.
The specification in Eq. (7) is equivalent with the hazard model specifications used in related bankruptcy studies, such as Hillegeist et al. (2004); Agarwal et al. (2008);Bauer et al. (2014) etc. Specifically, Shumway (2001) argues that a panel logit model, like the one in Eq. (7), is equivalent with a hazard rate model and therefore standard log-likelihood maximization procedures can be used to estimate the logit model in Eq. (7), with a minor adjustment that we explain below.
The model in Eq. (7) represents a multi-period logit model as it includes observations for each firm across time. However, the inclusion of multiple firm-year observations per firm yields understated standard errors because the log-likelihood objective function, which is maximized to estimate the multi-period logit model, assumes that each observation is independent from each other. This is a wrong assumption since firm observations at time t + 1 cannot be independent from firm observation at time t. Failing to address this econometric issue, could lead to wrong inference regarding the significance of the individual coefficients. Similar with Filipe et al. (2016), we use clustered-robust standard errors to adjust for the number of firms in the sample but also for heteroskedasticity (Huber 1967 andWhite 1980).
t n ln p n x n , + 1 − t n ln 1 − p n x n , = e a * Rate t + * prob i,t 1 + e a * Rate t + * prob i,t

Economic analysis of bankruptcy models
The analysis so far addressed the forecasting accuracy of the bankruptcy models. But how accuracy is economically beneficial for banks? In particular, Bauer et al. (2014) show that even small differences in the AUROCs between the models affect the profitability of a bank. Similar findings are found in Charalambous et al. (2020). Therefore, it would be interesting to investigate the effect of using models with maximized AUROC, on bank economic performance. Here, we follow the approach of Agarwal et al. (2008) and Bauer et al. (2014) to examine it by assuming a loan market worth $100 billion and banks compete to grant loans to individual firms. Each bank uses a bankruptcy model to evaluate the credit worthiness of their customers.

Calculating credit spreads
We estimate the models using data spanning the years 1990-2006 (70% of the sample). We sort firm-customers from this sample in 10 groups of equal size and a credit spread is calculated according to the following rule; Firms in the first group, which are firms with the lowest bankruptcy risk, are given a credit spread, k and firms in the remaining groups are given a credit spread, CS i , obtained from Blochlinger et al. (2006) and it is defined as follows: where p(Y = 1|S = i) and p(Y = 0|S = i) is the average probability of bankruptcy and nonbankruptcy respectively, for the i-th group, with i = 2, 3, …,10 and LGD is the loan loss upon default. Following Agarwal et al. (2008), the average probability of bankruptcy for the i-th group is the actual bankruptcy rate for that group, defined as the number of firms that went bankrupt the following year divided by the number of firms in the group. Furthermore, k = 0.3% and LGD = 45%.

Granting loans and measuring economic performance
To evaluate economic performance, we assume that banks compete to grant loans to prospective firm-customers between the period 2007-2015. Each bank uses a bankruptcy model that has been estimated in the period 1990-2006. The bank sorts those customers according to their riskiness and rejects the bottom 5% with highest risk. The remaining firms are classified in 10 groups of equal size and firms from each group are charged a credit spread that has been obtained from the period 1990-2006. Finally, the bank that charges the lowest credit spread for the customer (i.e. for the firm-year observation) is granting the loan. Two measures of profitability are used. The first one, Return on Assets (ROA) is defined as Profits/Assets lent and the second one, Return on Risk-Weighted Assets (RORWA) takes into consideration the riskiness of the assets, defined as Profits/ Risk-Weighted Assets. Risk-Weighted Assets are obtained from formulas provided by the Basel Committee on Banking Supervision (2006).

Results
In this section, we present the out-of-sample comparisons between models with maximized AUROC and models with maximized log-likelihood (traditional models). We use the bankruptcy years 1990-2006 for training and keep the bankruptcy years 2007-2015 as the testing set. We start our analysis by comparing their performance in terms of discriminatory power, information content and economic benefits, when forecasting bankruptcy one year ahead. 9 Next, using the same tests, we compare their performance by forecasting bankruptcy two years ahead and then when our sample consists of financially distressed firms. An additional analysis is performed in this section where we compare our methodology with other methods proposed in the literature to maximize AUROC using the same analysis as before. Finally, we provide out-of-sample discriminatory power comparisons, when we use quarterly data. 1 3 Table 4 shows the out-of-sample performance (2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014)(2015) of the models with maximized AUROC versus the models with maximized log-likelihood, in terms of discriminatory power.

AUROC results
Overall, models that are trained to maximize AUROC perform better out-of-sample compared to models trained to maximize the log-likelihood function, indicating that the function we introduced performs well out-of-sample in discriminating firms that will go bankrupt the next year. The effect by maximizing the AUROC, as expected, is more pronounced in the case of "private firms model" where only limited information is available (i.e. financial information), hence there is more space to improve the performance. For instance, the AUROC of the logistic and neural network model trained to maximize the log-likelihood function (LL) are 0.8991 and 0.9138 respectively where those trained to maximize AUROC are 0.9221 and 0.9332 respectively. DeLong tests indicate that AUROC differences are statistically significant at the 5% level (2.43 and 2.04 respectively). In contrast, the effect is less pronounced in the case of "public firms model", since the inclusion of market data in addition to financial data, further increases the forecasting power of the models. In fact, the AUROCs of models trained to maximize the log-likelihood function are quite high and specifically, 0.9425 for the logistic 10 model and 0.9440 for the neural network. This improves to 0.9470 and 0.9508 respectively when maximizing AUROC. Differences are statistically significant at the 10% level.
From the results in this section, we suggest using the neural network model trained to maximize AUROC since it is the best-performing model which is consistent with the notion that neural networks outperform simpler modeling approaches (Zhang et al. 1999;Kumar et al. 2007;Lessmann et al. 2015). For the user interested in simpler models, we suggest the implementation of the logistic model but trained to maximize AUROC.

Information content results
In this section we report the results from information contest tests. We compare the information contained in out-of-sample bankruptcy probabilities produced by models where the AUROC is maximized versus where the log-likelihood is maximized. Models 1 and 2 include the out-of-sample (2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014)(2015) bankruptcy probabilities produced by a neural network (Prob1) and a logistic model (Prob2) respectively, obtained by maximizing the AUROC. Models 3 and 4 include the bankruptcy probabilities produced by a neural network (Prob3) and a logistic model (Prob4) respectively, obtained by maximizing the loglikelihood function. Table 5 reports the results of logit models that include the out-of-sample bankruptcy probabilities as explanatory variables but also the annual bankruptcy rate (Rate) as the baseline hazard rate.
Panel A reports results from four logit regression models. Models 1-4 in the first four columns refer to the models where the corresponding bankruptcy probability which is included as predictor (Prob1-Prob4) is generated with financial data only. We re-estimate the four logit regressions, of which their results are presented in the next four columns of Table 4 but this time, the corresponding bankruptcy probability is generated with financial and market data.
According to the results, the bankruptcy probabilities in all cases are highly statistically significant, indicating that they carry significant information in predicting bankruptcy one year ahead, (coefficient estimates are significant at the 1% significance Table 5 Information content test results-Out of sample (2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014)(2015) This table reports results from information content tests. Panel A shows estimation of four logit models. Models 1 and 2 include out-of-sample (2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014)(2015) bankruptcy probabilities produced by a neural network and a logistic model respectively, whose coefficients are obtained by maximizing AUROC, with financial data as inputs. Models 3 and 4 include bankruptcy probabilities produced by a neural network and a logistic model respectively, whose coefficients are obtained by maximizing the log-likelihood function, with financial data as inputs (private firm models). We re-estimate the four logit models 1-4 as explained before and report results in the next four columns but this time, the out-of-sample bankruptcy probabilities are generated with financial and market data as inputs to the neural network and logistic models (public firm models) All logit regression models include the Rate as proxy for the baseline hazard rate, which is the prior year bankruptcy rate in our sample. The last two rows of the panel reports log-likelihood and pseudo-R 2 for each model. level). More importantly, bankruptcy probabilities produced by models with maximized AUROC (Prob1 and Prob2) contain significantly more information than bankruptcy probabilities produced by models with the log-likelihood being maximized (Prob3 and Prob4). This is especially evident by the substantially higher pseudo-R 2 of models 1 and 2 (29.05% and 26.36% respectively) compared to models 3 and 4 (18.84% and 8.56% respectively) for the private firms case. Similarly, pseudo-R 2 of models 1 and 2 (34.53% and 33.77% respectively) is substantially higher than models 3 and 4 (19.74% and 14.07% respectively) in the case of public firms models.
In panel B, we use the Vuong (1989) test-statistic to test for differences in model log-likelihood values between various (non-nested) models. Results show that, in the case of private firms, the log-likelihoods of models 1 and 2 are significantly different than models 3 and 4 (test-statistics are 5.75 and 7.37 respectively) whereas Vuong teststatistics are 4.37 and 8.07 respectively for public firm models.
Overall, our results suggest that models with maximized AUROC provide probability estimates that contain significantly more information about bankruptcies over the next year compared to models which are trained to maximize traditional functions, such as the loglikelihood function, even when the increase in AUROC is relatively small (as in the case Table 6 Economic performance results This table reports economic results for four banks in a competitive loan market worth $100 billion. Bank 1 uses a neural network model with maximized AUROC for estimating the bankruptcy score of its customers. Bank 2 uses a logistic model with maximized AUROC. Banks 3 and 4 are competitors to banks 1 and 2 respectively, using a neural network and a logistic model respectively, but by maximizing the log-likelihood function. The table reports results in two cases; When the inputs to the models include financial data (private firm models) and when the inputs include financial and market data (public firm models) The banks sort prospective customers (2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014)(2015) and reject the 5% of firms with the highest risk. The remaining firms are classified in 10 groups of equal size and for each group, a credit spread is calculated, as described in the main text. The bank that classifies the firm to the group with the lowest spread is finally granting the loan. Market share is the number of loans given divided by the number of firmyears, Revenues = (market size)*(market share)*(average spread), Loss = (market size)*(prior probability of bankruptcy)*(share of bankruptcies)*(loss given default). Profit = Revenues-Loss. Return on Assets is profits divided by market size*market share and Return on Risk-Weighted-Assets is profits divided by Risk-Weighted Assets, obtained from formulas provided by the Basel Accord (2006 of our public firm models as shown in Table 4). From all models, however, our proposed neural network model is the best performing one.

Economic performance results
So far, we have considered discriminatory power and information contest tests to assess model performance. However, a bank is generally interested in the economic benefits arising by using bankruptcy forecasting models in the decision-making process of granting loans to individual firms. Following Agarwal et al. (2008) and Bauer et al. (2014), we consider a loan market worth $100 billion and four banks are competing to grant loans to prospective firm customers. We hypothesize that banks 1 and 2 are more sophisticated banks, training models by maximizing AUROC (a neural network and a logistic model respectively) whereas banks 3 and 4 are more "naïve", training models by maximizing the log-likelihood function (a neural network and a logistic model respectively). In Table 6 we report the results, for both private and public firm models. Clearly, banks 1 and 2 which use models with maximized AUROC, respectively, manage loan portfolios with higher quality relative to banks 3 and 4. This is evident by the lower concentration of bankruptcies they attract. In particular, the bankruptcy rate of bank's 1 portfolio, which uses a neural network trained to maximize AUROC is 0.047% and 0.074% when using the private and public firms model respectively. In contrast, the bankruptcy rate of bank's 3 portfolio, which uses a neural network trained to maximize the log-likelihood is 0.17% when using the private and public firms model respectively. Similarly, bank 2 which uses a logistic model trained to maximize AUROC manages a credit portfolio with bankruptcy rate equal to 0.38% and 0.066% when using the private and public firms model whereas for the bank using models to maximize the log-likelihood function (bank 4) the rates are 0.56% and 0.53% respectively.
Consequently, banks 1 and 2 achieve superior economic performance 11 compared to banks 3 and 4 respectively, on a risk-adjusted basis. For example, considering the private firms model, bank 1 which uses a neural network model with maximized AUROC, earns 2.24% relative to bank 3 which uses a traditional neural network model (1.52%). Also, bank 2 which uses a logistic model trained to maximize AUROC, earns 1.27% relative to bank 4 which uses a traditional logistic model (0.91%). In the case of public firms model, similar insights are obtained; bank 1 earns higher risk-adjusted returns relative to bank 3 (1.98% and 1.68% respectively). This is also the case for banks 2 and 4 (1.53% and 1.20% respectively). Again, the neural network trained to maximize AUROC, provides the higher economic benefits for the bank using it (bank 1).

Forecasting bankruptcy two years ahead
In this section, we evaluate and compare the performance of the models by increasing the forecasting horizon to two years. This is a more challenging problem because the characteristics of bankrupt firms are less pronounced relative to one year prior to bankruptcy and therefore is more difficult to forecast bankruptcy. Also, identifying the signs of the crisis Table 7 Forecasting bankruptcy two years ahead; Maximizing AUROC vs maximizing LL This table reports out-of-sample performance in terms of discriminatory power, information content and economic impact, of models trained to maximize AUROC versus models trained to maximize the log-likelihood function, when forecasting bankruptcy two years ahead. Panel A reports the performance when financial variables are used in the models (private firms models) and when financial and market variables are used in the models (public firms models). Panel B reports DeLong (1988) test statistics for differences in AUROCs between the models of interest whereas panel C reports Vuong (1989) test statistics for differences information content between the models of interest earlier, although difficult, can help the management of the firm to take remedial actions to correct the adverse situation and avoid bankruptcy. We perform the same analysis as before and we provide the results in Table 7. As expected, performance has dropped relative to when forecasting bankruptcy one year ahead due to the increased difficulty of the problem. However, models which are trained to maximize AUROC outperform those trained to maximize the log-likelihood function. Starting from the AUROC, which is the focus of this paper to improve, we document outof-sample that it is significantly higher, especially for the neural networks case. Two years prior to bankruptcy, the neural network with maximized AUROC achieves an AUROC equal to 0.8678 and 0.8864 for the private and public firms models respectively whereas for the neural network trained to maximize the log-likelihood function, these are 0.8441 and 0.8571 respectively. DeLong tests indicate that the differences are significant at the 5% and 1% level respectively (test statistics are 2.18 and 3.17 respectively). A logistic model with accounting data as predictors (private firms model) trained to maximize AUROC, has significantly higher discriminatory power than a competing logistic model trained to maximize the log-likelihood function (0.8503 vs 0.8113 respectively). Differences in AUROCs are statistically significant at the 1% level according to the DeLong test (test statistic is 3.54). For the public firms case, an enhancement is achieved (0.8664 vs 0.8558) but the difference is not statistically significant (DeLong test statistic is 1.28).
The summary of the remaining tests is that the neural network with maximized AUROC performs significantly better in terms of information content and provides higher economic benefits relative to a competing neural network which is trained by maximizing the loglikelihood function. The former is also the best performing model in all tests. Finally, a logistic model with maximized AUROC provides significantly more information than a logistic with maximized log-likelihood function, albeit no economic benefits are achieved in this case.

Forecasting financial distress
In this section, we change the event and instead of forecasting bankruptcy, we forecast financial distress. There are several reasons as to why firm stakeholders should be interested in models forecasting financial distress more accurately. First, financial distress is a state prior to bankruptcy filing and therefore forecasting the early signs of the crisis may help, for instance the management of the firm, to take corrective measures in order to avoid further deterioration that may ultimately lead to bankruptcy in which case the firm loses most of its value (see for instance Asquith et al. 1994;Glover 2016). Second, forecasting the early signs of the crisis is more challenging to accomplish not only because it is the starting point of the crisis but also because financial distress is not a formal event like bankruptcy, thus we need to construct a financial distress indicator. In this study, we follow Keasey et al. (2015) and Gupta et al. (2018) and we consider a firm as financially distressed if the following conditions are satisfied; 1) Earnings Before Interest, Tax and Depreciation and Amortization (EBITDA) is less than financial expenses (i.e. interest payments) for two consecutive years 2) Total Debt is higher than the Net Worth of the firm for two consecutive years and 3) The firm experiences negative Net Worth growth between two consecutive years. The firm is classified as financially distressed in the year immediately following these three events. For forecasting purposes, we use the data two years before financial distress. For example, when the conditions are satisfied for the years t and t-1, then the firm is considered as financially distressed in the year t and we construct the variables at t-2 to Table 8 Forecasting financial distress; Maximizing AUROC vs maximizing LL This table reports out-of-sample performance in terms of discriminatory power, information content and economic impact, of models trained to maximize AUROC versus models trained to maximize the log-likelihood function, when forecasting financial distress. Panel A reports the performance when financial variables are used in the models (private firms models) and when financial and market variables are used in the models (public firms models). Panel B reports DeLong (1988) test statistics for differences in AUROCs between the models of interest whereas panel C reports Vuong (1989) test statistics for differences information content between the models of interest predict financial distress. Following these conditions, we generate an extensive database with 1,929 financially distressed firms. In Table 8 we report the out-of-sample results from this exercise. Starting from the private firms model, logistic and neural network models trained to maximize AUROC, significantly outperform models trained to maximize the log-likelihood function (0.9175 vs 0.8956 for neural networks and 0.8982 vs 0.8897 for logistic models). Differences in AUROCs are statistically significant at the 1% level (DeLong test statistics are 5.42 and 3.29 for neural networks and logistic models respectively). Similar results are found with respect to the public firms model (0.9000 vs 0.8870 for neural networks and 0.8824 vs 0.8753 for logistic models). Differences in AUROCs are statistically significant at the 1% and 5% level respectively (DeLong test statistics are 5.99 and 2.24 for neural network and logistic models respectively).
Regarding the remaining tests, we find that the both models which are trained to maximize AUROC provide significantly more information and there is more gain by banks using them relative to using models trained to maximize the log-likelihood function. Overall, from this test we conclude that our methodology can help the interested parties to improve bankruptcy forecasts, considering the harder nature of the problem, either using a neural network or a logistic model. Once again, the neural network constructed to maximize AUROC is the best performing one.

Comparing our methodology with other methodologies
In this section we use the same analysis as before to compare our proposed methodology with other methods of AUROC maximization proposed in the bankruptcy literature and the advantages (shortcomings) of our method (other methods) are discussed.
We consider two other approaches proposed by Miura et al. (2010) and Kraus et al. (2014), to maximize AUROC of credit scoring models. Miura et al. (2010) suggest a sigmoid function as an approximation of Eq. (1). Specifically, they maximize the following objective function 12 : However, unlike the function that we introduced previously, it treats all d i,j 's in the same way, whereas our function, give more emphasis on the "bad" cases, for example, when a healthy firm has higher bankruptcy score than a bankrupt firm. Further, the authors consider only a linear response function (the output is a linear score) and unlike our method, it cannot be used by models which employ probabilistic response functions such as logistic models and highly non-linear models such as neural networks. Instead, our methodology works well with logistic and neural networks, which are the among the most popular bankruptcy models, because they allow for a probabilistic response function. Finally, Kraus et al. (2014) suggest using directly Eq. (1) as the objective function and implementing derivative-free methods (such as Nelder et al. 1965) to optimize the coefficients. The optimization algorithm that is used, however, assumes that ∕  Table 9 Forecasting bankruptcy one year ahead; Our methodology vs alternative methodologies This table reports out-of-sample performance in terms of discriminatory power, information content and economic impact, using our methodology to maximize AUROC on a neural network versus alternative methods to maximize AUROC as proposed by Miura et al. (2010) and Kraus et al. (2014), denoted as KK (2014), when forecasting bankruptcy one year ahead. Panel A reports the performance when financial variables are used in the models (private firms models) and when financial and market variables are used in the models (public firms models). Panel B reports DeLong (1988) test statistics for differences in AUROCs between the models of interest whereas panel C reports Vuong (1989)

Table 10
Forecasting bankruptcy two years ahead; Our methodology vs alternative methodologies This table reports out-of-sample performance in terms of discriminatory power, information content and economic impact, using our method to maximize AUROC on a neural network versus alternative methods to maximize AUROC as proposed by Miura et al. (2010) and Kraus and Kuchenhoff (2014), denoted as KK (2014), when forecasting bankruptcy two years ahead. Panel A reports the performance when financial variables are used in the models (private firms models) and when financial and market variables are used in the models (public firms models). Panel B reports DeLong (1988) test statistics for differences in AUROCs between the models of interest whereas panel C reports Vuong (1989)

Table 11
Forecasting financial distress; Our methodology vs alternative methodologies This table reports out-of-sample performance in terms of discriminatory power, information content and economic impact, using our method to maximize AUROC on a neural network versus alternative methods to maximize AUROC as proposed by Miura et al. (2010) and Kraus and Kuchenhoff (2014), denoted as KK (2014), when forecasting financial distress. Panel A reports the performance when financial variables are used in the models (private firms models) and when financial and market variables are used in the models (public firms models). Panel B reports DeLong (1988) test statistics for differences in AUROCs between the models of interest whereas panel C reports Vuong (1989) test statistics for differences information content between the models of interest the objective function is continuous, which is not the case for Eq. (1). Also, this approach while is easy to implement, ignores information provided by the gradient which could increase the accuracy of the coefficients after the optimization process and thus we believe that using specifications with differentiable functions is a better choice. 13 Table 9 presents the out-of-sample results by forecasting bankruptcy one year ahead (we use our neural network model which consistently outperformed the competing models).
Overall, we find that the neural network trained with our method has higher discriminatory power than the other models but the differences is AUROCs are not statistically significant. Despite that, the neural network provides significantly more information relative to Miura et al. (2010) and KK (2014) which is also economically beneficial.
When we increase the complexity of the problem, however, by forecasting bankruptcy two years ahead and forecasting financial distress, our neural network model consistently outperforms the competing methods and differences in performance are statistically significant. Results regarding forecasting bankruptcy two years ahead, reported in Table 10, show that the neural network model more accurately discriminates bankrupt from healthy firms Table 12 Forecasting bankruptcy and financial distress using quarterly data This table reports AUROC results for logistic and neural network models, when their coefficients are estimated by maximizing the AUROC (reported in the second and third columns respectively) and when their coefficients are estimated by maximizing the log-likelihood function (denoted as LL and reported in the fourth and fifth columns respectively) The models are trained in the period 1990-2006 and the table reports the prediction performance, in terms of AUROC, in the out-of-sample period when predicting bankruptcy one quarter ahead (panel A), four quarters ahead (panel B), eight quarters ahead (panel C) but also financial distress (panel D). In the models, we incorporate either financial variables (Private Firms Model) or financial and market data (Public Firms Model) constructed using quarterly data . Differences are statistically significant at the 10% level using the private firms model and at the 1% level using the public firms model. In the remaining tests, we document that the neural network provides significantly more information about future bankruptcies than the Miura et al. (2010) and KK (2014) methods and the better performance is associated with higher economic benefits for the bank which uses our proposed method.
Finally, in Table 11 we report our results when forecasting financial distress. Consistent with previous results, our method achieves significantly higher discriminatory power than the other methods (0.9175 vs 0.8982 for the Miura et al. 2010 and0.8964 for KK 2014 in the private firm model case and 0.9000 vs 0.8825 and 0.8827 in the public firm model case). Differences are statistically significant at the 1% level.
In the remaining tests we confirm previous findings; our method provides significantly more information and the better performance overall is economically beneficial for the bank using our method as opposed to the competition.

Forecasting using quarterly data
Public firms issue financial information each quarter, thus investors can update their risk assessments more frequently as new information becomes available. In this section, we perform the same analysis as before but this time, the input variables to the models are constructed using quarterly data and we make predictions one, four and eight quarters ahead. Overall, the AUROC results reported in Table 12, are qualitatively similar with the results reported in the case where yearly data are used. More specifically, training the bankruptcy models with our proposed method, in all cases, improves the out-ofsample performance in terms of discriminatory power as opposed to training the models to maximize the traditional log-likelihood function.

Conclusions
The goal of this paper is to propose an alternative method to estimate the coefficients of bankruptcy forecasting models and specifically logistic and neural network models which are the most popular bankruptcy models used in prior research. In particular, we suggest those interested in forecasting bankruptcy, to obtain the coefficients by maximizing the discriminatory power as measured by the Area Under ROC curve (AUROC). In this study, a method is introduced and we highlight the benefits arising, out-of-sample, by using models which are trained to maximize AUROC over models trained with traditional methods, such as optimizing the log-likelihood function. Overall, we find that models trained to maximize AUROC outperform traditional methods, out-of-sample, in terms of discriminatory power, information content and economic impact. Our results hold when we test the method in different settings, such as forecasting bankruptcy one year ahead which is the most common horizon, forecasting bankruptcy two years ahead and forecasting financial distress which make forecasts more difficult (using yearly and quarterly data). Thus, forecasting bankruptcy accurately well in advance, which would be beneficial for the firm to take corrective measures, requires a more sophisticated estimation method, such as maximizing the AUROC function by using our method. From all models, the neural network trained with our method is the best performing one.
Next, we compare our method with alternative methods proposed in the literature and we provide both theoretical as well as empirical justifications as to why our method should be preferred. As expected, the results are more pronounced when we increase the forecasting difficulty, such as forecasting financial distress.
Our results have implications to the way bankruptcy forecasting is performed. Our proposed estimation approach provides, to those interested to forecast bankruptcy, a significant advancement over traditional methods that can be used by logistic and neural networks for better bankruptcy analysis and possibly can be extended to areas such as in credit risk analysis.

Conflict of interest he authors declare no conflicts of interest.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.