Macro-Financial Linkages and Heterogeneous Non-Performing Loans Projections
An Application to Ecuador1

Contributor Notes

Authors’ E-Mail Addresses: fgrigoli@imf.org; mmansilla@imf.org; msaldias@imf.org

We propose a stress testing framework of credit risk, which analyzes macro-financial linkages, generates consistent forecasts of macro-financial variables, and projects non-performing loans (NPL) on the basis of such forecasts. Economic contractions are generally associated with increases in NPL. However, despite the common assumption used in the empirical literature of homogeneous impact across banks, the strength of this relationship is often bank-specific, and imposing homogeneity may lead to over or underestimating the resilience of the financial system to macroeconomic woes. Our approach accounts for banks’ heterogeneous reaction to macro-financial shocks in a dynamic context and potential cross-sectional dependence across banks caused by common shocks. An application to Ecuador suggests that substantial heterogeneity is present and that this should be taken into account when trying to anticipate inflections in the quality of portfolio.

Abstract

We propose a stress testing framework of credit risk, which analyzes macro-financial linkages, generates consistent forecasts of macro-financial variables, and projects non-performing loans (NPL) on the basis of such forecasts. Economic contractions are generally associated with increases in NPL. However, despite the common assumption used in the empirical literature of homogeneous impact across banks, the strength of this relationship is often bank-specific, and imposing homogeneity may lead to over or underestimating the resilience of the financial system to macroeconomic woes. Our approach accounts for banks’ heterogeneous reaction to macro-financial shocks in a dynamic context and potential cross-sectional dependence across banks caused by common shocks. An application to Ecuador suggests that substantial heterogeneity is present and that this should be taken into account when trying to anticipate inflections in the quality of portfolio.

I. Introduction

One of the legacies of economic crises is the surge in nonperforming loans (NPL) in the financial system. The empirical literature unambiguously uncovered the negative relationship between economic activity and asset quality, both for advanced and developing economies.2 Even with appropriate macro-prudential regulation, the impact of a deteriorating economic cycle on the accumulation of NPLs is almost unavoidable. In turn, persistently high impaired assets pose significant challenges for policy makers as they threaten financial stability by jeopardizing the solvency of the system and slowing credit growth or affecting its quality. Thus, from the perspective of policy makers, regulators, and supervisors, accurately stress testing the resilience of the financial system’s loan portfolio to a deterioration in economic growth, and assessing whether shocks would have a systemic impact is critical, as it provides information to anticipate recapitalization needs, elaborate financial regulation changes, and develop crisis preparedness tools.

There is a number of arguments that can explain NPL behavior via economic activity. Economic downturns are often accompanied by higher unemployment, which affects the ability of debtors to service their debt, ultimately leading to an increase in NPL. In contrast, economic growth allows households and corporations to keep their finances buoyant and stay current in their debt payments, which is reflected in higher quality of banks’ loan portfolio. At any rate, changes in the macroeconomic context may take time to affect credit quality owing to many factors, including financial sector dynamics. For instance, NPL may increase sometime after the slowdown in economic activity as borrowers draw down savings to face payments. In some cases, financial crises were preceded by credit booms which, due to the accounting effect in the ratio, were accompanied by unusually low NPL ratios that upon the cycle reversal resulted in a delayed NPL response to the downturn. Also, the observed rise of NPL ratios may occur with a lag because some banks may concentrate lending in sectors which typically enjoy a grace period.

While the literature generally assumes that the impact of economic activity on NPL is homogenous across banks, in reality these are likely to be affected heterogeneously by lower economic activity (and other shocks). For example, if some banks concentrate their lending into sectors whose performance is not correlated with GDP, the relationship may not be as strong. Alternatively, if real GDP growth is concentrated in a specific sector and a given bank has a well-diversified portfolio, it may not suffer the impact of economic contractions as much. In other words, the homogeneity assumption could severely bias the results, possibly leading to an over or underestimation of the extent to which a negative shock to economic activity could translate into a systemic financial crisis. Similarly, looking at the aggregate series of NPL for evaluating recapitalization needs can be problematic, as capital cannot be reallocated across banks.3 At the same time, the presence of common shocks (e.g., the global financial crisis or financial sector regulations affecting all banks) could generate correlation across banks, making the cross-section independence hypothesis assumed in the literature unrealistic.

This paper proposes a three-stage stress testing framework to conduct credit risk assessments based on an econometric approach that takes into account dynamics, heterogeneity, and cross-sectional dependence (CSD). We present an application to Ecuador, which case is particularly interesting as the country is currently suffering from a worsening of its economic conditions due to mix of external and domestic shocks. The results identify the presence of macro-financial linkages, with interrelations and feedback effects between real and financial variables. Under July 2016 oil projections, the forecasts for GDP growth are expected to negatively affect the NPL ratio.4 For the average bank, the short-run increase in NPL associated with one percentage point (pp) fall in real GDP growth is 0.15 pp, while the long-run effect reaches 0.55 pp. This would more than double the weighted average NPL ratio for the financial system over the next two years. However, such weighted average hides considerable heterogeneity, with banks reaching an NPL ratio between 2 and 15 percent over the next two years. Such segmentation could help increase effectiveness and accuracy to crises prevention efforts by helping supervisors focus on institutions that are more sensitive to cyclical fluctuations, especially those with relatively lower balance sheet buffers.

The rest of the paper is organized as follows. Section II places the paper in the literature on solvency stress testing, reviewing the main contributions in the literature. Section III presents the recent economic developments in Ecuador and shows some stylized facts. Section IV discusses the econometric approach and presents the results. Section V concludes.

II. Literature Review

Top-down credit stress testing has become increasingly popular over the past twenty years. Financial institutions as well as supervisory authorities periodically assess vulnerabilities to adverse macro-financial scenarios to evaluate the resilience of the financial system. Similarly, this stress testing approach has become common in multilateral surveillance, and has been widely used in the aftermath of the global financial crisis in macro-prudential assessments and to boost market confidence in advanced economies at the peak of the crisis.5 It generally comprises four components. 6 The first is the macro-financial scenario design; the second looks at the impacts of the designed scenarios on banks’ performance; the third consists of a solvency calculation module; and the fourth analyzes the contagion and feedback effects. This paper focuses on the first two elements, carrying out a solvency stress test along the lines of Vazquez et al. (2012) and Wezel et al. (2014).7, 8

