Open access peer-reviewed chapter

Bayesian Graphical Model Application for Monetary Policy and Macroeconomic Performance in Nigeria

By David Oluseun Olayungbo

Submitted: June 28th 2018Reviewed: June 12th 2019Published: July 17th 2019

DOI: 10.5772/intechopen.87994

Downloaded: 105

Abstract

This study applies Bayesian graphical networks (BGN) using Bayesian graphical vector autoregressive (BGVAR) model with efficient Markov chain Monte Carlo (MCMC) Metropolis-Hastings (M-H) sampling algorithm in a dynamic interaction among monetary policies and macroeconomic performances in Nigeria for the period of 1986Q1–2017Q4. The motivation stems from the instability in the movement of exchange rate, inflation rate and interest rate in Nigeria over the past years as a result of the structure of the economy. In this way, the monetary authority periodically applies the various policy instruments to stabilize the economy using reserve and money supply as at when due. This study adapts VAR and SVAR structure to examine the dynamic interaction among variables of interest, using BN, to provide a better understanding of the monetary policy dynamics and fit the changing structure of the Nigeria’s economy as regards the dynamics in her economic structure. Our results show that inflation is the strong predictor of interest rate in Nigeria. A monetary policy of broad inflation targeting is recommended for the country.

Keywords

  • Bayesian graphical networks
  • SVAR
  • MCMC
  • M-H
  • Granger-causal inference
  • Nigeria

1. Introduction

A network can be described as a set of items, with nodes or vertices that are related by edges or links for a specific purpose. There are different types of networks. These networks can be social, economic, informational, technological, biological and so on. However, in this paper, we are interested in Bayesian graphical network (BGN) to investigate causal inferences among monetary policies and macroeconomic performances in Nigeria. Causal effects have been used in economic literature starting from Granger [1], Engle and Granger [2] and Sims [3] using vector autoregression (VAR) to Amisano and Giannini [4] and Blanchard and Quah [5] using structural vector autoregression (SVAR). The progression from VAR to SVAR has been a result of over-parameterization and identification problem associated with VAR that have limited its use for forecasting [6]. SVAR, in a way, has been able to overcome these problems through the application of the recursive and the non-recursive structural model [5]. The advantage of the BGN is that there are directed acyclic graphs (DAG) with nodes or vertices called variables having edges that indicate structural dependence among the variables of interest. Studies that have used graphical models to causal relationships are Pearl [7], Spirites et al. [8], Demiralp and Hoover [9] and Moneta [10], among others. Our intention is to combine BGN with SVAR, in order to examine interrelationship among monetary policies and macroeconomic performances in Nigeria. However, many previous studies have worked on causal inference using BGN-SVAR in developed countries, studies like Swanson and Granger [11], Demiralp and Hoover [9], Moneta and Spirites [12], Corander and Villani [13] and more recently Ahelegbey et al. [6]. The advantage of the BGN-SVAR is that causal influences are tracked down among variables of interest either instantaneously or with time lags using conditional probabilistic inference on the structural model.

Importantly, in this present study, we adopted the BGN-SVAR approach used by Ahelegbey et al. [6] and deal with the identification structure to derive the Bayesian networks using DAG. The Bayesian structural model is then simulated using the multi-move Markov chain Monte Carlo Metropolis-Hastings sampling algorithm. The analyses are done with a view to examine the causal relationship among monetary policy actions and macroeconomic performances in Nigeria. The BGN-SVAR method is more superior to the usual standard Granger causality test, thereby providing an important tool for policy implications, especially for an emerging country like Nigeria where monetary policy stance is a major factor to the macroeconomic performance of the economy. The paper is as follows. Section 2 gives the data source, variable definition and the descriptive statistics, and Section 3 provides an overview of the BGN-SVAR model. While Section 4 highlights the empirical discussion, Section 5 concludes and states the policy recommendation.

2. Data source and variable definition

All the data used for this study were sourced from the Statistical Bulletin [14] published by the Central Bank of Nigeria. The data range from the first quarter of 1986 to the fourth quarter of 2017. The choice of the scope of study is strictly informed by data availability. Monetary policy variable is measured using broad money supply normally referred to as M2, which is defined as the sum of currency in circulation, demand deposit and time deposit in the banking sector. The industrial output (Ind) comprises of crude oil petroleum, natural gas, solid minerals, coal mining, metal ores, quarrying, mining, manufacturing, oil refining and cement production. In addition, interest rate is measured with lending rate (Intr) on banks’ credit to the public. Inflation rate (Infl) is measured by the average price index of consumer goods over the period of study. Exchange rate (Exch), on the other hand, is the rate of change of the local currency, naira, to the United State (US) dollar. The industrial output and money supply are measured in the local currency, naira, while the exchange rate is measured in US dollars.

