Acta Universitatis Danubius. Œconomica, Vol 10, No 1 (2014)

The Kalman Filter Approach for Estimating the Natural Unemployment Rate in Romania



Mihaela Simionescu1



Abstract: The aim of this research is to determine the monthly natural rate of unemployment during the third quarter of 2013 in Romania. The Phillips curve approach is not valid for the Romanian economy, but Kalman filter is a suitable approach for computing the natural rate of unemployment. We make the assumption that the cyclical component follows a random walk. Predictions were made for the unemployment rate in Romania using Kalman approach during July-September 2013 and on this horizon an insignificant decrease was observed from a month to another. A value of 5.85% is expected for unemployment rate in Romania in September 2013.

Keywords: Kalman filter; natural rate of unemployment; forecasts; random walk

JEL Classification: E21; E27; C51; C53



  1. Introduction

This Kalman approach is usually applied in determining the natural unemployment rate, the value for each we have a reasonable level or a stability of inflation rate and wages. The Phillips curve used to describe the relationship between inflation and unemployment rate is not checked in Romania, but the state space models are valid.

The objective of this research is to determine the monthly natural unemployment rate in Romania and to make predictions using Kalman filter. There are not relevant studies till now for the Romanian economy.

The organisation of this research is clear: after a brief literature presentation of the quantitative methods used in predicting unemployment rate, we explained the used methodology. One-step-ahead predictions are made for unemployment rate in Romania during the third quarter of 2013 using Kalman filter.



  1. Recent Results in Literature

A complete study related to the Measurement of the natural rates, gaps, and deviation cycles is provided by Murasawa (2013). Claar (2005) estimated the natural rate of unemployment using the Kalman filter for the civilian unemployment rate in USA during 1977-2002. The author also studies the relationship between the natural rate of unemployment and other macroeconomic variables of the labour market. Moreover, Groenewold and Hagger (2002) pointed out before that the natural rate of unemployment is model dependent. Garlach-Kristen (2004) estimated the natural unemployment rate assuming that it follows a random walk, being a determinant of Beveridge curve. Valletta (2006) used the same approach of Beveridge curve, but utilizing regional data. Basistha and Startz (2008) reduced the uncertainty that affects the NAIRU natural rate of unemployment by using multiple indicators.

King and Morley (2003) estimated the natural rate of unemployment without the utilization of the Phillips curve, considering that the natural rate that varies in time is endogenous. Schreiber (2011) estimated the natural rate of unemployment for euro countries by using the integrated systems. Greenslade, Pierse, and Saleheen (2003) applied Kalman filter technique to England Phillips curve models for the NAIRU unemployment during 1973-2000. Meļihovs and Zasova (2009) determined the natural unemployment rate for Latvia using Phillips curve for quarterly data.

Two parallel disturbances are presented for unemployment: a permanent effect and a temporary one. The permanent component is represented by supply shocks that modify the full-employment level while the temporary effect does not modify this full-employment level of output as in the approach of King, Stock and Watson (1995), Staiger, Stock and Watson (1997) and Gordon (1998). According to Apel and Jansson (1999) the cyclical component of unemployment presents serial correlation. Proietti (2003) compared the accuracy of several predictions based on linear unobserved components models for monthly US unemployment rate, drawing the conclusion that the shocks are not persistent during the business cycle.

Camba-Mendez (2012) built conditional forecasts for unemployment rate using VAR models and Kalman filter techniques. Sermpinis, Stasinakis and Karathanasopoulos (2013) made predictions for US unemployment rate, using Neural Networks and compared the utility of Support Vector Regression (SVR) and Kalman Filter in combining these forecasts.



  1. Methodology

The Kalman filter is an econometric method for predicting the endogenous variables and for adjusting the estimated parameters in forecast equations. There are two systems of equations: a system of prediction equations and a system of update equations.

The stages for applying the Kalman filter are:

  1. The estimation of endogenous variables values using available prior information;

  2. The adjustment of estimated parameters using adjustment equations and the computation of prediction errors.

A state space model includes two equations:

Measurement equation (the relationship between the observed and the unobserved variables): yt = Htβt + Azt + et

Transition equation (the dynamic of state (unobserved)): βt = μ + Fβt-1 + vt

yt – data series

t –observed explanatory variables

Ht – variable coefficients of unobserved series

βt, A, F and F’ – constant coefficients

R and Q- state space parameters (matrix of covariance)

et and vt – shocks

Assumptions

et iid. N(0, R)

vt iid. N(0, Q)

E(et, vt) = 0

The objectives are:

1. The estimation of state space model parameters;

yt = Htβt + Azt + et

βt = μ + Fβt-1 + vt

et iid. N(0, R)

vt iid. N(0, Q)

2. Restoration of the unobserved state;

yt = Htβt + Azt + et

βt = μ + Fβt-1 + vt

et iid. N(0, R)

vt iid. N(0, Q)

βt/t-1 – the estimation of βt latent state according to the information till t-1 moment

βt/t – the estimation of βt state according to the information till t moment

Pt/t-1 - the βt covariance according to the information till t-1 moment

Pt/t C- the βt covariance according to the information till t moment

yt/t-1 P- the prediction of y using the information till t-1¬moment

ηt/t-1¬ = yt – yt/t-1 - error prediction

ft/t-1 -the variance of prediction error

The Kalman filter offers an optimal estimation for βt, conditioned by the information related to the Ht state space parameters: A, μ, F, R, Q.

We suppose that μ, F, R, Q are known. The recursive Kalman filters implies 3 stages:

1. We start with the supposed values at the initial moment 0: β0/0 si P0/0;

2. The prediction: the optimal prediction y1/0 at moment 1, using β1/0;

3. The update: the calculation of the prediction error, using the observed value for y at moment 1.

η1/0 = y1 – y1/0

The information included in the prediction error has data that can be recovered for redefining our assumption regarding the value that β could have

β1/1 = β1/0 + Kt η1/0

Kt - the Kalman gain (the importance accorded to the new information).

The predicted values

βt/t-1 = μ + Fβt-1/t-1

Pt/t-1 = FPt-1/t-1¬F' + Q

The prognosis for y and the error prediction

ηt/t-1 = yt – yt/t-1 = yt - ztβt/t-1

ft/t-1¬ = xtPt/t-1z't + R

The update

βt/t = βt/t-1 + Kt ηt/t-1

Pt/t = Pt/t-1 – KtZtPt/t-1

Kalman gain: Kt = Pt/t-1 z't (ft/t-1)-1.

The actual observed unemployment rate is the sum of two components: the natural unemployment rate quantifying the persistent shocks from the supply side (we assume it follows a random walk) and the cyclical unemployment that refers to the shocks from the demand side which are limited as persistence (this component exhibits the serial correlation).

 

 

 = 

 ~ N(0; )

 ~ N(0; )

E( ) = 0

A state space model for the natural unemployment can have the following form:

  , t=1,2,…,T (measurement equation)

Z=[1 1],  

  (transition equation)

T= ,  

 ~ N(0; )

 ~ N(0; )

E( ) = 0

Under these conditions the Kalman filter generates optimal predictions and updates of the state variables. The Kalman filter determines the estimator of the minimum square error of the state variables vector. There are two approaches in literature regarding the estimation of a variable using this filter. The first one assumes that the initial value of the non-stationary state variable can be fixed and unknown. On the other hand, the second approach considers that the initial value is random. The diffuse prior is specified. If we analyse the first observations, the approach is better even if it can generate numerical instability. If m is the number of state variables we utilize the approach with diffuse prior of Koopman, Shepard and Doornik (1998) and m predictions are provided. The unknown parameters that will be estimated are   and  . However, some authors give these parameters some reasonable values from the start. For we have to establish the value from the start and the log-likelihood function is computed. The variance of the shocks coming from the demand side ( ) is always greater than the variance of supply shocks ( ).



  1. The Computation of Natural Unemployment Rate and of the Predicted Unemployment