The macro-financial scenario design is usually based on statistical modeling. Jones et al. (2004) discusses advantages and disadvantages of structural models, including assumptions and restrictions. An alternative is the use of the vector autoregressive (VAR) models (or their vector error correction representation), which are sometimes preferred for their flexibility, smaller set of requirements, and ease of interpretation, as suggested by Åsberg and Shahnazarian (2008). Examples of these models are Hoggart et al. (2005), Van den End et al. (2006), Bank of Japan (2007), Jiménez and Mencía (2007), Vazquez et al. (2012), and Beck et al. (2013), where the selection of variables depends on the country’s characteristics. A similar approach is the global VAR (GVAR) used by Castrén et al. (2008) and Haldane et al. (2007), in which domestic and foreign variables interact simultaneously. However, data requirements associated with this technique often turn out to be prohibitive and weak erogeneity assumptions are required.9 A third approach to modeling macro-financial linkages is a pure statistical one based on simulations (Boss et al., 2006), which has the advantage of allowing for differences between marginal and multivariate distributions and changes in the correlation in stress scenarios. However, by nature it has limitations for policy analysis.

The models to map the macro-financial scenarios into stress scenarios at bank level can be divided into models based on borrower-level data and models assessing loan performance (Cihák, 2007). The former demands extensive information on the default risk at the household and corporate sector levels, and as a result is less common. The latter is often relying on NPL data (or loan loss provisions and default rates), and can be run at different degrees of aggregation (economy-, sector-, region-, or bank-level data, depending on availability). Our approach belongs to the latter group of models, given data availability for Ecuador.

Generally, the econometric approaches for loan performance rely heavily on panel data estimations of the determinants of loan impairment. The estimation technique ranges from ordinary least squares (OLS) with or without fixed effects, static and dynamic, and employing instrumental variables and cointegration techniques, to GMM-based methods (Arellano and Bond, 1991; Blundell and Bond, 1998). Applications of these sorts have the advantage of allowing for an easy mapping of the macro-financial forecasts and accommodating short samples and a large number of units. Recent examples are Vazquez et al. (2012) for Brazil and Wezel et al. (2014) for a small open economy with a large banking sector.10

The methods used in the literature only accommodate heterogeneity across banks by using fixed effects and do not commonly deal with issues of CSD caused by common shocks. However, both issues are likely to be present in bank-level data and impose a bias in the estimates if not accounted for. As an example, Memmel et al. (2014) show that common factors in Germany can represent a significant portion of credit losses. To address this issue, Pesaran et al. (2006) and Henry et al. (2013) propose using GVAR models, which require both the number of financial institutions and the timespan to be large. While being more flexible, these method requires large cross-section samples that are not generally available for emerging markets.

In the case of Ecuador, the dataset consists of a relatively small sample of banks and a somewhat long sample. Banks have different sizes in terms of assets, and they are heterogeneous in terms of their financial conditions and their ability to withstand shocks. Also, a high degree of cross-section dependence is likely to be present in the data, as banks are concentrated in certain credit segments (e.g., corporate and consumption) and at times geographically (e.g., coastal areas, possibly specialized in attending clients operating in the fishing or tourism sector).11 An appropriate statistical method would then require to address CSD while keeping the model suitable for a relatively smaller cross-section dimension, and address heterogeneity across banks.

III. Stylized Facts

In recent years Ecuador’s economy and its financial system have come under pressure due to the combined fallout of two significant shocks: (i) the oil price fall, and (ii) a real appreciation of the US dollar. Figure 1 presents a series of charts that depict this deterioration. As shown in the top left panel, real GDP growth has been highly correlated with the real oil price over the past years, even if oil price movements affect the Ecuadorian economy with some delay. As of the end of 2015, real GDP growth fell to about zero percent.

Figure 1.
Figure 1.

Variables’ Correlation and Dispersion

(Percent, yoy, unless otherwise specified)

Citation: IMF Working Papers 2016, 236; 10.5089/9781475559347.001.A001

Source: Authors’ calculations.

Since 2004, average NPL declined substantially from about 10 percent to less than 4 percent before the global financial crisis. When the crisis hit, NPL rebounded to above 5 percent, but declined again thereafter. With the end of the commodity supercycle, NPL started to rise again and are now above the levels observed in the aftermath of the global financial crisis. During the first quarter of 2016, NPL reached 3.9 percent. The top right panel contrasts the NPL ratio with real GDP growth, and clearly depicts a negative correlation.12

The Ecuadorian banking system is highly concentrated and composed by small, medium, and large banks, and the asset quality varies considerably across them. The middle left panel shows the median NPL ratio, along with the inter-quartile range and the 10th and 90th percentiles. At the beginning of the sample, only 10 percent of the banks had NPL ratios above 12 percent, but the median was as low as 6 percent. The deterioration started in 2014 presents much narrower bands around the median compared to previous critical periods, suggesting that NPL increased for all banks. The median NPL in the first quarter of 2016 reached 5.8 percent, close to the historical maximum in the sample.

Starting in early 2015, commercial banks experienced steady deposit withdrawals that lowered significantly their liquidity buffers. As suggested by the co-movement of real deposit and credit growth rates in the middle right panel, banks reacted by rationing credit. Interestingly, the real credit growth series seems to follow the deposit growth path with some delay, suggesting that banks have been able to lower lending in light of falling deposit levels. In this regard, the recent economic deterioration is not an exception.

At the same time, heterogeneity across banks in terms of deposit and credit growth is substantial. The lower left and right panels present the real credit growth and real deposit growth dispersion, respectively. These panels indicate that changes in real deposits and credit have become more homogenous over time. However, while the recent fall in deposits has been a common factor to all banks, real credit growth varied, and for more than 10 percent of banks it was still above 10 percent.

IV. Results

In this section we first describe the empirical strategy, and then present the results of the estimations.

A. Econometric Approach

We frame the empirical strategy in three stages. The first stage models the macro-financial linkages for Ecuador and generates forecasts that account for feedback effects between the real and the financial sectors. The second stage consists of a bank-level panel data estimation of the impact of real and financial variables on NPL. The third stage builds on the first two and simulates bank-specific NPL responses to the forecasted dynamics.

