Predicting Firms’ Financial Distress: An Empirical Analysis Using the F-Score Model

Financial performance of firms is very important to bankers, shareholders, potential investors, and creditors. The inability of firms to meet their liabilities will affect all its stakeholders and will result in negative consequences in the wider economy. The objective of the study is to explore the applicability of a distress prediction model which uses the F-Score and its components to identify firms which are at high risk of going into default. The study incorporates a prediction model and vast literature to address the research questions. The sample of the study is collected from publicly listed firms of the United States. In total, 81 financially distressed firms wereextracted from the UCLA-LoPucki Bankruptcy Research Database during 2009–2017. This study found that the relationship of the F-Score and probability of firms going into financial distress is significant. This study also demonstrated that firms which are at risk of distress tend to record a negative cash flow from operations (CFO) and showed a greater decline in return on assets (ROA) in the year prior to default. This study extends the existing literature by supporting a model which has not been widely used in the area of financial distress predictions.


Introduction
When a firm goes into bankruptcy or when it defaults on its financial obligations, it will have negative impacts on its stakeholders such as lenders, creditors, customers, employees, and shareholders. Lenders and creditors will lose their capital loaned. The operations of customers' businesses might be halted due to sudden interruptions in the supply of raw materials, which will adversely affect revenues, margins, and reputations. Employees will lose income, whereas shareholders might suffer losses from capital invested. In addition, these stakeholders are often debtors themselves and this affects their ability to meet financial obligations. Therefore, the bankruptcy of one firm not only affects its direct stakeholders, it will also have spillover effects on the economy as a whole, and its severity is dependent on the size of the company. Based on the above, it is evident that business failures will result in serious consequences to many parties. As such, in the past few decades, the area of bankruptcy and financial distress prediction has attracted the interests of many researchers. In the United States (US), businesses can declare bankruptcy by filing for Chapter 7 or Chapter 11 with the United States Courts. Chapter 7 is also known as liquidation bankruptcy, whereas Chapter 11 is known as reorganization bankruptcy (United States Courts 2017a, 2017b. Chapter 7 filings are for companies where reorganizations are not feasible and all of their assets will be liquidated to pay off debts. In contrast, Chapter 11 filings allow companies to restructure debts and revive their businesses. After a few years of decline in the number of bankruptcy filings in the United States, the number of commercial bankruptcy filings increased again in 2016 (Brown 2017). The to- J. Risk Financial Manag. 2021, 14, 199 2 of 16 tal assets involved in Chapter 7 and Chapter 11 filings in 2016 were in excess of $100 billion (Brown 2017). Based on the quantum of assets involved, it is evident that the cost associated with bankruptcies is high. These costs not only have to be borne by the distressed companies; they will also result in negative consequences to their creditors due to the inability of the distressed companies to meet their obligations. Furthermore, as the severity of economic crises is increasing, it is expected that the losses associated with bankruptcies will also be much higher. Therefore, it is very important for businesses to take proactive actions to minimize such losses. In order to do that, managers need to understand the factors which will lead to financial distress and have the tools available to enable them to make timely decisions. To date, many studies have been conducted to address the issue of financial distress predictions. Since the 1960s, many credit risk prediction models have been introduced. However, the performance of existing models is inconsistent (Tinoco and Wilson 2013) and thus, there is a need to continue to search for a better model that fits the current and future market environment. This study has primarily referenced the study conducted by Agrawal (2015), which demonstrated that the F-Score can be used to predict financial distress. The study tested the usability of the aggregate F-Score (Model 1) and the discrete variables that made up the F-Score (Model 2). The study by Agrawal (2015) was done using a sample of 270 Indian firms, which comprised of 135 defaulting firms and 135 non-defaulting firms.
However, the Indian market alone is small in size as compared to other markets such as the US market, and evidence from the Indian market is insufficient to be used as a representation of the global market. In addition, the operating environment for companies in India is also different from the US market. In India, access to capital is limited, which hinders the ability of firms to raise funds in times of necessity and in turn affects the survival of businesses. According to Sivabalan (2017), the market capitalization of US is close to $30.0 trillion, whereas the market capitalization of the Indian market is much lower at $2.3 trillion. With global market capitalization at close to $100 trillion (Edwards 2017), the US market represents about one third of the world's equity market. The bankruptcies of any large public companies in the US will impact many stakeholders, including both domestic and international stakeholders. In addition, US is also a major import and export partner to many countries around the world and the negative consequences of bankruptcies in the US will spillover to other countries, creating a domino effect in the world economy. Thus, given its size and its significance to the global economy, this study has been conducted using data from companies listed in the US market.
The global economy has evolved over the years and the current business environment is very different as compared to the 1980s.Based on this background, the drivers of corporate failure might have changed as well and thus, there is a need to keep on exploring and find the best default prediction models by incorporating the most relevant variables for the current economic environment. The F-Score developed by (Piotroski 2000) is one potential model.
In addition to the commonly used ratios in existing bankruptcy prediction models, the F-Score also takes into account the quality of earnings of firms, which is represented by the accruals variable in the model. This is an important differentiating factor, as research by Janes (2005) showed that the addition of accruals can significantly improve the results of distress predictions. The study by Agrawal (2015) also provided significant evidence that the F-Score has default predictive ability for firms in India. Piotroski (2000) first developed the F-Score as a measure to select good stocks for investments. The F-Score ranges between zero to nine and is the summation of scores based on nine financial variables. Piotroski (2000) demonstrated that companies with a high F-Score give better stock returns as compared to those with a low F-Score. As such, this implies that companies with a low F-Score will give low stock returns. According to Campbell et al. (2011), companies with a high probability of failure give low average stock returns as well.Given this similarity, we hypothesize that firms with a low F-Score have a higher tendency to be financially distressed as compared to the firms with a high F-Score. This study contributes by providing evidence that the relationship of the F-Score and probability of firms going into financial distress is significant. This study also demonstrates that firms which are at risk of distress tend to record a negative cash flow from operations (CFO) and show a greater decline in return on assets (ROA) in the year prior to default. Finally, this study contributes to existing literature by examining the usefulness of the Piotroski's F-Score and its components in predicting financial distress among firms in the United States.