In this research the data set is represented by the unemployment rate in Romania (denoted by u) registered in the period 1992: January- 2013: June. The unemployment rate is an indicator used to measure the unemployment intensity, being computed as a ratio of number of registered unemployed people and the active population. One-step-ahead predictions are made on the horizon 2013: July- 2013: September. The data series are provided by the National Institute of Statistics.

The natural unemployment rate is determined for diffuse prior and different values of    represents the starting value of the state space model.

 = , where   is the error term of the model that explains the evolution of the unemployment rate using the natural unemployment rate

The estimations based on Kalman filter are made in EViews:

@ signal ur= sv1+ sv2

@ state sv1= sv1(-1) + [var=exp(c(2))]

@ state sv2= c(4)* sv2(-1) + [var=exp(c(3))]

The state space models for different values of starting value of   are presented in Appendix 1. The proposed models in literature are also valid for Romania.



Table 1. The Natural Unemployment rate for Different Values of Starting Values (July 2013-September 2013)

Month

Unemployment rate (%)

(dynamic forecasts)


  =1

  =0.9

  =0.8

  =0.7

  =0.5

  =0.3

July 2013

5.52

5.516

5.516

5.517

5.5177

5.518

August 2013

5.517

5.515

5.515

5.515

5.518

5.517

September 2013

5.518

5.515

5.516

5.5166

5.517

5.517

Dynamic forecasts are made for different values of   (July 2013-September 2013). These values include not only the natural unemployment rate, but also the cyclical component. For July 2013 the Kalman filter approach predicts a rate of 5.88% for the unemployment rate, followed by an insignificant decrease till 5.87% in August 2013 and 5.85% in September 2013.

Table 2. Dynamic Forecasts of the Unemployment Rate for Different Values of Starting Values   (July 2013-September 2013)

Month

Unemployment rate (%)

(dynamic forecasts)


  =1

  =0.9

  =0.8

  =0.7

  =0.5

  =0.3

July 2013

5.8862

5.8862

5.88621

5.886239

5.886226

5.886235

August 2013

5.87249

5.87253

5.87246

5.87251

5.87248

5.87250

September 2013

5.85878

5.85885

5.85874

5.85881

5.85877

5.85880

The differences between the forecasts corresponding to a certain month are insignificant. The increase in the value of   does not imply necessary an increase in the value of the unemployment rate. For July 2013, the most accurate unemployment rate forecast was registered for the case of   =0.5 (with an absolute error of 0.59622 percentage points).

The one-step-ahead forecasts based on Kalman filter and the actual values of unemployment rate are represented in the following graph.



Figure 1. The Actual and Predicted Values of Monthly Unemployment rate in Romania (1992: January- June: 2013)

As we can observe, the differences between the actual values and the predicted ones are low. In 2002 the greatest unemployment rates were registered.



  1. Conclusions

An important conclusion is that the classical state space model used in literature to determine the natural unemployment rate provided expected results for the Romanian economy. A very slow decrease in the monthly unemployment rate is observed during the third quarter of 2013 when Kalman approach is used. A value of 5.85% is predicted for September 2013.

This research provides pertinent results regarding the prediction of unemployment rate in Romania, but the study could be improved by assessing the forecasts accuracy and making the comparison with other predictive quantitative techniques.



  1. References

Apel, M. & Jansson, P. (1999). System of Estimates of Potential Output and the NAIRU. Empirical Economics, No. 23, pp. 378-388.

Basistha, A. & Startz, R. (2008). Measuring the NAIRU with Reduced Uncertainty: A Multiple-Indicator Common-Cycle Approach. The Review of Economics and Statistics, MIT Press, Vol. 90, No. 4, November, pp. 805-811.

Camba-Mendez, G. (2012). Conditional Forecasts on SVAR Models Using the Kalman Filter. Economics Letters, Vol. 115, No. 3, June, pp. 376-378.

Claar, V. (2005). The Kalman Approach to Estimating the Natural Rate of Unemployment. Proceedings of Rijeka School of economics and Business, No. 23, pp. 1-24.