In the first stage, we estimate a VAR (p) representation of the macro-financial linkages by ordinary least squares (OLS) method for the aggregate financial system:

zt=b+Σj=jpDjzt1+Σj=1pFjxtj+et(1)

where the vector zt includes real GDP, real credit provided by private banks, and real deposits of private banks (all in logs); xt denotes the exogenous variable (log of) real price of oil; b is a vector of intercepts; D and F are a matrix and a vector of coefficients, respectively; and et is a zero mean white noise vector of errors. Such parsimonious specification captures the linkages between real and financial variables, allowing for exogenous shocks coming from changes in real oil prices. In the case of Ecuador, oil prices are a proxy of liquidity in the financial system through exports and through fiscal revenues.13 We embed inflation in the system by expressing deposits and credit in real terms.

With variables integrated of order one, and in absence of cointegration, the VAR can still be estimated provided that the eigenvalues lie within the unit circle, i.e. the dynamic system is stable. We can recover the structural VAR (SVAR) representation from the reduced form in equation (1):

A0zt=α+Σj=1pAjztj+Σj=1pFjxtj+ɛt(2)

where ɛt is a vector of shocks. As shown in Fernández-Villaverde et al. (2005), since Aj = A0D and ɛt = A0et, the mapping between et and ɛt is given by A01ɛt.

Given that A01 contains nine parameters and we only have six distinct covariances in Σe, we need to impose restrictions on contemporaneous relations between variables. Thus, we assume the following Cholesky decomposition, with the variables set in the aforementioned order:

A01=[a1100a21a220a31a23a33](3)

Such restrictions are plausible and have a sound economic meaning. Consistent with the stylized facts described above, the ordering of the variable assumes that real shocks affect the financial sector within the same period, and in particular, that increases (withdrawal) in deposits are reflected in credit expansion (rationing) during the same period. Thus, the impulse response functions help depicting the relationships among the endogenous variables of the VAR model and therefore the macro-financial linkages at work in the economy.

Finally, we generate forecasts of the endogenous variables, conditional on the evolution of the exogenous one. In other words, we obtain consistent projections for real GDP, real deposit, and real credit (and their growth rates) for a two-year horizon based on the IMF-World Economic Outlook (WEO) projection for oil prices.

The second stage of the empirical approach relies on a bank-level panel dataset of NPL to quantify the sensitivity of NPL to changes in macro-financial conditions. In particular, we model the growth rate of the logistic transformation of NPL for bank i at time t, yit, with the following Autoregressive Distributed Lag (ARDL) specification:14

yi,t=α+Σk=1nβkxtx+Σk=1nδkzi,tk+ui,t(4)

where xt is the real GDP growth and zit denotes the bank-specific variables real deposit growth and real credit growth; β and δ are the relative coefficients; a is the constant term; η is the maximum number of lags; and ɛit is the idiosyncratic disturbance term assumed to be independent across banks and serially uncorrelated. Other variables, including real salaries, real GDP growth of trading partners, and global interest rates resulted non-significant and were dropped from the specification. The parsimonious model allows to create a direct mapping with the variables used in the first stage of the analysis.

We estimate equation (4) using pooled ordinary least squares (POLS) applied to the panel sample of quarterly observations, correcting standard errors for heteroskedasticity and autocorrelation. Such specification also deals with endogeneity by excluding the contemporaneous terms of the independent variables. However, this estimation suffers from other econometric issues: lack of dynamics, omitted variable bias, parameter heterogeneity across banks, and CSD.

Dynamics of the dependent variable are likely to be an important factor in the estimation because NPL is generally persistent. Thus, we specify a target adjustment model yi,tyi,t-1 = (1—γ)(Yi,tyi,t-1), where γ is the adjustment parameter. Thus, if γ = 0, then yit = Yit, meaning that the adjustment takes place immediately. We then introduce dynamics in the following equation:

yi,t=α+γyi,t1+Σk=1nβkxtx+Σk=1nδkzi,tk+ui,t(5)

where an incomplete adjustment for which γ ≠ 0 leads to a form of state dependence where yi,t-1 determines yi,t. The omission of the lag of the dependent variable would result in unobserved heterogeneity from the correlation between yi,t-1 and yi,t To further address the omitted variable bias, we modify equation (4) to include bank-specific fixed effects (OLSFE).

Most of the empirical literature in this area addresses endogeneity concerns by estimating some version of the Generalized Method of Moments (GMM) estimator (Vazquez et al., 2012; Klein, 2013; Wezel et al., 2014), instead of allowing for lagged effects. The difference GMM estimator assumes that the idiosyncratic error is serially uncorrelated and that past values of the endogenous variables are not correlated with the current error. These conditions allow the use of the second lag (and higher) of the dependent variable as instruments for its first lag, and second (and higher) lags of endogenous variables as instruments for the endogenous variables. Blundell et al. (2000) and Bond et al. (2001) show that the difference GMM estimator has poor finite sample properties and that the estimator performs weakly in samples with limited time dimension and when the dependent variable is persistent. Thus, Arellano and Bond (1991) and Blundell and Bond (1998) propose the system GMM (SGMM) estimator, which increases efficiency by estimating a system of two simultaneous equations, one in levels (with lagged first differences as instruments) and the other in first differences (with lagged levels as instruments). This estimator requires the additional identifying assumption that the instruments are exogenous to the fixed effects. Thus, we estimate the following equation with the asymptotically more efficient two-step SGMM:

yi,t=α+γyi,t1+Σk=0nβkxtk+Σk=0nδkzi,tk+ci+ui,t(6)

where ci represents the unobserved bank-specific heterogeneity.15

All models discussed so far (and generally employed in the literature on NPL determinants) neglect parameter heterogeneity, which is likely to be a major issue in the data at hand. Pesaran and Smith (1995) and Haque et al. (1999) show that in cross-country panel data the assumption of slope homogeneity may not hold, and this leads to inconsistency of the estimates. Thus, they propose the mean group (MG) estimator for stationary panels by Pesaran and Smith (1995), which accounts for parameter heterogeneity and constructs simple mean estimates across the estimates derived from separate bank regression.