2.1. Descriptive statistics of variables

The descriptive statistics as presented in Table 1 shows industrial output, exchange rate, inflation rate, interest rate and money supply variables in their unit form with 124 observations for the period of study. The difference between the average mean value and the maximum value of exchange rate, inflation rate and interest shows the high volatility of these variables. The average values for exchange rate, inflation rate and interest rate are 93.86, 19.25 and 14.1%, respectively, while the maximum values are 304.72, 73.1 and 26.7%. The differences between the average values and the minimum values also support the volatility behaviour of the variables. In addition, the volatility movement in these variables as shown in Figure 1 equally confirms their fluctuation. The skewness of all the variables implies positive skewness of the data distribution. Finally, the significance of the probability of the Jarque-Bera at 5% indicates that all the variables reject the acceptance of the null hypothesis of normal distribution except for the exchange rate variable.

Figure 1.

The volatility behaviour of exchange rate, inflation rate and interest rate.

2.2. Unit root test

In order to ascertain the order of integration of our variables, we used both augmented [15, 16] tests. The results from Table 2 show that all the variables are nonstationary using ADF test. The PP, on the other hand, shows all the variables to be nonstationary except industrial output. Given the stationarity of industrial output with PP, we further conducted structural break following Perron [17] in that its presence can bias the unit root result. The structural break results show that all the variables are nonstationary and support the ADF test. The structural break result is not presented here but available upon request. Given the nonstationarity of all the variables, we proceeded to examine their cointegrating relationship over the period of study.

2.3. Cointegration test

We followed Johansen [18] cointegration technique that compares the trace and the eigenvalue with their critical values for the rejection or acceptance of the null hypothesis of no cointegration. The optimal lag lenght of 1 was chosen following Schwarz Criterion’s (SC) result in Table A4 at the Appendix. The cointegration results presented in Table 3 show that the trace statistics is greater than the critical value at 5% significance level at r=0. This hypothesis testing takes us to the next cointegrating vector, r1, where the trace statistics is less than the critical value. We therefore conclude that there is a long-run relationship among the variables.

3. Bayesian graphical VAR model

The starting point of the BGVAR is from VAR process proposed by Sim [3] as dynamic endogenous variables specified as

Yt=B0+B1Yt1+B2Yt2++BρYtρ+εtE1

Eq. (1) is a vector autoregressive process of order ρ, and it can be respecified in a reduced form as

Yt=B0+i=1ρBiYti+εtE2

Eq. (2) can further be specified in SVAR process following Amisano and Giannini [4] and Blanchard and Quah [5] with a lower triangular matrix Aas

AYt=B0+i=1ρBiYti+εtE3

Eq. (3) can be written in an inverted form as

Yt=β0+i=1ρβiYti+utE4

where β0=A1B0, βi=A1Biand 1iρare lags of the parameter matrices and ut=A1Σεtwhere Σis a diagonal matrix of variance and covariance matrix. Assume in econometric term that Yti=Xti'; then Eq.(4) can be expressed as

Yt=β0+i=1ρβiXti+utE5

Eq. (5) can be written in matrix form as

Y=β'X+UE6

The solution to the SVAR model can be achieved through the parameter identification by placing restrictions on the lower triangular matrix Aor the Bdiagonal matrix following the relevant economic theories for the underlying variables. The general method of solving the structural dynamics of the SVAR model after placing the necessary restrictions to attain identification is to determine the effects of the shocks on the contemporaneous variables through the impulse response functions. The impulse response function can be represented through the diagonal matrix of the covariance matrix as

Σ=A1ΣεtA1E7