Garlach-Kristen, P. (2004). Estimating the Natural Rate of Unemployment in Hong Kong, Hong Kong Institute of Economics and Business Strategy. Working Paper, pp. 1-13.

Greenslade, J. V.; Pierse, R. G. & Saleheen, J. (2003). A Kalman Filter Approach to Estimating the UK NAIRU. Bank of England, Working Papers 179. Bank of England.

Groenewold, N. & Hagger, A. J. (2000). The Natural Rate of Unemployment in Australia: Estimates from a Structural Var. Australian Economic Papers, 39, pp. 121-137.

King, T. B. & Morley, J. (2003). In Search of the Natural Rate of Unemployment. Working Paper. Washington University in St. Louis.

King, R. G.; Stock J. H. & Watson M. W. (1995). Temporalinstability of the Unemployment-Inflation Relationship. Economic Perspectives, Vol. 14, No. 3, pp. 2-12.

Murasawa, Y (2013). Measuring the Natural Rates, Gaps, and Deviation Cycles. Empirical Economics, Online publication date: 14-Sep-2013.

Proietti, T. (2003). Forecasting the US Unemployment Rate. Computational Statistics & Data Analysis, Vol. 42, No. 3, 28 March, pp. 451–476.

Schreiber, S. (2011). Estimating the Natural Rate of Unemployment in Euro-Area Countries with Co-Integrated Systems. Applied Economics, 1-21, Online publication date: 15-Feb-2011.

Sermpinis, G.; Stasinakis, C. & Karathanasopoulos, A. (2013). Kalman Filter and SVR Combinations in Forecasting US Unemployment. Artificial Intelligence Applications and Innovations IFIP Advances in Information and Communication Technology, Vol. 412, pp. 506-515.

Staiger, D.; Stock, J. H. & Watson, M. W. (1997). The NAIRU, Unemployment and Monetary Policy. Journal of Economic Perspectives, Vol. 11, No. 1, pp. 33-49.

Gordon, R. J. (1998). The Time-Varying NAIRU and Its Implications for Economic Policy. Journal of Economic Perspectives, Vol. 11, No. 1, pp. 11-32.

Valletta, R. G. (2006). Why Has The U.S. Beveridge Curve Shifted Back? New Evidence Using Regional Data. Working Paper Series 2005-25, Federal Reserve Bank of San Francisco.

Meļihovs, A. & Zasova, A. (2009). Assessment of the Natural Rate of Unemployment And Capacity Utilisation in Latvia. Baltic Journal of Economics, Baltic International Centre for Economic Policy Studies, Vol. 9, No. 2, December, pp. 25-46.



APPENDIX 1

  =1

Sspace: SS01

Method: Maximum likelihood (Marquardt)


Included observations: 258

Convergence achieved after 1 iteration


Coefficient

Std. Error

z-Statistic

Prob.

C(1)

-1.694572

0.025524

-66.39032

0.0000

C(2)

0.997666

0.003013

331.1242

0.0000


Final State

Root MSE

z-Statistic

Prob.

SV1

5.886231

0.428577

13.73437

0.0000

Log likelihood

-150.2963

Akaike info criterion

1.180591

Parameters

2

Schwarz criterion

1.208134

Diffuse priors

0

Hannan-Quinn criter.

1.191666

Unknown  

Sspace: SS01

Method: Maximum likelihood (Marquardt)


Included observations: 258

Convergence achieved after 15 iterations


Coefficient

Std. Error

z-Statistic

Prob.

C(1)

-1.695059

0.025506

-66.45642

0.0000

C(2)

0.997660

0.003009

331.5723

0.0000


Final State

Root MSE

z-Statistic

Prob.

SV1

5.886193

0.428472

13.73763

0.0000

Log likelihood

-150.2964

Akaike info criterion

1.180592

Parameters

2

Schwarz criterion

1.208134

Diffuse priors

0

Hannan-Quinn criter.

1.191667



  =0.9

Sspace: SS01

Method: Maximum likelihood (Marquardt)


Included observations: 258

Convergence achieved after 1 iteration