In addition, Phillips and Sul (2007) show that if CSD exists across units, then estimated parameters may be significantly biased and identification problems may be present. To deal with this, Pesaran (2006) present the common correlated effects MG (CCEMG) estimator, and Eberhardt and Bond (2009) and Eberhardt and Teal (2010) propose the augmented MG (AMG) estimator, both of which account for parameter heterogeneity and CSD, albeit in a different way. The CCEMG estimator augments the regression equation with cross-section averages of dependent and independent variables (to be interpreted as nuisance), as a way to control for the unobserved common factors, while the AMG regards the unobserved common factors as the common dynamic process and estimates it in two steps. First, it estimates a pooled difference OLS model with time dummy variables and saves the estimated coefficients as the common dynamic process. Second, the common dynamic process is added to the regression equation either by subtracting it from the dependent variable (I-AMG) or by including it in each of the bank-specific regressions. Finally, both for the CCEMG and the AMG estimators, bank-specific estimates are averaged as in Pesaran and Smith (1995). In this paper, we estimate the following equation using the AMG estimator:

yi,t=αi+γiyi,t1+Σk=1nβi,kxtk+Σk=1nδi,kzi,tk+diμ^t+ui,t(7)

where the vector μ^t contains the quarter dummy coefficients extracted from the pooled regression in first differences.

In the third and last stage of the empirical strategy, we follow Vazquez et al. (2012) and project NPL two years out of sample. While the AMG estimation of equation (7) provides the average coefficients α, βk, δk, and d, the VAR provides a consistent set of forecasts for xt-k and zt-k that account for macro-financial linkages. Furthermore, we assume the common dynamic factor μ^t to maintain the same value as in the last observed period. Bank-specific NPL projections are then averaged using the stocks of credit as weights.

B. Analysis

We start by analyzing the stationarity properties of the series used in the VAR model.16 Our sample starts in the third quarter of 2002 and extends to the fourth quarter of 2015. Figure 2 plots the logs of real GDP, real deposits, real credit, and real oil price, as well as the first differences in percent. A visual inspection suggests that first differences are stationary, while levels appear to be integrated of order one. We formally test for the presence of unit root with the augmented Dickey Fuller (ADF) test. In particular, we perform the tests excluding deterministic components, including the intercept, and including the intercept and a trend. The results corroborate the findings of the visual inspection. 17 Namely, variables in levels are non-stationary, while the first differences are stationary. Despite the presence of unit roots, however, the series in levels do not cointegrate.18

Figure 2.
Figure 2.

Stationarity of VAR Series

Citation: IMF Working Papers 2016, 236; 10.5089/9781475559347.001.A001

Source: Authors’ calculations.

Conventional tests reveal that the VAR model is well specified. Standard information criteria are used to select the lag length of the VAR. The likelihood-ratio test statistic, the final prediction error, the Akaike’s and the Hannan and Quinn information criteria suggest two lags, while the Schwarz’s Bayesian information criterion and suggests one lag.19 Therefore, we opt for two lags. Despite the absence of cointegration, a VAR can still be specified in levels as long as it is stable. The test for VAR stability reveals that the eigenvalues lie in the unit circle.20 While the normality of the residuals is not required for the generation of the impulse response functions (IRFs), these have still to be serially uncorrelated. Thus, we perform a Lagrange multiplier test and the Portmanteau test, which do not reject the null of no autocorrelation.21

We then estimate the IRFs, which describe the interactions among real GDP and financial variables over three years, thereby providing a picture of the macro-financial linkages. We first analyze the response of real GDP. As shown in Figure 3, real GDP reacts contemporaneusly to shocks in its level, and the effect tends to be persistent. A positive shock to deposits has a significant and short-lived impact on real GDP, as the effect dies out after six quarters. Finally, real GDP does not show a significant reaction to real credit.22 We now turn to responses of real deposits. A one-standard deviation shock to real GDP produced an increase in real reposits that lasts about two years. The impact of real deposits on itself is strongly significant and takes longer than three years to set in. Shocks to real credit generate an insignificant effect on real deposits, which becomes insignificant after eight months. Finally, we analyze responses of real credit. A one-standard deviation shock to real GDP has a significant effect on real credit only after three quarters, and it lasts for about three years. Real credit also reacts positively to shocks in real deposits, however in this case the effects is immediate, and also sets in after three years. Finally, a shock in real credit has an immediate and short-lived effect on itself, as it phases out in three quarters.

Figure 3.
Figure 3.

Macro-Financial Linkages

(Quarterly responses to one-standard deviation shocks with two-standard-deviation confidence intervals)

Citation: IMF Working Papers 2016, 236; 10.5089/9781475559347.001.A001

Source: Authors’ calculations.

Finally, we proceed to generate forecasts for the endogenous variables of the VAR model. The forecasts are conditional on the expected oil price projections of the IMF-WEO over 2016-17, and provide a consistent set of estimates that take into account the feedback effects among real and financial variables. Figure 4 presents the forecasts in growth rates. Under current oil price assumptions, real GDP growth is projected to remain in negative territory through the end of 2017, albeit recovering above current levels from the trough in the last quarter of 2016 due to higher projected oil prices. Real deposit growth is projected to continue falling until mid-2016 and recover thereafter. Real credit growth presents a somewhat delayed fall with respect to real deposits, hitting the bottom at the end of 2016.

Figure 4.
Figure 4.

SVAR Conditional Forecasts

(Forecast in percent, yoy, with one-standard-deviation confidence intervals)

Citation: IMF Working Papers 2016, 236; 10.5089/9781475559347.001.A001

Notes: Forecasts are conditional on the IMF-WEO oil price projections.Source: Authors’ calculations.

We now move to the second step of the analysis, which analyzes the relationship between the macro-financial variables and NPL using a bank-level panel dataset.23 For this purpose, we rely on a sample that covers 22 banks over the decade 2005-15.

Common shocks in this context can be both external (e.g., the global financial crisis) and domestic (e.g., financial regulation changes affecting all banks) and, if not appropriately accounted for, can cause CSD, which can bias the parameters of interest (Philips and Sul, 2003; Andrews, 2005). Non-linearity can be observed if the deterioration in NPL accelerates when real GDP growth surpasses certain thresholds. For example, one could argue that NPL would increase faster when real GDP growth is in negative territory, or that the relationship would flatten out for high levels of real GDP growth. Similarly, it is very likely that banks react differently to shocks in real GDP growth. For example, as banks concentrate their lending into sectors whose performance is not necessarily correlated with GDP, the relationship may not be strong. Also, if real GDP growth is concentrated in, say, oil-related sectors and a given bank has a well-diversified portfolio, it may not suffer real GDP growth decelerations.