where Σεt=εtεt'and the covariance matrix are assumed to be an identity matrix. The Cholesky decomposition where all the elements above the diagonal are zero (lower triangular matrix) is the usual way of solving the identification problem in SVAR to identify the structural shocks. The system thus becomes exactly identified by comparing the known elements, n2+n2, and the unknown elements, n2n+n, of the covariance matrix [5]. Interestingly, the graphical model can be represented in the form of SVAR following Ahelegbey et al. [6] exposition. A graphical model can be described as the representation of the conditional relationships among random variables. Graphical models are generated in nodes and edges. The nodes house the variables, while the edges point to their relationships. A bivariate graphical model can be written as XY,meaning Xcauses Y,where Y(child) is the dependent variable and X(parent) is the independent variable. In a multivariate setting, the graphical model can be expressed as XYZand can be interpreted to mean that the relationship between Xand Zis conditional on variable Y.Assume Yt=Yt11Yt22Ytnnwhere Ytiis a realization of the ith variable at time t. The graphical network of a DAG can be written in the form of Eq. (1) as YtsJYtiwhere Bρ,ij0and 0sρ. This implies that past value of YtsJat time lag s, with s<t,causes the future value of Yti. This explains the notion that cause precedes effect in time. Therefore, following past studies such as Corander and Villani [13] and Ahelegbey et al. [6] among others, a network structure of a DAG is described as G=VA, where Vis a finite set of nodes symbolizing Yt=Yt11Yt22Ytnnin this case, while Ais also a finite set of directed edges denoting YtsJYtistated earlier. In other words, the graphical model can be specified in VAR representation as Bs=GSφs,where sis a period lag, Bis the structural parameters of the interdependent variable Yt, Gsis the binary connectivity matrix and φsis a matrix of coefficient of lag s.At period sfor 0sρ, Gs,ij=1implies a causal effect of YtsJYti, and Gs,ij=0implies no causal relationship between YtsJand Ytiwith φs,ijindicating the quantity of the causal effects of YtsJYti.The graphical model can further be written in a VAR form from Eq. (2) as.

Yt=G0φ0+i=1ρGiφiYti+εtE8

Following Koop [19], Geweke [20] and Olayungbo [21], the posterior distribution in Bayes’ theorem can be written in continuous form as

PθY=PYθPθPYθPθE9

where PYθis the likelihood function and Pθis the prior. In proportionality form, Eq. (9) becomes

PθYPYθPθE10

Given the parameters to be estimated in our models, the posterior becomes

PμΣy1GYPYμΣy1GPΣy1GE11

where μ=0and the likelihood function is generated with respect to Eq. (8) as

