Abstract: This study employs the bootstrap autoregressive distributed lag (ARDL) approach alongside the dynamic ARDL simulations technique to investigate the non-linear effect of public debt on public expenditure in Nigeria during the 1981–2020 period. The result of the bootstrap bounds test illustrates the presence of a long-term relationship between public expenditure and public debt (along with oil rents, output growth and urbanisation). Further, the estimation results indicate that the effect of public debt on public expenditure is non-linear. In particular, public expenditure increases at early stages of rising public debt but declines at latter phases when public debt grows beyond specific threshold. This empirical outcome is further validated by the dynamic ARDL simulations approach which shows a significant decline in predicted public expenditure after short-term expansion due to counterfactual shock in public debt. Thus, policies which diversify public revenue from oil production and a reversal of the rising trend in public debt are recommended to avert the adverse welfare implications of declining public expenditure.
Keywords: public expenditure, public debt, bootstrap ARDL, dynamic ARDL simulations, Nigeria.
Articles
Non-Linear Effect of Government Debt on Public Expenditure in Nigeria: Insight from Bootstrap ARDL Procedure

Recepción: 15 Enero 2022
Aprobación: 24 Abril 2022
The role of public debt in bridging financial resource gap in developed and developing economies cannot be over-stressed (Ncanywa & Masoga, 2018). However, there seems to be disagreement amongst scholars on the precise effect of public debt on public expenditure. For instance, Buchanan and Wagner (1977) claimed that debt-financed deficit tends to produce higher levels of public expenditure. The authors’ argument rests on the premise that rising debt (and/or deficit) creates a counterfactual perception amongst citizens that the cost of public goods and services have become cheaper, thus, they will increase their demand for such goods and services (including transfers), leading to an expansion in public expenditure (Buchanan & Wagner, 1977; Iiyambo & Kaulihowa, 2020). A number of empirical studies appear to support this view (Abu et al., 2022; Del Monte & Pennacchio, 2020; Fagge & Ibrahim, 2018; Iiyambo & Kaulihowa, 2020; Ncanywa & Masoga, 2018; Ouattara, 2006; Ukwueze, 2015).
On the contrary, recent debates on the implication of growing public debt underline the deleterious effect of government debt on public spending. For instance, certain researchers have argued that huge debt burden tends to generate large diminution in public expenditure and/or adjustments in components of public expenditure such as spending on education and health, among others (Cashel-Cordo & Craig, 1990; Fosu, 2007, 2010; Mahdavi, 2004; Ncanywa & Masoga, 2018; Ndour, 2017; Picarelli et al., 2019; Shabbir & Yasin, 2015; Quattri & Fosu, 2012). To them, the fall in public expenditure is attributable, in part, to debt-induced liquidity constraints, the need to ensure continual servicing of accumulated debts, and/or because expenditure cuts are more quickly applicable rather than revenue generation when implementing fiscal austerity as a deficit reducing strategy (Fosu, 2007; Mahdavi, 2004; Picarelli et al., 2019; Shabbir & Yasin, 2015; Shonchoy, 2010; Taylor, 1993). Moreover, Lora and Olivera (2006) ascribed the decline in public expenditure to the disincentive which higher debt stock creates for further indebtedness.
Besides the conflicting views on the impact of debt on public expenditure, certain authors including Shonchoy (2010) and Ndour (2017) emphasised the potential effect of public debt on public expenditure to be non-linear. Whereas they seem to agree on the increasing effect of public debt on public expenditure, the authors believed that after a specific limit or threshold the continuous increase in the stock of public debt would lead to diminution in both productive and non-productive expenditure. This argument aligns with the “debt overhang” hypothesis of Krugman (1988) which asserts that public (investment) expenditure tends to decline following an increase in debt burden beyond certain threshold. Mahdavi (2004) attributed the decline in public expenditure, following a continuous increase in debt stock over and above a specific threshold, to shift in budgetary allocations away from some (or all) components of public expenditure in favour of the relative share of interest payments.
Like in most countries, public debt stock has been on the increase in Nigeria since she gained independence from the British due in part to the growing size of the country’s population, increased participation of the government in the economy and, most importantly, the dwindling oil revenue which accounts for over 70 percent of the government’s income (Central Bank of Nigeria [CBN], 2020). Available statistics demonstrates that Nigeria’s public debt stock maintained an expansionary stance between 1981 and 2020, except in 2006, when the country’s debt stock declined significantly by about 47.77 percent, following the implementation of the Paris Club debt relief (CBN, 2020). However, the sustained increase in public debt saw Nigeria emerging as the most indebted country in West Africa, and second in Sub-Saharan Africa, only behind South Africa (World Bank, 2021). In the same vein, Nigeria’s public expenditure has remained relatively high (Aladejare, 2019; CBN, 2020).