Figure 5 illustrates the potential for the distorting effects arising from common shocks, non-linearity, and CSD in the dataset. The upper panel shows the maximum NPL ratio by bank over time. There is a clear clustering around specific dates which can be associated to common shocks, e.g. the global financial crisis in 2009 and the recent worsening of economic conditions in 2016. While this is only a prima facie evidence, it suggests a remarkable presence of common shocks. The mid panel depicts a fractional polynomial regression line (along with a 95 percent confidence interval) for NPL against real GDP growth. While this descriptive analysis is highly stylized and there are other factors beyond real GDP growth that affect NPL, it still suggests that non-linearity is less of a concern for our dataset. Finally, the lower panel plots the same fractional polynomial regression for every bank, as a way to test whether departing from the assumption of homogenous parameters is justified. The chart illustrates well the potential for misspecification in the NPL-real GDP growth relationship, and suggests that heterogeneous parameters need to be introduced.

Figure 5.
Figure 5.

Relationship between NPL and Real GDP Growth

(Percent)

Citation: IMF Working Papers 2016, 236; 10.5089/9781475559347.001.A001

Notes: In the mid panel the maximum of the horizontal axis has been set to 10 o ease the visual inspection of the data.Source: Authors’ calculations.

Table 1 presents the results of the panel estimations. Column 1 shows the results of the static POLS. In Column 2 we control for the persistence of NPL by adding its lag, and we also introduce bank-specific fixed effects to deal with omitted variable bias. In column 3, we report the results of the SGMM, which addresses endogeneity using a different instrumentation technique.24 In Column 4, we allow for parameter heterogeneity using the MG estimator. Finally, in Columns 5 and 6 we also account for CSD. Column 5 presents the results using the AMG-I estimator, and Columns 6 using the AMG estimator. Our preferred model is presented in Column 6 for which the results from the empirical tests (Pesaran, 2004) suggest that CSD is effectively removed. Three of the 22 banks representing the private financial sector in Ecuador are dropped from the dataset due to their very short time series. Depending on the estimator used, the number of observations ranges between 704 and 718. In the case of the AMG estimators, a fourth bank is dropped due to the estimation technique.

Table 1.

Panel Regressions Results

(Dependent variable: logistic transformation of NPL)

article image
Notes: Standard errors in parentheses are corrected for heteroskedasticity and autocorrelation of the error term. AMG estimations correct for cross-sectional dependence. The SGMM estimation uses a collapsed instrument matrix and performs the Windmeijer (2005) correction of the covariance matrix. The null hypothesis for the Hansen J -test is that the full set of instruments is valid. ***, **, * next to a number indicate statistical significance at 1, 5 and 10 percent, respectively. Source: Authors’ calculations.

Unsurprisingly, the results present some commonalities and some differences across estimators, consistent with the idea that they address progressively a larger set of econometric issues. The lag of the dependent variable is always significant, suggesting persistence in NPL. Real GDP growth is a strong determinant of NPL. In all estimations, the first lag of real GDP growth is negative and significant. The second lag is generally not significant, but in our preferred model that controls for CSD, it still has a negative and significant impact. The third lag is always positive but not always significant, suggesting some quick reversal following a shock. The results for other variables are not consistent across estimators, but the coefficients present much smaller magnitudes.25

The ARDL specification allows retrieving the short-run and the long-run impact. More formally, the former is calculated as Σk=1nxi,k, while the latter is equal to Σk=1nxi,k/(1α). Figure 6 presents an overview of the average impact of a one pp shock in real GDP growth on NPL, calculated at the NPL ratio observed in 2016Q1 for every bank and using bank-specific coefficients. It is clear that the impact can be different by bank and that, in some cases, it contravenes theory. In our case, for example, the short-run and the long-run effects are counterintuitively positive in four banks, albeit close to zero. Bank-specific coefficients should be interpreted as merely indicative as individual estimates may provide weak signals, while averages represent very plausible estimates (Boyd and Smith, 2002; Baltagi et al., 2003).26 Relying on average coefficients, we can calculate a rule of thumb in the case of Ecuador, for which in the short run NPL will increase by 0.15 pp for every pp fall in real GDP growth, while in the long run NPL would increase by 0.55 pp.

Figure 6.
Figure 6.

Effect of Real GDP Growth Shock on NPL

(Impact of one pp real GDP growth shock evaluated at 2016Q1 NPL ratio)

Citation: IMF Working Papers 2016, 236; 10.5089/9781475559347.001.A001

Notes: Banks are ordered by the magnitude of the short-run effect. Source: Authors’ calculations.

We finally move to project the path of NPL weighted average over the next two years. Figure 7 depicts the NPL projection through the end of 2017 using the variables forecasted in the first stage and the coefficients obtained in the second stage. The confidence intervals generated from the VAR estimation are used to project a pessimistic and an optimistic scenario. The results suggest that in the baseline scenario, under the July 2016 oil price projections, NPL may increase from 3.8 percent in March 2015 to 6.4 percent at end-2016 and 9.3 percent at end-2017, peaking up to 9.7 percent during the second quarter of 2017. These results can help authorities to enhance their preparedness against anticipated critical capital needs (or upcoming credit expansions) at the aggregate level.

Weighted (or simple) averages hide heterogeneity across banks. Figure 8 shows the deterioration in the financial system for the three scenarios derived from the VAR forecasts. As shown, there is great heterogeneity across banks even under the baseline scenario, in which some banks reach an NPL ratio of almost 15 percent and others do not pass 2 percent. This more granular outcome would allow supervisors to better plan their activities at the micro level and map banks’ strategies according to their expected NPLs paths.

Figure 7.
Figure 7.

Weighted Average of NPL Projections

(Percent)

Citation: IMF Working Papers 2016, 236; 10.5089/9781475559347.001.A001

Notes: Confidence bands of the VAR forecasts are used to build optimistic and pessimistic scenarios. Source: Authors’ calculations.
Figure 8.
Figure 8.