A Review of Literature on Development of Default Prediction Models
The subject of bankruptcies prediction is a well-researched topic among finance researchers. Some of the earliest studies such as those done by Beaver (1966) and  found that analysis of financial ratios calculated from financial statements can be used to predict financial failure and the risk of bankruptcy.
The first notable bankruptcy prediction model was developed by Altman (1968) which is widely known as the Z-Score. In developing the Z-score, Altman (1968) used multivariate discriminant analysis and identified five variables which are most effective in predicting corporate bankruptcies. The five variables are the working capital/total assets ratio, retained earnings/total assets ratio, earnings before interest and taxes/total assets ratio, market value of equity/book value of total debt ratio, and sales/total assets ratio. The Z-score is made up of the weighted average combination of these five ratios. In the study by , 95% of all firms in the sample was correctly categorized into two groupings, i.e., bankrupt and non-bankrupt. Based on the findings in the study, the Z-Score's predictive ability decreases as the number of years before bankruptcy increases.
Subsequently, in response to the changes in accounting standards in the 1970s, Altman et al. (1977) further researched and modified the Z-score to develop a new model for bankruptcy prediction named the Zeta model. The Zeta model was also developed by using multivariate discriminant analysis, which was the same technique used by . The Zeta model consisted of seven variables, namely return on assets, stability of earnings, interest coverage, cumulative profitability, liquidity, capitalization, and size. It was found that, in the year prior to bankruptcies, the accuracies of both the Z-score and Zeta models are about the same at 95% and 96%, respectively. However, when financial data between two to five years prior to bankruptcies were used, the Zeta model is more accurate in predicting bankruptcies than the Z-score (Altman et al. 1977).
The O-score is another prediction model which was developed by Ohlson (1980) to predict corporate bankruptcies. The O-score has nine variables which are mainly derived from the balance sheet and income statement. Ohlson (1980) is one of the first researchers to use conditional logistic regression analysis in studying bankruptcy risks. The O-Score was derived from a study of more than 2000 companies (Ohlson 1980), whereas the Z-Score was developed based on a study of 66 companies . Until today, the Z-score, Zeta model, and O-score have remained as popular reference points for research in the area of financial distress and bankruptcy predictions. Even in recent years, researchers such as Le Maux and Morin (2011), Tinoco and Wilson (2013), and Agrawal (2015) have continued to show that the information contained in a company's balance sheet, profit and loss statements, as well as cash flow statements can be used to predict the risk of payment defaults.
In the past decade, due to rapid advancements in technologies, many researchers had developed bankruptcy prediction models using advanced techniques such as neural networks (Odom and Sharda 1990;Salchenberger et al. 1992;Coats and Fant 1993;Yang et al. 1999;Kim and Kang 2010), genetic algorithms (Varetto 1998;Shin and Lee 2002), and decision trees (Pompe and Feelders 1997). These models are newer and were developed using modern statistical computation techniques. However, Adnan Aziz and Dar (2006) performed a comparison of available prediction models and found that there is no marked difference, in terms of prediction accuracy, between traditional statistical models and more advanced models, as mentioned above. This is consistent with the argument by Bandyopadhyay (2006), that there are no unifying agreements between researchers on which models have highest predictive accuracies. Adnan Aziz and Dar (2006) also found that accounting-based statistical models such as multivariate discriminant analysis and logistic models are most reliable.
The Z-Score and O-Score were developed more than 30 years ago. The fact that today's businesses are still using these models in predicting financial distress in their day to day operations reflects the strong appeal of these prediction models to society, due to their relatively simplistic nature of their applications (Kumar and Kumar 2012). However, Kumar and Kumar (2012) revisited the Z-Score and O-Score by testing these models using the financial data of Texmo Industries for a five year period between 2005-2006 and 2009-2010. The study found that the O-Score performed slightly better than the Z-Score but the predictive abilities of both models are unsatisfactory.

