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

Forecast Intervals for Inflation Rate

and Unemployment Rate in Romania



Mihaela Simionescu1



Abstract: The main objective of this research is to construct forecast intervals for inflation and unemployment rate in Romania. Two types of techniques were employed: bootstrap technique (t-percentile method) and historical error technique (root mean square error method- RMSE). The forecast intervals based on point forecasts of National Bank of Romania (NBR) include more actual values of quarterly inflation rate during Q1:2011-Q4:2013. The proposed prediction intervals for quarterly inflation and unemployment rate contain the registered values. Considering as constant the error from previous year, we will build forecast intervals for annual inflation and unemployment rate based predictions provided by two anonymous experts on the horizon 2004-2015.All the forecast intervals for inflation rate based on first expert expectations included the actual values during 2004-2013.

Keywords: forecast intervals; point prediction; inflation rate; unemployment rate

JEL Classification: C51; C53



1 Introduction

The point forecasts did not provide any information regarding the degree of accuracy. On the other hand, the forecast intervals allow the evaluation of future uncertainty and the comparison between the forecasting methods, indicating the strategies to be applied for desired results.

The main aim of this paper is to construct forecast intervals for inflation and unemployment rate in Romania. Excepting some forecast intervals proposed by (Bratu, 2012, p. 146), in Romania prediction intervals for macroeconomic variables were not proposed. The most frequently used method for constructing forecast intervals is the historical errors method that supposes keeping constant an accuracy measure. The bootstrap method is also used when the distribution type of the sample is unknown.

A grid bootstrap was used to compute the median without bias by (Gospodinov, 2002, p. 86), but the evolution of the events should be characterized by a high persistence. The main disadvantage is the high volume of computations. (Guan, 2003, p. 79).

The sieve bootstrap technique allows for consistent estimators of conditional repartition in the case of non-parametric prediction intervals (Alonso, Pena and Romo, 2003, p. 182). In Romania there is a strong correlation between inflation and money, making us to believe if the money earning went too far (Croitoru, 2013, p. 6). Therefore, the inflation forecasting should be taken under control.

In this paper we used as forecasting methods to build prediction intervals the historical errors method and the bootstrap technique.

It would be necessary to continue the research and build some Bayesian forecast intervals and to compare the results with those obtained using usual methods. The paper continues with the methodological framework, providing different types of quarterly and annual forecast intervals for the two variables. Moreover, the intervals are constructed using as point forecasts the anticipations of two forecasters during 2004-2015. It seems that first expert generated better inflation forecast intervals than the second one.





2. Methodology

The prediction interval that uses the historical errors method considers that errors follow a normal distribution of zero mean and standard deviation that equals the root mean squared error (RMSE) of the histrorical forecasts. Given a certain level of significance (α), the forecast intervals are built as it follows:

  (1)

  the point prediction of the variable Y given at time t for the period (t+k)

 - quantile α/2 of normal distribution of zero mean and standard deviation equalled to1

The following multiple linear regression is considered:

  (1)

Y- vector (length: nx1)

X- matrix (length: nxp)

 - vector of parameters (length: px1)

u- vector of error terms (length: nx1)

 - residuals ( =Y-X )

 - parameter estimator ( )

The form of bootstrap model is:

  (2)

Y*- vector (length: nx1)

X*- matrix (lenth: nxp)

 - parameter estimator ( )

 random element

The selected sample is:  . The random term from theoretical bootstrap process uses modified residuals:

  (3)

The theoretical process is computed as:

  (4)

i=1,2,..,n

b-order of iteration

 - resampled from  

Given the random variable   ( ), the interval for   considers that   has Student distribution (n-p degrees of freedom). For a level of confidence of (1-2 ) the interval is:

  (5)

  (6)

The percentile-t bootstrap method is based on  estimation. We build a bootstrap table, the values of   are:

  (7)

The percentile-t forecast interval for   is:

  (8)

