1

MODELLING BAHRAIN’S ECONOMY

A VECTOR AUTOREGRESSION (VAR) APPROACH

RIZWAN TAHIR

AHMED ABDUL GHANI

ABSTRACT: The last decade or so has witnessed a significant growth in research studies applying Vector Autoregression (VAR) technique for macroeconomic modelling. A VAR technique pioneered by Sims and popularized by researchers such as Litterman and Doan is useful particularly when knowledge about “true” structural relations is absent. This study represents the first attempt to apply such a technique to Bahraini yearly data (1971 – 2002) for five key macroeconomic variables. The results of the study indicated that all key macroeconomic variables are interlinked and influence each other. Oil exports is pure exogenous and is unaffected by any other macroeconomic variable. There is an evidence of bi-directional money-income causality and uni-directional causality from Government expenditure to CPI.. Impulse response and variance decomposition analysis suggest that fiscal policy is relatively more effective in short-run and monetary policy is more effective in the long run. Both oil exports and money are the sources of variation in GDP in the long-run and Government expenditures is the source of variation in the short-run. The result also suggests that money and oil exports are the important sources of variations in CPI in the long-run and Government expenditures is important source in the short-run. Inflation is found to be a fiscal phenomenon in the short-run but monetary phenomenon in the long-run. This result supports the view of monetarists.

KEYWORDS: VAR MODELLING, MACROECONOMIC POLICIES, OIL EXPORTS, BAHRAIN

1. INTRODUCTION

Being pioneer of oil producer in the Arabian Gulf region, Bahrain witnessed the prospects of potential economic prosperity in 1932 with the discovery of oil. Although oil exports contributed significantly in achieving higher levels of GDP over past few decades, its volatile nature (because of oil prices) and gradually decreasing share in GDP provided a challenge of maintaining higher levels of GDP. As a result, export base was diversified to non-oil products like Petrochemicals and Aluminium whose share in GDP has gradually increased. Inspite of diversifying sources of GDP, the rates of real GDP growth have showed wide fluctuations of more than 8 percent to negative 2 percent, over the period of last ten to fifteen years.

In order to analyze the sources of fluctuations in GDP growth, a standard complete structural macro model is probably desirable. However, such a model is derived on the basis of economic theory. Thereby two major problems arise: the theory must be exact enough to identify the endogenous and exogenous variables and the functional form connecting them. The second problem concerns the identification problem of recovering structural parameters from estimated reduced form. Out of these problems another class of nonstructural models: Vector autoregressive (VAR models) have been evolved, pioneered by Sims (1980) and popularized by researchers such as Litterman (1984) and Doan(1984). VAR model does not require any explicit economic theory to estimate a model. It uses only the observed time series properties of the data to forecast economic variables.

The VAR models have many applications (see Cooley and Leory, 1985). They are used to determine how each endogenous variable responds over time to a shock in that variable and in every other endogenous variable. VAR models are useful for analysis of the effect of alternative monetary or fiscal policies (Sims, 1982). The VAR models also provide a straightforward way of predicting the values of set of economic variables at any given point in time.

Our study represents the first attempt to apply such an approach in the case of Bahrain. In this paper, we develop and estimate an annual macroeconometric model for the economy of Bahrain over the period 1971 to 2002 using VAR technique proposed by Litterman (1984) and Sims (1980, 1982 & 1986). The main focus of this study is to analyze empirically the strength of short-run and long-run impacts of anticipated and unanticipated macroeconomic policies and oil exports shock (or innovations) on Bahrain’s macroeconomy.

The paper is divided into five parts. VAR approach is outlined next, followed by a discussion on the data used to estimate the model. The next part discusses the empirical results, concluding with the summary and conclusions.

2. THE VAR METHODOLOGY[i]

The methodology of the VAR is briefly described here. A k-equation VAR can be represented in a matrix form as follows:

A(L)Yt = A + Ut (1)

and

A(L) = I – H1L1 – H2L2 - …..HkLn (2)

Yt is an kx1 vector of variables, A is an kx1 vector of constants, and Ut is an kx1 vector of random variables. Equation (2) is an kxk matrix of normalized polynomial in lag operator L (Lm Yt=Yt-1) with the first entry of each polynomial on A’s being unity.