Source: Authors’ calculation based on data collected from CBN (2020)
Nonetheless, a cursory look at statistics from the CBN (2020) illustrates that the pattern of movement of the (unadjusted) public expenditure and public debt stock of the central government is not very clear, as they moved in the same direction in some years and opposite in other periods (Figure 1). For example, whereas public debt stock rose from N23.83 billion in 1982 to N32.80 billion in 1983, overall public expenditure declined from N11.9 billion to N9.60 billion during the same period. Similarly, public expenditure declined following increases in public debt stock in 1994, 2000, 2012 and 2014. Whereas rising public debt in the 1984–1993, 1995, 1997–1999, 2001–2004, 2007–2011, 2013, and 2015–2020 periods coincided with an expansion in public expenditure, the reduction in government debt in 1996, 2005 and 2006 saw public expenditure increase in the same period.
Given the unclear pattern in the movements of public debt and expenditure in Nigeria, it is possible that the public expenditure–debt relationship is non-linear. Therefore, the main objective of this study is to evaluate whether the public expenditure–debt relationship in Nigeria is non-linear. The present study is relevant and contributes to the literature in a number of ways. This study is the first attempt (to our knowledge) to investigate the non-linear effect of public debt on public expenditure in Nigeria. Second, the study is also the first to identify the threshold point in the public expenditure and public debt relationship in Nigeria. Third, the current study adopts the bootstrap ARDL-bounds testing method advanced by McNown et al. (2018) alongside the dynamic ARDL simulations procedure of Jordan and Philips (2018) to estimate the public expenditure–debt relation in Nigeria. These techniques are robust given that they address major shortcomings in the conventional ARDL approach of Pesaran et al. (2001). The various advantages of these techniques are discussed in subsequent section. Lastly, by exploring the effect of debt on public expenditure and identifying the threshold, this study is expected to advance practical insight into fiscal policy-making in Nigeria.
The remainder of the study is structured as follows. The second section discusses relevant empirical studies, while section three consists of theoretical framework, model formulation and econometric techniques. Section four contains presentation and discussion of results, and section five is for conclusion and recommendations.
Several studies have assessed the effect of public debt on public expenditure in both developing and developed countries. However, the empirical literature is typically polarised along two dimensions, with studies illustrating the positive effect of debt on expenditure on the one hand, and those that found a deleterious effect of debt on public expenditure on the other hand. For instance, using a sample of 34 Organization of Economic Cooperation and Development (OECD) countries, Del Monte and Pennacchio (2020) indicated that the impact of public debt burden on public expenditure is positive. In the same vein, Ouattara (2006) reported a positive effect of public debt (servicing) on non-developmental expenditure (i.e., wages, salaries and subsidies) in developing countries.
In contrast, Picarelli et al. (2019) confirmed that debt burden is detrimental to public investment in selected European Union (EU) countries. Similarly, a number of studies illustrate that public debt (servicing and burden) has an adverse impact on aggregate and/or components of public expenditure in developing countries in Sub-Saharan Africa, East Asia and Pacific, Europe and Central Asia, Latin America and Caribbean, Middle East and North Africa, and South Asia, including highly indebted poor countries (Cashel-Cordo & Craig, 1990; Fosu, 2007, 2008, 2010; Lora & Olivera, 2006; Mahdavi, 2004; Quattri & Fosu, 2012; Shabbir & Yasin, 2015; Shonchoy, 2010; Stephens, 2001).
At country-specific level, researchers have also assessed the impact of public debt on public expenditure, particularly in developing countries such as Bangladesh, Kenya, Namibia, Nigeria, South Africa, including Uganda, and found that public debt (stock, servicing) is positively related to aggregate and/or components of public expenditure (Abu et al., 2022; Ahaisibwe, 2018; Fagge & Ibrahim, 2018; Firoj et al., 2018; Iiyambo & Kaulihowa, 2020; Kiminyei, 2014, 2019; Ncanywa & Masoga, 2018; Omodero, 2019; Ukwueze, 2015). On the contrary, some authors discovered that public debt (stock and servicing) decelerates public expenditure (Adamu, 2016; Aregbeyen & Akpan, 2013). Moreover, other studies confirmed an adverse impact of external debt on public expenditure and a positive influence of domestic debt on public expenditure (Aladejana et al., 2021; Kanano, 2006; Omodero, 2019; Orhewere & Ogbeide-Osaretin, 2020).
Besides these conflicting empirical outcomes, it is apparent that researchers did not deem it important to investigate the non-linear effect of public debt on public expenditure. In addition, despite the marked increase in public expenditure and debt, and the unclear pattern in their movements in Nigeria, there is obvious dearth of empirical studies which explore whether the public expenditure–debt relationship in the country is non-linear. Therefore, the present study contributes to the literature by assessing the non-linear impact of public debt on public expenditure in Nigeria during the 1981–2020 period. The study also considers the threshold of the public expenditure–debt relationship. Further, robust estimation techniques such as the bootstrap ARDL-bounds testing approach and the dynamic ARDL simulations method are adopted to estimate the relationship between the variables.
Theoretical Framework and Model Specification
The literature on the effect of public debt on public expenditure is characterised by conflicting opinions. On the one hand, scholars including Buchanan and Wagner (1977) argued that public debt generates higher levels of public expenditure. On the contrary, a number of studies illustrate that the impact of public debt on public expenditure is deleterious (Fosu, 2007; Lora & Olivera, 2006; Mahdavi, 2004; Picarelli et al., 2019; Shonchoy, 2010; Taylor, 1993).
Furthermore, certain authors suggested that the effect of public debt on public expenditure can be non-linear (Ndour, 2017; Shonchoy, 2010). In line with this view, it is argued that rather than a direct/linear effect, the influence of public debt on public expenditure is positive at the initial stage of public debt growth and then negative after public debt grows beyond a specific threshold. Interestingly, Ndour (2017) included the interaction term of public debt in the public expenditure model to ascertain the non-linearity and threshold of the relationship. In other words, public expenditure can be modelled to incorporate both public debt and public debt squared.
In addition, evidence has shown that public expenditure can be influenced by a number of macroeconomic factors. For instance, since proceeds from the sales of oil accounts for over 70 percent of income of the government, the role of oil rents (OIL) in impacting public expenditure level in Nigeria cannot be over-stressed. Besides, authors have advocated the inclusion of oil rents when modelling economic variables because the Nigerian economy depends on events in the oil and gas industry (Abu, 2017; Abu & Gamal, 2020; Abu & Karim, 2021; Abu & Staniewski, 2019). Equally, it is well documented in the literature that output growth (Y) is a potential determinant of public expenditure. Interestingly, recent studies on public expenditure dynamics in Nigeria illustrated that real output growth rate is a relevant determinant of public expenditure in the country (Aladejare, 2013, 2019, 2020; Ukwueze, 2015). Additionally, the role of urbanisation (URBAN) in public expenditure has been advanced in the literature (Akanbi, 2014; Akanbi & Schoeman, 2010; Aregbeyen & Akpan, 2013; Firoj et al., 2018; Shonchoy, 2010). Typically, the greater the concentration of a country’s population in urban centres, the higher the demand for public goods and services (such as education, roads and transportation), leading to rising public spending on infrastructure and public utilities (Shonchoy, 2010).
Therefore, a simple model in which public expenditure (PEX) is dependent on public debt (DEBT), its quadratic term (DEBT2), alongside oil rents (OIL), real output growth rate (Y), and the rate of urbanisation (URBAN) is specified as follows:

where c, αi, and μ represent the intercept, coefficient of the regressors, and the stochastic disturbance term with zero mean and constant variance, respectively.
To explore the co-integrating relationship between the variables, the bootstrap ARDL-bounds testing to co-integration proposed by McNown et al. (2018) is employed. This approach is adopted due to the number of advantages which it has over the traditional ARDL framework of Pesaran et al. (2001). For example, the approach addresses the issue of weak size and power properties that characterise the traditional ARDL bounds-testing. Based on Monte Carlo simulations, McNown et al. (2018) demonstrated that the bootstrap ARDL procedure performs better than the asymptotic test in the traditional ARDL technique with regard to size and power properties. In other words, besides allowing for endogeneity and feedback that may exist in the system, the bootstrap co-integration procedure is shown to be more robust and less biased in comparison with the asymptotic test in the traditional ARDL bounds testing. Moreover, by incorporating an additional test on the lagged level(s) of the independent variable(s) to complement the existing F- and t-tests presented in the traditional ARDL bounds-testing framework, the approach is superior to the conventional ARDL because it increases the power of the F-test, and by extension provides a better insight on the co-integration status of the system. The bootstrap ARDL method is also preferred to the traditional ARDL co-integration approach due to its ability to eliminate the issue of inconclusive inferences which may arise when using the traditional ARDL procedure (Goh et al., 2017; McNown et al., 2018).
To perform the bootstrap ARDL bounds-testing, an unrestricted error correction ARDL(p, q1, q2, q3, q4, q5) model to be estimated is specified as follows:

where ∆ represents the difference operator; c is the constant term; δ and β are the coefficients of the lagged dependent variables and the regressors; λi is the vector of coefficients of the lagged levels of dependent and independent variables; t = 1…T denotes time; ε is the stochastic error term with zero mean and constant variance.
According to McNown et al. (2018), it is required that the following three null hypotheses are rejected to establish co-integration amongst the series:

If all the three null hypotheses are rejected (i.e., if the three statistics exceed the respective critical values at specific level of significance), it is concluded that co-integration exists. Else, if all tests except the test on the lagged independent variable(s) are significant, the degenerate case #1 (degenerate lagged independent variable(s) case) occurs. Here the lag of the dependent variable is solely responsible for the joint significance of the error correction terms. However, if the t-test on the lagged level of the dependent variable is insignificant whereas the other two are significant, degenerate case #2 (degenerate lagged dependent variable case) is established. Hence, all three null hypotheses must be rejected for valid conclusion on co-integration to be reached as either of the cases indicates the case of no co-integration. The critical values for the three tests are generated through the bootstrap procedure as outlined in McNown et al. (2018). The major benefit of the bootstrap generated critical values is its elimination of the likelihood of inconclusive or indeterminate inferences from the ARDL bounds-testing, which is a likely occurrence (McNown et al., 2018).
In addition to the bootstrap ARDL cointegration method, the dynamic ARDL simulations approach proposed by Jordan and Philips (2018) is employed to examine the response of public expenditure following a counterfactual change/shock in public debt. One of the main advantages of the dynamic ARDL simulations procedure over the traditional ARDL framework is its ability to better interpret the significance of short-term and long-term effects of ARDL models, which becomes difficult as the lag structure gets more complex. The unique function of the procedure is to simulate and automatically visualise the impact of a counterfactual change in one weakly exogenous regressor at a single point in time using stochastic simulation techniques, while holding all else equal (Khan et al., 2020). The approach has been utilised in recent studies to predict future outcomes in social, economic and climate variables following a counterfactual shock (Ali et al., 2021; Khan et al., 2020; Khan et al., 2021; Olasehinde-Williams & Oshodi, 2021; Sarkodie et al., 2020; Sarkodie & Owusu, 2020, 2021; Shabbir et al., 2020). For the purpose of the present study, the simulation graph would be generated using the error correction algorithm of the dynamic ARDL procedure through 5,000 simulations of the public debt parameter from the unrestricted ARDL equation specified in equation (2).
The present study uses Nigeria’s annual time-series dataset covering 1981 to 2020. The choice of the period is guided by data availability. The data are measured as follows. PEX is measured as a ratio of central government aggregate (recurrent and capital) expenditure to the GDP. Whereas researchers are often faced with a dilemma of either using the absolute value of public expenditure or its ratio to the GDP, the measure of public expenditure as a ratio of the GDP has come to enjoy widespread acceptance because it reflects the participation of the government in the economy (García & Martos, 2004). Besides, this measure of public expenditure was employed in a number of empirical studies (Akanbi & Schoeman, 2010; Cockx & Francken, 2015; Jibir & Aluthge, 2019; Mehmood & Sadiq, 2010; Picarelli et al., 2019; Shonchoy, 2010; Ukwueze, 2015). In the same vein, DEBT is measured as the share of total public debt of the federal government to the GDP (as used in Del Monte & Pennacchio, 2020; Picarelli et al., 2019; Shabbir & Yasin, 2015). We used this broader measure of public debt, rather than total debt of the central government, because it reflects both the extent of indebtedness of the government and the capacity of a country to pay back its outstanding debts.
Moreover, OIL is captured by the ratio of oil rents to the GDP. Oil rents (the difference between the value of crude oil production at regional prices and total costs of production) indicate the actual amount that the central government realises annually from the sales of crude oil after making provision for the cost of production (World Bank, 2021). Y is the annual percentage growth rate of real GDP, and it reflects the degree of expansion and/or contraction of economic activities in a country. Lastly, URBAN is measured as the percentage growth rate of urban population. The data on the variables were collected from various sources. Specifically, the data on aggregate public expenditure and public debt of the central government, and GDP were sourced from the Central Bank of Nigeria (CBN) annual statistical bulletin, while oil rents, real GDP growth rate and urban population growth rate were collected from the World Bank’s World Development Indicators (WDI).
The computed summary statistics of the variables are presented in Table 1. The results demonstrate that average public expenditure (as share of the GDP) between 1981 and 2020 is 8.44 percent, while the ratio of public debt to the GDP is 31.85 percent during the same period. Additionally, the results indicate that average oil rents (as share of GDP), real output growth rate and urban population growth rate are 11.83 percent, 3.03 percent and 4.73 percent, respectively.