Heterogeneity in NPL Projections

(Percent)

Citation: IMF Working Papers 2016, 236; 10.5089/9781475559347.001.A001

Source: Authors’ calculations.

V. Conclusions and Policy Implications

It is common to observe high NPL ratios during downturns, owing to the contraction in economic activity and the consequent reduced ability of borrowers to service their debt. In turn, high and persistent NPL ratios may eventually affect the stability of the financial system by raising solvency concerns and slowing down credit growth or affecting its quality. It is therefore crucial to measure to the best extent possible the resilience of banks to a worsening in economic conditions.

Stress testing to conduct solvency assessment became increasingly popular over the past twenty years. However, most econometric analyses assume homogeneity and independence across banks. In reality such assumptions are unrealistic due to many factors, including banks’ varying degree of specialization and diversification across borrowers and sectors of activity. Such heterogeneity becomes particularly relevant when assessing recapitalization needs, as capital cannot be reallocated across banks, warranting a more granular approach. At the same time, shocks such as global crises or changes in the financial regulation are common to all banks, inducing CSD in the sample.

In this paper, we present an application of stress testing to Ecuador in line with the literature, but accounting for banks’ heterogeneous reaction to shocks and CSD. After assessing the dynamics among real and financial variables, we generate forecasts for the same variables. Under July 2016 oil projections, the forecasts for macro-financial variables are expected to negatively affect the NPL ratio. Our results suggest that for the average bank, the short-run effect of a one pp fall in real GDP growth is 0.15 pp, while the long-run effect reaches 0.55 pp. This would more than double the weighted average NPL ratio for the financial system over the next two years. However, such weighted average hides considerable heterogeneity, with banks reaching an NPL ratio between about 2 and 15 percent.

From a policy perspective, this exercise is designed to supplement regular surveillance and provide more granularity to stress testing tools. Being able to appropriately anticipate inflections in the quality of portfolio allows a more efficient surveillance of the financial system, by focusing on the supervisory agenda to protect financial stability, engaging with the banks’ management, and requiring mitigating measures in a timely manner, which could comprise additional provisions, new capital, better risk management, selective lending, among others. Also, the methodology presented would allow supervisors to target their efforts across institutions and possibly mitigate systemic risks via regulatory measures (e.g., by imposing generic provisions or enhancing macro prudential measures). Similarly, the methodology could serve to target efforts across different sub-segments of activity. For instance, to the extent information is available, the framework is adaptable to analysis by sector, business line, or region, which can help design and fine tune policy responses. In sum, the framework presented in this paper allows for a better quantification of NPL projections and therefore avoid over or underestimation of systemic risk.

Appendix I. VAR Specification Tests

Table A1.

Augmented Dickey-Fuller Unit Root Test for VAR Series

article image
Notes: The null hypothesis is that the series has a unit root. The lagged differences are included in the specifications to obtain white noise residuals. The Schwartz Information Criterion is used to select the optimal lag length. ***, **, * next to a number indicate statistical significance at 1, 5 and 10 percent, respectively. Source: Authors’ calculations.
Table A2.

Johansen Cointegration Test

article image
Notes: The null hypothesis is that the series do not have a cointegration relationship. The Schwartz Information Criterion is used to select the optimal lag length. The critical values are from MacKinnon (1991). ***, **, * next to a number indicate statistical significance at 1, 5 and 10 percent, respectively. Source: Authors’ calculations.
Table A3.

VAR Lag Selection

article image
Notes: * next to a number indicate the preferred lag lenght. Source: Authors’ calculations.
Table A4.

VAR Stability

article image
Source: Authors’ calculations.
Table A5.

VAR Residual Autocorrelation

article image
Notes: The null hypotheses of the Lagrange multiplier (LM) and Portmanteau Test is that there is no autocorrelation at the indicated lag. ***, **, * next to a number indicate statistical significance at 1, 5 and 10 percent, respectively. Source: Authors’ calculations.