Since the right-hand side of the equations in the system contains only the predetermined variables, the error terms are assumed to be serially uncorrelated with constant variance and zero mean.

Hence, each equation in the system can be estimated using OLS. Moreover, OLS estimates are consistent and asymptotically efficient. Even though the error terms are correlated across equations, Seemingly Unrelated Regression (SUR) do not add to the efficiency of the estimation procedure since all regressions have identical right-hand-side variables. However, before estimating the model, the lag length must be chosen. If L is the lag length, number of coefficients to be estimated is k(kL + c), where c is the number of constants. The VAR model presented above indicates that the current innovations (Ut) are unanticipated but becomes part of the information set in the next period. The implication is that the anticipated impact of a variable is captured in the coefficients of lagged polynomials while the residuals capture unforeseen contemporaneously events. A joint F-test on the lagged polynomials provides information regarding the impact of the anticipated portion of the right-hand side variables. The impact of the unanticipated policy shocks (i.e. the policy variables such as changes in money supply and government expenditures) on other economic variables can be analyzed by employing the “impulse response functions” (IRFs) and “variance decompositions” (VDCs) that are obtained from a moving average representation of the VAR model given below[equations (3) & (4)]:

Yt = Constant + H(L)U (3)

and

H(L) = I + H1L + H2L + ….. (4)

where H is the coefficient matrix of the moving average representation, which can be obtained by successive substitution in equations (1) and (2). The elements of the H matrix trace the response over time of a variable i due to a unit shock given to a variable j.

The impulse response functions make it possible to analyze the dynamic behavior of the target variables due to unanticipated shocks in the policy variables. Variance decompositions show the portion of variance in the prediction for each variable in the system that is attributable to its own innovations and to shocks to other variables in the system.

3. DATA ANALYSIS

The data for the study is obtained from various publications of Government of Bahrain. Because of the unavailability of quarterly data, the model is estimated using yearly data from 1971 to 2002. The variables included in the VAR model are gross domestic product (GDP), consumer price index (CPI), government expenditures (GEXP), value of oil exports (XOIL) and money supply (M1). The objective here is to study the dynamics of the variables or the inter-relationship between these key macroeconomic variables, in particular, the influence of policy variables, such as Government expenditures or money stock and oil exports on economic activities. All the variables are in logarithms. It is important to note that in order to capture the finer details of the economy, a larger model with more variables would be desirable. However, with VAR models, one runs into serious degrees of freedom problems when the variables are many, especially with yearly data. For example, with 5 variables (k = 5) and 2 lags (L = 2), one would need to estimate 10 (kL) parameters (excluding the intercept) for each equation. Also, if the number of lags is increased by one, then with the same model, the estimated parameters increase to 15.

Unit root and cointegration tests

As a pre-requisite certain properties of the variables in the model must be checked in order to determine the appropriate specification for VAR estimation. The order of integration for each variable is determined using Augmented Dickey and Fuller (1979) and Phillips and Perron (1988) tests. The results of these tests are reported in table I. With the exception of ADF test with constant for GDP & CPI and with constant ant trend for CPI and PP test with constant for GDP, CPI & GEXP, all other ADP and PP tests for variables in log levels indicate that they are non-stationary. When first differenced in log, we find the evidence that the variables are stationary. ADF test with constant & trend and PP test with constant indicate the presence of two unit roots in CPI and GEXP. Since the results, overall, tend to suggest non-stationarity in log levels of the variables but stationarity in their log first differences, we proceed by contending that the variables belong to the I(1) process.

TABLE I

TESTS FOR UNIT ROOTS