Source: Authors’ computation using Eviews 12.
PEX, DEBT, OIL, Y and URBAN denote public expenditure, public debt, oil rents, real output growth, and urbanisation, respectively.Prior to modelling the relationship between the variables, unit root test was conducted to ascertain the stationarity property of the variables. The traditional Augmented Dickey-Fuller (ADF), Philips-Perron (PP), and Kwiatkowski-Phillips-Schmidt-Shin (KPSS) unit root tests, alongside the Zivot-Andrews (ZA) test of Zivot and Andrews (1992) are employed for this purpose. Since ZA test explicitly accounts for the possibility of structural break in series, the test is adopted as a complement to the conventional tests which are weak in the presence of structural breaks.
The results of the tests summarised in Table 2 present mixed outcomes. However, the test statistics generally confirm that OIL, Y and URBAN are stationary at level (i.e., they have no unit root), while PE X, and DEBT are only stationary after taking their first difference. Thus, the results demonstrate that the variables are a mixture of I(0) and I(1). Since the bounds-testing procedure allows for series to have different orders of integration (at least not greater than order one – I(1)), these findings provide the justification for adopting the bootstrap ARDL bounds-testing approach (Pradhan & Hiremath, 2020; Pradhan et al., 2021).

Source: Authors’ computation using Eviews 12.
ADF’s critical values for intercept are given as: -3.62, -2.94 and -2.61, at 1%, 5% and 10% levels, respectively. PP’s critical values are -3.61 (1%), -2.94 (5%), and -2.61 (10%).KPSS’s critical values are 0.74, 0.46, and 0.35, at 1%, 5% and 10% levels, respectively. Zivot and Andrews’ (1992) critical values for structural change in the level shift are -5.34 (1%), -4.93 (5%) and -4.58 (10%). ADF, PP andZA est the null of a unit root against the (trend-) stationary (with a one-time break, in the case of ZA) alternative. In contrast, KPSS tests the null of stationary against unit root alternative. The optimal lag length selection in ADF and ZA is based on Akaike’s (1979) Information Criteria (AIC), while Bartlett kernel and the Newey-West methods are followed using PP and KPSS process. Asterisks (***), (**) and (*) indicate significance at 1%, 5% and 10% level, respectively.To determine the co-integrating relationship between the series, the bootstrap ARDL bounds testing proposed by McNown et al. (2018) is employed, and the result reported in Table 3. The result reveal that the values of the overall F-statistic on lagged level variables (4.65), t-test on the lagged level dependent variable (-4.42) and F-statistic on lagged level independent variables (5.34) all exceed the critical values generated through the bootstrap procedure proposed by McNown et al. (2018) at the 5 percent level. Therefore, there is sufficient evidence to reject the null hypothesis of no co-integrating relationship between public expenditure and public debt (including public debt squared, oil rents, output growth and urbanisation).

