Econometrics for Policy Analysis A Train The Trainer Workshop Oct 22-28, 2016
|
|
- Juliet Bruce
- 5 years ago
- Views:
Transcription
1 Econometrics for Policy Analysis A Train The Trainer Workshop Delivered by Dr. Nathaniel E. Urama Department of Economics, University of Nigeria, Nsukka
2 Loading Time Series data in E-views: Review For the purpose of this training, we will make use of E-views. The first thing therefore is to load the time series data in E-views. Steps: Save the data in Excel (xlsx or csv comma delimited). Open the E-view software; For Version 8 and above, click on Open Foreign files, select the file, open, next, next, finish. Double click any of the series to cross check. If you see NA in place of data values, your loading was wrong, if you see the data values, you are good to go!! Delivered by Dr. Nathaniel E. Urama 2
3 Loading Time Series data in E-views: Review cont. For Versions below 8, Click on file, new, workfile On the Workfile create Dialogue box below, select the frequency, Enter start date e.g.(1970) oor 1970Q1 and the end date, type the workfile name e.g. (TTT1), ok. Delivered by Dr. Nathaniel E. Urama 3
4 Loading Time Series data in E-views: Review cont. On the workfile TTT1 user interface, click; Proc, Import, Read text lotus excel. Select the file containing the data saved in csv comma delimited, Enter the number of series, ok. Delivered by Dr. Nathaniel E. Urama 4
5 Class Activity 1 1. Load the data set of Nigerian Macro Variables provided. 2. Follow the above steps to load the data set in EViews depending your version. 3. Revision on how to: i. manipulate the data like; logging, differencing and lagging a variable; generating another variable that is a mathematical combination of other variables etc. ii. Plotting graph/ Chats of the Series to visualise there statistical properties. iii. Test the variables for Structural breaks iv. Test the variables for Stationarity or level of Integration. v. Test for stationarity Under structural Breaks. Delivered by Dr. Nathaniel E. Urama vi. Cointegration tests 5
6 Univariate Time Series models AR, MA, ARMA, ARIMA, ARMAX, ARIMAX or ARFIMA? What are these Models? Compare to structural Models a-theoretical Why use them? Serial correlation: presence of correlation between residuals and their lagged values OLS is no longer efficient among linear estimators. Standard errors computed using the textbook OLS formula are not correct, and are generally understated. If there are lagged dependent variables on the right-hand side of the equation specification, OLS estimates are biased and inconsistent. Lack of data date for structural models ( absence, incorrect frequency, immeasurable) Better for Forecasting Delivered by Dr. Nathaniel E. Urama 6
7 AR(p) Model: A review An autoregressive model is one where the current value of a variable depends upon only its previous values and a white noise error term. AR p Y t = α 1 Y t 1 + α 2 Y t α p Y t p + ε t p = j=1 α j Y t j + ε t (1) Using a lag operator L such that L k Y t = Y t k p The AR(p) is given as: Y t = j=1 α j L j Y t + ε t p Y t j=1 p The term: 1 j=1 the AR model. α j L j p Y t = ε t or 1 j=1 α j L j Y t = ε t (2) α j L j is the characteristics polynomial of AR 1 is given as 1 αl Y t = ε t (3) Delivered by Dr. Nathaniel E. Urama 7
8 MA(p) Model : A review An MA(q) model a linear combination of white noise processes, so that yt depends on the current and previous values of a white noise disturbance term. MA q Y t = ε t + θ 1 ε t 1 + θ 2 ε t θ p ε t q q = ε t + j=1 = θ(l)ε t θ j ε t j (4) q For θ L = 1 + j=1 θ j L j (5) q The term: 1 + j=1 θ j L j is the characteristics polynomial of the MA model. 8 where ε t are the independent and identically distributed innovations for the process Delivered by Dr. Nathaniel E. Urama
9 MA(p) Model : A review The distinguishing properties of the moving average process of order q given above are: 1. E y t = μ 2. var y t = γ 0 = 1 + θ θ θ q 2 σ 2 3. cov y t = γ s = 1 + θ θ θ 2 q σ 2 for s = 1,2,, q 0 for s > q 9 Delivered by Dr. Nathaniel E. Urama
10 ARMA model: A review ARMA p, q is a combination of Ar(p) amd MA(q) as follows: p ARMA p, q Y t = j=1 q α j Y t j + ε t + j=1 θ j ε t j (6) ARMA(1,1) is given as: 1 αl Y t = 1 + θl ε t (7) Seasonal AR and MA Terms: Due to seasonal patterns in most monthly and quarterly data, Box and Jenkins (1976) recommend the use of seasonal autoregressive (SAR) and seasonal moving average (SMA) terms in the ARMA process. SAR p is a seasonal AR term with lag p and it adds to an existing AR, a polynomial with lag p given as 1 p L p : A second order AR process for quarterly data can be written as; 1 α 1 L 1 α 2 L L 4 Delivered by Dr. Nathaniel E. Urama Y t = ε t (8) 10
11 AR, MA and ARMA, : A review (8) on expansion will give: Y t = α 1 Y t 1 + α 2 Y t 2 4 Y t 4 α 1 4 Y t 5 α 2 4 Y t 6 + ε t (9) For seasonal moving average with lag q, the resulting MA lag structure is obtained from the product of the lag polynomial specified by the MA terms and the one specified by any SMA terms. For a second order MA without seasonality, the process is written as : Y t = ε t + θ 1 ε t 1 + θ 2 ε t 2 2 = ε t + j=1 θ j ε t j (10) This in the lag form is given as: Y t = 1 + θ 1 L 1 θ 2 L 2 ε t (11) Delivered by Dr. Nathaniel E. Urama 11
12 AR, MA and ARMA, : A review If the data for (11) is quarterly for example, we introduce the SMA(4) given as 1 + φ 4 L 4 in the MA term. This will give: Y t = 1 + θ 1 L 1 θ 2 L φ 4 L 4 ε t (12) (12) on Expansion will give: Y t = ε t + θ 1 ε t 1 + θ 2 ε t 2 + φ 4 ε t 4 + θ 1 φ 4 ε t 5 + θ 2 φ 4 ε t 6 (13) The parameter φ is associated with the seasonal part of the MA process. Delivered by Dr. Nathaniel E. Urama 12
13 ARIMA and ARIMAX models The AR, MA and ARMA models discussed before assumes that the series in question is at least weakly stationary. Why? (see Gujarati, 2004, pp. 840) Since most time series are not stationary, there is need to account for this in our ARMA model. Hence, the need for ARIMA model. In our previous class, a series that must be differenced d times for it to become stationary is said to integrated of order d i.e. I(d) ARIMA (p,d,q) is an ARMA(p,q) model of non-stationary series differenced d times to make it stationary. Delivered by Dr. Nathaniel E. Urama 13
14 Estimating ARIMA models:the BJ [Box Jenkins] Methods Revisited How do we identify the value of P, d and q for an ARIMA(p, d, q) models? The BJ methodology has an answer and consists of the following steps: Differencing to achieve Stationarity Identification Estimation Diagnostic Checking Forecasting Delivered by Dr. Nathaniel E. Urama 14
15 The BJ [Box Jenkins] Methods Differencing the Series to achieve stationarity Identify the model to be tentatively estimated Estimate the parameters of the tentative model Diagnostic Checking of the adequacy of the model Is the model adequate? No Yes Use the model for Forecasting and control Delivered by Dr. Nathaniel E. Urama 15
16 The BJ [Box Jenkins] Methods cont. 1. Differencing to achieve stationarity: If a series is integrated of order d, it has to be differenced d times to make it stationary. 1. Identification: This involves determining the order of the model required to capture the dynamic features of the data. 2. The chief tools are autocorrelations (AC), the and partial autocorrelations (PAC), and the correlogram of the series. Covariance at lag K AC at lag k is given as ρ k = With γ k = t=1 n k Y t Y n Y t+k Y, variance = γ k γ 0 γ 0 = Y t Y 2 n A plot of ρk against k is known as the sample correlogram. Note if n is small, use n-1 as the divisor. Asymptotically however, both give same result 16 Delivered by Dr. Nathaniel E. Urama
17 Finding the value of d in ARIMA (p, d, q) Model The autocorrelation is the correlation coefficient of the current value of the series with the series lagged a certain number of periods. For the purely white noise(stationary) process the autocorrelations bars/values are not significantly different from zero at various lags. For pure RWM, the autocorrelation coefficients at various lags are statistically different from zero. Generating AC and PAC in Eviews: Click Quick, series statistics, correlogram, on the dialogue box that will appear, select level, ok Or Double click the series, View, correlogram, Level, Enter the lag length OK. Delivered by Dr. Nathaniel E. Urama 17
18 Correlogram of a white noise random process Delivered by Dr. Nathaniel E. Urama 18
19 Correlogram of a random walk time series Delivered by Dr. Nathaniel E. Urama 19
20 Finding the value of d in ARIMA (p, d, q) Model cont. If a time series exhibits white noise, the sample autocorrelation coefficients ρ k are approximately ρ k ~N 0, 1 n Hence, following the properties of the standard normal distribution, the 95% confidence interval for any (population) ρ k is: ρ k ± n If the preceding interval includes the value of zero, we do not reject the hypothesis that the true ρk is zero. If this interval does not include 0, we reject the hypothesis that the true ρk is zero. Alternatively, we can use the asymptotic confidence limit for the ρ k which is ± Delivered by Dr. Nathaniel E. Urama
21 Finding the value of d in ARIMA (p, d, q) Model cont. If the autocorrelation coefficients ρ k are statistically zero, the value of d=0. Otherwise, it may not be zero. Run the correlogram for the 1 st, 2 nd, until the dth higher order difference of the series where autocorrelation coefficients ρ k becomes statistically zero. Steps in Eviews: Double click the series, View, correlogram, select the order and enter the lag length, ok For series needing more than 2 nd differencing, Quick, series statistics, correlogram, type D(Y, n) as the series name where y is the series name and d is the number of differencing, select level, enter the lag length, ok. 21 Delivered by Dr. Nathaniel E. Urama
22 Finding the value of d in ARIMA (p, d, q) Model cont. For series having seasonal terms, click Quick, series statistics, correlogram type d(y, n, s) as the series name where y is the series, with seasonal difference at lag s, and d the number of differencing, level, enter the lag length, ok. d y, n, s = d 1 L 2 (1 L s )y More accurately, conduct the structural break and the stationarity test to determine the order of the integration of the series which is the d. Why Structural Break? The presence of structural breaks in a series may lead to wrong conclusion that the series is not stationary. See the example below for the annual series of OER in our data set. Delivered by Dr. Nathaniel E. Urama 22
23 The BJ [Box Jenkins] Methods cont. Correlogram of OER below show that it is not stationary Delivered by Dr. Nathaniel E. Urama 23
24 The BJ [Box Jenkins] Methods cont. Correlogram of OER below show that it is stationary after the first difference i.e., the D may be 1 Delivered by Dr. Nathaniel E. Urama 24
25 The BJ [Box Jenkins] Methods cont. The ADF unit root test confirm that ORE is I(1) Delivered by Dr. Nathaniel E. Urama 25
26 Finding the value of d in ARIMA (p, d, q) Model cont. The ADF results are confirmed by the KPSS results below. Delivered by Dr. Nathaniel E. Urama 26
27 Finding the value of d in ARIMA (p, d, q) Model cont. A test of the series OER for structural breaks shows the rejection of null of no structural break as shown in the next 2 Slides below. Delivered by Dr. Nathaniel E. Urama 27
28 Finding the value of d in ARIMA (p, d, q) Model cont. A test of the series OER for structural Using the multiple breakpoint test. Delivered by Dr. Nathaniel E. Urama 28
29 Finding the value of d in ARIMA (p, d, q) Model cont. A Look at the series graph reveals the breaks. Delivered by Dr. Nathaniel E. Urama 29
30 Finding the value of d in ARIMA (p, d, q) Model cont. Zivo-Andrews unit root test shows a rejection of the null hypothesis that OER has a unit root with a structural break. This shows that OER may be I(0). Delivered by Dr. Nathaniel E. Urama 30
31 Finding the value of d in ARIMA (p, d, q) Model cont. In E-views 9 and above, one can conduct unit root with structural break test directly. The result as shown below also show that considering the structural breaks, OER may be I(0). Delivered by Dr. Nathaniel E. Urama 31
32 Finding the values of p and q in ARIMA (p, d, q) Model The PAC measures the correlation of Y values that are K periods apart after removing the correlation from the intermediate lags. This means that it measures the added predictive power of the kth lag in a model containing Y t-1 to Y t-k terms as the regressors. The last two columns reported in the correlogram are the Ljung-Box Q-statistics and their p-values. 2 m ρ Ljung-Box (LB) = n(n + 2) k k=1 n k ~χ2 m The Q-statistics test the joint hypothesis that all the ρk up to certain lags are simultaneously equal to zero. if the computed Q exceeds the critical Q value, in the Chi Square distribution, reject the Oct 22-28, null 2016 hypothesis that all the (true) ρks are zero; at least some of them must be nonzero. 32 Delivered by Dr. Nathaniel E. Urama
33 The correlogram of the ARIMA Model The correlogram view compares the autocorrelation pattern of the structural residuals and that of the estimated model The structural residuals are the residuals after removing the effect of the fitted exogenous regressors but not the ARIMA terms For a properly specified model, the residual and theoretical (estimated) autocorrelations and partial autocorrelations should be close. Steps in E-view: From the ARIMA Diagnostic Views dialog box, select correlogram, enter number of lags, graph, OK, to display the graph or, table, OK, to displays the table Delivered by Dr. Nathaniel E. Urama 33
34 Autocorrelation Function (ACF) and Correlogram How does a sample correlogram enable us to find out if a particular series is AR(p) MA(p) or ARMA(p,q)? Delivered by Dr. Nathaniel E. Urama 34
35 BJ Identification Methods: An e.g. AC and PAC of AR(1), MA(1) and ARMA(1, 1) and ARIMA(1, 1,1) for the Annual OER Septs for AR(1) in Eviews: Quick, estimate equation, enter the equation as oer c ar 1, ok. Or d(oer) c ar 1, Ok if the series is I(1). On the equation dialogue box (DB), click on View, ARMA structure. On the ARMA Diagnostic views DB, click Correlogram, graph, ok Delivered by Dr. Nathaniel E. Urama
36 AC and PAC of AR(1) process of Annual OER From E-views 9 Delivered by Dr. Nathaniel E. Urama 36
37 AC and PAC of MA(1) process of OER Delivered by Dr. Nathaniel E. Urama 37
38 AC and PAC of ARMA(1, 1) process of OER From E-views 9 Delivered by Dr. Nathaniel E. Urama 38
39 AC and PAC of ARMA(1, 1, 1) process of Annual OER Delivered by Dr. Nathaniel E. Urama 39
40 AC and PAC of AR(1) process of Quarterly OER From E-views 8 From E-views 9 Delivered by Dr. Nathaniel E. Urama 40
41 AC and PAC of MA(1) process of Quarterly OER Delivered by Dr. Nathaniel E. Urama 41
42 AC and PAC of ARMA(1,1,1) process of Quarterly OER From E-views 8 Delivered by Dr. Nathaniel E. Urama 42
43 ARIMA (p, d, q) Model Selection: Information Criteria Since the actual data ACF and PACF rarely exhibits the theoretical pattern shown, identification will be difficult using them An alternative is to use the information criteria namely: Akaike (AIC), Schwarz (SBIC) and Hannan-Queen (HQIC). AIC = ln σ 2 SBIC = ln σ 2 HQIC = ln σ 2 + 2k T + k ln T T + 2k T ln(ln( T)) Where σ 2 is the residual variance = RSS T The objective is the select the model with the least IC 43 Delivered by Dr. Nathaniel E. Urama
44 ARIMA (p, d, q) Selection: Information Criteria Cont. Which criterion should be preferred if they suggest different model order? SBIC embodies much stiffer penalty term than the AIC. SBIC is strongly consistent but inefficient. AIC is not consistent but generally efficient. Hence, SBIC asymptotically gives the correct model order than the AIC. However, the average variation in selected model orders from different samples within a given population is greater with SBIC than AIC In conclusion then, no criterion is absolutely superior to others. 44 Delivered by Dr. Nathaniel E. Urama
45 BJ Identification Methods: An e.g. Cont. Since the ACFs and PACFs didn t show the exact resemblance of of any of the theoretical AR(p), MA(q), nor ARMA(p, d, q), though closer to ARIMA(1, 1, 1) we resort to the Information criteria method. The Information criteria involves estimating different models and comparing their IC values. Tedious!!! testing up to ARMA(p, q) means running (P + 1)X(Q + 1) Models for non seasonal series and(p + 1)X(Q + 1)X4 models for Quarterly series with seasonal pattern. Thanks to E-views for simplifying the case. Steps in Eview 8 and above: Double click the series, Proc. Add-ins, Automatic ARIMA selection, Select the max AR,and MA terms and the IC, OK. For this particular series, the result is as shown below 45 This shows ARIMA( 1, 1, 1) with seasonal pattern in the 4 th lag of the moving average Delivered by Dr. Nathaniel E. Urama
46 SBIC selection for ARMA(P,Q) for p and q up to 3 each Model Selection Criterion Model Selection Criterion Model Selection Criterion Model Selection Criterion Model Selection Criterion 1 (0,0)(0,0) (0,3)(0,4) (1,2)(0,0) (2,1)(4,0) (3,0)(4,4) (0,0)(0,0) (0,3)(0,0) (1,2)(0,4) (2,1)(4,0) (3,1)(0,0) (0,0)(0,4) (0,3)(0,0) (1,2)(4,0) (2,1)(4,4) (3,1)(0,0) (0,0)(0,0) (0,3)(0,4) (1,2)(4,0) (2,2)(0,0) (3,1)(0,4) (0,0)(0,0) (0,3)(4,0) (1,2)(4,4) (2,2)(0,0) (3,1)(0,0) (0,0)(0,4) (0,3)(4,0) (1,3)(0,0) (2,2)(0,4) (3,1)(0,0) (0,0)(4,0) (0,3)(4,4) (1,3)(0,0) (2,2)(0,0) (3,1)(0,4) (0,0)(4,0) (1,0)(0,0) (1,3)(0,4) (2,2)(0,0) (3,1)(4,0) (0,0)(4,4) (1,0)(0,0) (1,3)(0,0) (2,2)(0,4) (3,1)(4,0) (0,1)(0,0) (1,0)(0,4) (1,3)(0,0) (2,2)(4,0) (3,1)(4,4) (0,1)(0,0) (1,0)(0,0) (1,3)(0,4) (2,2)(4,0) (3,2)(0,0) (0,1)(0,4) (1,0)(0,0) (1,3)(4,0) (2,2)(4,4) (3,2)(0,0) (0,1)(0,0) (1,0)(0,4) (1,3)(4,0) (2,3)(0,0) (3,2)(0,4) (0,1)(0,0) (1,0)(4,0) (1,3)(4,4) (2,3)(0,0) (3,2)(0,0) (0,1)(0,4) (1,0)(4,0) (2,0)(0,0) (2,3)(0,4) (3,2)(0,0) (0,1)(4,0) (1,0)(4,4) (2,0)(0,0) (2,3)(0,0) (3,2)(0,4) (0,1)(4,0) (1,1)(0,0) (2,0)(0,4) (2,3)(0,0) (3,2)(4,0) (0,1)(4,4) (1,1)(0,0) (2,0)(0,0) (2,3)(0,4) (3,2)(4,0) (0,2)(0,0) (1,1)(0,4) (2,0)(0,0) (2,3)(4,0) (3,2)(4,4) (0,2)(0,0) (1,1)(0,0) (2,0)(0,4) (2,3)(4,0) (3,3)(0,0) (0,2)(0,4) (1,1)(0,0) (2,0)(4,0) (2,3)(4,4) (3,3)(0,0) (0,2)(0,0) (1,1)(0,4) (2,0)(4,0) (3,0)(0,0) (3,3)(0,4) (0,2)(0,0) (1,1)(4,0) (2,0)(4,4) (3,0)(0,0) (3,3)(0,0) (0,2)(0,4) (1,1)(4,0) (2,1)(0,0) (3,0)(0,4) (3,3)(0,0) (0,2)(4,0) (1,1)(4,4) (2,1)(0,0) (3,0)(0,0) (3,3)(0,4) (0,2)(4,0) (1,2)(0,0) (2,1)(0,4) (3,0)(0,0) (3,3)(4,0) (0,2)(4,4) (1,2)(0,0) (2,1)(0,0) (3,0)(0,4) (3,3)(4,0) Delivered by Dr. Nathaniel E. Urama (0,3)(0,0) (1,2)(0,4) (2,1)(0,0) (3,0)(4,0) (3,3)(4,4) (0,3)(0,0) (1,2)(0,0) (2,1)(0,4) (3,0)(4,0)
47 Estimating the ARIMA Model Quick Estimate equation Select LS- Least Square(NLS and ARMA) Type D(OER) ar(1) ma(1) Sma(4) Delivered by Dr. Nathaniel E. Urama 47
48 Result of ARIMA(1, 1, 1) SMA(4) Model Delivered by Dr. Nathaniel E. Urama 48
49 AC and PAC for ARIMA(1, 1, 1) SMA(4) Model Delivered by Dr. Nathaniel E. Urama 49
50 Class Activity Determine the values of p, d and q ARIMA (p, d, q) model of the following Nigerian series: CPI, INF, RIR., Ycul, using, 1. Correlagram 2. Information Criteria 7.2 Estimate the ARIMA (p, d, q) model of the above Nigerian series and use the AC and PAC to check the correctness f the Model Order. Delivered by Dr. Nathaniel E. Urama 50
51 ARMA Diagnostic Checking. Over fitting : Deliberately fitting a larger model than that identified to see if all other extra terms would be insignificant. Need for Parsimonious Model Inclusion of irrelevant variables will increase the coefficient standard error and decrease the parameter significance. Whether this happens of not depends on how much the RSS falls and the size of T(number of observations) relative to K(number of Parameters). If T is large, relative to K, the decrease in RSS will outweigh the decrease in Degrees of freedom (T-K) and hence, coefficient standard error will increase. In line with Brooks, Models that are profligate might be inclined to fit the 51 data specific features, which would not be replicated out-of-sample. Delivered by Dr. Nathaniel E. Urama
52 ARIMA Diagnostic Check: Overfitting The added ar(2) and ar(3) are all insignificant. Delivered by Dr. Nathaniel E. Urama 52
53 ARMA Diagnostic Checking Cont. Residual diagnostics: Checks the residual for evidence of linear dependence which, if present, suggest that the model is inadequate. The correlogram ( ACF and PACF) If the acf and pacf confirms stationarity of the residual, the model is correctly specified, Steps in E-view: View/Residual Diagnostics/Correlogram-Qstatistic Delivered by Dr. Nathaniel E. Urama 53
54 ARIMA Diagnostic Check: Correlogram Delivered by Dr. Nathaniel E. Urama 54
55 ARMA Diagnostic Checking Cont. Residual diagnostics through the serial Correlation LM test If the LM test fail to reject the null hypothesis of no serial Correlation, the model is correctly specified. Steps in E-view: Ensure that ML is bot the Used ARMA method. In Estimate DB, click option and ensure that ML is not the selected method View/Residual Diagnostics/Serial Correlation LM Test 55 Delivered by Dr. Nathaniel E. Urama
56 ARIMA Diagnostic Check: LM Test Delivered by Dr. Nathaniel E. Urama 56
57 Notes on BJ diagnostic tests It essentially involves only autocorrelation tests I can only find under and not over parameterized model(brooks. 2008, P.231) 57 Delivered by Dr. Nathaniel E. Urama
58 Further ARIMA Equation Diagnostics The Stationarity and the Invertibility condition The Stationarity Condition Needed to avoid the previous values of the error term having nondeclining effect on the current values of Y t as the time progresses. Given a general ARMA model in the lag polynomial form α L and θ L as α L Y t = θ L ε t The reported roots are the inverse roots of the characteristics polynomials: α x 1 = 0, and θ x 1 = 0 58 Which may be imaginary, but should have modulus not greater than 1 Delivered by Dr. Nathaniel E. Urama
59 Further ARIMA Diagnostic Tests Cont. The above is same as the roots of the characteristics equation, 1 α 1 L 1 α 2 L 2 α p L p = 0 all being outside the unit circle. (see Brooks, 2008, p.216) The roots view displays the inverse roots of the AR and/or MA characteristic polynomial, either as a graph or as a table. The graph view plots the inverse roots in the complex plane the horizontal axis is the real part and the vertical axis is the imaginary part of each root. Delivered by Dr. Nathaniel E. Urama 59
60 Further ARIMA Diagnostic Tests Cont. If the estimated ARMA process is (covariance) stationary, all the AR inverse roots should be inside the unit circle. If the estimated ARMA process is invertible, all MA inverse roots should be inside the unit circle. If the AR has a real Inverse roots with r > 1, a pair of complex reciprocal roots with modulus greater than one, it means that the autoregressive process is explosive. If the MA has reciprocal roots outside the unit circle, we say that the MA process is noninvertible, which makes interpreting and using the MA results difficult. Delivered by Dr. Nathaniel E. Urama 60
61 Further ARIMA Diagnostic Tests Cont. If the MA has roots with Modulus close to 1, It is a sign that we have over differenced the series. It makes estimation and forecasting difficult. According to Hamilton (1994a, p. 65) however, noninvertibility is not a serious problem since there is always an equivalent representation of the MA model with the inverse roots inside the unit circle. The Iteration should continue with different staring values until invertibility is achieved. Delivered by Dr. Nathaniel E. Urama 61
62 Further ARIMA Diagnostic Tests Cont. Steps in E-view: On the ARIMA Equation DB, click view, ARMA structure and the ARMA diagnostic views below will appear Select Roots, Graph, ok and the graph in the next slide will appear. Delivered by Dr. Nathaniel E. Urama 62
63 The Inverse Roots Diagnostic test for ARIMA Delivered by Dr. Nathaniel E. Urama 63 All the Inverse roots of both the AR and MA lies inside the unit circle. The AR is stationary and the MA is ivertible
64 The Inverse Roots of AR/MA in a Tabular form Both the AR and the MA has no root outside the unit circle. Hence, AR is stationary and MA is invertible. However, the MA has roots close to 1, a sign that we may have over differenced the Delivered by Dr. Nathaniel E. Urama 64 series.
65 The Inverse Roots Diagnostic test for ARMA Reducing the number of differencing by 1 gives us: an ARIMA(1, 0, 1) model or an ARMA(1, 1) model with the result as below: Delivered by Dr. Nathaniel E. Urama 65
66 Inverse roots of the AR/MA Polynomial for the ARMA(1,1) All the Inverse roots of both the AR and MA lies inside the unit circle. Delivered by Dr. Nathaniel E. Urama 66
67 Inverse roots of the AR/MA Polynomial for the ARMA(1,1) All the Inverse roots of both the AR and MA still lies inside the unit circle but the roots of the MA are now far from How ever, since the value of the IC criteria are all smaller for the differenced series, we will adopt it. Delivered by Dr. Nathaniel E. Urama
68 Impulse Response of an ARIMA model Traces the response of the ARIMA part of the estimated equation to one time shock in the innovation. The accumulated response is the response to step impulse where the same shock occurs in every period from the first. If the estimated ARMA model is stationary, the impulse responses will asymptote to zero, while the accumulated responses will asymptote to its longrun value. Steps in E-view: From the ARIMA Diagnostic Views dialog box, select Impulse Response, enter number of periods, graph, OK, to display the graph or, table, 68 OK, to display the table. Delivered by Dr. Nathaniel E. Urama
69 Impulse Response of our ARIMA model Delivered by Dr. Nathaniel E. Urama 69
70 Class Activity 8 Carry out a diagnostic check on the ARIMA models you estimated during Class activity 7 and interpret them. Conduct an impulse response on all of them. End of day 6 Assignment: Read Brooks section 5.9 and adapt the model there to test if Uncovered Interest rate Parity held in Nigeria between 1980 and Oct 22-28, Delivered by Dr. Nathaniel E. Urama
71 Forecasting with ARIMA and ARIMAX models In-sample vs Out-of Sample Dynamic vs static Forecast Accuracy Decomposition of the Forecast Errors Class Activity 9: Day 7: Course Details Multivariate models VAR, Parsimonious VAR (PVAR), VEC, ARDL Series stationarity and the VAR model Stability and Stationarity in VAR VAR lag length Selection, Impulse Response and Variance Decomposition Forecasting with VAR In-sample vs Out-of Sample Dynamic vs static VAR diagnostic Checking Class activity 10: Delivered by Dr. Nathaniel E. Urama 71
72 Forecasting in Econometrics What is Forecasting? The role of forecasting in Economics and Policy Decision Forecasting Accuracy Vs Statistical Significance of a model Structural Forecasting: Forecasting from Parametric/theoretical models return predictions derived from arbitrage pricing models, long-term exchange rate prediction based on purchasing power parity or uncovered interest parity theory. 72 Delivered by Dr. Nathaniel E. Urama
73 Forecasting in Econometrics Time series Forecasting: Forecasting from a-theoretic/non-parametric model. forecast the future values of a series given its previous values and/or previous values of an error term. Is there a clear distinction? Which one is better? Point and Interval Forecasts. In-sample vs Out-of Sample forecasts Dynamic Vs static forecasts Rolling Vs Recursive forecasts 73 Delivered by Dr. Nathaniel E. Urama
74 Dynamic Vs Static forecasts The dynamic and static is all about how the lagged value of Y that appears on the right-hand side of the equation should be evaluated Dynamic Calculates, multi-step forecasts from the first period in the forecast sample. From an AR p model Y t = α 0 + α 1 Y t 1 + α 2 Y t α p Y t p + ε t The Forecast for Y t is gotten from the actual values of the lagged Y: p Y t = α 0 + i=1 α i Y t i Then forecast of Y t+s uses previously forecasted values for the lagged dependent variables. s Y t+s = α 0 + i=1 p α i Y t+s i + i=s+1 α i Y t+s i 74 Delivered by Dr. Nathaniel E. Urama
75 Dynamic Vs Static forecasts Dynamic Forecast : Is true multi-step and recursive forecasts. It uses the recursively computed forecast of the lagged value of the dependent variable. requires that data for the exogenous variables be available for every observation in the forecast sample, only available when the estimated equation contains dynamic components that values for any lagged dependent variables be observed at the start of the forecast sample In Dynamic Forecast, any missing values for the explanatory variables will generate an NA for that observation and in all subsequent observation (E-views User guide) 75 Delivered by Dr. Nathaniel E. Urama
76 Static Forecast The Static Forecast calculates a sequence of one-step ahead forecasts, using the actual, rather than forecasted values for lagged dependent variables, if available. From an AR p model Y t = α 0 + α 1 Y t 1 + α 2 Y t α p Y t p + ε t The Forecast for both Y t and Y t+s are gotten from the actual values of the lagged Y if available: p Y t+s = α 0 + i=1 α i Y t+s i 76 Delivered by Dr. Nathaniel E. Urama
77 Rolling Vs Recursive forecasts Rolling Fix the length of the in-sample period used to estimate the model so that the start and the end date increases successively by one observation. Recursive Fix the initial estimation date but additional observation are added one at a time to the estimation period (Brooks, 2008, p. 246) 77 Delivered by Dr. Nathaniel E. Urama
78 Forecasting with MA(q) model A moving average process has a memory only of length q, and this limits the sensible forecasting horizon Given an MA q : Y t = θ 0 + ε t + θ 1 ε t 1 + θ 2 ε t θ q ε t q * With the assumption of constant Parameter, What hold at time t will also hold at time t + s. Hence, Y t+1 = θ 0 + ε t+1 + θ 1 ε t + θ 2 ε t θ q ε t q+1 ** Y t+2 = θ 0 + ε t+2 + θ 1 ε t+1 + θ 2 ε t + + θ q ε t q+2 *** 78 Y t+q = θ 0 + ε t+q + θ 1 ε t+q 1 + θ 2 ε t+q θ q ε t q+q Delivered by Dr. Nathaniel E. Urama
79 Forecasting with MA(q) model Cont. The forecasts of Y t+1 given as f t,1 is gotten thus: f t,1 = E(Y t+1 Ωt ) = E θ 0 + ε t+1 + θ 1 ε t + θ 2 ε t θ q ε t q+1 Ω t Where E(Y t+1 Ωt ) means the expected value of Y at time t + 1 given all the information available up to and including time t. Since E(ε t+s Ωt ) = 0 s > 0, f t,1 = E(Y t+1 Ωt ) = θ 0 + θ 1 ε t + θ 2 ε t θ q ε t q+1 f t,q = E(Y t+q Ωt ) = θ 0 + θ q ε t f t,s = E(Y t+s Ωt ) = θ 0 s q + 1 What the last two equations show is that the MA(q) process has a memory of only q periods and all forecasts above q periods ahead collapse to the intercept if there is intercept, 0 otherwise. Delivered by Dr. Nathaniel E. Urama 79
80 Forecasting with AR(p) model Unlike a moving average process, an autoregressive process has infinite memory (see Brooks, 2008, p ). (Explain this) Given an AR p : Y t = α 0 + α 1 Y t 1 +α 2 Y t α p Y t p + ε t With the assumption of constant Parameter, What hold at time t will also hold at time t + s. Hence, Y t+1 = α 0 + ε t+1 + α 1 Y t +α 2 Y t α q Y t p+1 Y t+2 = α 0 + ε t+2 + α 1 Y t+1 +α 2 Y t + + α q Y t p+2 Y t+p = α 0 + ε t+p + α 1 Y t+p 1 +α 2 Y t+p α p Y t p+p 80 Delivered by Dr. Nathaniel E. Urama
81 Forecasting with AR(p) model Cont. The forecasts of Y t+1 given as f t,1 is : f t,1 = E(Y t+1 Ωt ) = E α 0 + ε t+1 + α 1 Y t + α 2 Y t α p Y t p+1 Ω t = α 0 +α 1 Y t + α 2 Y t α p Y t p+1 f t,2 = E(Y t+2 Ωt ) = E α 0 + ε t+2 + α 1 Y t+1 + α 2 Y t + + α p Y t p+2 Ω t = α 0 + α 1 E Y t+1 t + α 2 E Y t + + α p E Y t p+2 = α 0 + α 1 f t,1 + α 2 y t + + α p Y t p+2 f t,p = E(Y t+p Ωt ) = E α 0 + ε t+p + α 1 Y t+p 1 + α 2 Y t+p α p Y t+p p Ω t = α 0 + α 1 E Y t+p 1 t + α 2 E Y t+p α p E Y t+p p = α 0 + α 1 f t,p 1 + α 2 f t,p α p Y t Etc. These confirms that the AR(p) process has infinite memory 81 Delivered by Dr. Nathaniel E. Urama
82 Forecasting with ARIMA models Let f t,s be a forecast from ARMA(P, q) model at time t for s steps into the future. p Then, f t,s = i=1 q α i f t,s i + j=1 θ j ε t+s j The above is known as Forecast function Where f t,s = y t+s,s 0 ; and ε t+s = 0, s > 0 = ε t+s, S 0 α i and θ j are the autoregressive and moving average coefficients, respectively (Brooks, 2008, p. 249) The value of s (the forecast horizon) is determined by the ARIMA structure. Delivered by Dr. Nathaniel E. Urama 82
83 Forecasting OER with AR(1) process in E-view Estimate D(OER) with Ar(1) process: In the Command window, Click Quick, estimate equation In the equation estimation DB, type: D(oer) c ar(1) Since our workfile data rage is 1972Q3 to 2012Q4, for us to be able to obtain out-ofsample forecast, we reduce the sample range in the equation estimation DB to : 1972Q3 to 2011Q4, click ok. In the result window, click forecast. The Forecast Db will appear. In the Forecast DB, Select the series to forecast (will this be in all cases?), type the name for the forecast series, Select the method, and the forecast sample, ok Delivered by Dr. Nathaniel E. Urama 83
84 Dynamic Forecasting OER with AR(1) process in E-view Forecast: OER_IN_DF Actual: OER Forecast sample: 1972Q3 2012Q4 Included observations: 162 Root Mean Squared Error Mean Absolute Error Mean Abs. Percent Error Theil Inequality Coefficient Bias Proportion Variance Proportion Covariance Proportion OER_IN_DF ± 2 S.E. 84 Delivered by Dr. Nathaniel E. Urama
85 Static Forecasting OER with AR(1) process in E-view Forecast: OER_OUT_SF Actual: OER Forecast sample: 2005Q3 2012Q4 Included observations: 30 Root Mean Squared Error Mean Absolute Error Mean Abs. Percent Error Theil Inequality Coefficient Bias Proportion Variance Proportion Covariance Proportion OER_OUT_SF Delivered by Dr. Nathaniel E. Urama ± 2 S.E.
86 Forecasting OER with AR(1) process in E-view To see the difference between dynamic and Static static forecast, increase the workfile Dynamic range to 2013Q4 but having data only up to 2012Q4 and try to forecast to 2013Q4. Delivered by Dr. Nathaniel E. Urama 86
87 162 Forecasting OER with AR(1) process in E-view IV I II III IV I II III IV I II III IV Delivered by Dr. Nathaniel E. Urama 87 OER OER_AR_DF OER_AR_SF
88 Forecasting OER with MA(1) process in E-view Estimate D(OER) with MA(1) process: In the Command window, Click Quick, estimate equation In the equation estimation DB, type: D(oer) c ma(1) sma(4) reduce the sample range in the equation estimation DB to 1972Q3 to 2011Q4, click ok. Delivered by Dr. Nathaniel E. Urama 88
89 Forecasting OER with MA(1) process in E-view The dynamic Dynamic forecast covers all the forecast horizon Staticwhile the static forecast only one period above the available data 89 Delivered by Dr. Nathaniel E. Urama
90 162 Forecasting OER with MA(1) process in E-view IV I II III IV I II III IV I II III IV OER OER_MA_DF OER_MA_SF 90 Delivered by Dr. Nathaniel E. Urama
91 162 Forecasting OER with MA(1) process in E-view IV I II III IV I II III IV I II III IV OER OER_AR_SF OER_MA_SF OER_AR_DF OER_MA_DF Delivered by Dr. Nathaniel E. Urama 91
92 Forecasting OER with ARIMA(1,1,1) process in E-view Table of ARIAM forecast result 92 Delivered by Dr. Nathaniel E. Urama
93 62 Forecasting OER with ARIMA(1,1,1) process in E-view Gaph of ARIAM forecast result IV I II III IV I II III IV I II III IV OER OER_ARIMA_DF OER_ARIMA_SF Delivered by Dr. Nathaniel E. Urama 93
94 Forecast Accuracy The level of accuracy of a forecast is determined using the forecast error (FE). The FE for a given observation is the difference between the actual value and its forecast. It will therefore be positive if the forecast is too low and negative if the forecast is too high. Thus, before the forecast errors are aggregated, they are usually squared or the absolute value taken (Brooks, 2008 p ) as in the next slide. Delivered by Dr. Nathaniel E. Urama 94
95 Forecast Accuracy criteria Mean Square Error MSE = 1 T T1 1 T t=t1 Y t+s f t,s 2 Mean Absolute Error MAE = 1 T T1 1 T t=t1 Y t+s f t,s Mean Absolute % Error MAPE = 100 T T1 1 T t=t1 Y t+s f t,s Y t+s Adjusted MAPE AMAPE = 100 T T1 1 T t=t1 Y t+s f t,s Y t+s +f t,s T is total sample size (in-sample + out-of-sample) T 1 is the first out-of-sample forecast observation. 95 Delivered by Dr. Nathaniel E. Urama
96 Decision with the forecast accuracy result MSE, MAE and MAPE can not be used to judge a single forecast: They are unbounded from above. The only way is to compare results from different models for same data and forecast period and choose the model with the lowest value Which of the criterion is better? Use MSE when large forecast errors are disproportionately more serious than smaller errors Authors like Makridakis (1993, p. 528) see MAPE as the best. Brooks (2008, p. 253) argues that compared to MAPE, AMAPE is better(except when the series and the forecasts can take on opposite signs) as, it corrects for the problem of asymmetry between the actual and forecast values If they Y t+s and the f t,s are equal and opposite, the denominator in AMAPE turns to 0 Delivered by Dr. Nathaniel E. Urama 96
97 Theil s U-statistic Theil s U-statistic is another criterion for judging forecast accuracy. It is given as: U = t=t T 1 T t=t1 Y t+s f 2 t,s Y t+s Y t+s fb 2 t,s Y t+s fb t,s is the forecast obtained from a benchmark model (typically a simple model such as a naive or random walk). U = 1 implies that both the model and the benchmark are equally in(accurate) U < 1 implies that the model is superior to the benchmark U > 1 implies that the model is inferior to the benchmark Delivered by Dr. Nathaniel E. Urama 97
98 Decomposition of the forecast errors The mean squared forecast error can be decomposed into: a bias proportion, a variance proportion and a covariance proportion. The bias component measures the extent to which the mean of the forecasts is different from the mean of the actual data the variance component measures the difference between the variation of the forecasts and the variation of the actual data, while the covariance component captures any remaining unsystematic part of the forecast errors. Accurate forecasts would be unbiased and also have a small variance proportion, so that most of the forecast error should be attributable to the covariance, see Pindyck and Rubinfeld (1998, p ). Delivered by Dr. Nathaniel E. Urama 98
99 Statistical versus financial or economic loss functions. The statistical loss functions vs practical usefulness. Models with lower forecast error may not necessarily be more profitable when used in business forecasting Rather, Leitch and Tanner, (1991) as cited in Brooks (2008) argue that models that can accurately forecast the sign, or predict turning points in a series are more profitable. These are computed as: % correct sign prediction = For z t+s = 1 if Y t,s f t,s > 0, 0 otherwise 1 T T T 1 1 t=t 1 z t+s % correct direction of change predictions = For z t+s = 1 if Y t,s Y t, f t,s Y t, > 0, 0 otherwise 1 T T T 1 1 t=t 1 z t+s 99 Delivered by Dr. Nathaniel E. Urama
100 Class activity 9 1. Carry out dynamic and static forecast of 8 period ahead values of the series modeled in class activity Model the same variable with ARMA(1,1) process and compare their statistical loss function 3. Compute both the % correct sign prediction and the % correct direction of change prediction for the two models and compare them. 100 Delivered by Dr. Nathaniel E. Urama
Econometrics for Policy Analysis A Train The Trainer Workshop Oct 22-28, 2016 Organized by African Heritage Institution
Econometrics for Policy Analysis A Train The Trainer Workshop Oct 22-28, 2016 Organized by African Heritage Institution Delivered by Dr. Nathaniel E. Urama Department of Economics, University of Nigeria,
More informationTIME SERIES ANALYSIS AND FORECASTING USING THE STATISTICAL MODEL ARIMA
CHAPTER 6 TIME SERIES ANALYSIS AND FORECASTING USING THE STATISTICAL MODEL ARIMA 6.1. Introduction A time series is a sequence of observations ordered in time. A basic assumption in the time series analysis
More informationUnivariate ARIMA Models
Univariate ARIMA Models ARIMA Model Building Steps: Identification: Using graphs, statistics, ACFs and PACFs, transformations, etc. to achieve stationary and tentatively identify patterns and model components.
More informationat least 50 and preferably 100 observations should be available to build a proper model
III Box-Jenkins Methods 1. Pros and Cons of ARIMA Forecasting a) need for data at least 50 and preferably 100 observations should be available to build a proper model used most frequently for hourly or
More informationBox-Jenkins ARIMA Advanced Time Series
Box-Jenkins ARIMA Advanced Time Series www.realoptionsvaluation.com ROV Technical Papers Series: Volume 25 Theory In This Issue 1. Learn about Risk Simulator s ARIMA and Auto ARIMA modules. 2. Find out
More informationLecture 2: Univariate Time Series
Lecture 2: Univariate Time Series Analysis: Conditional and Unconditional Densities, Stationarity, ARMA Processes Prof. Massimo Guidolin 20192 Financial Econometrics Spring/Winter 2017 Overview Motivation:
More informationProblem Set 2: Box-Jenkins methodology
Problem Set : Box-Jenkins methodology 1) For an AR1) process we have: γ0) = σ ε 1 φ σ ε γ0) = 1 φ Hence, For a MA1) process, p lim R = φ γ0) = 1 + θ )σ ε σ ε 1 = γ0) 1 + θ Therefore, p lim R = 1 1 1 +
More informationEmpirical Market Microstructure Analysis (EMMA)
Empirical Market Microstructure Analysis (EMMA) Lecture 3: Statistical Building Blocks and Econometric Basics Prof. Dr. Michael Stein michael.stein@vwl.uni-freiburg.de Albert-Ludwigs-University of Freiburg
More informationARIMA Modelling and Forecasting
ARIMA Modelling and Forecasting Economic time series often appear nonstationary, because of trends, seasonal patterns, cycles, etc. However, the differences may appear stationary. Δx t x t x t 1 (first
More informationChapter 12: An introduction to Time Series Analysis. Chapter 12: An introduction to Time Series Analysis
Chapter 12: An introduction to Time Series Analysis Introduction In this chapter, we will discuss forecasting with single-series (univariate) Box-Jenkins models. The common name of the models is Auto-Regressive
More informationAutoregressive Moving Average (ARMA) Models and their Practical Applications
Autoregressive Moving Average (ARMA) Models and their Practical Applications Massimo Guidolin February 2018 1 Essential Concepts in Time Series Analysis 1.1 Time Series and Their Properties Time series:
More informationMODELING INFLATION RATES IN NIGERIA: BOX-JENKINS APPROACH. I. U. Moffat and A. E. David Department of Mathematics & Statistics, University of Uyo, Uyo
Vol.4, No.2, pp.2-27, April 216 MODELING INFLATION RATES IN NIGERIA: BOX-JENKINS APPROACH I. U. Moffat and A. E. David Department of Mathematics & Statistics, University of Uyo, Uyo ABSTRACT: This study
More informationCh 6. Model Specification. Time Series Analysis
We start to build ARIMA(p,d,q) models. The subjects include: 1 how to determine p, d, q for a given series (Chapter 6); 2 how to estimate the parameters (φ s and θ s) of a specific ARIMA(p,d,q) model (Chapter
More informationChapter 8: Model Diagnostics
Chapter 8: Model Diagnostics Model diagnostics involve checking how well the model fits. If the model fits poorly, we consider changing the specification of the model. A major tool of model diagnostics
More informationCh 8. MODEL DIAGNOSTICS. Time Series Analysis
Model diagnostics is concerned with testing the goodness of fit of a model and, if the fit is poor, suggesting appropriate modifications. We shall present two complementary approaches: analysis of residuals
More informationProf. Dr. Roland Füss Lecture Series in Applied Econometrics Summer Term Introduction to Time Series Analysis
Introduction to Time Series Analysis 1 Contents: I. Basics of Time Series Analysis... 4 I.1 Stationarity... 5 I.2 Autocorrelation Function... 9 I.3 Partial Autocorrelation Function (PACF)... 14 I.4 Transformation
More informationLecture 3: Autoregressive Moving Average (ARMA) Models and their Practical Applications
Lecture 3: Autoregressive Moving Average (ARMA) Models and their Practical Applications Prof. Massimo Guidolin 20192 Financial Econometrics Winter/Spring 2018 Overview Moving average processes Autoregressive
More informationFE570 Financial Markets and Trading. Stevens Institute of Technology
FE570 Financial Markets and Trading Lecture 5. Linear Time Series Analysis and Its Applications (Ref. Joel Hasbrouck - Empirical Market Microstructure ) Steve Yang Stevens Institute of Technology 9/25/2012
More informationLecture 4a: ARMA Model
Lecture 4a: ARMA Model 1 2 Big Picture Most often our goal is to find a statistical model to describe real time series (estimation), and then predict the future (forecasting) One particularly popular model
More informationMultivariate Time Series Analysis and Its Applications [Tsay (2005), chapter 8]
1 Multivariate Time Series Analysis and Its Applications [Tsay (2005), chapter 8] Insights: Price movements in one market can spread easily and instantly to another market [economic globalization and internet
More informationEASTERN MEDITERRANEAN UNIVERSITY ECON 604, FALL 2007 DEPARTMENT OF ECONOMICS MEHMET BALCILAR ARIMA MODELS: IDENTIFICATION
ARIMA MODELS: IDENTIFICATION A. Autocorrelations and Partial Autocorrelations 1. Summary of What We Know So Far: a) Series y t is to be modeled by Box-Jenkins methods. The first step was to convert y t
More informationAustrian Inflation Rate
Austrian Inflation Rate Course of Econometric Forecasting Nadir Shahzad Virkun Tomas Sedliacik Goal and Data Selection Our goal is to find a relatively accurate procedure in order to forecast the Austrian
More informationUnivariate Time Series Analysis; ARIMA Models
Econometrics 2 Fall 24 Univariate Time Series Analysis; ARIMA Models Heino Bohn Nielsen of4 Outline of the Lecture () Introduction to univariate time series analysis. (2) Stationarity. (3) Characterizing
More informationEconometric Forecasting
Robert M. Kunst robert.kunst@univie.ac.at University of Vienna and Institute for Advanced Studies Vienna October 1, 2014 Outline Introduction Model-free extrapolation Univariate time-series models Trend
More informationClassic Time Series Analysis
Classic Time Series Analysis Concepts and Definitions Let Y be a random number with PDF f Y t ~f,t Define t =E[Y t ] m(t) is known as the trend Define the autocovariance t, s =COV [Y t,y s ] =E[ Y t t
More informationFORECASTING SUGARCANE PRODUCTION IN INDIA WITH ARIMA MODEL
FORECASTING SUGARCANE PRODUCTION IN INDIA WITH ARIMA MODEL B. N. MANDAL Abstract: Yearly sugarcane production data for the period of - to - of India were analyzed by time-series methods. Autocorrelation
More informationA SARIMAX coupled modelling applied to individual load curves intraday forecasting
A SARIMAX coupled modelling applied to individual load curves intraday forecasting Frédéric Proïa Workshop EDF Institut Henri Poincaré - Paris 05 avril 2012 INRIA Bordeaux Sud-Ouest Institut de Mathématiques
More informationECONOMETRIA II. CURSO 2009/2010 LAB # 3
ECONOMETRIA II. CURSO 2009/2010 LAB # 3 BOX-JENKINS METHODOLOGY The Box Jenkins approach combines the moving average and the autorregresive models. Although both models were already known, the contribution
More informationAdvanced Econometrics
Advanced Econometrics Marco Sunder Nov 04 2010 Marco Sunder Advanced Econometrics 1/ 25 Contents 1 2 3 Marco Sunder Advanced Econometrics 2/ 25 Music Marco Sunder Advanced Econometrics 3/ 25 Music Marco
More informationEconometrics II Heij et al. Chapter 7.1
Chapter 7.1 p. 1/2 Econometrics II Heij et al. Chapter 7.1 Linear Time Series Models for Stationary data Marius Ooms Tinbergen Institute Amsterdam Chapter 7.1 p. 2/2 Program Introduction Modelling philosophy
More information{ } Stochastic processes. Models for time series. Specification of a process. Specification of a process. , X t3. ,...X tn }
Stochastic processes Time series are an example of a stochastic or random process Models for time series A stochastic process is 'a statistical phenomenon that evolves in time according to probabilistic
More informationDynamic Time Series Regression: A Panacea for Spurious Correlations
International Journal of Scientific and Research Publications, Volume 6, Issue 10, October 2016 337 Dynamic Time Series Regression: A Panacea for Spurious Correlations Emmanuel Alphonsus Akpan *, Imoh
More informationFinancial Econometrics Review Session Notes 3
Financial Econometrics Review Session Notes 3 Nina Boyarchenko January 22, 2010 Contents 1 k-step ahead forecast and forecast errors 2 1.1 Example 1: stationary series.......................... 2 1.2 Example
More informationThe ARIMA Procedure: The ARIMA Procedure
Page 1 of 120 Overview: ARIMA Procedure Getting Started: ARIMA Procedure The Three Stages of ARIMA Modeling Identification Stage Estimation and Diagnostic Checking Stage Forecasting Stage Using ARIMA Procedure
More informationAutomatic Forecasting
Automatic Forecasting Summary The Automatic Forecasting procedure is designed to forecast future values of time series data. A time series consists of a set of sequential numeric data taken at equally
More informationCircle a single answer for each multiple choice question. Your choice should be made clearly.
TEST #1 STA 4853 March 4, 215 Name: Please read the following directions. DO NOT TURN THE PAGE UNTIL INSTRUCTED TO DO SO Directions This exam is closed book and closed notes. There are 31 questions. Circle
More information5 Autoregressive-Moving-Average Modeling
5 Autoregressive-Moving-Average Modeling 5. Purpose. Autoregressive-moving-average (ARMA models are mathematical models of the persistence, or autocorrelation, in a time series. ARMA models are widely
More informationSection 2 NABE ASTEF 65
Section 2 NABE ASTEF 65 Econometric (Structural) Models 66 67 The Multiple Regression Model 68 69 Assumptions 70 Components of Model Endogenous variables -- Dependent variables, values of which are determined
More informationAPPLIED ECONOMETRIC TIME SERIES 4TH EDITION
APPLIED ECONOMETRIC TIME SERIES 4TH EDITION Chapter 2: STATIONARY TIME-SERIES MODELS WALTER ENDERS, UNIVERSITY OF ALABAMA Copyright 2015 John Wiley & Sons, Inc. Section 1 STOCHASTIC DIFFERENCE EQUATION
More informationModelling Monthly Rainfall Data of Port Harcourt, Nigeria by Seasonal Box-Jenkins Methods
International Journal of Sciences Research Article (ISSN 2305-3925) Volume 2, Issue July 2013 http://www.ijsciences.com Modelling Monthly Rainfall Data of Port Harcourt, Nigeria by Seasonal Box-Jenkins
More informationTESTING FOR CO-INTEGRATION
Bo Sjö 2010-12-05 TESTING FOR CO-INTEGRATION To be used in combination with Sjö (2008) Testing for Unit Roots and Cointegration A Guide. Instructions: Use the Johansen method to test for Purchasing Power
More informationA SEASONAL TIME SERIES MODEL FOR NIGERIAN MONTHLY AIR TRAFFIC DATA
www.arpapress.com/volumes/vol14issue3/ijrras_14_3_14.pdf A SEASONAL TIME SERIES MODEL FOR NIGERIAN MONTHLY AIR TRAFFIC DATA Ette Harrison Etuk Department of Mathematics/Computer Science, Rivers State University
More informationVector Autoregressive Model. Vector Autoregressions II. Estimation of Vector Autoregressions II. Estimation of Vector Autoregressions I.
Vector Autoregressive Model Vector Autoregressions II Empirical Macroeconomics - Lect 2 Dr. Ana Beatriz Galvao Queen Mary University of London January 2012 A VAR(p) model of the m 1 vector of time series
More informationEconometrics of financial markets, -solutions to seminar 1. Problem 1
Econometrics of financial markets, -solutions to seminar 1. Problem 1 a) Estimate with OLS. For any regression y i α + βx i + u i for OLS to be unbiased we need cov (u i,x j )0 i, j. For the autoregressive
More informationRomanian Economic and Business Review Vol. 3, No. 3 THE EVOLUTION OF SNP PETROM STOCK LIST - STUDY THROUGH AUTOREGRESSIVE MODELS
THE EVOLUTION OF SNP PETROM STOCK LIST - STUDY THROUGH AUTOREGRESSIVE MODELS Marian Zaharia, Ioana Zaheu, and Elena Roxana Stan Abstract Stock exchange market is one of the most dynamic and unpredictable
More informationNon-Stationary Time Series and Unit Root Testing
Econometrics II Non-Stationary Time Series and Unit Root Testing Morten Nyboe Tabor Course Outline: Non-Stationary Time Series and Unit Root Testing 1 Stationarity and Deviation from Stationarity Trend-Stationarity
More informationTopic 4 Unit Roots. Gerald P. Dwyer. February Clemson University
Topic 4 Unit Roots Gerald P. Dwyer Clemson University February 2016 Outline 1 Unit Roots Introduction Trend and Difference Stationary Autocorrelations of Series That Have Deterministic or Stochastic Trends
More informationModule 3. Descriptive Time Series Statistics and Introduction to Time Series Models
Module 3 Descriptive Time Series Statistics and Introduction to Time Series Models Class notes for Statistics 451: Applied Time Series Iowa State University Copyright 2015 W Q Meeker November 11, 2015
More informationMidterm Suggested Solutions
CUHK Dept. of Economics Spring 2011 ECON 4120 Sung Y. Park Midterm Suggested Solutions Q1 (a) In time series, autocorrelation measures the correlation between y t and its lag y t τ. It is defined as. ρ(τ)
More informationLab: Box-Jenkins Methodology - US Wholesale Price Indicator
Lab: Box-Jenkins Methodology - US Wholesale Price Indicator In this lab we explore the Box-Jenkins methodology by applying it to a time-series data set comprising quarterly observations of the US Wholesale
More informationApplied time-series analysis
Robert M. Kunst robert.kunst@univie.ac.at University of Vienna and Institute for Advanced Studies Vienna October 18, 2011 Outline Introduction and overview Econometric Time-Series Analysis In principle,
More informationEconometría 2: Análisis de series de Tiempo
Econometría 2: Análisis de series de Tiempo Karoll GOMEZ kgomezp@unal.edu.co http://karollgomez.wordpress.com Segundo semestre 2016 IX. Vector Time Series Models VARMA Models A. 1. Motivation: The vector
More informationCHAPTER 8 FORECASTING PRACTICE I
CHAPTER 8 FORECASTING PRACTICE I Sometimes we find time series with mixed AR and MA properties (ACF and PACF) We then can use mixed models: ARMA(p,q) These slides are based on: González-Rivera: Forecasting
More informationCHAPTER 8 MODEL DIAGNOSTICS. 8.1 Residual Analysis
CHAPTER 8 MODEL DIAGNOSTICS We have now discussed methods for specifying models and for efficiently estimating the parameters in those models. Model diagnostics, or model criticism, is concerned with testing
More informationForecasting using R. Rob J Hyndman. 2.4 Non-seasonal ARIMA models. Forecasting using R 1
Forecasting using R Rob J Hyndman 2.4 Non-seasonal ARIMA models Forecasting using R 1 Outline 1 Autoregressive models 2 Moving average models 3 Non-seasonal ARIMA models 4 Partial autocorrelations 5 Estimation
More informationTime Series Forecasting: A Tool for Out - Sample Model Selection and Evaluation
AMERICAN JOURNAL OF SCIENTIFIC AND INDUSTRIAL RESEARCH 214, Science Huβ, http://www.scihub.org/ajsir ISSN: 2153-649X, doi:1.5251/ajsir.214.5.6.185.194 Time Series Forecasting: A Tool for Out - Sample Model
More informationUnivariate linear models
Univariate linear models The specification process of an univariate ARIMA model is based on the theoretical properties of the different processes and it is also important the observation and interpretation
More informationLecture 5: Unit Roots, Cointegration and Error Correction Models The Spurious Regression Problem
Lecture 5: Unit Roots, Cointegration and Error Correction Models The Spurious Regression Problem Prof. Massimo Guidolin 20192 Financial Econometrics Winter/Spring 2018 Overview Stochastic vs. deterministic
More informationARIMA Models. Richard G. Pierse
ARIMA Models Richard G. Pierse 1 Introduction Time Series Analysis looks at the properties of time series from a purely statistical point of view. No attempt is made to relate variables using a priori
More informationStationarity and Cointegration analysis. Tinashe Bvirindi
Stationarity and Cointegration analysis By Tinashe Bvirindi tbvirindi@gmail.com layout Unit root testing Cointegration Vector Auto-regressions Cointegration in Multivariate systems Introduction Stationarity
More informationCircle the single best answer for each multiple choice question. Your choice should be made clearly.
TEST #1 STA 4853 March 6, 2017 Name: Please read the following directions. DO NOT TURN THE PAGE UNTIL INSTRUCTED TO DO SO Directions This exam is closed book and closed notes. There are 32 multiple choice
More informationA lecture on time series
A lecture on time series Bernt Arne Ødegaard 15 November 2018 Contents 1 Survey of time series issues 1 1.1 Definition.............................. 2 1.2 Examples.............................. 2 1.3 Typical
More informationLATVIAN GDP: TIME SERIES FORECASTING USING VECTOR AUTO REGRESSION
LATVIAN GDP: TIME SERIES FORECASTING USING VECTOR AUTO REGRESSION BEZRUCKO Aleksandrs, (LV) Abstract: The target goal of this work is to develop a methodology of forecasting Latvian GDP using ARMA (AutoRegressive-Moving-Average)
More informationArma-Arch Modeling Of The Returns Of First Bank Of Nigeria
Arma-Arch Modeling Of The Returns Of First Bank Of Nigeria Emmanuel Alphonsus Akpan Imoh Udo Moffat Department of Mathematics and Statistics University of Uyo, Nigeria Ntiedo Bassey Ekpo Department of
More informationAR, MA and ARMA models
AR, MA and AR by Hedibert Lopes P Based on Tsay s Analysis of Financial Time Series (3rd edition) P 1 Stationarity 2 3 4 5 6 7 P 8 9 10 11 Outline P Linear Time Series Analysis and Its Applications For
More informationQuantitative Finance I
Quantitative Finance I Linear AR and MA Models (Lecture 4) Winter Semester 01/013 by Lukas Vacha * If viewed in.pdf format - for full functionality use Mathematica 7 (or higher) notebook (.nb) version
More informationPart 1. Multiple Choice (50 questions, 1 point each) Part 2. Problems/Short Answer (10 questions, 5 points each)
GROUND RULES: This exam contains two parts: Part 1. Multiple Choice (50 questions, 1 point each) Part 2. Problems/Short Answer (10 questions, 5 points each) The maximum number of points on this exam is
More informationEconometrics I: Univariate Time Series Econometrics (1)
Econometrics I: Dipartimento di Economia Politica e Metodi Quantitativi University of Pavia Overview of the Lecture 1 st EViews Session VI: Some Theoretical Premises 2 Overview of the Lecture 1 st EViews
More informationTime Series Econometrics 4 Vijayamohanan Pillai N
Time Series Econometrics 4 Vijayamohanan Pillai N Vijayamohan: CDS MPhil: Time Series 5 1 Autoregressive Moving Average Process: ARMA(p, q) Vijayamohan: CDS MPhil: Time Series 5 2 1 Autoregressive Moving
More informationStochastic Modelling Solutions to Exercises on Time Series
Stochastic Modelling Solutions to Exercises on Time Series Dr. Iqbal Owadally March 3, 2003 Solutions to Elementary Problems Q1. (i) (1 0.5B)X t = Z t. The characteristic equation 1 0.5z = 0 does not have
More informationCovariance Stationary Time Series. Example: Independent White Noise (IWN(0,σ 2 )) Y t = ε t, ε t iid N(0,σ 2 )
Covariance Stationary Time Series Stochastic Process: sequence of rv s ordered by time {Y t } {...,Y 1,Y 0,Y 1,...} Defn: {Y t } is covariance stationary if E[Y t ]μ for all t cov(y t,y t j )E[(Y t μ)(y
More informationTime Series Analysis. James D. Hamilton PRINCETON UNIVERSITY PRESS PRINCETON, NEW JERSEY
Time Series Analysis James D. Hamilton PRINCETON UNIVERSITY PRESS PRINCETON, NEW JERSEY & Contents PREFACE xiii 1 1.1. 1.2. Difference Equations First-Order Difference Equations 1 /?th-order Difference
More informationVector autoregressions, VAR
1 / 45 Vector autoregressions, VAR Chapter 2 Financial Econometrics Michael Hauser WS17/18 2 / 45 Content Cross-correlations VAR model in standard/reduced form Properties of VAR(1), VAR(p) Structural VAR,
More informationForecasting of the Austrian Inflation Rate
Forecasting of the Austrian Inflation Rate Case Study for the Course of Econometric Forecasting Winter Semester 2007 by Nadir Shahzad Virkun Tomas Sedliacik Goal setting and Data selection The goal of
More informationThe Identification of ARIMA Models
APPENDIX 4 The Identification of ARIMA Models As we have established in a previous lecture, there is a one-to-one correspondence between the parameters of an ARMA(p, q) model, including the variance of
More informationMODELLING TIME SERIES WITH CONDITIONAL HETEROSCEDASTICITY
MODELLING TIME SERIES WITH CONDITIONAL HETEROSCEDASTICITY The simple ARCH Model Eva Rubliková Ekonomická univerzita Bratislava Manuela Magalhães Hill Department of Quantitative Methods, INSTITUTO SUPERIOR
More informationNon-Stationary Time Series and Unit Root Testing
Econometrics II Non-Stationary Time Series and Unit Root Testing Morten Nyboe Tabor Course Outline: Non-Stationary Time Series and Unit Root Testing 1 Stationarity and Deviation from Stationarity Trend-Stationarity
More informationAuthor: Yesuf M. Awel 1c. Affiliation: 1 PhD, Economist-Consultant; P.O Box , Addis Ababa, Ethiopia. c.
ISSN: 2415-0304 (Print) ISSN: 2522-2465 (Online) Indexing/Abstracting Forecasting GDP Growth: Application of Autoregressive Integrated Moving Average Model Author: Yesuf M. Awel 1c Affiliation: 1 PhD,
More information1 Linear Difference Equations
ARMA Handout Jialin Yu 1 Linear Difference Equations First order systems Let {ε t } t=1 denote an input sequence and {y t} t=1 sequence generated by denote an output y t = φy t 1 + ε t t = 1, 2,... with
More informationRoss Bettinger, Analytical Consultant, Seattle, WA
ABSTRACT DYNAMIC REGRESSION IN ARIMA MODELING Ross Bettinger, Analytical Consultant, Seattle, WA Box-Jenkins time series models that contain exogenous predictor variables are called dynamic regression
More informationMultivariate Time Series: VAR(p) Processes and Models
Multivariate Time Series: VAR(p) Processes and Models A VAR(p) model, for p > 0 is X t = φ 0 + Φ 1 X t 1 + + Φ p X t p + A t, where X t, φ 0, and X t i are k-vectors, Φ 1,..., Φ p are k k matrices, with
More information7. Integrated Processes
7. Integrated Processes Up to now: Analysis of stationary processes (stationary ARMA(p, q) processes) Problem: Many economic time series exhibit non-stationary patterns over time 226 Example: We consider
More informationNon-Stationary Time Series and Unit Root Testing
Econometrics II Non-Stationary Time Series and Unit Root Testing Morten Nyboe Tabor Course Outline: Non-Stationary Time Series and Unit Root Testing 1 Stationarity and Deviation from Stationarity Trend-Stationarity
More information11.1 Gujarati(2003): Chapter 12
11.1 Gujarati(2003): Chapter 12 Time Series Data 11.2 Time series process of economic variables e.g., GDP, M1, interest rate, echange rate, imports, eports, inflation rate, etc. Realization An observed
More informationForecasting. Simon Shaw 2005/06 Semester II
Forecasting Simon Shaw s.c.shaw@maths.bath.ac.uk 2005/06 Semester II 1 Introduction A critical aspect of managing any business is planning for the future. events is called forecasting. Predicting future
More informationChapter 6: Model Specification for Time Series
Chapter 6: Model Specification for Time Series The ARIMA(p, d, q) class of models as a broad class can describe many real time series. Model specification for ARIMA(p, d, q) models involves 1. Choosing
More informationFöreläsning /31
1/31 Föreläsning 10 090420 Chapter 13 Econometric Modeling: Model Speci cation and Diagnostic testing 2/31 Types of speci cation errors Consider the following models: Y i = β 1 + β 2 X i + β 3 X 2 i +
More informationApplied Econometrics. Applied Econometrics. Applied Econometrics. Applied Econometrics. What is Autocorrelation. Applied Econometrics
Autocorrelation 1. What is 2. What causes 3. First and higher orders 4. Consequences of 5. Detecting 6. Resolving Learning Objectives 1. Understand meaning of in the CLRM 2. What causes 3. Distinguish
More informationEconometrics. 9) Heteroscedasticity and autocorrelation
30C00200 Econometrics 9) Heteroscedasticity and autocorrelation Timo Kuosmanen Professor, Ph.D. http://nomepre.net/index.php/timokuosmanen Today s topics Heteroscedasticity Possible causes Testing for
More informationECON/FIN 250: Forecasting in Finance and Economics: Section 8: Forecast Examples: Part 1
ECON/FIN 250: Forecasting in Finance and Economics: Section 8: Forecast Examples: Part 1 Patrick Herb Brandeis University Spring 2016 Patrick Herb (Brandeis University) Forecast Examples: Part 1 ECON/FIN
More information10. Time series regression and forecasting
10. Time series regression and forecasting Key feature of this section: Analysis of data on a single entity observed at multiple points in time (time series data) Typical research questions: What is the
More informationWORKSHOP. Introductory Econometrics with EViews. Asst. Prof. Dr. Kemal Bağzıbağlı Department of Economic
WORKSHOP on Introductory Econometrics with EViews Asst. Prof. Dr. Kemal Bağzıbağlı Department of Economic Res. Asst. Pejman Bahramian PhD Candidate, Department of Economic Res. Asst. Gizem Uzuner MSc Student,
More informationARIMA modeling to forecast area and production of rice in West Bengal
Journal of Crop and Weed, 9(2):26-31(2013) ARIMA modeling to forecast area and production of rice in West Bengal R. BISWAS AND B. BHATTACHARYYA Department of Agricultural Statistics Bidhan Chandra Krishi
More informationIntroduction to ARMA and GARCH processes
Introduction to ARMA and GARCH processes Fulvio Corsi SNS Pisa 3 March 2010 Fulvio Corsi Introduction to ARMA () and GARCH processes SNS Pisa 3 March 2010 1 / 24 Stationarity Strict stationarity: (X 1,
More informationLECTURE 11. Introduction to Econometrics. Autocorrelation
LECTURE 11 Introduction to Econometrics Autocorrelation November 29, 2016 1 / 24 ON PREVIOUS LECTURES We discussed the specification of a regression equation Specification consists of choosing: 1. correct
More informationECON/FIN 250: Forecasting in Finance and Economics: Section 6: Standard Univariate Models
ECON/FIN 250: Forecasting in Finance and Economics: Section 6: Standard Univariate Models Patrick Herb Brandeis University Spring 2016 Patrick Herb (Brandeis University) Standard Univariate Models ECON/FIN
More informationARIMA Models. Jamie Monogan. January 25, University of Georgia. Jamie Monogan (UGA) ARIMA Models January 25, / 38
ARIMA Models Jamie Monogan University of Georgia January 25, 2012 Jamie Monogan (UGA) ARIMA Models January 25, 2012 1 / 38 Objectives By the end of this meeting, participants should be able to: Describe
More information7. Integrated Processes
7. Integrated Processes Up to now: Analysis of stationary processes (stationary ARMA(p, q) processes) Problem: Many economic time series exhibit non-stationary patterns over time 226 Example: We consider
More informationForecasting Bangladesh's Inflation through Econometric Models
American Journal of Economics and Business Administration Original Research Paper Forecasting Bangladesh's Inflation through Econometric Models 1,2 Nazmul Islam 1 Department of Humanities, Bangladesh University
More informationFrequency Forecasting using Time Series ARIMA model
Frequency Forecasting using Time Series ARIMA model Manish Kumar Tikariha DGM(O) NSPCL Bhilai Abstract In view of stringent regulatory stance and recent tariff guidelines, Deviation Settlement mechanism
More information