Forecasting using R. Rob J Hyndman. 2.5 Seasonal ARIMA models. Forecasting using R 1
|
|
- Neal Greene
- 5 years ago
- Views:
Transcription
1 Forecasting using R Rob J Hyndman 2.5 Seasonal ARIMA models Forecasting using R 1
2 Outline 1 Backshift notation reviewed 2 Seasonal ARIMA models 3 ARIMA vs ETS 4 Lab session 12 Forecasting using R Backshift notation reviewed 2
3 Backshift notation A very useful notational device is the backward shift operator, B, which is used as follows: By t = y t 1. In other words, B, operating on y t, has the effect of shifting the data back one period. Two applications of B to y t shifts the data back two periods: B(By t ) = B 2 y t = y t 2. For monthly data, if we wish to shift attention to the same month last year, then B 12 is used, and the notation is B 12 y t = y t 12. Forecasting using R Backshift notation reviewed 3
4 Backshift notation A very useful notational device is the backward shift operator, B, which is used as follows: By t = y t 1. In other words, B, operating on y t, has the effect of shifting the data back one period. Two applications of B to y t shifts the data back two periods: B(By t ) = B 2 y t = y t 2. For monthly data, if we wish to shift attention to the same month last year, then B 12 is used, and the notation is B 12 y t = y t 12. Forecasting using R Backshift notation reviewed 3
5 Backshift notation A very useful notational device is the backward shift operator, B, which is used as follows: By t = y t 1. In other words, B, operating on y t, has the effect of shifting the data back one period. Two applications of B to y t shifts the data back two periods: B(By t ) = B 2 y t = y t 2. For monthly data, if we wish to shift attention to the same month last year, then B 12 is used, and the notation is B 12 y t = y t 12. Forecasting using R Backshift notation reviewed 3
6 Backshift notation A very useful notational device is the backward shift operator, B, which is used as follows: By t = y t 1. In other words, B, operating on y t, has the effect of shifting the data back one period. Two applications of B to y t shifts the data back two periods: B(By t ) = B 2 y t = y t 2. For monthly data, if we wish to shift attention to the same month last year, then B 12 is used, and the notation is B 12 y t = y t 12. Forecasting using R Backshift notation reviewed 3
7 Backshift notation First difference: 1 B. Double difference: (1 B) 2. dth-order difference: (1 B) d y t. Seasonal difference: 1 B m. Seasonal difference followed by a first difference: (1 B)(1 B m ). Multiply terms together together to see the combined effect: (1 B)(1 B m )y t = (1 B B m + B m+1 )y t = y t y t 1 y t m + y t m 1. Forecasting using R Backshift notation reviewed 4
8 Backshift notation for ARIMA ARMA model: y t = c + φ 1 y t φ p y t p + e t + θ 1 e t θ q e t q = c + φ 1 By t + + φ p B p y t + e t + θ 1 Be t + + θ q B q e t φ(b)y t = c + θ(b)e t where φ(b) = 1 φ 1 B φ p B p and θ(b) = 1 + θ 1 B + + θ q B q. ARIMA(1,1,1) model: (1 φ 1 B) (1 B)y t = c + (1 + θ 1 B)e t AR(1) First MA(1) difference Forecasting using R Backshift notation reviewed 5
9 Backshift notation for ARIMA ARMA model: y t = c + φ 1 y t φ p y t p + e t + θ 1 e t θ q e t q = c + φ 1 By t + + φ p B p y t + e t + θ 1 Be t + + θ q B q e t φ(b)y t = c + θ(b)e t where φ(b) = 1 φ 1 B φ p B p and θ(b) = 1 + θ 1 B + + θ q B q. ARIMA(1,1,1) model: (1 φ 1 B) (1 B)y t = c + (1 + θ 1 B)e t AR(1) First MA(1) difference Forecasting using R Backshift notation reviewed 5
10 Backshift notation for ARIMA ARIMA(p, d, q) model: (1 φ 1 B φ p B p ) (1 B) d y t = c + (1 + θ 1 B + + θ q B q )e t AR(p) d differences MA(q) Forecasting using R Backshift notation reviewed 6
11 Outline 1 Backshift notation reviewed 2 Seasonal ARIMA models 3 ARIMA vs ETS 4 Lab session 12 Forecasting using R Seasonal ARIMA models 7
12 Seasonal ARIMA models ARIMA (p, d, q) }{{} Non-seasonal part of the model (P, D, Q) m }{{} Seasonal part of of the model where m = number of observations per year. Forecasting using R Seasonal ARIMA models 8
13 Seasonal ARIMA models E.g., ARIMA(1, 1, 1)(1, 1, 1) 4 model (without constant) (1 φ 1 B)(1 Φ 1 B 4 )(1 B)(1 B 4 )y t = (1+θ 1 B)(1+Θ 1 B 4 )e t. ) ( ) ( ) Non-seasonal Non-seasonal Non-seasonal AR(1) difference MA(1) ( ) Seasonal AR(1) ( ) Seasonal difference ( ) Seasonal MA(1) Forecasting using R Seasonal ARIMA models 9
14 Seasonal ARIMA models E.g., ARIMA(1, 1, 1)(1, 1, 1) 4 model (without constant) (1 φ 1 B)(1 Φ 1 B 4 )(1 B)(1 B 4 )y t = (1+θ 1 B)(1+Θ 1 B 4 )e t. ) ( ) ( ) Non-seasonal Non-seasonal Non-seasonal AR(1) difference MA(1) ( ) Seasonal AR(1) ( ) Seasonal difference ( ) Seasonal MA(1) Forecasting using R Seasonal ARIMA models 9
15 Seasonal ARIMA models E.g., ARIMA(1, 1, 1)(1, 1, 1) 4 model (without constant) (1 φ 1 B)(1 Φ 1 B 4 )(1 B)(1 B 4 )y t = (1+θ 1 B)(1+Θ 1 B 4 )e t. ) ( ) ( ) Non-seasonal Non-seasonal Non-seasonal AR(1) difference MA(1) ( ) Seasonal AR(1) ( ) Seasonal difference ( ) Seasonal MA(1) Forecasting using R Seasonal ARIMA models 9
16 Seasonal ARIMA models E.g., ARIMA(1, 1, 1)(1, 1, 1) 4 model (without constant) (1 φ 1 B)(1 Φ 1 B 4 )(1 B)(1 B 4 )y t = (1+θ 1 B)(1+Θ 1 B 4 )e t. All the factors can be multiplied out and the general model written as follows: y t = (1 + φ 1 )y t 1 φ 1 y t 2 + (1 + Φ 1 )y t 4 (1 + φ 1 + Φ 1 + φ 1 Φ 1 )y t 5 + (φ 1 + φ 1 Φ 1 )y t 6 Φ 1 y t 8 + (Φ 1 + φ 1 Φ 1 )y t 9 φ 1 Φ 1 y t 10 + e t + θ 1 e t 1 + Θ 1 e t 4 + θ 1 Θ 1 e t 5. Forecasting using R Seasonal ARIMA models 10
17 Common ARIMA models In the US Census Bureau uses the following models most often: ARIMA(0,1,1)(0,1,1) m ARIMA(0,1,2)(0,1,1) m ARIMA(2,1,0)(0,1,1) m ARIMA(0,2,2)(0,1,1) m ARIMA(2,1,2)(0,1,1) m with log transformation with log transformation with log transformation with log transformation with no transformation Forecasting using R Seasonal ARIMA models 11
18 Seasonal ARIMA models The seasonal part of an AR or MA model will be seen in the seasonal lags of the PACF and ACF. ARIMA(0,0,0)(0,0,1) 12 will show: a spike at lag 12 in the ACF but no other significant spikes. The PACF will show exponential decay in the seasonal lags; that is, at lags 12, 24, 36,... ARIMA(0,0,0)(1,0,0) 12 will show: exponential decay in the seasonal lags of the ACF a single significant spike at lag 12 in the PACF. Forecasting using R Seasonal ARIMA models 12
19 European quarterly retail trade autoplot(euretail) + xlab("year") + ylab("retai 100 Retail index Year Forecasting using R Seasonal ARIMA models 13
20 European quarterly retail trade ggtsdisplay(diff(euretail,4)) diff(euretail, 4) 2 x Time ACF PACF Lag Lag Forecasting using R Seasonal ARIMA models 14
21 European quarterly retail trade ggtsdisplay(diff(diff(euretail,4))) x diff(diff(euretail, 4)) Time ACF PACF Lag Lag Forecasting using R Seasonal ARIMA models 15
22 European quarterly retail trade d = 1 and D = 1 seems necessary. Significant spike at lag 1 in ACF suggests non-seasonal MA(1) component. Significant spike at lag 4 in ACF suggests seasonal MA(1) component. Initial candidate model: ARIMA(0,1,1)(0,1,1) 4. We could also have started with ARIMA(1,1,0)(1,1,0) 4. Forecasting using R Seasonal ARIMA models 16
23 European quarterly retail trade fit <- Arima(euretail, order=c(0,1,1), seasonal=c(0,1,1)) ggtsdisplay(residuals(fit)) x residuals(fit) Time ACF 0.0 PACF Forecasting using LagR Seasonal ARIMA models Lag 17
24 European quarterly retail trade ACF and PACF of residuals show significant spikes at lag 2, and maybe lag 3. AICc of ARIMA(0,1,2)(0,1,1) 4 model is AICc of ARIMA(0,1,3)(0,1,1) 4 model is fit <- Arima(euretail, order=c(0,1,3), seasonal=c(0,1,1)) ggtsdisplay(residuals(fit)) Forecasting using R Seasonal ARIMA models 18
25 European quarterly retail trade ACF and PACF of residuals show significant spikes at lag 2, and maybe lag 3. AICc of ARIMA(0,1,2)(0,1,1) 4 model is AICc of ARIMA(0,1,3)(0,1,1) 4 model is fit <- Arima(euretail, order=c(0,1,3), seasonal=c(0,1,1)) ggtsdisplay(residuals(fit)) Forecasting using R Seasonal ARIMA models 18
26 European quarterly retail trade residuals(fit) x Time ACF PACF Lag Lag Forecasting using R Seasonal ARIMA models 19
27 European quarterly retail trade res <- residuals(fit) Box.test(res, lag=16, fitdf=4, type="ljung") ## ## Box-Ljung test ## ## data: res ## X-squared = , df = 12, p-value = Forecasting using R Seasonal ARIMA models 20
28 European quarterly retail trade autoplot(forecast(fit, h=12)) Forecasts from ARIMA(0,1,3)(0,1,1)[4] 100 y 95 level Forecasting using R Time Seasonal ARIMA models 21
29 European quarterly retail trade auto.arima(euretail) ## Series: euretail ## ARIMA(1,1,2)(0,1,1)[4] ## ## Coefficients: ## ar1 ma1 ma2 sma1 ## ## s.e ## ## sigma^2 estimated as : log likelihood= ## AIC=69.37 AICc=70.51 BIC=79.76 Forecasting using R Seasonal ARIMA models 22
30 European quarterly retail trade auto.arima(euretail, stepwise=false, approximation=false) ## Series: euretail ## ARIMA(0,1,3)(0,1,1)[4] ## ## Coefficients: ## ma1 ma2 ma3 sma1 ## ## s.e ## ## sigma^2 estimated as : log likelihood=-28.7 ## AIC=67.4 AICc=68.53 BIC=77.78 Forecasting using R Seasonal ARIMA models 23
31 Cortecosteroid drug sales H02 sales (million scripts) Log H02 sales Year Forecasting using R Seasonal ARIMA models 24
32 Cortecosteroid drug sales 0.4 Seasonally differenced H02 scripts 0.2 x Year ACF 0.2 PACF Lag Lag Forecasting using R Seasonal ARIMA models 25
33 Cortecosteroid drug sales Choose D = 1 and d = 0. Spikes in PACF at lags 12 and 24 suggest seasonal AR(2) term. Spikes in PACF suggests possible non-seasonal AR(3) term. Initial candidate model: ARIMA(3,0,0)(2,1,0) 12. Forecasting using R Seasonal ARIMA models 26
34 Cortecosteroid drug sales Model AICc ARIMA(3,0,0)(2,1,0) ARIMA(3,0,1)(2,1,0) ARIMA(3,0,2)(2,1,0) ARIMA(3,0,1)(1,1,0) ARIMA(3,0,1)(0,1,1) ARIMA(3,0,1)(0,1,2) ARIMA(3,0,1)(1,1,1) Forecasting using R Seasonal ARIMA models 27
35 Cortecosteroid drug sales (fit <- Arima(h02, order=c(3,0,1), seasonal=c(0,1,2), lambda=0)) ## Series: h02 ## ARIMA(3,0,1)(0,1,2)[12] ## Box Cox transformation: lambda= 0 ## ## Coefficients: ## ar1 ar2 ar3 ma1 sma1 sma2 ## ## s.e ## ## sigma^2 estimated as : log likelihood= ## AIC= AICc= BIC= Forecasting using R Seasonal ARIMA models 28
36 Cortecosteroid drug sales ggtsdisplay(residuals(fit)) 0.2 residuals(fit) x Time ACF 0.0 PACF Lag Lag Forecasting using R Seasonal ARIMA models 29
37 Cortecosteroid drug sales Box.test(residuals(fit), lag=36, fitdf=6, type="ljung") ## ## Box-Ljung test ## ## data: residuals(fit) ## X-squared = , df = 30, p-value = Forecasting using R Seasonal ARIMA models 30
38 Cortecosteroid drug sales fit <- auto.arima(h02, lambda=0, d=0, D=1, max.order=9, stepwise=false, approximation=false) ggtsdisplay(residuals(fit)) residuals(fit) x ACF Time PACF Lag Lag Forecasting using R Seasonal ARIMA models 31
39 Cortecosteroid drug sales Box.test(residuals(fit), lag=36, fitdf=8, type="ljung") ## ## Box-Ljung test ## ## data: residuals(fit) ## X-squared = , df = 28, p-value = Forecasting using R Seasonal ARIMA models 32
40 Cortecosteroid drug sales Model RMSE ARIMA(3,0,0)(2,1,0)[12] ARIMA(3,0,1)(2,1,0)[12] ARIMA(3,0,2)(2,1,0)[12] ARIMA(3,0,1)(1,1,0)[12] ARIMA(3,0,1)(0,1,1)[12] ARIMA(3,0,1)(0,1,2)[12] ARIMA(3,0,1)(1,1,1)[12] ARIMA(4,0,3)(0,1,1)[12] ARIMA(3,0,3)(0,1,1)[12] ARIMA(4,0,2)(0,1,1)[12] ARIMA(3,0,2)(0,1,1)[12] ARIMA(2,1,3)(0,1,1)[12] ARIMA(2,1,4)(0,1,1)[12] ARIMA(2,1,5)(0,1,1)[12] Forecasting using R Seasonal ARIMA models 33
41 Cortecosteroid drug sales Models with lowest AICc values tend to give slightly better results than the other models. AICc comparisons must have the same orders of differencing. But RMSE test set comparisons can involve any models. No model passes all the residual tests. Use the best model available, even if it does not pass all tests. In this case, the ARIMA(3,0,1)(0,1,2) 12 has the lowest RMSE value and the best AICc value for models with fewer than 6 parameters. Forecasting using R Seasonal ARIMA models 34
42 Cortecosteroid drug sales fit <- Arima(h02, order=c(3,0,1), seasonal=c(0,1,2), lambda=0) autoplot(forecast(fit)) + ylab("h02 sales (million scripts)") + xlab("year") Forecasts from ARIMA(3,0,1)(0,1,2)[12] H02 sales (million scripts) level Year Forecasting using R Seasonal ARIMA models 35
43 Outline 1 Backshift notation reviewed 2 Seasonal ARIMA models 3 ARIMA vs ETS 4 Lab session 12 Forecasting using R ARIMA vs ETS 36
44 ARIMA vs ETS Myth that ARIMA models are more general than exponential smoothing. Linear exponential smoothing models all special cases of ARIMA models. Non-linear exponential smoothing models have no equivalent ARIMA counterparts. Many ARIMA models have no exponential smoothing counterparts. ETS models all non-stationary. Models with seasonality or non-damped trend (or both) have two unit roots; all other models have one unit root. Forecasting using R ARIMA vs ETS 37
45 Equivalences ETS model ARIMA model Parameters ETS(A,N,N) ARIMA(0,1,1) θ 1 = α 1 ETS(A,A,N) ARIMA(0,2,2) θ 1 = α + β 2 θ 2 = 1 α ETS(A,A,N) ARIMA(1,1,2) φ 1 = φ θ 1 = α + φβ 1 φ θ 2 = (1 α)φ ETS(A,N,A) ARIMA(0,0,m)(0,1,0) m ETS(A,A,A) ARIMA(0,1,m + 1)(0,1,0) m ETS(A,A,A) ARIMA(1,0,m + 1)(0,1,0) m Forecasting using R ARIMA vs ETS 38
46 Outline 1 Backshift notation reviewed 2 Seasonal ARIMA models 3 ARIMA vs ETS 4 Lab session 12 Forecasting using R Lab session 12 39
47 Lab Session 12 Forecasting using R Lab session 12 40
Forecasting 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 informationForecasting using R. Rob J Hyndman. 3.2 Dynamic regression. Forecasting using R 1
Forecasting using R Rob J Hyndman 3.2 Dynamic regression Forecasting using R 1 Outline 1 Regression with ARIMA errors 2 Stochastic and deterministic trends 3 Periodic seasonality 4 Lab session 14 5 Dynamic
More informationTransformations for variance stabilization
FORECASTING USING R Transformations for variance stabilization Rob Hyndman Author, forecast Variance stabilization If the data show increasing variation as the level of the series increases, then a transformation
More informationForecasting using R. Rob J Hyndman. 2.3 Stationarity and differencing. Forecasting using R 1
Forecasting using R Rob J Hyndman 2.3 Stationarity and differencing Forecasting using R 1 Outline 1 Stationarity 2 Differencing 3 Unit root tests 4 Lab session 10 5 Backshift notation Forecasting using
More informationForecasting: principles and practice 1
Forecasting: principles and practice Rob J Hyndman 2.3 Stationarity and differencing Forecasting: principles and practice 1 Outline 1 Stationarity 2 Differencing 3 Unit root tests 4 Lab session 10 5 Backshift
More informationExponentially weighted forecasts
FORECASTING USING R Exponentially weighted forecasts Rob Hyndman Author, forecast Simple exponential smoothing Forecasting Notation: ŷ t+h t = y t+h y 1,...,y t point forecast of Forecast Equation: given
More information477/577 In-class Exercise 5 : Fitting Wine Sales
477/577 In-class Exercise 5 : Fitting Wine Sales (due Fri 4/06/2017) Name: Use this file as a template for your report. Submit your code and comments together with (selected) output from R console. Your
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 informationTime Series I Time Domain Methods
Astrostatistics Summer School Penn State University University Park, PA 16802 May 21, 2007 Overview Filtering and the Likelihood Function Time series is the study of data consisting of a sequence of DEPENDENT
More informationAR(p) + I(d) + MA(q) = ARIMA(p, d, q)
AR(p) + I(d) + MA(q) = ARIMA(p, d, q) Outline 1 4.1: Nonstationarity in the Mean 2 ARIMA Arthur Berg AR(p) + I(d)+ MA(q) = ARIMA(p, d, q) 2/ 19 Deterministic Trend Models Polynomial Trend Consider the
More informationFORECASTING USING R. Dynamic regression. Rob Hyndman Author, forecast
FORECASTING USING R Dynamic regression Rob Hyndman Author, forecast Dynamic regression Regression model with ARIMA errors y t = β 0 + β 1 x 1,t + + β r x r,t + e t y t modeled as function of r explanatory
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 informationSTAT 436 / Lecture 16: Key
STAT 436 / 536 - Lecture 16: Key Modeling Non-Stationary Time Series Many time series models are non-stationary. Recall a time series is stationary if the mean and variance are constant in time and the
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 informationSTA 6857 ARIMA and SARIMA Models ( 3.8 and 3.9)
STA 6857 ARIMA and SARIMA Models ( 3.8 and 3.9) Outline 1 Building ARIMA Models 2 SARIMA 3 Homework 4c Arthur Berg STA 6857 ARIMA and SARIMA Models ( 3.8 and 3.9) 2/ 34 Outline 1 Building ARIMA Models
More informationSTA 6857 ARIMA and SARIMA Models ( 3.8 and 3.9) Outline. Return Rate. US Gross National Product
STA 6857 ARIMA and SARIMA Models ( 3.8 and 3.9) Outline 1 Building ARIMA Models 2 SARIMA 3 Homework 4c Arthur Berg STA 6857 ARIMA and SARIMA Models ( 3.8 and 3.9) 2/ 34 Return Rate Suppose x t is the value
More informationWe will only present the general ideas on how to obtain. follow closely the AR(1) and AR(2) cases presented before.
ACF and PACF of an AR(p) We will only present the general ideas on how to obtain the ACF and PACF of an AR(p) model since the details follow closely the AR(1) and AR(2) cases presented before. Recall that
More informationFirstly, the dataset is cleaned and the years and months are separated to provide better distinction (sample below).
Project: Forecasting Sales Step 1: Plan Your Analysis Answer the following questions to help you plan out your analysis: 1. Does the dataset meet the criteria of a time series dataset? Make sure to explore
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 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 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 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 informationMinitab Project Report - Assignment 6
.. Sunspot data Minitab Project Report - Assignment Time Series Plot of y Time Series Plot of X y X 7 9 7 9 The data have a wavy pattern. However, they do not show any seasonality. There seem to be an
More informationSeasonal Autoregressive Integrated Moving Average Model for Precipitation Time Series
Journal of Mathematics and Statistics 8 (4): 500-505, 2012 ISSN 1549-3644 2012 doi:10.3844/jmssp.2012.500.505 Published Online 8 (4) 2012 (http://www.thescipub.com/jmss.toc) Seasonal Autoregressive Integrated
More informationThe log transformation produces a time series whose variance can be treated as constant over time.
TAT 520 Homework 6 Fall 2017 Note: Problem 5 is mandatory for graduate students and extra credit for undergraduates. 1) The quarterly earnings per share for 1960-1980 are in the object in the TA package.
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 informationEstimation and application of best ARIMA model for forecasting the uranium price.
Estimation and application of best ARIMA model for forecasting the uranium price. Medeu Amangeldi May 13, 2018 Capstone Project Superviser: Dongming Wei Second reader: Zhenisbek Assylbekov Abstract This
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 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 informationForecasting: principles and practice. Rob J Hyndman 3.2 Forecasting with multiple seasonality
Forecasting: principles and practice Rob J Hyndman 3.2 Forecasting with multiple seasonality 1 Outline 1 Time series with complex seasonality 2 STL with multiple seasonal periods 3 Dynamic harmonic regression
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 informationTime Series Outlier Detection
Time Series Outlier Detection Tingyi Zhu July 28, 2016 Tingyi Zhu Time Series Outlier Detection July 28, 2016 1 / 42 Outline Time Series Basics Outliers Detection in Single Time Series Outlier Series Detection
More informationHomework 4. 1 Data analysis problems
Homework 4 1 Data analysis problems This week we will be analyzing a number of data sets. We are going to build ARIMA models using the steps outlined in class. It is also a good idea to read section 3.8
More informationSuan Sunandha Rajabhat University
Forecasting Exchange Rate between Thai Baht and the US Dollar Using Time Series Analysis Kunya Bowornchockchai Suan Sunandha Rajabhat University INTRODUCTION The objective of this research is to forecast
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 informationBasics: Definitions and Notation. Stationarity. A More Formal Definition
Basics: Definitions and Notation A Univariate is a sequence of measurements of the same variable collected over (usually regular intervals of) time. Usual assumption in many time series techniques is that
More informationPart 1. Multiple Choice (40 questions, 1 point each) Part 2. Problems/Short Answer (10 questions, 6 points each)
GROUND RULES: This exam contains two parts: Part 1. Multiple Choice (40 questions, 1 point each) Part 2. Problems/Short Answer (10 questions, 6 points each) The maximum number of points on this exam is
More informationSTAT Financial Time Series
STAT 6104 - Financial Time Series Chapter 4 - Estimation in the time Domain Chun Yip Yau (CUHK) STAT 6104:Financial Time Series 1 / 46 Agenda 1 Introduction 2 Moment Estimates 3 Autoregressive Models (AR
More informationEmpirical Approach to Modelling and Forecasting Inflation in Ghana
Current Research Journal of Economic Theory 4(3): 83-87, 2012 ISSN: 2042-485X Maxwell Scientific Organization, 2012 Submitted: April 13, 2012 Accepted: May 06, 2012 Published: June 30, 2012 Empirical Approach
More informationTime Series Analysis -- An Introduction -- AMS 586
Time Series Analysis -- An Introduction -- AMS 586 1 Objectives of time series analysis Data description Data interpretation Modeling Control Prediction & Forecasting 2 Time-Series Data Numerical data
More informationUsing Analysis of Time Series to Forecast numbers of The Patients with Malignant Tumors in Anbar Provinc
Using Analysis of Time Series to Forecast numbers of The Patients with Malignant Tumors in Anbar Provinc /. ) ( ) / (Box & Jenkins).(.(2010-2006) ARIMA(2,1,0). Abstract: The aim of this research is to
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 informationStat 5100 Handout #12.e Notes: ARIMA Models (Unit 7) Key here: after stationary, identify dependence structure (and use for forecasting)
Stat 5100 Handout #12.e Notes: ARIMA Models (Unit 7) Key here: after stationary, identify dependence structure (and use for forecasting) (overshort example) White noise H 0 : Let Z t be the stationary
More informationApplied Forecasting (LECTURENOTES) Prof. Rozenn Dahyot
Applied Forecasting (LECTURENOTES) Prof. Rozenn Dahyot SCHOOL OF COMPUTER SCIENCE AND STATISTICS TRINITY COLLEGE DUBLIN IRELAND https://www.scss.tcd.ie/rozenn.dahyot Michaelmas Term 2017 Contents 1 Introduction
More informationLecture 19 Box-Jenkins Seasonal Models
Lecture 19 Box-Jenkins Seasonal Models If the time series is nonstationary with respect to its variance, then we can stabilize the variance of the time series by using a pre-differencing transformation.
More informationARMA (and ARIMA) models are often expressed in backshift notation.
Backshift Notation ARMA (and ARIMA) models are often expressed in backshift notation. B is the backshift operator (also called the lag operator ). It operates on time series, and means back up by one time
More informationLecture 5: Estimation of time series
Lecture 5, page 1 Lecture 5: Estimation of time series Outline of lesson 5 (chapter 4) (Extended version of the book): a.) Model formulation Explorative analyses Model formulation b.) Model estimation
More informationData %>% Power %>% R You ready?
Data %>% Power %>% R You ready? R a powerful open source programming language available in the database as Oracle R Enterprise that lets you do virtually anything Especially (but not only) Statistics Machine
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 informationAutoregressive Integrated Moving Average Model to Predict Graduate Unemployment in Indonesia
DOI 10.1515/ptse-2017-0005 PTSE 12 (1): 43-50 Autoregressive Integrated Moving Average Model to Predict Graduate Unemployment in Indonesia Umi MAHMUDAH u_mudah@yahoo.com (State Islamic University of Pekalongan,
More informationLesson 13: Box-Jenkins Modeling Strategy for building ARMA models
Lesson 13: Box-Jenkins Modeling Strategy for building ARMA models Facoltà di Economia Università dell Aquila umberto.triacca@gmail.com Introduction In this lesson we present a method to construct an ARMA(p,
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 informationHomework 5, Problem 1 Andrii Baryshpolets 6 April 2017
Homework 5, Problem 1 Andrii Baryshpolets 6 April 2017 Total Private Residential Construction Spending library(quandl) Warning: package 'Quandl' was built under R version 3.3.3 Loading required package:
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 informationSTAT 443 Final Exam Review. 1 Basic Definitions. 2 Statistical Tests. L A TEXer: W. Kong
STAT 443 Final Exam Review L A TEXer: W Kong 1 Basic Definitions Definition 11 The time series {X t } with E[X 2 t ] < is said to be weakly stationary if: 1 µ X (t) = E[X t ] is independent of t 2 γ X
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 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 informationDevelopment of Demand Forecasting Models for Improved Customer Service in Nigeria Soft Drink Industry_ Case of Coca-Cola Company Enugu
International Journal of Scientific Research Engineering & Technology (IJSRET), ISSN 2278 882 Volume 5, Issue 4, April 26 259 Development of Demand Forecasting Models for Improved Customer Service in Nigeria
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 information3 Theory of stationary random processes
3 Theory of stationary random processes 3.1 Linear filters and the General linear process A filter is a transformation of one random sequence {U t } into another, {Y t }. A linear filter is a transformation
More informationTime Series Analysis
Time Series Analysis hm@imm.dtu.dk Informatics and Mathematical Modelling Technical University of Denmark DK-2800 Kgs. Lyngby 1 Outline of the lecture Chapter 9 Multivariate time series 2 Transfer function
More informationStat 153 Time Series. Problem Set 4
Stat Time Series Problem Set Problem I generated 000 observations from the MA() model with parameter 0.7 using the following R function and then fitted the ARMA(, ) model to the data. order = c(, 0, )),
More informationSome Time-Series Models
Some Time-Series Models Outline 1. Stochastic processes and their properties 2. Stationary processes 3. Some properties of the autocorrelation function 4. Some useful models Purely random processes, random
More informationBooth School of Business, University of Chicago Business 41914, Spring Quarter 2017, Mr. Ruey S. Tsay. Solutions to Midterm
Booth School of Business, University of Chicago Business 41914, Spring Quarter 017, Mr Ruey S Tsay Solutions to Midterm Problem A: (51 points; 3 points per question) Answer briefly the following questions
More informationARIMA Models. Jamie Monogan. January 16, University of Georgia. Jamie Monogan (UGA) ARIMA Models January 16, / 27
ARIMA Models Jamie Monogan University of Georgia January 16, 2018 Jamie Monogan (UGA) ARIMA Models January 16, 2018 1 / 27 Objectives By the end of this meeting, participants should be able to: Argue why
More informationUniversity of Oxford. Statistical Methods Autocorrelation. Identification and Estimation
University of Oxford Statistical Methods Autocorrelation Identification and Estimation Dr. Órlaith Burke Michaelmas Term, 2011 Department of Statistics, 1 South Parks Road, Oxford OX1 3TG Contents 1 Model
More informationSeasonal Adjustment using X-13ARIMA-SEATS
Seasonal Adjustment using X-13ARIMA-SEATS Revised: 10/9/2017 Summary... 1 Data Input... 3 Limitations... 4 Analysis Options... 5 Tables and Graphs... 6 Analysis Summary... 7 Data Table... 9 Trend-Cycle
More informationA Comparison of the Forecast Performance of. Double Seasonal ARIMA and Double Seasonal. ARFIMA Models of Electricity Load Demand
Applied Mathematical Sciences, Vol. 6, 0, no. 35, 6705-67 A Comparison of the Forecast Performance of Double Seasonal ARIMA and Double Seasonal ARFIMA Models of Electricity Load Demand Siti Normah Hassan
More informationA Data-Driven Model for Software Reliability Prediction
A Data-Driven Model for Software Reliability Prediction Author: Jung-Hua Lo IEEE International Conference on Granular Computing (2012) Young Taek Kim KAIST SE Lab. 9/4/2013 Contents Introduction Background
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 informationTMA4285 December 2015 Time series models, solution.
Norwegian University of Science and Technology Department of Mathematical Sciences Page of 5 TMA4285 December 205 Time series models, solution. Problem a) (i) The slow decay of the ACF of z t suggest that
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 informationDaily sales forecasting in foodservice: Developing a model for application in an expert system
Daily sales forecasting in foodservice: Developing a model for application in an expert system BPJ 420 Final Report 28 September 2016 H.F. Versluis 11300389 Final Year B.Eng. (Industrial) ii Executive
More informationUnit root problem, solution of difference equations Simple deterministic model, question of unit root
Unit root problem, solution of difference equations Simple deterministic model, question of unit root (1 φ 1 L)X t = µ, Solution X t φ 1 X t 1 = µ X t = A + Bz t with unknown z and unknown A (clearly X
More informationSugarcane Productivity in Bihar- A Forecast through ARIMA Model
Available online at www.ijpab.com Kumar et al Int. J. Pure App. Biosci. 5 (6): 1042-1051 (2017) ISSN: 2320 7051 DOI: http://dx.doi.org/10.18782/2320-7051.5838 ISSN: 2320 7051 Int. J. Pure App. Biosci.
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 informationFinal Examination 7/6/2011
The Islamic University of Gaza Faculty of Commerce Department of Economics & Applied Statistics Time Series Analysis - Dr. Samir Safi Spring Semester 211 Final Examination 7/6/211 Name: ID: INSTRUCTIONS:
More informationStat 565. (S)Arima & Forecasting. Charlotte Wickham. stat565.cwick.co.nz. Feb
Stat 565 (S)Arima & Forecasting Feb 2 2016 Charlotte Wickham stat565.cwick.co.nz Today A note from HW #3 Pick up with ARIMA processes Introduction to forecasting HW #3 The sample autocorrelation coefficients
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 informationElectricity Consumption Forecasting in the Khan Younis Province Using Exponential Smoothing and Box - Jenkins Methods: A Modeling Viewpoint.
Electricity Consumption Forecasting in the Khan Younis Province Using Exponential Smoothing and Box - Jenkins Methods: A Modeling Viewpoint. August 26, 2015 The Islamic University of Gaza Faculty of Science
More informationAnalysis of Violent Crime in Los Angeles County
Analysis of Violent Crime in Los Angeles County Xiaohong Huang UID: 004693375 March 20, 2017 Abstract Violent crime can have a negative impact to the victims and the neighborhoods. It can affect people
More information2. An Introduction to Moving Average Models and ARMA Models
. An Introduction to Moving Average Models and ARMA Models.1 White Noise. The MA(1) model.3 The MA(q) model..4 Estimation and forecasting of MA models..5 ARMA(p,q) models. The Moving Average (MA) models
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 informationForecasting Gold Price. A Comparative Study
Course of Financial Econometrics FIRM Forecasting Gold Price A Comparative Study Alessio Azzutti, University of Florence Abstract: This paper seeks to evaluate the appropriateness of a variety of existing
More informationAutomatic seasonal auto regressive moving average models and unit root test detection
ISSN 1750-9653, England, UK International Journal of Management Science and Engineering Management Vol. 3 (2008) No. 4, pp. 266-274 Automatic seasonal auto regressive moving average models and unit root
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 informationTime Series Analysis of United States of America Crude Oil and Petroleum Products Importations from Saudi Arabia
International Journal of Applied Science and Technology Vol. 5, No. 5; October 2015 Time Series Analysis of United States of America Crude Oil and Petroleum Products Importations from Saudi Arabia Olayan
More informationMore About Exponential Smoothing
More About Exponential Smoothing There is a wealth of methods and literature for the task of time series forecasting. A very good overview and access to methods that are simple to handle can be found here:
More informationMarcel Dettling. Applied Time Series Analysis SS 2013 Week 05. ETH Zürich, March 18, Institute for Data Analysis and Process Design
Marcel Dettling Institute for Data Analysis and Process Design Zurich University of Applied Sciences marcel.dettling@zhaw.ch http://stat.ethz.ch/~dettling ETH Zürich, March 18, 2013 1 Basics of Modeling
More informationModelling using ARMA processes
Modelling using ARMA processes Step 1. ARMA model identification; Step 2. ARMA parameter estimation Step 3. ARMA model selection ; Step 4. ARMA model checking; Step 5. forecasting from ARMA models. 33
More informationEXAMINATIONS OF THE ROYAL STATISTICAL SOCIETY
EXAMINATIONS OF THE ROYAL STATISTICAL SOCIETY GRADUATE DIPLOMA, 011 MODULE 3 : Stochastic processes and time series Time allowed: Three Hours Candidates should answer FIVE questions. All questions carry
More informationDesign of Time Series Model for Road Accident Fatal Death in Tamilnadu
Volume 109 No. 8 2016, 225-232 ISSN: 1311-8080 (printed version); ISSN: 1314-3395 (on-line version) url: http://www.ijpam.eu ijpam.eu Design of Time Series Model for Road Accident Fatal Death in Tamilnadu
More informationApplied Time Series Topics
Applied Time Series Topics Ivan Medovikov Brock University April 16, 2013 Ivan Medovikov, Brock University Applied Time Series Topics 1/34 Overview 1. Non-stationary data and consequences 2. Trends and
More informationTime-Series analysis for wind speed forecasting
Malaya Journal of Matematik, Vol. S, No. 1, 55-61, 2018 https://doi.org/10.26637/mjm0s01/11 Time-Series analysis for wind speed forecasting Garima Jain1 Abstract In this paper, an enormous amount of study
More information03 Time series with trend and seasonality components. Andrius Buteikis,
03 Time series with trend and seasonality components Andrius Buteikis, andrius.buteikis@mif.vu.lt http://web.vu.lt/mif/a.buteikis/ Time series with deterministic components Up until now we assumed our
More informationModule 4. Stationary Time Series Models Part 1 MA Models and Their Properties
Module 4 Stationary Time Series Models Part 1 MA Models and Their Properties Class notes for Statistics 451: Applied Time Series Iowa State University Copyright 2015 W. Q. Meeker. February 14, 2016 20h
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 informationExercises - Time series analysis
Descriptive analysis of a time series (1) Estimate the trend of the series of gasoline consumption in Spain using a straight line in the period from 1945 to 1995 and generate forecasts for 24 months. Compare
More informationSTAT 520 FORECASTING AND TIME SERIES 2013 FALL Homework 05
STAT 520 FORECASTING AND TIME SERIES 2013 FALL Homework 05 1. ibm data: The random walk model of first differences is chosen to be the suggest model of ibm data. That is (1 B)Y t = e t where e t is a mean
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 information