Financial Ratios and Prediction of Financial Distress
According to Le Maux and Morin (2011), bankrupt companies usually exhibit similar financial attributes a few years prior to the event. Since the early 1960s, many studies such as those done by Altman (1968) 2018) had demonstrated that financial statements analysis can be used as a proactive approach in identifying and managing companies with a high tendency to go into financial distress. This is consistent with the findings by Jones (2017), which revealed that ratios calculated from financial statements are some of the strongest predictors of corporate bankruptcies, whereas external factors such as macroeconomic indicators, analyst recommendations, as well as industry variables are the weakest.
The financial ratios used in both the Z-Score and Zeta model are calculated from data in the balance sheet and profit and loss statements only and do not contain information from cash flows statements. In more recent years, more and more researchers argued that relying on data from balance sheet and profit and loss statements alone to predict financial distress is insufficient. John (1976) argued that the data in the balance sheet and profit and loss statements are susceptible to manipulation. According to Kordestani et al. (2011), companies are unable to make many distortions in the cash flow statements and thus, this increases the reliability of the data. Gentry et al. (1985) demonstrated that when cash flows variables were used in combination with accounting financial ratios, the accuracy of bankruptcies prediction is higher. Gentry et al. (1985) also found that dividend, investment, and receivables cash flows have significant importance in predicting bankruptcies. This is consistent with the study done by Almamy et al. (2016) which demonstrated that the predictive ability of the Z-Score increased when operating cash flow is added into the model.
In addition, Sloan (2005) demonstrated that improvements in earnings, which were driven by alteration of accruals, is not sustainable. Agrawal (2015) argued that other than the commonly used ratios, the quality of earnings is also a vital signal of firm's financial position. This is consistent with the results by Janes (2005), which demonstrated that the use of accruals as variable can significantly improve the results of distress prediction models.
One such model which incorporates the above-mentioned variables is the F-Score by Piotroski (2000). Agrawal (2015) had shown that this model is useful in the prediction of defaults in Indian firms. Kolte et al. (2017) conducted a study to assess the applicability of the F-Score in predicting the bankruptcy of Kingfishers Airline. The results from the study showed that Kingfishers Airline's F-Scores in the three years prior to its bankruptcy were low.Thus, this study contributes to existing literature by investigating the usefulness of the F-Score in predicting financial distress for public companies in United States.

A Review of Statistical Techniques in Financial Distress and Bankruptcy Prediction Models
A review of literature revealed that the commonly used statistical techniques in developing accounting-based bankruptcy prediction models are multiple discriminant analysis, logistic regression, the probabilistic model, and linear probability (Altman and Saunders 1998). Of these four approaches, multiple discriminant analysis is the most widely used technique, followed by logistic regressions (Altman and Saunders 1998;Bandyopadhyay 2006). The Z-Score and Zeta models were developed using multiple discriminant analysis, whereas the O-Score is a logistic regression model.
However, the usage of discriminant analysis requires the independent variables to have a multivariate normal distribution (Angelini et al. 2008). According to the study done by Lennox (1999), cash flow and leverage have non-linear effects on the probability of bankruptcy of a firm. Thus, for bankruptcy prediction models incorporating cash flow and leverage variables, it is more suitable to adopt logistic regression analysis instead of discriminant analysis, as the logistic model does not require the variables to satisfy the assumption on normality. Begley et al. (1996), Lennox (1999), and Mihalovic (2016) also showed that logistic and probabilistic models can predict bankruptcy more accurately than discriminant analysis. Begley et al. (1996), Lennox (1999), and Pongsatat et al. (2004) concluded that generally, the O-Score, which was developed using logistic analysis, has higher predictive accuracy as compared to the Z-Score. The strength of logistic models is further supported by Adnan Aziz and Dar (2006), who concluded that logistic regression models are one of the most reliable models in the area of bankruptcy prediction due to their stability and low Type 1 and Type 2 errors. Based on the above evidence, the most appropriate statistical technique for this study is the logistic regression model.

