Acta Universitatis Danubius. Œconomica, Vol 12, No 5 (2016)

A Bayesian Approach for the Analysis of Macroeconomic Dynamic in Case of Emerging Countries-Monetary and Fiscal Policy Model



Georgiana-Alina Ionita1



Abstract:The paper proposes the analysis of the main drivers of the economic growth in Central and Eastern Europe, in three emerging countries: Czech Republic, Hungary and Poland, with a development stage similar with that of Romania. Given the vulnerabilities of the Central and Eastern Europe region at the beginning and during the recent global economic and financial crisis, there is an increased interest to identify the models that can describe the principalcharacteristics of the Central and Eastern Europe macroeconomic variables: gross domestic product, investment, wages and salaries, inflation, hours worked, consumption and themonetary variable- interest rate. Moreover, another scope is to analyze the frictions that describe the evolution of the seven data series, as the stochastic dynamic of the macroeconomic model is driven by orthogonal structural shocks.

Keywords:monetary; policy; frictions; shocks

JEL Classification: C01; D50; B22; C4.



1. Introduction

I have proposed to analyze in the current thesis the degree in which the model responses to the requests of the Central and Eastern Europe economies, based on a study completed for the following countries: Poland, Hungary and Czech Republic. In this purpose I have used a Bayesian approach for the estimation of this forward-looking model, in a general equilibrium framework.

The current model is not a simple monetary policy model, involving also the fiscal policy, so it can be analyzed taking into account the policy shocks involved: price and interest rate shocks (as monetary policy shocks) and also the exogenous spending shock (which includes together with net exports also the government spending, instrument of the fiscal policy).

So, the model which is an improvement and a simplified version of the one proposed by the authors Frank Smets and RafWouters in 2003 in the article “An Estimated Stochastic Dynamic General Equilibrium Model for the Euro Area has the purpose to analyze the main drivers of the economic growth, putting an eye on the frictions of these four economies (seven frictions- reduced to the number of the seven observed variables).

As a result, the research proposes an overview of the dynamic of monetary policy that should be conducted in the face of multiple sources of uncertainty, including model and parameter uncertainty as well as uncertainty about future shocks.



2. The Non-Linear Model

2.1. Producers

The product with destination of final consumption, ,is composed of good for intermediate consumption, that are bought by final good producers, grouped into and sold in a perfectly competitive market.

They maximize the profit obtained, as per maximization function:

Max , s.t. ( ),

, are the prices of goods for final consumption and for intermediate consumption and G is a function having the characteristic of being strictly concave on one hand and being an increasing function on the other hand, with the property: G(1)=1.

istheprice mark-up shock, and follows and exogenous ARMA processln .

2.2. Intermediate Goods Producers

The producers of goods with the destination of intermediate consumption follow the technology equation:

Y (i)= - , represents capital in services form, represents the input of labor, while is a fixed cost.

represents the labor growth rate and represents the shock of the productivity factor and is defined as:

ln .

The profit of any firm from the economy follows the following equation: - , represents nominal wage or salary rate and represents the capital rental rate.

The model assumes that firms are able to adjust prices used with probability 1- in each period.

The resulting price that is optimal is obtained from the following maximization function:

Max (

s.t. .

represents firms probability of price adjustment, is inflation defined as = is the discount factor for firms,

= and

Theindex of prices obtained has the following equation form:

P =(1- .

2.3. Households

Households in order to maximize utility function are able to choose final consumption , the number of working hours , acquisition of bonds ,capital investment and capital utilization .The utility function is as follows: , related to the following equation of the budget:

Capital at moment “t” has the following form: ,

represents the stochastic premiumresulted frominvestment in bonds, follows the stochastic process: ln = ,

where, has the significance ofrate of depreciation, S(.) is a function for quantifying the cost adjusting, with S( )=0, S’( )=0, S’’(.)>0. illustrates the shock of investment component and is described by the equation: ln = . represents taxes, while represents the dividends for distribution.

2.4. Intermediate labor and labor packers

Households provide their labor forthe intermediate labor union, the labor provided having the following form: . (16)

Labor packers are those who acquire labor from unions and distribute the labor to producers of intermediary goods. Their maximization profit function has the following form:

Max , s.t. ( )

,where and represent the price of total and for intermediary labor services , while H is a function strictly concave , which follows H(1)=1 and is increasing.

is an exogenous shocks of wages mark-up, and has the form of an ARMA process:

From the first order conditions of the labor packers, we obtain:

.

Combining this condition with the zero profit condition we obtain an expression for the wage cost for the intermediate goods producers:

, where is defined as an ARMA process:

ln = .

The dividends that are received by households from labor unions are included in the constraint of budget for households:

In case of unions there are also nominal rigidities as proposed by Calvo (1983), more precisely wages being adjusted with the probability of 1- .

In case of unions who readjust wages, the optimization rule consists of choosing a wage in order to maximize subsequent wage income in case when unions keep this fixed wage.

The expression of aggregate wageobtained is:

.

2.5. Monetary policy and government budget constraint

The interest rate established by central bank by taking into account the deviation of output and inflation from the targeted levels is as follows:

.

is the value of nominal interest rate at steady-state and is the natural gross domestic product.

The parameter represents the interest ratesmoothness, while the definition of the shock of monetary policy, , is: .

The constraint of government budget is described as follow : . represents the nominal lump-sum taxes and the government spending in relation with the steady-state output follows the process:

.



3. The Linearized Model

The aggregate constraint for the linearized model is described as follows: (1) ,

The interpretation of the resource constraint is that gross domestic product ( ) is absorbed by investment ( ), consumption ( ), capital utilization costs (expressed in relation to the capital utilization rate ( ) and the exogenous spending shock ( )). In addition, is the state-state share of consumption in output and is equal to , where and are the steady-state exogenous spending-output ratio and investment-output ratio.

In addition, , where is the steady-state rental rate of capital and is the steady-state capital-output ratio.

The dynamics of consumption function is described as follows:

(2) , where , and .

Consumption ( ) is defined in a direct relation with the past and expected consumption, with workinghours expected increase , with ex-ante real interest-rate ( ) and theshock factor- .

The investment dynamic is described by the function below:

(3) ,where , , represents represent the steady state elasticity of the cost function for capital adjustment, while represents the factor of discount.

The corresponding arbitrage equation for the capital is given by:

(4): in direct relation with expected value and expectedcapital rental rate ( ) and inverse relation with therisk premium shock and ex-ante real interest rate.

The production function is described by the following equation: (5) .

The production of output by firms is based on two components: capital ( ) and labor ( ).

Current capital services used ( ) depends onpreviously capital and capital utilization rate ( ): (6) .

Theaccumulated capital ( ) is depends on investment and its technology disturbance:

(8) , where and .

In case of goods market, price mark-up, ( ), is defined by differentiating the marginal product of labor ( ) and the real wage ( ):

As firms that are price adjusters have the objective of profit maximization, this leads to the New-Keynesian Philips curve described below: (10) , where , and .

Rental rate of capital is in an inverse relation with the ratio capital-labor and in a direct relation with the real wages: (11) .

In case oflabor market the wage mark-up is described by the difference between real salaries and marginal substitution rate working-consumption : (12) , where represents elasticity of labor supply in relation with real salaries and illustrates the habit in consumption of households.

In terms of real wages it is assumed that are adjusted progressive in order to reach the target level of wage mark-up, as follows:

(13) , where

, , and .

The monetary policy equation is described as follows:

(14): .

The authorities responsible with the monetary describe a generalized Taylor rule and adjust the interest rate ( ) based on information regarding inflation and output gap. The parameter captures interest rate smoothing degree.

The linearized system is composed of 14 equations and 14 endogenous variables: , , , , , , , , , , , , , .

The stochastic behavior of the system of linear rational expectations equations is driven by seven exogenous disturbances which follow a first-order autoregressive dynamic: productivityfactor , risk premium , investment technology ,wage mark-up , price mark-up , exogenous spending and monetary policy .



4. Solution and Analysis of the Model

4.1. Econometric Estimation Methodology

In the current research I have proposed to analyze the way in which the proposed DSGE model responses to the economic matters of the three analyzed countries: Poland, Hungary and Czech Republic.

I have chosen to estimate the DSGE model through a Bayesian approach, using Matlab and Dynare tool as this approach is able to give estimates based on the a-priori distributions assumed of the parameters and the information brought by the seven observed variables (extracted with the Kalman filter maximization of the likelihood function), in line with F.Canova, 2007, pp. 26-45, inthe article “Methods for Applied Macroeconomic Research”, and Johannes Pfeifer, 2014, pp. 1-150 in the articleA Guide to Specifying Observation Equations for the Estimation of DSGE Models.

Moreover, using the Bayesian approach I was able to include in the estimation the seven shocks of the model (each one for each observed variable) and to estimate their standard deviation.

I have used for the estimation of the model 350,000 Metropolis-Hastings draws, with a scale of 0.3, in case of each of the three countries, in line with Geweke, J. (1998), “Using simulation methods for Bayesian econometric models: inference, development and communication”, pp. 1-75 and with Kimball, M. (1995), “The quantitative analytics of the basic neo-monetarist model”, pp. 1241-1277.

4.2. Calibration and prior distribution of the parameters and parameter estimates

One important step of estimation of a DSGE model consists of calibration of the model’s parameters.

Regarding the choice of the prior distribution, a number of parameters were kept fixed from the start of the exercise.

The standard errors of shocks are defined as inverse-gamma distributions, having the characteristics described in the table below.

Table 1. Parameters calibrated -shock factors

Symbol

Description

Value

 

the t-1 term coeficient, in the AR process

0.968

 

the t-1 term coefficient of the AR process

0.27

the t-1 term coefficient of the AR process

0.9928

the t-1 term coefficient of the AR process

0.7165

the t-1 term coefficient of the AR process

0

the t-1 term coefficient of the process

0

the t-1 term coefficient of the process

0

the coefficient of the IID-Normal price mark-up shock (error term)

0

the coefficient of the IID-Normal wage mark-up shock (error term)

0

 

price mark-up disturbance

10

 

wage-markup disturbance

10

Other parameters, such as the discount factor, , is calibrated to be 0.9994, which is the mean of the sample of the quarterly real interest rate, while the exogenous spending-GDP ratio is set at 18%. The demand elasticity for labor, , is equal with 1.5, the degree of wage stickiness, , is 0.7937, the degree of price stickiness, , is 0.75, which is in accordance with the assumptions with the calibrations proposed by the author Mihai Copaciu (2012), in the article “Estimation of an open economy DSGE model with financial and employment frictions for Romania”.

Table 2. Parameters calibrated

Symbol

Description

Value

 

share of capital in production

0.4

 

is the discount factor applied by households.

0.999

 

elasticity of intertemporal substitution between consumption and leisure

1.5

 

reflects the presence of fixed costs in production (one plus the share of fixed costs in production)

1.5

 

the degree in which the productivity shock impacts the exogenous spending

0.51

 

the steady-state elasticity of the capital adjustment cost function

6.014

 

the habit parameter in consumption (external consumption habit )

0.636

 

the degree of wages stickiness

0.794

 

the elasticity of labour supply with respect to the real wage

1.942

 

the degree of price stickiness

0.75

 

wages indexation to past inflation indicator

0.324

 

the degree of price indexation to past inflation

0.329

 

the degree of capital utilization

0.27

 

depreciation rate of capital

0.0246

 

demand elasticity for labor

1.5

 

steady-state exogenous spending-output ratio

0.18



Finally, parameters used for the quantification of monetary policy rule follow a standard Taylor rule and are calibrated as described in the table below:

Table 3. Parameters calibrated-Taylor rule

Symbol

Description

Value

 

the coeficient of inflation in the interest rate Taylor rule

1.488

 

the degree of interest rate smoothing

0.826

 

the coeficient of the output-gap in the interest rate Taylor rule

0.059

 

the coeficient of the output gap variation (t, t-1)

0.224

I have also analyzed the posterior distribution of parameters and standard deviation of shocks, before and after the optimization using Metropolis-Hastings algorithm. The results of the posterior distribution of the parameters obtained through the Metropolis- Hastings sampling algorithm was based on 350,000 draws and 2 Markov Chains.

4.3. Data Set

The dataset used in this paper were collected from Eurostat Database, Federal Reserve Economic Data -FRED –St. Louis Fed and European Central Bank-ECB, National Bank of Czech Republic, Hungary and Poland for the period: 2001q1 and 2015q4, for a quarterly frequency.

The seven observed quarterly macroeconomic variables are the following: real GDP, consumption, investment and wages, hours worked, GDP deflator and short-term nominal interest rate.

The GDP, investment (gross fixed capital formation-namq_gdp_c) and consumption series (final consumption expenditure in current prices (namq_fcs_c) were collected from Eurostat site, being expressed in Eur millions, in current prices. Regarding the real wages, I have collected from Eurostat the quarterly wages series ([namq_nace10_c]), representing the gross wages and salaries in Eur millions. All these four variables: investment, consumption, wages and gross domestic product are deflated by gross domestic product deflator for the use of real variables and expressed per capita, being divided by the number of employees (16 years to 64 years, part-time and full-time employment-Eurostat series: [lfsq_epgais]).

All the real variables (real gross domestic product, investment, consumption andsalaries) have been expressed in 100 times log, in order to get an evidence of the real growth rates.

In terms of monetary policy, given the high volatility of ON interest rate, I have collected the daily rates for PRIBOR 3M, BUBOR 3M and, WIBOR 3M, from 2001q1 and I have computed an average interest rate for each quarter during 2001 q1 to 2015 q4.

The series used in Dynare soft for the observable variables were: the log difference of real GDP, the log difference of real consumption, the log difference of real investment, the log difference of real wages, the log of the index of hours worked, the log difference of GDP deflator and the quarterly interest rate.

Once all the observed variables are seasonally adjusted and tested for stationarity in Eviews 7, imported in Matlab for a further processing with Dynare 4.4.0 tool.



5. Econometric Analysis of the Model

5.1. Interpretation of Apriori and Aposteriori Distribution Graphs

For each parameter the graphs illustrates the prior and the posterior distribution in one figure (the grey line representing the prior, the black line the posterior). The distributions are different from each other, so the parameters are identified.

The dotted green line represents the value at the posterior mode, which in case of all parameters is approximately the center of the posterior distribution.

The posterior distribution is also distributed around the mode. Moreover, in almost all cases the posterior distribution has the form of a normal distribution and is close, but still different from the a-priori distributions.

This means that the observable variables bring new information than the assumptions made regarding the a-priori distribution.

Each graph contains the a-priori assumptions that were made regarding each parameter’s mean and, respectively, the a-posteriori results with the mode estimated by numerical methods, posterior approximations of the standard deviations, all obtained after the maximization algorithm, the posterior mean obtained from MH algorithm for the 350,000 draws on the 2 chain that I chose.

Overall, according with the a-priori and a-posteriori distribution graphs and t-statistics, all estimated parameters are significantly different from zero. This is true for all the standard errors of all the shocks.



Poland

Graph 9. A-priori and a-posteriori distribution

Graph 10. A-priori and a-posteriori distribution



Graph 1. A-priori and a-posteriori distribution

Graph 2. A-priori and a-posteriori distribution



Hungary

Graph 3. A-priori and a-posteriori distribution

Graph 4. A-priori and a-posteriori distribution



Graph 5. A-priori and a-posteriori distribution

Graph 6. A-priori and a-posteriori distribution

The tables described above contains in the first part the a-priori assumptions for the mean of each parameter and in the second part the a-posteriori results, estimated by numerical methods, the posteriori estimation of the standard deviation (obtained after the maximization algorithm), the posterior mean, the 10% and 90% percentiles (confidence intervals) obtained from the MH algorithm, for the 350,000 draws on the 2 chains that I have chosen.

5.2. The Model’s Stability Analysis –Blanchard-Kahn Condition

The system is stable according with the analysis of the eigenvalues of the system.

In order to meet the Blanchard-Kahn condition, there must be as many roots larger than one in modulus as there are forward-looking variables in the model.

According with the output of Matlab, in the model there are 12 eigenvalue(s) larger than 1 in modulus for 12 forward-looking variable(s), which means that the Blanchard-Kahn condition is met. The rank condition is also verified.

5.3. Analysis of Smoothed Variables and Shocks

The posterior distribution of smoothed endogenous variables and shocks, i.e. infers about the unobserved state variables using all available information up to T: .

Using the smoothed shocks as a check for the model, the shock realizations are around zero (being white noises).



Hungary

Poland

Graph 7. Smoothed shocks

Graph 8-Smoothed shocks



Hungary

Poland

Graph 9. Smoothed Variables

Graph 10. Smoothed Variables

According with the estimates all the seven shocks are significant, taking into account the 90% HPD (highest probability density) confidence interval, that indicates the fact that with a 90% probability, the parameter is in the calculated interval.

5.4. Brooks and Gelman’s Markov Chain Monte Carlo (MCMC) Convergence

As resulted from the convergence graphs analysis for 2 Monte Carlo Chains, from the three economies analyzed, I concluded that in case of Czech Republic the convergence condition of the series is not respected, taking into consideration a number of 350,000 MCMC draws. The acceptance ratios of the MCMC are lower than the acceptance limit of 25% (chain 1: 9.452% and chain 2: 8.3891%). As a result, analysis of the steady-state does not lead to accurate result (the reason why this country was excluded from the further analysis).

The current acceptance ratio per chain were the following, at the level of each country: for Romania- chain 1: 26.5536% and chain 2: 23.9725%, for Poland-chain 1:21.1708% and chain 2:24.6608% and for Hungary-chain 1:36.2359% and chain 2:37.0039%.



Poland

Hungary

Graph 11

Graph 12



Czech Republic





6. Conclusion

The proposal of this working paper is to analyze whether the model proposed for the economy of United Statesby Frank Smets and RafaelWouters in 2002, pp. 1-70 and 2003, pp. 1-57 and also in articles for the economy of Euro Area published in 2005, pp. 1-52 and pp. 2-36, and in 2003, pp. 1123-1175 and pp. 1-57 is suitable for the economies from Central and Eastern Europe.

Taking into consideration the tests performed I conclude through the current working paper that the results of the model are effective for the economy of Poland and Hungary. This is confirmed also by the analysis of the estimations performed: convergence analysis, Blanchard-Kahn stability of the model, the analysis of the a-priori and a-posteriori distributions and also by the analysis of parameters and variables identification. For the three economies analyzed, I concluded that in case of Czech Republic the convergence condition of the series is not respected, taking into consideration a number of 350,000 draws. As a result, analysis of the model for this country is not accurate, so it was excluded from the further analysis.

To conclude, the model has achieved the purpose of analyzing of the main drivers of the economic growth, putting an eye on the frictions of these three economies (seven shocks- reduced to the number of the seven observed variables).

The model can be analyzed also taking into account the policy shocks: exogenous spending shock (that includes government spending-fiscal policy instrument), price and interest rate shocks (monetary policy shocks).

Regarding the monetary policy, I can conclude that the Taylor rule proposed by the model approximates the behavior of the Central banks of the two Central and Eastern Europe countries: Poland and Hungary.



7. References

Canova, F. (2007). Methods for Applied Macroeconomic Research. New Jersey: Princeton University Press, pp. 26-45.

Frank Smets, RafWouters, (2002). An estimated dynamic stochastic general equilibrium model. Research series, pp. 1-70.

Geweke, J. (1998). Using simulation methods for Bayesian econometric models: inference, development and communication, mimeo. University of Minnesota and Federal Reserve Bank of Minneapolis, pp. 1-75.

Pfeifer, Johannes (2014). A Guide to Specifying Observation Equations for the Estimation of DSGE Models. Research series, pp. 1-150.

Kimball, M. (1995). The quantitative analytics of the basic neo-monetarist model. Journal of Money, Credit and Banking, 27(4), pp. 1241-1277.

Smets, F. & R. Wouters (2003a). An estimated Dynamic Stochastic General Equilibrium Model of the Euro Area. Journal of the European Economic Association 1(5), pp. 1123-1175.

Smets, F. & Wouters, R. (2003b). Shocks and Frictions in US business cycles: a Bayesian DSGE Approach, European Central Bank, pp. 1-57.

Smets, F. & Wouters, R. (2005). Comparing Shocks and Frictions in US and Euro Business Cycles: a Bayesian DSGE Approach. Journal of Applied Econometrics, pp. 2-36.

Smets, F. & Wouters, R. (2005). An Estimated Two-Country DSGE Model for the Euro Area and the US economy. Research series, pp. 1-52.

Smets, F. & Wouters, R. (2007). Shocks and Frictions in US and Euro Business Cycles: a Bayesian DSGE Approach. Journal of Applied Econometrics, pp. 1-57.



1PhD Student, Academy of Economic Studies, Doctoral School of Finance and Law, Bucharest, Romania, Address: Piata Romana, Nr. 6, Cladirea Ion N. Angelescu, Camera 0327, Romania, Tel.: +4021 319 1900, Corresponding author: Georgiana.alina.ionita@gmail.com.

AUDŒ, Vol. 12, no. 5, pp. 51-69

Refbacks

  • There are currently no refbacks.
Creative Commons License
This work is licensed under a Creative Commons Attribution 4.0 International License.