Coefficient

Std. Error

z-Statistic

Prob.

C(1)

-1.694995

0.025518

-66.42405

0.0000

C(2)

0.997670

0.003014

330.9585

0.0000


Final State

Root MSE

z-Statistic

Prob.

SV1

5.886252

0.428486

13.73733

0.0000

Log likelihood

-150.2964

Akaike info criterion

1.180592

Parameters

2

Schwarz criterion

1.208134

Diffuse priors

0

Hannan-Quinn criter.

1.191667


  =0.8

Sspace: SS01

Method: Maximum likelihood (Marquardt)


Included observations: 258

Convergence achieved after 1 iteration


Coefficient

Std. Error

z-Statistic

Prob.

C(1)

-1.694879

0.025515

-66.42797

0.0000

C(2)

0.997664

0.003011

331.3113

0.0000


Final State

Root MSE

z-Statistic

Prob.

SV1

5.886217

0.428511

13.73645

0.0000

Log likelihood

-150.2963

Akaike info criterion

1.180592

Parameters

2

Schwarz criterion

1.208134

Diffuse priors

0

Hannan-Quinn criter.

1.191667

  =0.7

Sspace: SS01

Method: Maximum likelihood (Marquardt)


Included observations: 258

Convergence achieved after 1 iteration


Coefficient

Std. Error

z-Statistic

Prob.

C(1)

-1.694806

0.025520

-66.41028

0.0000

C(2)

0.997668

0.003014

331.0630

0.0000


Final State

Root MSE

z-Statistic

Prob.

SV1

5.886240

0.428526

13.73600

0.0000

Log likelihood

-150.2963

Akaike info criterion

1.180592

Parameters

2

Schwarz criterion

1.208134

Diffuse priors

0

Hannan-Quinn criter.

1.191667

  =0.5

Sspace: SS01

Method: Maximum likelihood (Marquardt)


Included observations: 258

Convergence achieved after 1 iteration


Coefficient

Std. Error

z-Statistic

Prob.

C(1)

-1.694716

0.025520

-66.40697

0.0000

C(2)

0.997666

0.003012

331.1881

0.0000


Final State

Root MSE

z-Statistic

Prob.

SV1

5.886227

0.428546

13.73536

0.0000

Log likelihood

-150.2963

Akaike info criterion

1.180592

Parameters

2

Schwarz criterion

1.208134

Diffuse priors

0

Hannan-Quinn criter.

1.191666

  =0.3

Sspace: SS01

Method: Maximum likelihood (Marquardt)


Included observations: 258

Convergence achieved after 1 iteration


Coefficient

Std. Error

z-Statistic

Prob.

C(1)

-1.694646

0.025523

-66.39575

0.0000

C(2)

0.997667

0.003013

331.0847

0.0000


Final State

Root MSE

z-Statistic

Prob.

SV1

5.886236

0.428561

13.73489

0.0000

Log likelihood

-150.2963

Akaike info criterion

1.180591

Parameters

2

Schwarz criterion

1.208134

Diffuse priors

0

Hannan-Quinn criter.

1.191666

  =0

Sspace: SS01

Method: Maximum likelihood (Marquardt)


Included observations: 258

Convergence achieved after 1 iteration


Coefficient

Std. Error

z-Statistic

Prob.

C(1)

-1.694508

0.025527

-66.38215

0.0000

C(2)

0.997667

0.003013

331.0771

0.0000


Final State

Root MSE

z-Statistic

Prob.

SV1

5.886235

0.428590

13.73394

0.0000

Log likelihood

-150.2963

Akaike info criterion

1.180591

Parameters

2

Schwarz criterion

1.208134

Diffuse priors

0

Hannan-Quinn criter.

1.191666



1 PhD, Researcher, Romanian Academy- Institute for Economic Forecasting, Address: No. 13, Calea 13 Septembrie, District 5, 76-117, Bucharest, Romania, Corresponding author: mihaela_mb1@yahoo.com.

AUDŒ, Vol. 10, no. 1, pp. 148-159

Refbacks

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