Sample Description and Statistical Technique
The sample of financially distressed firms for this study comprise of companies which had filed for Chapter 11 bankruptcy protection in the United States between January 2009 and September 2017. The list of companies was extracted from UCLA-LoPucki Bankruptcy Research Database and the original list consists of 306 large public companies.
It was deemed that companies in financial sector are unsuitable for this study, as the performance metrics for companies in these industries differ from other industries. On top of that, the variables in the F-Score do not take into account variables which are deemed important in measuring the performance of financial institutions such as the bank capital, loan ratio (Cleary and Hebb 2016), and net interest margin (Goyal and Bhatia 2016). At times of distressed, firms in the financial sector tend to receive government interventions in the form of bailouts or capital support as well (Cabrera et al. 2016;Dieckmann and Plank 2011;Kirshner 2015;Berger et al. 2016). Thus, companies in financial industries such as banks and insurance companies are excluded from this study. Subsequently, the financial data for remaining companies was extracted from Bloomberg and companies with incomplete financial information were excluded from this study. This includes i.
Companies whereby financial data is unavailable for the period between one to three years prior to bankruptcy. ii.
Financial data available is insufficient to calculate all the ratios included in the F-Score.
As a result, the sample of financially distressed firms ended up with a total of 81 firms. Table 1 shows the year-wise distribution of the sample for financially distressed firms. The financial data for the variables used in this study was mainly collected from Bloomberg and the EDGAR database on the website of the Securities and Exchange Commission of the United States. The industry-wise distribution of the sample is as per Table 2. In the past 3 years, many firms in the Crude Petroleum and Natural Gas industry were negatively affected by the decline in global oil prices. Therefore, there is a high concentration of firms in this industry in our sample of financially distressed firms. The high number of bankrupt companies in year 2016 mainly comprised of companies in the crude petroleum and natural gas industries. Moreover, the worldwide oil price issue gave this industry more potential for being financially distressed and thus became an interesting case to examine (Ward 1994). The sample of non-distressed firms was formed by using matched-pair sampling technique. This technique has been used in many studies relating to financial distress, such as the studies done by Agrawal (2015), Beaver (1966), , and Begley et al. (1996). Matching pairs of non-distressed firms were chosen based on matching industries, the nearest total assets, and the financial year of reporting. The matching of industries was done using the Standard Industry Classification (SIC) Code in the United States. As such, the overall sample size became162 firms, which wasmade up of 81 distressed firms and 81 non-distressed firms.
In order for us to have a meaningful study, it is vital to ensure that the differences of the total assets of the 2 samples of companies are not statistically significant (Agrawal 2015). To confirm this, an independent samples t-test was performed. The results from the independent samples t-test and descriptive statistics for the total assets of the both distressed firms sample and non-distressed firms sample are presented in Table 3. As the p-value is at 0.944, which is greater than 0.05, this shows that the difference in total assets for the distressed firms sample and non-distressed firms sample are not statistically significant. In addition, the overall sample was split into an estimation sample and a hold-out sample. Firms in the period from 2015 to 2017 were taken to represent the estimation sample, whereas the remaining firms from 2009 to 2014 were taken to represent the hold-out sample. Therefore, the estimation sample comprised of 57 distressed firms and 57 non-distressed firms. On the other hand, the hold-out sample comprised of 24 distressed firms and 24 non-distressed firms. As per a similar study done by Agrawal (2015), this study also used logistic regression to predict the probability of default for a firm based on the respective aggregate Piotroski's F-Score for each firm and also the nine individual components that made up the aggregate F-Score. Based on logistic regression, the probability of an event happening was established as follows: where P (Y) = Probability of event Y happening z = linear combination of independent variables as represented by: β 0 + β 1×1 + β 2 X 2 + β 3 X 3 + . . . + β n X n + µ.