Source: Authors’ computation using Eviews 12.
Asterisk (**) denotes significance at the 5% level based on critical values generated from the bootstrap procedure (with 2000 replications) of McNown et al. (2018). F1 represents the F-statistic for the lagged level variables, F2 denotes the F-statistic for the lagged level of the independent variables, and t is the t-statistic for the lagged level of the dependent variable. The optimal lag-length is suggested by AIC.Following the confirmation of cointegrating relationship amongst the series, the results of the estimated short-run and long-run model of the selected ARDL model with optimal lag length (2,1,1,1,1,1) suggested by AIC, are summarised in Panel A and Panel B of Table 4, respectively.
The results demonstrate that public debt has a significant positive effect on public expenditure both in the short- and long-run at the 1 percent level. A percent increase in public debt raises public expenditure by 0.179 percent and by 0.355 percent in the short- and long-run, respectively. The positive linear effect of public debt on public expenditure is consistent with previous findings (Abu et al., 2022; Ahaisibwe, 2018; Del Monte & Pennacchio, 2020; Fagge & Ibrahim, 2018; Iiyambo & Kaulihowa, 2020; Ncanywa & Masoga, 2018; Ouattara, 2006; Ukwueze, 2015). In contrast, the coefficient of squared public debt is statistically significant and negatively related to public expenditure both in the short- and long-term, at the 5 percent level. A percent increase in public debt squared causes public expenditure to decline by 0.002 percent in the long- and short-term, respectively, at the 5 percent level.
The implication of this outcome is that the effect of public debt on public expenditure in Nigeria is non-linear. In other words, there is the tendency of public expenditure rising in tandem with increases in public debt until it reaches a specific limit. However, further increase in public debt beyond this specific threshold leads to a decline in public expenditure. Interestingly, the non-linear impact of debt on public expenditure validates the postulations of Shonchoy (2010) and Ndour (2017). Although public expenditure in Nigeria has maintained an expansionary stance, its share in productive sectors of the economy (such as investment, education, health, and roads) has continued to decline, with a larger proportion of budgetary allocation devoted to servicing accumulated debt.