PYΣy1G=2πnT2Σy1T2exp12(Σy1t=1TYtYt'E12

The prior density, PG, is chosen from a uniform prior as PG1; the inverse of the variance-covariance matrix of the error term follows a Wishart distribution, i.e., Σy1W(ν,t=1TY0Y0'1; and βis set to 0.5.The prior density is then written as

PΣy1G=1Kn(ν,Y0Y0'Σy1νn12exp12(Σy1t=1TY0Y0'E13

The posterior distribution is written with the likelihood function and prior density in Eqs. (12) and (13) as.

PμΣy1Y2πnT2Σy1T2exp12(Σy1t=1TYtYt'.1Kn(ν,Y0Y0'Σy1νn12exp12(Σy1t=1TY0Y0'E14

Eq. (14) which is the posterior distribution with the likelihood function and the prior density in Eqs. (12) and (13) is estimated with Markov chain Monte Carlo sampling methods and Metropolis-Hastings recursively to obtain the posterior means. Samples are drawn from the posterior distribution, PΣy1GY, given PYΣy1Gand PΣy1Gby using the MCMC and M-H algorithm (see Ahelegbey et al. [6] for more exposition).

3.1. Granger causality test

Following Granger [1], a pairwise Granger causal (P-GC) relationship that conditions a variable on other variables and their time lags is investigated. The reality about the P-GC causality in the graphical network analysis is that it is a directed and forward causal relationship among the dependence variable structures. The P-GC VAR ρmodel is stated as

Yti=s=1ραsYtsi+s=0ρβsYtsjE15

4. Empirical analysis and discussion

Following Ahelegbey et al. [6], we sampled and derive separately the multivariate autoregressive (MAR) and the multivariate instantaneous (MIN) system using M-H algorithm specified in Appendix 2. The MAR network is a contemporaneous interaction of the variables from their past realizations to current realizations, while the MIN network is a contemporaneous interaction of the variables in their current realizations. At the implementation level, we ordered our variables in the Ytvector as Yt=(M2, Ind, Infl, Intr, Exch) from Eq. (2). It should be noted that our results are not sensitive to the variable ordering, in which case, any order can be taken.

After choosing an optimal lag length of 1 using Akaike information criterion (AIC) and Schwarz criterion (SC) (see Appendix 1), the BGN-VAR model ran 20,000 Gibbs iteration each for the MAR and MIN, making a total number of iteration to be 40,000, out of which 20,000 was set as burn-ins to achieve convergence. The results of the P-GC VAR for the MIN and MAR are presented in Figures 2 and 3. The dark green (light) color implies strong (weak) dependence between the dependent and independent variables. The row variables are the independent or explanatory variables, while the column variables are the dependent or response variables. The result of the MIN in Figure 2 shows strong evidence of causal relationship with M2Infl, IndInfl, InflIntr. The results can be interpreted to mean a stronger evidence of effects of both money supply and industrial output on inflation rate in Nigeria. This implies that the increase of money in circulation increases industrial demand and production; however, the resultant effects lead to inflationary rate due to high cost of production in Nigeria. High cost of production for firms and companies generates increase in prices of goods and services due to the use of generators to produce. The country generates 5000 megawatts (MW) of electricity which is not enough in residential needs let alone industrial needs; hence, most firms result to the use of generating set which results to high cost of production and high prices. We also found stronger effects of inflation rate behaviour on interest rate movement in Nigeria. This implies that instantaneous change in money supply and industrial output are major determinants of inflation rate in Nigeria. We also found strong effects of inflation rate on interest rate. This implies that resultant inflationary effects compel monetary authority and commercial banks to choose high interest rate to stabilize and ensure reasonable returns on investment. Furthermore, Figure 3 gives the contemporaneous autoregressive structure of the variable of interest. The MAR result shows the causal edges as, firstly, M2t1M2, Ind and Infl; secondly, Indt1Exch; thirdly, Inflt1Intr, Infl and Ind; fourthly, Intrt1Intr; and, lastly, Excht1Exch. The causal edges can be interpreted to mean that inflation rate, industrial output and money supply respond strongly to immediate past lag of money supply. Furthermore, we found past value of industrial output to Granger-cause exchange rate. In addition, past value of inflation Granger causes interest rate, current inflation and industrial output. This outcome corroborates the MIN result that inflation rate is a strong predictor of interest rate in Nigeria. The fourth and last causal edges can be explained to mean both past interest rate and exchange rate are strong predictors of their current states.

Figure 2.

Multivariate instantaneous structure (MIN).

Figure 3.

Multivariate autoregressive structure (MAR).

5. Conclusion and policy recommendations

This study examines the dynamic interactions among monetary policies and macroeconomic performances in Nigeria over the period of 1986Q1–2017Q4 with the application of BNG-SVAR. The use of the BNG-SVAR comes from the dynamic response of monetary policy to macroeconomic indicators in Nigeria. The nonstationary process of the data used led to test for their cointegration. The cointegration results show the existence of long-run relationship among the variables of interest. The P-GC results from the BNG-SVAR with the MCMC and M-H sampling techniques show that inflation is a strong predictor of interest rate in Nigeria given both the contemporaneous instantaneous (MIN) and the contemporaneous autoregressive (MAR) results. This is more reason why the Central Bank of Nigeria (CBN) in its period monetary policy committee (MPC) targets the inflation rate by choosing appropriate monetary policy rate (MPR) in response to the inflation rate in the country. This study, therefore, recommends that the inflation targeting in Nigeria should be broad and not limited to changing the MPR only. The fiscal discipline should be ensured from the Ministry of finance and the executive arm of the government. In any case, any fiscal policy should be directed towards the productive sectors of the economy. A major determinant of inflationary pressure in Nigeria is incessant supply of electricity. There should be massive investment in power generation and transmission in the country to eliminate the additional cost of production that drives prices up through the use of generating sets for production. Finally, the exchange rate policy should be to ensure domestic production to drive prices down rather than reliance on imported goods that promote imported inflation.

A.1. Appendix 1

See the Table A4.

VAR lag order selection criteria
Endogenous variables: M2, Ind, Infl, Intr, Exch
LagLogLLRFPEAICSC
0−4973.01NA7.41e + 2982.9668383.08297
1−4220.511429.7534.02e + 24*70.84180*71.53868*
2−4206.6425.189284.85e + 2471.0273872.30498
3−4190.3828.180715.64e + 2471.1730773.03140
4−4168.0736.812225.96e + 2471.2179073.65696
5−4122.2571.78438*4.28e + 2470.8709073.89069
6−4100.8931.686854.66e + 2470.9315474.53205
7−4081.6426.949655.32e + 2471.0273875.20861
8−4056.9732.488245.61e + 2471.0328075.79476

Table A4.

Optimal lag selection results.

Indicates the optimal lag length, where LogL, LR, FPE, AIC and SC indicate log likelihood, likelihood ratio, final prediction error, Akaike information criterion and Schwarz criterion.


Industrial outputExchange rateInflation rateInterest rateMoney supply
Mean32,497.8493.8619.2514.14,916,951
Median30,995.49117.1611.2513.661,294,950
Maximum52,931.79304.7273.126.723,388,300
Minimum18,998.233.762.14622,730.8
Std. dev7778.2569.4918.146,665,696
Skewness0.530.241.50.391.244417
Kurtosis2.652.533.883.873.17
Jarque-Bera6.532.3650.677.1732.15
Prob.0.030.310.000.030.00
Sum403011,638.862386.881748.66.10E + 08
Sum sq. dev7.44E + 09593,982.64.03E+041977.45.47E + 15
Observation124124124124124

Table 1.

Descriptive statistics for selected variables.

Augmented Dickey-FullerPhillips-Peron
VariablesLevelsFirst diff.VariablesLevelFirst diff.
Ind−1.7046−6.2244Ind−4.2882
Exch0.8368−8.4777Exch1.2655−8.3739
Infl1.4409−6.4756Infl−2.8021−10.695
Intr−2.8834−10.4781Intr−2.7495−11.2467
M25.1489−11.3326Mss5.5596−11.7113

Table 2.

Unit root tests.

Note: The critical values at 1, 5 and 10%, respectively, are −3.4846, −2.8853 and −2.5749 for both ADF and PP. Ind, Exch, Infl, Intr and M2 indicate industry output, exchange rate, inflation rate, interest rate and money supply, respectively.

Coint. rankEigenvalueTrace stat.Critical valueProb.
r=00.44108.6469.820.00
r10.1337.2947.860.33
r20.0919.8929.80.43
r30.058.2815.50.44
r40.011.793.840.18
Coint. rankEigenvalueEigen stat.Critical valueProb.
r=00.4471.3533.880.00
r10.1317.3927.580.54
r20.0911.6121.130.58
r30.056.4814.260.55
r40.011.83.840.18

Table 3.

Johansen cointegration results.

Note: The null hypothesis, H0, of no cointegration is rejected when the value of the trace and maximal eigen statistics is greater than the critical values at 5% significance level.

A.2. Appendix 2

The MCMC and the M-H algorithms are a proposal distribution to sample a new graph Gconditioned on a graph Gwith acceptance probability given as.

MGG=min PYGPYGPGPGZGGZGG1where PYGis the likelihood function, PGis the prior density and, finally, ZGGis the proposal distribution.

A.3. Appendix 3

Inverse Wishart Prior Posterior MCMC

The procedure for Gibbs sampling for the independent-normal Wishart prior is as follows:

  1. Draw Gikfrom the normal pGiYΣ.

  2. Draw Σi1kfrom the Wishart pΣi1YGi.

  3. Repeat steps 1 and 2N(20,000) times, and discard the first Nburn(10,000) iterations as burn-ins.

A.4. Appendix 4

Identification structure of the BGN-VAR results.

Bρ=0.900000.600.5000.700.500000.50.7000.6000.6B_0=000.800000000000000.5000.500000

© 2019 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

David Oluseun Olayungbo (July 17th 2019). Bayesian Graphical Model Application for Monetary Policy and Macroeconomic Performance in Nigeria, Bayesian Networks - Advances and Novel Applications, Douglas McNair, IntechOpen, DOI: 10.5772/intechopen.87994. Available from:

chapter statistics

105total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

Introductory Chapter: Timeliness of Advantages of Bayesian Networks

By Douglas S. McNair

Related Book

Advances in Statistical Methodologies and Their Application to Real Problems

Edited by Tsukasa Hokimoto

First chapter

Why the Decision‐Theoretic Perspective Misrepresents Frequentist Inference: Revisiting Stein’s Paradox and Admissibility

By Aris Spanos

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us