Variable Descriptions
Model 1 uses the aggregate Piotroski's F-Score (Piotroski 2000) calculated as the total of the following independent variables. Though the F-Score comprised of nine components, only seven components werebeing used in Model 2 as the variable ROA and Accruals were removed due to high multicollinearity. In addition, the linearity of all the remaining seven components of F-Score with respect of the logit of the outcome was assessed via the Box-Tidwell procedure (Box and Tidwell 1962). A Bonferroni correction was applied using all 15terms in the model, resulting in statistical significance being accepted when p< 0.00333 (Tabachnick et al. 2007). Based on this assessment, all the components were found to be linearly related to the logit of the outcome except CFO. Therefore, the CFO variable was recoded into a nominal variable. Thus, the components of the Piotroski's F-Score used in Model 2 are ∆ROA, CFO, ∆Margin, ∆Turnover, ∆Leverage, ∆Liquidity and Eq_Offer.

Results
Independent t-tests were run for Model 1 and Model 2 and the results together with their respective descriptive statistics arepresented in Tables 4 and 5, respectively. It is observed that on average, distressed firms tend to have a lower Piotroski's F-Score as compared to non-distressed firms. On the other hand, for Model 2, the decline in ROA for distressed firms is significantly greater than non-distressed firms. It is also observed that on average, distressed firms tend to record negative CFO as compared to non-distressed firms.
Logistic regression was run to determine the effects of Piotroski's F-Score on the likelihood that a firm will be in financial distress. Both models exhibited a statistically significant chi-square value at a 0.01 level of significance. It was observed that the Nagelkerke R 2 value of Model 2 is higher than Model 1, which implies that Model 2 explained 20.7% of the variations in the dependent variable in comparison to 11.1% in Model 1.  This study also found that the F-Score has a significant negative relationship with the probability of financial distress, which means that firms with a lower F-Score have ahigher probability to go into financial distress (See Table 6). These findings are consistent with the study done by Agrawal (2015). This finding provides a clear indication for bankers, shareholders, and creditors that they should avoid investing in or giving credit to those firms that have a low F-score. However, the distressed firms can also be identified by looking at individual measures such as ROA and CFO, as this study's findings show that the decline in ROA for distressed firms is higher than non-distressed firms and they tend to record a negative CFO. Note: ***, ** and * denotes statistical significance at 0.01, 0.05, and 0.10 levels, respectively.
As it was found that there are collinearity issues with ROA and Accruals, a new aggregate F-Score (F-Score 2) was calculated and logistic regression was run for F-Score 2. It was found that the F-Score 2 also exhibited a statistically significant chi-square at the 0.00 significance level (please refer to Appendix A for results for SPSS). As this result is similar to the results for Model 1, these two variables were retained in the calculations of the aggregate F-Score in Model 1.
On top of that, in Model 2, ∆ROA was found to have a significant inverse relationship with probability of distress. This signifies that a greater decline in ROA is linked to a higher probability of distress. On the other hand, CFO has asignificant positive relationship with the probability of distress, which means that distressed firms tend to record a negative CFO. However, in the study by Agrawal (2015), ∆ROA and CFO were insignificant in relation to the probability of distress.
Tables 7 and 8 shows the classification matrix of Model 1 and Model 2 and both models are shown to have the same predictive ability in the estimation sample. However, Model 1 has a higher overall accuracy of 79.2% in the hold-out sample as compared to 72.9% in Model 2. This shows that Model 1 is a better model in the prediction of financial distress.  Notes: a indicates the percentage of distressed firms correctly classified as distressed, b indicates the percentage of non-distressed firms correctly classified as non-distressed, whereas c indicates the overall accuracy estimated as the average of a and b .