Source: Authors’ computation using Eviews 12.
Note The model is estimated by setting the maximum lag to 4, and the optimum lag-length is suggested by AIC. ∆is the first difference operator. Asterisks (***), (**) and (*) denote significance at 1%, 5%, and 10% respectively. Values in parentheses (.) in panel A and B are the t-ratio, and values in square parentheses [.] in panel C are the probability values of the LM test statistics. Importar imagen denote the Breusch-Godfrey serial correlation, Breusch-Pagan-Godfrey heteroscedasticity, Jarque-Bera normality and Ramsey RESET functional form LM tests, respectively.Moreover, the findings suggest that public expenditure and public debt relationship in Nigeria has an inverted U-shaped feature. Based on the quadratic specification in Equation 1, the estimated turning points/threshold of public debt and public expenditure are 44.75 percent and 11.43 percent, respectively. A graphical representation of the quadratic function is provided in Figure 2. The plots illustrate the increase in public expenditure following an increase in public debt prior to its threshold of 11.43 percent. However, over and above this threshold, further expansion in public debt (as a percentage of the GDP) leads to diminution in public expenditure (as a share of GDP).
Regarding other variables, the results indicate that oil rents is positively related to public expenditure in the short- and long-term, and significant at 1 percent level. A percent increase in oil rents raises public expenditure (as a share of the GDP) by 0.313 percent and 0.162 percent in the long- and short-term, respectively. The positive relationship between oil rents and public expenditure is consistent with the outcome of past studies (Abu et al., 2022; Adedokun, 2018; Aladejare, 2020; Aregbeyen & Fasanya, 2017; Jibir & Aluthge, 2019). Thus, an increase in oil rents raises government’s income, enhances its capacity to embark on developmental projects and investments, leading to an expansion in the level of public expenditure.

Source: Authors’ computation based on long-run estimates in Table 4
Further, the results reveal that output growth and urbanisation are positively related to public expenditure in the long term at the 10 percent and 5 percent level, respectively. A percent increase in output growth raises the level of public expenditure in the long term by 0.054 percent. The positive effect of output growth on public expenditure validates Wagner’s law of increasing state activities in Nigeria and recent empirical studies in the country (Aladejare, 2013, 2019; Ukwueze, 2015). The outcome suggests that the greater the growth of the economy, the larger would be the expansion in public expenditure in the country.
In addition, the results demonstrate that an increase in urban population growth rate by a percent leads to an increase in public expenditure by 1.33 percent. This outcome aligns with the discovery of a number of studies (Akanbi, 2014; Akanbi & Schoeman, 2010; Aregbeyen & Akpan, 2013; Firoj et al., 2018; Shonchoy, 2010). Since rising urbanisation is often accompanied by higher demand for public goods and services such as infrastructures, healthcare and education, etc., it will push up public expenditure.
Furthermore, the coefficient of error correction term lagged by one period (ECTt–1) is significant, correctly signed, and suggests that about 87 percent of the deviation in the short term is corrected within a year.
The results of the diagnostic tests presented in panel C of Table 4 demonstrate that the ARDL model is free from the issues of serial-correlation, misspecification error or bias, heteroscedasticity and normality. More so, the cumulative sum of recursive residuals (CUSUM) and cumulative sum of squares of recursive residuals (CUSUMQ) plots of Brown, Durbin and Evans (1975) presented in Figure 3 and Figure 4 indicate the stability of the parameters in the long term, and hence the suitability of the empirical outcome for policy making.

Following the evidence of cointegration and the estimation results using the bootstrap ARDL method, the dynamic ARDL simulations method was employed to ascertain the response of public expenditure following a counterfactual change (or shock) in public debt over a 50-year period (i.e. from 2020 to 2070).
Considering a 5 percent counterfactual change in public debt on public expenditure over the 50-year period, the corresponding spikeplot of the dynamically simulated ARDL is presented in Figure 5. The spikeplot illustrates that the immediate impact of a counterfactual 5 percentage point increase in public debt in the tenth period is a significant increase in public expenditure by over 10 percent point in 2030 from about 8 percent point in pre-shock periods between 2020 and 2029. Over time and in subsequent periods, especially immediately – a period after the hypothetical shock in the tenth period, the predicted public expenditure declines to about 9 percent point. However, between 2032 and 2070, public expenditure stabilises within the 9 percent point region with little changes. In essence, this suggests the decline in the ratio of public expenditure to the GDP in the long term following significant increases in public debt (as a share of the GDP), and thus in consonance with the ARDL estimates.