References

  • Andrews, Donald W.K., 2005, “Cross-section regression with common shocks,Econometrica, Vol. 73, pp. 15511585.

  • Arellano, Manuel, and Stephen Bond, 1991, “Some Tests of Specification for Panel Data: Monte Carlo Evidence and an Application to Employment Equations,Review of Economic

    • Search Google Scholar
    • Export Citation
  • Studies, Vol. 58, pp. 277297.

  • Åsberg, Per, and Hovick Shahnazarian, 2008, “Macroeconomic Impact on Expected Default Frequency,Sveriges Riksbank Working Paper No. 219.

    • Search Google Scholar
    • Export Citation
  • Baltagi, Badi H., Georges Bresson, James M. Griffin, and Alain Pirotte, 2003, “Homogeneous, Heterogeneous or Shrinkage Estimators? Some Empirical Evidence from French Regional Gasoline Consumption,Empirical Economics, Vol. 28, No. 4, pp. 795811.

    • Search Google Scholar
    • Export Citation
  • Bank of Japan, 2007, “The Framework for Macro Stress-testing of Credit Risk: Incorporating Transition in Borrower Classifications,Financial System Report, September.

    • Search Google Scholar
    • Export Citation
  • Beck, Roland, Petr Jakubik, and Anamaria Piloiu, 2013, “Non-performing Loans: What Matters in Addition to the Economic Cycle?,Working Paper Series 1515, European Central Bank.

    • Search Google Scholar
    • Export Citation
  • Blundell, Richard, and Stephen Bond, 1998, “Initial Conditions and Moment Restrictions in Dynamic Panel Data Models,Journal of Econometrics, Vol. 87, pp. 115143.

    • Search Google Scholar
    • Export Citation
  • Blundell, Richard, and Stephen Bond, and Frank Windmeijer, 2000, “Estimation in Dynamic Panel Data Models: Improving on the Performance of the Standard GMM Estimators,The Institute of Fiscal Studies Working Paper, No. 00/12.

    • Search Google Scholar
    • Export Citation
  • Bond, Stephen, Anke Hoeffler, and Jonathan Temple, 2001, GMM Estimation of Empirical Growth Models, mimeo, September.

  • Boyd, Derick A. C., and Ron P. Smith, 2002, “Some Econometric Issues in Measuring the Monetary Transmission Mechanism, with an Application to Developing Countries.In Monetary Transmission in Diverse Economies, ed. Lavan Mahadeva, and Peter Sinclair, Cambridge University Press.

    • Search Google Scholar
    • Export Citation
  • Boss, Michael, Gerald Krenn, Claus Puhr, and Martin Summer, 2006, “Systemic Risk Monitor: A Model for Systemic Risk Analysis and Stress Testing of Banking Systems,Internal Technical Document, Oesterreichische Nationalbank.

    • Search Google Scholar
    • Export Citation
  • Camacho, María Isabel, Ivonne Cilio, and Juan Pablo Erráez, 2015, Nota Técnica No. 78. Banco Central del Ecuador.

  • Castren, Olli, Trevor Fitzpatrick, and Matthias Sydow, 2008, “Assessing Portfolio Credit Risk Changes in a Sample of EU Large and Complex Banking Groups in Reaction to Macroeconomic Shocks.Mimeo

    • Search Google Scholar
    • Export Citation
  • Cihák, Martin, 2007, “Introduction to Applied Stress Testing,IMF Working Papers 07/59, International Monetary Fund.

  • Drehmann, Mathias, Steffen Sorensen, and Marco Stringa, 2010, “The Integrated Impact of Credit and Interest Rate Risk on Banks: A Dynamic Framework and Stress Testing Application,Journal of Banking and Finance, Vol. 34, No. 4, pp. 713729.

    • Search Google Scholar
    • Export Citation
  • Eberhardt, Markus, and Stephen Bond, 2009, “Cross-Section Dependence in Nonstationary Panel Models: A Novel Estimator,MPRA Paper 17692, University Library of Munich, Germany.

    • Search Google Scholar
    • Export Citation
  • Eberhardt, Markus, and Stephen Bond, and Francis Teal, 2011, “Econometrics For Grumblers: A New Look At The Literature On Cross-Country Growth Empirics”, Journal of Economic Surveys, Vol. 25, No. 1, pp. 109155.

    • Search Google Scholar
    • Export Citation
  • Espinoza, Raphael, and Ananthakrishnan Prasad, 2010, “Nonperforming Loans in the GCC Banking System and their Macroeconomic Effects,IMF Working Papers 10/224, International Monetary Fund.

    • Search Google Scholar
    • Export Citation
  • Foglia, Antonella, 2009, ““Stress Testing Credit Risk: A Survey of Authorities’ Approaches,International Journal of Central Banking, Vol. 5, No. 3, pp. 945.

    • Search Google Scholar
    • Export Citation
  • Haldane, Andrew, Simon Hall, and Silvia Pezzini, 2007, “A New Approach to Assessing Risks to Financial Stability,Bank of England Financial Stability Working Paper No. 2, April.

    • Search Google Scholar
    • Export Citation
  • Haque, Nadeem U., M. Hashem Pesaran, Sharma Sharma,Neglected Heterogeneity and Dynamics in Cross-Country Savings Regressions,IMF Working Papers 99/128, International Monetary Fund.

    • Search Google Scholar
    • Export Citation
  • Henry, Jérôme, Christoffer Kok, Adrien Amzallag, Patrizia Baudino, Inês Cabral, Maciej Grodzicki, Marco Gross, Grzegorz Halaj, Markus Kolb, Miha Leber, Cosimo Pancaro, and Matt Sydow, 2013, “A Macro Stress Testing Framework for Assessing Systemic Risks in the Banking Sector,Occasional Paper Series 152, European Central Bank.

    • Search Google Scholar
    • Export Citation
  • Hoggarth, Glenn, Steffen Sorensen, and Lea Zicchino, 2005, “Stress Tests of UK Banks using a VAR Approach,Bank of England working papers 282, Bank of England.

    • Search Google Scholar
    • Export Citation
  • International Monetary Fund, 2015. Staff Discussion Note. A Strategy to Resolve Europe’s Problem Loans. SDN 15/19. September.

  • Jiménez, Gabriel, and Javier Mencía, 2007, “Modelling the Distribution of Credit Losses with Observable and Latent Factors,Banco de España, Working Paper No. 0709.

    • Search Google Scholar
    • Export Citation
  • Klein, Nir, 2013, “Non-Performing Loans in CESEE: Determinants and Macroeconomic Performance,IMF Working Paper 13/72, International Monetary Fund.

    • Search Google Scholar
    • Export Citation
  • Laeven, Luc, and Fabián Valencia, 2012, “Systemic banking Crises Database: An Update,IMF Working Papers 12/163, International Monetary Fund.

    • Search Google Scholar
    • Export Citation
  • Memmel, Christoph, Yalin Gunduz, and Peter Raupach, 2014, “Common Drivers of Default Risk,Journal of Financial Stability, Vol. 16, pp.232247.

    • Search Google Scholar
    • Export Citation
  • Mwanza, Nkusu, 2011, “Nonperforming Loans and Macrofinancial Vulnerabilities in Advanced Economies,IMF Working Papers 11/161, International Monetary Fund.

    • Search Google Scholar
    • Export Citation
  • Ong, LilIan and Ceyla Pazarbasioglu, 2013, “Credibility and Crisis Stress Testing,IMF Working Papers 13/178, International Monetary Fund.

    • Search Google Scholar
    • Export Citation
  • Pesaran, M. Hashem, 2004, “General Diagnostic Tests for Cross Section Dependence in Panels,IZA Discussion Paper No. 1240.

  • Pesaran, M. Hashem, Ronald P. Smith, 1995, “Estimating long-run relationships from dynamic heterogeneous panels,Journal of Econometrics, Vol. 68, No. 1, pp. 79113.

    • Search Google Scholar
    • Export Citation
  • Pesaran, M. Hashem, Ronald P. Smith, Til Schuermann, Bjorn-Jakob Treutler, and Scott Weiner, 2006, “Macroeconomic Dynamics and Credit Risk: A Global Perspective,Journal of Money, Credit and Banking, Vol. 38, No. 5, pp. 12111261.

    • Search Google Scholar
    • Export Citation
  • Phillips, D. Sul, 2007, “Bias in dynamic panel estimation with fixed effects, incidental trends and cross section dependence,Journal of Econometrics, Vol. 137, No. 1, pp. 162188.

    • Search Google Scholar
    • Export Citation
  • van den End, Jan W., Marco Hoeberichts, and Mostafa Tabbae, 2006, “Modelling Scenario Analysis and Macro Stress-Testing,De Nederlandsche Bank Working Paper No. 119.

    • Search Google Scholar
    • Export Citation
  • Vazquez, Francisco, Benjamin Tabak, and Marcos Souto, 2012, “A Macro Stress Test Model of Credit Risk for the Brazilian Banking Sector,Journal of Financial Stability, Vol. 8, No. 2, pp. 6983.

    • Search Google Scholar
    • Export Citation
  • Wezel, Torsten, Michel Chanta, and Manuel Luy, 2014, “A Practical Example of the Nonperforming Loans Projection Approach to Stress Testing,” in Ong, Li., A Guide to IMF Stress Testing: Methods and Models, International Monetary Fund, Washington.

    • Search Google Scholar
    • Export Citation
  • Windmeijer, Frank, 2005, “A finite sample correction for the variance of linear efficient two-step GMM estimators,Journal of Econometrics, Vol. 126, No. 1, pp. 2551.

    • Search Google Scholar
    • Export Citation
