Arima Model.ppt u82

ARIMA Modelling and Forecasting

Introduction • Describe the stationarity of the AR process • Determine the mean and variance of the AR process • Assess the importance of Box-Jenkins methodology • Describe the various types of forecasts • Evaluate the measures of forecast importance

Stationarity of the AR process • If an AR model is not stationary, this implies that previous values of the error term will have a non-declining effect on the current value of the dependent variable. • This implies that the coefficients on the MA process would not converge to zero as the lag length increases. • For an AR model to be stationary, the coefficients on the corresponding MA process decline with lag length, converging on 0.

AR Process • The test for stationarity in an AR model (with p lags)is that the roots of the characteristic equation lie outside the unit circle (i.e. > 1), where the characteristic equation is:

1  1z  2 z  ...  p z 2

p

0

Unit Root • When testing for stationarity for any variable, we describe it as testing for a ‘unit root’, this is based on this same idea. • The most basic AR model is the AR(1) model, on which most tests for stationarity are based, such as the Dickey-Fuller test (covered later)

Unit Root Test yt  yt 1  ut yt  Ly t  ut (1  L) yt  ut (1  z )  0 (characteristic equation)

Unit Root Test • With the AR(1) model, the characteristic equation of (1-z)= 0, suggests that it has a root of z = 1. This lies on the unit circle, rather than outside it, so we conclude that it is nonstationary. • As we increase the lags in the AR model, so the potential number of roots increases, so for 2 lags, we have a quadratic equation producing 2 roots, for the model to be stationary, they both need to lie outside the unit circle.

AR Process • The (unconditional mean) for an AR(1) process, with a constant (μ) is given by:

yt    1 yt 1  ut  AR(1)

 E ( yt )  1  1

AR Process • The (unconditional) variance for an AR process of order 1 (excluding the constant) is:

var( yt ) 

1  1



2

2 (1  1 )

Box-Jenkins Methodology • This is a method for estimating ARIMA models, based on the ACF and PACF as a means of determining the stationarity of the variable in question and the lag lengths of the ARIMA model. • Although the ACF and PACF methods for determining the lag length in an ARIMA model are commonly used, there are other methods termed information criteria which can also be used (covered later)

Box-Jenkins • The Box-Jenkins approach typically comprises four parts: - Identification of the model - Estimation, usually OLS -Diagnostic checking (mostly for autocorrelation) - Forecasting

ACF and PACF • • • • • • • • • • • • • • • • • • • •

Autocorrelation

. |*. .*| . .|. .*| . .|. . |*. . |*. .*| . .|. .*| . .|. . |**

| | | | | | | | | | | |

Partial Correlation

. |*. .*| . .|. .*| . .|. . |*. . |*. .*| . .|. .*| . . |*. . |*.

| | | | | | | | | | | |

1 2 3 4 5 6 7 8 9 10 11 12

AC

PAC

Q-Stat

Prob

0.093 -0.073 -0.023 -0.106 -0.003 0.140 0.162 -0.103 -0.047 -0.124 0.028 0.217

0.093 -0.082 -0.008 -0.111 0.016 0.124 0.141 -0.128 -0.003 -0.117 0.088 0.166

0.7580 1.2307 1.2772 2.3010 2.3017 4.1334 6.6260 7.6356 7.8513 9.3624 9.4434 14.231

0.384 0.540 0.735 0.681 0.806 0.659 0.469 0.470 0.549 0.498 0.581 0.286

Identification • Identification of the most appropriate model is the most important part of the process, where it becomes as much ‘art’ as ‘science’. • The first step is to determine if the variable is stationary, this can be done with the correlogram. If it is not stationary it needs to be first-differenced. (it may need to be differenced again to induce stationarity) • The next stage is to determine the p and q in the ARIMA (p, I, q) model (the I refers to how many times the data needs to be differenced to produce a stationary series)

Identification • To determine the appropriate lag structure in the AR part of the model, the PACF or Partial correlogram is used, where the number of nonzero points of the PACF determine where the AR lags need to be included. • To determine the MA lag structure, the ACF or correlogram is used, again the non-zero points suggest where the lags should be included. • Seasonal dummy variables may also need to be included if the data exhibits seasonal effects.

Diagnostic Checks • With this approach we only test for autocorrelation usually, using the Q or LjungBox statistic. • If there is evidence of autocorrelation, we need to go back to the identification stage and respecify the model, by adding more lags. • A criticism of this approach is that it fails to identify if the model is too big or overparameterised, it only tells us if it is too small.