Variables / ADF / PP
No constant & no trend / Constant / Constant & trend / No constant & no trend / Constant / Constant & trend
Log Levels
GDP / 1.397471 / -5.0341801 / -2.104510 / 1.772711 / -4.1690351 / -2.870965
CPI / 0.250873 / -8.1962641 / -7.6826551 / 1.594314 / -5.4168561 / -2.304288
GEXP / 1.381006 / 1.070895 / -1.185481 / 1.611615 / -5.0597011 / -3.363113
XOIL / 0.673281 / -3.085783 / -2.674549 / 0.692473 / -3.183007 / -2.727758
M1 / 1.684001 / -1.543266 / -2.491304 / 2.624443 / -1.835356 / -2.068293
Log
first differences
GDP / -2.867525* / -3.539993** / -4.451694* / -2.771254* / -3.660785** / -4.603926*
CPI / -6.051877* / -5.121595* / -2.7078432 / -1.697754*** / -2.1844072 / -4.048458**
GEXP / -2.784017* / -2.608528*** / -1.4313552 / -2.072438** / -2.6085282 / -3.297932***
XOIL / -5.472073* / -5.537618* / -5.854153* / -5.472099* / -5.537852* / -5.905655*
M1 / -2.518638** / -3.252550** / -3.233277*** / -2.460782** / -3.255888** / -3.182406**8

Notes:

1 reject null hypothesis (series has no unit root)

2 cannot reject null hypothesis (series has a unit root)

* reject null hypothesis (unit root) at 1 percent level;

** reject null hypothesis (unit root) at 5 percent level;

*** reject null hypothesis (unit root) at 10 percent level;

TABLE II

MacKinnon(1996) CRITICAL VALUES FOR ADF & PP UNIT ROOT TESTS

Level of significance / ADF / PP
No constant & no trend / Constant / Constant & trend / No constant & no trend / Constant / Constant & trend
Log Levels
1 percent / -2.644302 / -3.670170 / -4.356068 / -2.641672 / -3.661661 / -4.284580
5 percent / -1.952473 / -2.963972 / -3.595026 / -1.952066 / -2.960411 / -3.562882
10 percent / -1.610211 / -2.621007 / -3.233456 / -1.610400 / -2.619160 / -3.215267
Log
first differences
1 percent / -2.644302 / -3.670170 / -4.296729 / -2.644302 / -3.670170 / -4.296729
5 percent / -1.952473 / -2.963972 / -3.56879 / -1.952473 / -2.963972 / -3.568379
10 percent / -1.610211 / -2.621007 / -3.218382 / -1.610211 / -2.621007 / -3.218382

Since the five variables are noted to be I(1), there exists the possibility that they share a long-run equilibrium relationship, as was pointed out by Engle and Granger (1987). To test this, we used Eviews[ii] which implements VAR-based cointegration tests using the methodology developed in Johansen (1991, 1995). In formulating the dynamic model for the test, the question of whether an intercept and trend should enter the short- and/or long-run model is raised (Harris, 1995, p.95). We used all five deterministic trend models[iii] considered by Johansen (1995, pp. 80,84). The number of cointegrating relations from all five models, on the basis of trace statistics and the maximal eigenvalue statistics using critical values from Osterwald-Lenum (1992) at 5 percent level, are summarized in table III.

TABLE III

SELECTED NUMBER OF COINTEGRATING RELATIONS BY MODEL

TEST TYPE / MODEL 1 / MODEL 2 / MODEL 3 / MODEL 4 / MODEL 5
TRACE / 4 / 4 / 2 / 3 / 3
MAXIMUM EIGENVALUE / 2 / 2 / 2 / 3 / 3

Notes: the selection of cointegrating relations is based on .05 level critical values from Osterwald-Lenum (1992)


There is an evidence of minimum two and maximum four cointegrating relations. Generally, there are two different ways of specifying a VAR when the time series under study are cointegrated - an unrestricted VAR in levels or a VECM. Which specification is more appropriate remains debatable. While the VECM conveniently combines the long-run behavior of the variables and their short-run relations and thus can better reflect the relationship among the variables, there is no guarantee that imposing restriction of cointegration can be a reliable basis for making structural inferences (Faust and Leeper, 1997). Moreover, current finding is still unclear on whether the VECM outperforms the unrestricted VAR at all forecasting horizons. Naka and Tufte (1997) found that the two methods have comparable performance at short horizons. The support for the use of the unrestricted VAR can also be found in Clements and Hendry (1995), Engle and Yoo (1987) and Hoffman and Rasche (1996). Accordingly, with low computational burden required by the VAR in levels, we implement the VAR using the variables in levels.