1

We thank Catherine Patillo, Peter Pedroni, Luca Ricci, Hui Tong, Fabián Valencia, and the participants to the WHD seminar held in July for their comments and suggestions. We also thank Adrian Robles Villamil for excellent research assistance.

2

See for instance Espinoza and Prasad (2010), Mwanza (2011), and Laeven and Valencia (2012) for empirical evidence across countries using alternative empirical approaches.

3

This is less of a concern when disaggregating across business lines (or loan types) for a given bank, as capital can be reallocated.

4

More favorable oil price projections would change the GDP forecasts and therefore the projected NPL ratios. With illustrative purposes, we use the most current oil price projections (July 2016) at the moment in which the paper is written.

5

For instance, the Fed’s Supervisory Capital Assessment Program in 2009 and the Committee of European Banking Supervisors in 2009 and 2010 conducted for the first time the so-called crisis stress tests (Ong and Pazarbasioglu, 2013) as an overall strategy to rebuild public confidence in a banking system in the aftermath of global financial crisis. They set the stage to use this type of stress tests as a tool in financial crisis management

6

See Foglia (2009) and Henry et al. (2013) for a comprehensive overview of stress testing techniques across countries and recent developments.

7

As for the third and four elements, we do not provide results in terms of capital adequacy ratios for this paper but they can be easily obtained by applying assumptions on loan loss provisions, and the analysis of contagion and feedback effects is beyond the scope of the paper.

8

See Drehmann et al. (2010) for an integrated approach for both solvency and interest rate risk.

9

While VAR models (and their variations) are extensively used, many authors expressed concerns about the inability of these models to capture non-linearities.

10

Vazquez et al. (2012) additionally applies the model to disaggregated credit loan portfolios by economic activity and then compute aggregate NPL ratios for each bank and the whole system.

11

For an account of concentration across credit lines, see Camacho et al., 2015.

12

Negative correlation is even stronger when looking at non-oil real GDP growth, however, real GDP growth is preferred here because it is generally the measure used for generating forecasts by the authorities.

13

Vazquez et al. (2012) include the yield curve slope in the model to account for relevant monetary policy shocks in Brazil, Bank of Japan (2007), Beck et al. (2013), and Hoggart et al. (2005) include stock prices, nominal or real exchange rates, and lending interest rates. Owing to Ecuador’s full dollarization, capped interest rates, and a bank-based financial system, these variables do not add information to the model.

14

In line with Vazquez et al. (2012) and Wezel et al. (2014), we apply the logistic transformation to the NPL ratio yi,t = ln[ni,t/(1–ni,t)], where ni,t is the NPL ratio, to create an unrestricted variable to be used in the regression and therefore avoid non-normality of the error term.

15

The two-step variant presents estimates of the standard errors that tend to be severely downward biased (Arellano and Bond, 1991; Blundell and Bond, 1998). However, we implement the finite-sample correction of the two-step covariance matrix derived by Windmeijer (2005), which produces unbiased standard errors.

16

All variables are seasonally adjusted using Census X-13.

17

See Table A1 in Appendix II.

18

See Table A2 in Appendix II.

19

See Table A3 in Appendix II.

20

See Table A4 in Appendix II.

21

See Table A5 in Appendix II.

22

This counterintuitive result might be associated to the fact that real deposits dynamics anticipate real credit ones, as shown in the previous section. Granger causality tests also confirm that changes in real deposits anticipate changes in real credit.

23

NPL, real deposit growth, and real credit growth are available at monthly frequency. We seasonally adjusted the monthly data using Census X-13 and calculated quarterly averages to match the frequency of the GDP series.

24

Ideally, the SGMM specification could include the contemporaneous lags of the independent variables as long as the set of instruments is valid and exogenous. In any case, the results of the SGMM estimation should be taken with caution as the numbers of instruments is higher than the groups.

25

Additional lags of the independent variables do not significantly change the results.

26

Boyd and Smith (2002) also note that if there are omitted variables in units’ estimates which are correlated with the observed covariates, these will lead to bias in these individuals’ estimates of the observed covariates. However, averaging estimates across units would cancel out biases. The issue, however, will persist if the omitted variable bias is structural across all units of the panel.

Macro-Financial Linkages and Heterogeneous Non-Performing Loans Projections: An Application to Ecuador
Author: Mr. Francesco Grigoli, Mr. Mario Mansilla, and Martín Saldías
  • View in gallery

    Variables’ Correlation and Dispersion

    (Percent, yoy, unless otherwise specified)

  • View in gallery

    Stationarity of VAR Series

  • View in gallery

    Macro-Financial Linkages

    (Quarterly responses to one-standard deviation shocks with two-standard-deviation confidence intervals)

  • View in gallery

    SVAR Conditional Forecasts

    (Forecast in percent, yoy, with one-standard-deviation confidence intervals)

  • View in gallery

    Relationship between NPL and Real GDP Growth

    (Percent)

  • View in gallery

    Effect of Real GDP Growth Shock on NPL

    (Impact of one pp real GDP growth shock evaluated at 2016Q1 NPL ratio)

  • View in gallery

    Weighted Average of NPL Projections

    (Percent)

  • View in gallery

    Heterogeneity in NPL Projections

    (Percent)