ARIMA Example • Following the Box-Jenkins methodology, the following ARIMA(2,1,1) model was produced:

yˆ t  0.7  0.6yt 1  0.3yt 2  0.1ut 1 (0.1) (0.2)

(0.1)

R  0.12, LB(2)  3.21 2

(0.01)

Parsimonious Model • The aim of this type of modelling is to produce a model that is parsimonious, or as small as possible, whilst ing the diagnostic checks. • A parsimonious model is desirable because including irrelevant lags in the model increases the coefficient standard errors and therefore reduces the t-statistics. • Models that incorporate large numbers of lags, tend not to forecast well, as they fit data specific features, explaining much of the noise or random features in the data.

Forecasting • One of the most important tests of how well a model performs is how well it forecasts. • One of the most useful models for forecasting is the ARIMA model. • To produce dynamic forecasts the model needs to include lags of either the variables or error .

Types of Forecast • Forecasts can be either in-sample or outof-sample forecasts. • In general the out-of sample forecasts are a better test of how well the model works, as the forecast uses data not included in the estimation of the model. • To conduct out-of-sample forecasts, we need to leave some observations at the end of our sample for this purpose

Types of Forecasts • A one-step-ahead is a forecast for the next observation only. • A multi-step-ahead forecast is for 1,2,3,…s steps ahead. • A recursive window for the forecast means that the initial estimation date is fixed but the additional observations are added one by one to the estimation time span. • A rolling window is where the estimation time period is fixed but the start and end dates successively increase by 1.

Conditional Forecasting • A conditional expectation is one that is taken for time t + 1, conditional upon or given all information available up to and including time t (this is important later). It is written as:

E ( yt 1 | t ) t  all information at time t

Measuring Forecast Accuracy • To determine how accurate a forecast is, the simplest method is to plot the forecast against the actual values as a direct comparison • In addition it may be worthwhile to compare the turning points, this is particularly important in finance. • There are a number of methods to determine accuracy of the forecast, often more than one is included in a set of results.

Tests of Forecast Accuracy • Tests of forecast accuracy are based on the difference between the forecast of the variables value at time t and the actual value at time t. The closer the two are together and the smaller the forecast error, the better the forecast. • There are a variety of statistics measuring this accuracy, mostly based on an average of the errors between the actual and forecast values at time t. • However these statistics provide little information in isolation, as they are unbounded from above and their value depends on the units of the variables being forecast.

Mean Squared Error (MSE) • The MSE statistic can be defined as: T 1 2 MSE  ( y  f )  t s t ,s T  (T1  1) t T1

T  total sample size T1  first out - of - sample forecastobservation f t ,s  s - step - ahead forecastsmade at time t yt  s  the actual value at time t

MSE Example Steps Ahead 1 2 3

Forecast

Actual

0.1 0.25 0.5

0.15 0.20 0.40

Squared Error 0.0025 0.0025 0.01

MSE  (0.0025  0.0025  0.01) / 3  0.005

Forecast Accuracy • There are a number of other measures used and these include: – Mean Absolute Error – Mean Average Prediction Error – Chows test for predictive failure – Theil’s U-statistic (where the forecast is compared to that of a benchmark model) – Root Mean Square Error ( the square root of the MSE)

Theil’s U-statistic

u

T



(

t T1 (

yt  s  f t ,s

)2

xt  s yt  s  fbt ,s xt  s

)

fbt ,s  forecastfrom a benchmarkmodel U  1  both models forecast with equal accuracy U  1  model forecastis best U  1  benchmarkmodel forecastsbest

Financial Loss Functions • When using financial data, it is not always the accuracy of the forecast that matters, very often it is the accuracy of forecasts relating to the sign of future returns or turning points that is important, as this decides if a profit is made. • It is possible to estimate statistics that give the % of correct sign predictions and correct direction change predictions • These statistics work on the principle that if the value of the forecast and actual value are the same, when multiplied together they produce a positive sign. Similarly with the change in forecast and actual values giving the correct direction change.

Correct Sign Predictions • The formula for the % of correct sign predictions is as follows: T

1 % correctsigns  zt  s  T - (T1  1) t T1 Z t  s  1if ( yt  s * f t ,s )  0 Z t  s  0 otherwise

Conclusion • When using AR models, whether the series is stationary of not determine how stable it is. • Box-Jenkins methodology is part art and part science. • Forecasting of time series is an important measure of how well a model works • There are many measures of how accurate a forecast is, usually a number are calculated to determine if the forecast is acceptable, although all have faults.