For the observation with number f of the exogenous variable X, the prediction is calculated using the model (Y- dependent variable):   . Having a normal distribution of the errors and the confidence interval (1-2  the standard prediction interval is:

  (9)

A prediction interval for   is based on forecast error  . The future value  

  (10)

It is based on a retrieval of an empirical distribution of the modified residuals. For replication b, the prediction error is:

  (11)

  (12)

The bootstrap forecast error is:

  (13)

Given the empirical distribution of   ( , the percentiles are employed to determine the bootstrap prediction intervals ( ). The percentile prediction interval is:

  (14)

For percentile-t prediction interval, standard deviation estimator ( ) is :

  (15)

  (16)

The statistic   is determined:

  (17)

The percentile-t prediction interval has the form:

  (18)

















3. Forecast Intervals

Using the quarter point forecasts and the prediction intervals provided by the National Bank of Romania, we built some forecast intervals based on historical errors methods by keeping constant the root mean square error (RMSE) of the previous 4 quarter. The horizon is 2011:Q1-2015:Q4.

Table 1. Forecast intervals for the inflation rate predicted by National Bank of Romania (2011:Q1-2015:Q4)

Quarter

Forecast interval

Forecast interval- historical error method

Point forecast

Actual values


Lower limit

Upper limit

Lower limit

Upper limit



T1:2011

7.48

7.95

-2.96

17.96

7.5

1.013

T2:2011

7.93

8.05

-0.73

16.59

7.93

1.005

T3:2011

3.45

3.58

-1.69

8.59

3.45

0.990

T4:2011

3.14

3.25

-0.99

7.27

3.14

1.012

T1:2012

1.43

2.52

-3.01

7.81

2.40

1.010

T2:2012

1.35

3.44

-5.02

9.10

2.04

1.002

T3:2012

2.46

5.20

-3.40

14.06

5.33

1.022

T4:2012

1.57

4.93

-4.42

14.32

4.95

1.009

T1:2013

1.34

5.33

-3.78

14.16

5.19

1.003

T2:2013

1.04

5.98

-2.91

14.69

5.89

1.003

T3:2013

0.62

4.77

-4.04

11.06

3.51

0.988

T4:2013

0.81

5.18

-2.18

9.20

3.51

1.006

T1:2014

1.02

7.93

-1.24

6.60

2.68


T2:2014

1.5

3.5

-1.14

6.70

2.78


T3:2014

1.5

3.5

-0.84

7.00

3.08


T4:2014

1.5

3.5

-0.73

7.11

3.19


T1:2015

1.5

3.5

-1.72

6.12

2.2


T2:2015

1.5

3.5

-2.12

5.72

1.8


T3:2015

1.5

3.5

-1.32

6.52

2.6


T4:2015

1.5

3.5

-1.12

6.72

2.8


Source: own computations

In the period from 2011 to 2013 only two forecast intervals of NBR include the actual values of inflation rate. The prediction intervals based on historical RMSE contain all the actual values during 2011-2013.

The variables with quarterly data that are used are: index of consumer prices that will be used in computing inflation rate, real exchange rate and unemployment rate on the period 2000:Q-2014:Q4. The quarterly forecasts will be made for 2011-2014, after the aggregation of data for obtaining annual values. The Tramo/Seats method was applied to get seasonally adjusted data. The logarithm was applied for the index of consumer prices. The data in first difference was computed for unemployment rate and exchange rate (d_ur and d_er).

The seasonally adjusted and stationarized index of consumer prices is denoted by log_ip. The following valid model was obtained:

  (19)

(std. error=0,08) (std. error=0,02)

(t-calc.=13.62) (t-calc.=-11.24)

According to Breusch-Godfrey test for the first lag, the errors are independent. The hypothesis of errors normal distrbution is checked using Jarque-Bera test and we do not have enought evidence to reject the normal repartition. According to White test, he errors are homoskedastic. The results of the application of these test are presented in Appendix 1.

For the seasonally adjusted and first differentiated quarterly unemployment rate (ur) an autoregressive model of order 1 is built, for which the errors are independent, homoskedastic and they follow a normal repartition (Appendix 2).

  (20)

Table 2. Point forecasts and bootstraped forecast intervals using the linear regression model for inflation rate (%) (percentile-t method) (horizon: 2011:Q1-2015:Q4)

Quarter

Point forecasts

Forecast intervals for inflation rate

Actual values



Intervals limits


Q1:2011

1.0114

0.0245

1.9983

1.013

Q2:2011

1.0088

0.0219

1.9957

1.005

Q3:2011

1.0059

0.0190

1.9928

0.990

Q4:2011

1.0075

0.0206

1.9944

1.012

Q1:2012

1.0044

0.0175

1.9913

1.010

Q2:2012

1.0032

0.0163

1.9901

1.002

Q3:2012

1.0012

0.0143

1.9881

1.022

Q4:2012

1.0047

0.0178

1.9916

1.009

Q1:2013

1.0035

0.0166

1.9904

1.003

Q2:2013

1.0029

0.0160

1.9898

1.003

Q3:2013

1.0031

0.0162

1.9900

0.988

Q4:2013

1.0032

0.0163

1.9901

1.006

Q3:2014

1.0034

0.0165

1.9903


Q4:2014

1.0033

0.0164

1.9902


Q3:2014

1.002

0.0151

1.9889


Q4:2014

1.0021

0.0152

1.9890


Q3:2015

1.002

0.0151

1.9889


Q4:2015

1.0013

0.0144

1.9882


Q3:2015

1.0012

0.0143

1.9881


Q4:2015

1.001

0.0141

1.9879


Source: authors’ computations

As we can see in the table above, the inferior and superior limits of the bootstrap intrvals have ranges with low variations. The results are close of the desired monetary policy in Romania, but the intervals are too narrow and the registered inflation rate for inflation is located out of these intervals. The reasons for this fact are related to the underestimated point forecasts for inflation based on linear regression model. All the forecast intervala based on percentile-t method include the actual values of inflation rate.



Table 3. Point forecasts and bootstraped forecast intervals using the linear regression model for unemployment rate (%) (percentile-t method) (horizon: 2011:Q1-2015:Q4)

Quarter

Point forecasts

Forecast intervals for inflation rate

Actual values



Intervals limits


Q1:2011

7.21

5.46

8.95

7.20

Q2:2011

7.27

5.52

9.01

7.40

Q3:2011

7.41

5.66

9.15

7.40

Q4:2011

7.41

5.66

9.15

7.40

Q1:2012

7.37

5.63

9.12

7.30

Q2:2012

7.21

5.47

8.96

7.00

Q3:2012

7.04

5.29

8.78

7.10

Q4:2012

7.07

5.33

8.82

7.00

Q1:2013

7.07

5.32

8.81

7.20

Q2:2013

7.24

5.49

8.98

7.30

Q3:2013

7.31

5.56

9.05

7.30

Q4:2013

7.31

5.56

9.05

7.30

Q1:2014

7.33

5.59

9.07

7.20

Q2:2014

7.4

5.66

9.14

7.20

Q3:2014

7.41

5.67

9.15


Q4:2014

7.43

5.69

9.17


Q1:2015

7.45

5.71

9.19


Q2:2015

7.45

5.71

9.19


Q3:2015

7.5

5.76

9.24


Q4:2015

7.53

5.79

9.27


Source: authors’ computations

Starting with 2013, the unemployment rate has a slow tendency of increase. The variations of range for forecast intervals for unemployment rate are rather small, because the differencies between predicted unemployment are low from a quarter to another. All the forecast intervals based on percentile-t method include the actual values of unemployment rate.

Table 4. Point forecasts and forecast intervals for qurterly inflation rate and unemployment rate (%) based on historical error methods (horizon: 2011:Q1-2015:Q4)

Quarter

Forecast intervals of inflation rate based on historical RMSE of the previous 4 quarters

Forecast intervals of unemployment rate based on historical RMSE of the previous 4 quarters


Intervals limits

Intervals limits

Q1:2011

-9.448

11.471

6.86

7.55

Q2:2011

-7.655

9.673

6.93

7.60

Q3:2011

-4.130

6.142

7.03

7.78

Q4:2011

-3.121

5.136

6.92

7.89

Q1:2012

-4.407

6.415

6.83

7.91

Q2:2012

-6.057

8.063

6.65

7.77

Q3:2012

-7.731

9.734

6.46

7.61

Q4:2012

-8.365

10.375

6.58

7.57

Q1:2013

-7.963

9.970

6.58

7.55

Q2:2013

-7.799

9.805

6.79

7.68

Q3:2013

-6.551

8.557

6.91

7.70

Q4:2013

-4.687

6.694

6.95

7.66

Q1:2014

-2.917

4.923

7.10

7.56

Q2:2014

-2.917

4.923

7.09

7.71

Q3:2014

-2.918

4.922

7.03

7.79

Q4:2014

-2.919

4.920

7.00

7.86

Q1:2015

-2.921

4.917

6.96

7.94

Q2:2015

-2.923

4.915

7.05

7.67

Q3:2015

-2.928

4.912

7.34

7.87

Q4:2015

-2.929

4.911

7.56

7.96

Source: authors’ computations

Forecasts of inflation and unemployment rate provided by this method seem reasonable,the lenght of intervals being rather big. However, if we go in time, these intervals become narrower. All the forecast intervala based on historical error method include the actual values of inflation and unemployment rate.

Considering constant the error from previous year, we will build forecast intervals for inflation and unemployment rate based on two experts’ predictions on the horizon 2004-2015. Some point forecasts are provided by (Dobrescu, 2013, p. 10).

Table 5. Prediction intervals for annual inflation rate (%) based on historical errors method (horizon: 2004-2015)

Year

Forecast intervals based on first expert forecasts

Forecast intervals based on second expert predictions

Actual inflation rate

2004

3.99

18.97

5.24

18.56

15.3

2005

10.13

17.35

3.32

14.68

11.9

2006

7.82

9.38

3.08

10.92

9

2007

3.90

7.42

9.4

8.06

6.56

2008

1.33

15.67

6.03

11.7

4.84

2009

1.19

10.01

7.93

11.07

7.85

2010

4.81

7.99

5.00

7.40

5.59

2011

3.17

7.04

8.29

9.931

6.09

2012

2.15

6.85

5.37

10.263

3.3

2013

-3.19

12.93

-4.58

2.13

3.98

2014

-3.194

4.806

-07.18

2.10


2015

-3.628

5.638

-8.18

2.2201


Source: authors’ computations

The intervals range for inflation rate is extremly variable in the period 2004-2012. The range is larger during 2013-2015. All the forecast intervals based on first expert anticipations include the actual values of inflation rate while only 5 out of 10 intervals on the horizon 2004-2013 contain the second expert prognosis.

Table 6.Forecast intervals for annual unemployment rate (%) based on historical errors method (horizon: 2004-2015)

Year

Forecast intervals based on first expert forecasts

Forecast intervals based on second expert predictions

Actual unemployment rate

2004

6.808

7.592

6.8240

9.1760

7.4

2005

4.754

11.066

4.7640

11.0360

6.3

2006

4.748

9.452

4.0760

11.5240

5.9

2007

1.638

11.282

0.5440

14.6560

4

2008

3.536

7.064

1.5200

13.2800

4.4

2009

3.400

13.200

3.3040

13.4960

5.8

2010

6.636

10.164

7.2040

7.5960

7.5

2011

6.604

7.812

6.3240

8.6760

6.9

2012

4.748

9.452

4.2680

10.9320

5.9

2013

3.136

6.664

3.4320

6.5680

7.3

2014

5.836

9.364

5.4320

8.5680


2015

5.945

9.567

5.4734

8.5834


Source: authors’ computations

In 2007 the highest range for prediction intervals was obtained for both experts. 9 out of 10 forecast intervals based on first expert anticipations and the second one predictions include the actual values of inflation rate during 2004-2013. For the last year in the horizon both forecasters anticipated lower unemployment rates.









4. Conclusion

The forecast intervals are a way to reflect the uncertainty that affects the forecasting process. For inflation rate and unemployment rate point predictions forecast intervals were built for Romania, providing a better framework for establishing the decision making process. The annual inflation rate forecasts of the first expert anticipation generated precise prediction intervals when bootstrapping and historical errors methods are applied during 2004-2013. However, the intervals are quite large. A future direction of research would be the construction of forecast intervals using Bayesian method.



5. Acknowledgement

This article is a result of the project POSDRU/159/1.5/S/137926, Routes of academic excellence in doctoral and post-doctoral research, being co-funded by the European Social Fund through The Sectorial Operational Programme for Human Resources Development 2007-2013, coordinated by The Romanian Academy.



6. References

Alonso, M., Pena, D. & Romo, J. (2000). Sieve Bootstrap Prediction Intervals. Proceedings in Computational Statistics 14th Symposium, pp. 181-186, Utrecht.

Bratu, M. (2012). Forecast Intervals for Inflation in Romania. Timisoara Journal of Economics, 5(1 (17)), pp. 145-152.

Croitoru, L. (2013). What Good is Higher Inflation? To Avoid or Escape the Liquidity Trap. Romanian Journal of Economic Forecast, Vol. 16, No. 3, pp. 5-25.

Dobrescu, E. (2013). Updating the Romanian Economic Macromodel. Journal for Economic Forecasting, Vol. 16, No. 4, pp. 5-31.

Gospodinov, N. (2002). Median unbiased forecasts for highly persistent autoregressive processes. Journal of Econometrics, Vol. 111, No. 1, pp. 85-101.

Guan, W. (2003). From the help desk: bootstrapped standard errors. The Stata Journal, Vol. 3, No. 1, pp. 71–80.



APPENDIX 1. Linear regression model for quarterly index of consumer prices



Variable

Coefficient

Std. Error

t-Statistic

Prob.

C

0.119341

0.008761

13.62264

0.0000

Curs_schimb__SA

-0.026202

0.002331

-11.24136

0.0000

R-squared

0.700613

Mean dependent var

0.022474

Adjusted R-squared

0.695068

S.D. dependent var

0.021395

S.E. of regression

0.011814

Akaike info criterion

-6.003922

Sum squared resid

0.007537

Schwarz criterion

-5.931588

Log likelihood

170.1098

F-statistic

126.3683

Durbin-Watson stat

1.032398

Prob(F-statistic)

0.000000



White Heteroskedasticity Test:

F-statistic

1.284795

Probability

0.285184

Obs*R-squared

2.589492

Probability

0.273967



Breusch-Godfrey Serial Correlation LM Test:

F-statistic

6.08290

Probability

0.191

Obs*R-squared

3.03713

Probability

0.305





APPENDIX 2. Autoregressive model for quarterly unemployment rate

Variable

Coefficient

Std. Error

t-Statistic

Prob.

C

0.005242

0.042195

0.124224

0.9016

U(-1)

0.309934

0.131379

2.359076

0.0221

R-squared

0.096677

Mean dependent var

0.009259

Adjusted R-squared

0.079305

S.D. dependent var

0.322881

S.E. of regression

0.309814

Akaike info criterion

0.530642

Sum squared resid

4.991194

Schwarz criterion

0.604308

Log likelihood

-12.32734

F-statistic

5.565238

Durbin-Watson stat

1.894308

Prob(F-statistic)

0.022112



Breusch-Godfrey Serial Correlation LM Test:

F-statistic

1.422747

Probability

0.238472

Obs*R-squared

1.465554

Probability

0.226049



White Heteroskedasticity Test:

F-statistic

0.352616

Probability

0.704547

Obs*R-squared

0.736532

Probability

0.691933



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

AUDŒ, Vol. 10, no. 5, pp. 39- 51

Refbacks

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