Note. The forecasting is executed using the Dynamic ARDL Simulations. The counterfactual shock in predicted public expenditure is achieved based on 5,000 simulations. Black dots show average predicted value, while shaded lines (from darkest to lightest) show the 75, 90, and 95 percentiles of the predictions from the simulations akin to a credible interval or confidence interval.
Source: Authors’ computation using ‘dynardl’ program in Stata
This study assesses the non-linear effect of public debt on public expenditure in Nigeria following an expansion in both variables during the 1981–2020 period. Employing the bootstrap ARDL bounds testing approach and the dynamic ARDL simulations procedure, the cointegration test result illustrates the presence of a long-term relationship between public expenditure and public debt (along with oil rents, real output growth, and urbanisation). The estimation results confirm the non-linear effect of public debt on public expenditure. More so, the results indicate that the public expenditure and public debt relationship has an inverted U-shaped feature. Public expenditure increases at early stages of rising public debt, but it declines at the latter phases of rising public debt. Based on the quadratic specification, a threshold of 44.75 percent and 11.43 percent is attributed to public debt and public expenditure, respectively. Further, the dynamic ARDL simulations results reveal a significant increase in public expenditure in the short term following a counterfactual debt shock and a gradual decline in public expenditure in subsequent periods.
In the light of these findings, the present study recommends that the Nigerian government reduce the debt stock to avert the long-term consequence of falling public expenditure. This measure is important as recent increase in public expenditure in Nigeria is on account of servicing accumulated debt rather than raising productive components of public expenditure in the country. It is also recommended that government’s revenue source be diversified away from oil due to the volatile nature of its price and its non-sustainability for long-term development plans.

Source: Authors’ calculation based on data collected from CBN (2020)

Source: Authors’ computation using Eviews 12.
PEX, DEBT, OIL, Y and URBAN denote public expenditure, public debt, oil rents, real output growth, and urbanisation, respectively.
Source: Authors’ computation using Eviews 12.
ADF’s critical values for intercept are given as: -3.62, -2.94 and -2.61, at 1%, 5% and 10% levels, respectively. PP’s critical values are -3.61 (1%), -2.94 (5%), and -2.61 (10%).KPSS’s critical values are 0.74, 0.46, and 0.35, at 1%, 5% and 10% levels, respectively. Zivot and Andrews’ (1992) critical values for structural change in the level shift are -5.34 (1%), -4.93 (5%) and -4.58 (10%). ADF, PP andZA est the null of a unit root against the (trend-) stationary (with a one-time break, in the case of ZA) alternative. In contrast, KPSS tests the null of stationary against unit root alternative. The optimal lag length selection in ADF and ZA is based on Akaike’s (1979) Information Criteria (AIC), while Bartlett kernel and the Newey-West methods are followed using PP and KPSS process. Asterisks (***), (**) and (*) indicate significance at 1%, 5% and 10% level, respectively.
Source: Authors’ computation using Eviews 12.
Asterisk (**) denotes significance at the 5% level based on critical values generated from the bootstrap procedure (with 2000 replications) of McNown et al. (2018). F1 represents the F-statistic for the lagged level variables, F2 denotes the F-statistic for the lagged level of the independent variables, and t is the t-statistic for the lagged level of the dependent variable. The optimal lag-length is suggested by AIC.
Source: Authors’ computation using Eviews 12.
Note The model is estimated by setting the maximum lag to 4, and the optimum lag-length is suggested by AIC. ∆is the first difference operator. Asterisks (***), (**) and (*) denote significance at 1%, 5%, and 10% respectively. Values in parentheses (.) in panel A and B are the t-ratio, and values in square parentheses [.] in panel C are the probability values of the LM test statistics. Importar imagen denote the Breusch-Godfrey serial correlation, Breusch-Pagan-Godfrey heteroscedasticity, Jarque-Bera normality and Ramsey RESET functional form LM tests, respectively.
Source: Authors’ computation based on long-run estimates in Table 4


Note. The forecasting is executed using the Dynamic ARDL Simulations. The counterfactual shock in predicted public expenditure is achieved based on 5,000 simulations. Black dots show average predicted value, while shaded lines (from darkest to lightest) show the 75, 90, and 95 percentiles of the predictions from the simulations akin to a credible interval or confidence interval.
Source: Authors’ computation using ‘dynardl’ program in Stata