Discussion
The finding in this study contradicts the results found by Agrawal (2015), which showed that the predictive ability of the model which uses the individual components of the F-Score (Model 2) is higher than when using the aggregate Piotroski's F-Score (Model 1). However, the variables in Model 2 used in this study comprised of only seven components of the F-Score, whereas the Model 2 investigated in the study by Agrawal (2015) used all nine components as per the original F-Score introduced by Piotroski (2000). In this study, the variables ROA and Accruals wereremoved due to high multicollinearity with the ∆ROA variable. Thus, it is unfeasible to compare the predictive ability of Model 2 in this study versus Model 2 from the study by Agrawal (2015) due to the different components used in both models.
As there was insufficient evidence that fundamental analysis is useful for the Indian market , it is difficult to explain the different findings between this study and the results found by Agrawal (2015). As Piotroski (2000) first modeled the F-Score using data from firms in the United States, it is possible that the aggregate F-Score is more suited for the US market.Nevertheless, the findings in this study and the study by Agrawal (2015) revealed that the F-Score can be useful in bankruptcy and financial distress prediction. The usage of the F-Score in predicting financial distress has not been widely researched and only Agrawal (2015) wasfound to have conducted a similar study for the Indian market. This study contributes to existing literature by providing evidence for the usefulness of the F-Score as a new financial distress prediction model.

Conclusions
Predicting a firm's financial distress hasalways been a topic of interest for researchers and practitioners, as it affects bankers, shareholders, potential investors, and creditors. On many occasions, a group of firms' inability to meet their liabilities negatively affects the economic stability of a country. Based on the current literature, researchers have mostly used accounting information based models to predict the financial distress of the firms. However, the findings using accounting information based models are not always consistent nor robust; as such, this study used alternative models that use Piotroski F scores forpredicting financial distress of a sample from the United States of America.
The results of this study present some significant implications for investment and lending decisions. It shall help bankers, shareholders, potential investors, and creditors in terms making investment decisions by knowing which companies might end up in financial distress in no time. However, it should also be highlighted that the current accountingbased models are simple and easy-to-use for predicting the likelihood of financial distress for any firm.

Limitation and Future Research
One of the limitations of this study is that the samples only include public listed companies in the US, with the average assets size in excess of $2 billion. In comparison with large corporations, small firms typically behave differently due to their size, lack of economies of scale, market power, and their ability to innovate. As such, the factors which can lead to bankruptcy and financial distress in smaller companies are different from larger companies. Therefore, it is uncertain if the F-Score is useful in predicting financial distress in smaller firms and more research can be done to explore this area.
In addition, as Piotroski (2000) first developed the F-Score to select best valued stocks for investments, the market value of equity is incorporated into the model as a means to capture market variable. Thus, the model is only useful for companies which are listed on stock exchanges and where the required information is available to the public. Moreover, the F-Score was developed using firms in the United States, whose financial market is one of the more transparent and robust markets in the world. In countries where there is poor transparency and regulations, the performance of the F-Score in distress predictions can be different.
Although the sample used in this study spanned across many industries, there is a heavy concentration in the oil and gas industry due to the collapse of oil price in the last few years. Additionally, the sample size used in this firm is relatively small. Moreover, the study did not consider the most recent year data in the analysis. Therefore, future studies should be conducted by including more companies from other industries with a larger sample size and various asset sizes. Plus, future studies can consider matching multiple healthy companies with a similar size but different F-scores instead of considering only one matching company. It is also proposed that future research can be conducted collecting data for sufficient number of years and running some time series analysis to see whether the study can have better conclusions about predicting firms' financial distress. Finally, advanced data analysis techniques such as machine learning models can be used for bankruptcy prediction. Future studies can use alternative methods in predicting defaults and compare with the existing findings from the F-score models.
Author Contributions: M.R. carried out the literature review and the statistical analysis, and drafted the manuscript. C.L.S. helped with the data collections, data analysis, and discussion. M.A.K.M. participated in the discussion, and communicated with the editor. All authors have read and agreed to the published version of the manuscript.
Funding: This research was partially funded by the Research Cell, Noakhali Science and Technology University.

Conflicts of Interest:
The authors declare that they have no competing interests.
Appendix A multiple healthy companies with a similar size but different F-scores instead of considering only one matching company. It is also proposed that future research can be conducted collecting data for sufficient number of years and running some time series analysis to see whether the study can have better conclusions about predicting firms' financial distress. Finally, advanced data analysis techniques such as machine learning models can be used for bankruptcy prediction. Future studies can use alternative methods in predicting defaults and compare with the existing findings from the F-score models.
Author Contributions: M.R. carried out the literature review and the statistical analysis, and drafted the manuscript. C.L.S. helped with the data collections, data analysis, and discussion. M.A.K.M. participated in the discussion, and communicated with the editor. All authors have read and agreed to the published version of the manuscript.
Funding: This research was partially funded by the Research Cell, Noakhali Science and Technology University.

Conflicts of Interest:
The authors declare that they have no competing interests.