Forecasting: principles and practice. Rob J Hyndman 3.2 Forecasting with multiple seasonality

Size: px
Start display at page:

Download "Forecasting: principles and practice. Rob J Hyndman 3.2 Forecasting with multiple seasonality"

Transcription

1 Forecasting: principles and practice Rob J Hyndman 3.2 Forecasting with multiple seasonality 1

2 Outline 1 Time series with complex seasonality 2 STL with multiple seasonal periods 3 Dynamic harmonic regression 4 TBATS model 5 Lab session 21 2

3 Examples autoplot(gasoline) + xlab("year") + ylab("thousands of barrels per day") + ggtitle("weekly US finished motor gasoline products") Weekly US finished motor gasoline products Thousands of barrels per day

4 Examples 5 minute call volume at North American bank 400 Call volume Weeks 400 Call volume Weeks 4

5 Examples Turkish daily electricity demand 25 Electricity Demand (GW) Year 5

6 Outline 1 Time series with complex seasonality 2 STL with multiple seasonal periods 3 Dynamic harmonic regression 4 TBATS model 5 Lab session 21 6

7 STL with multiple seasonal periods calls %>% mstl() %>% autoplot() + xlab("week") Week Data Trend Seasonal169 Seasonal845 Remainder 7

8 . STL with multiple seasonal periods calls %>% stlf() %>% autoplot() + xlab("week") Forecasts from STL + ETS(M,N,N) level Week 8

9 . STL with multiple seasonal periods calls %>% stlf() %>% autoplot(include=5*169) + xlab("week") Forecasts from STL + ETS(M,N,N) level Week 9

10 Outline 1 Time series with complex seasonality 2 STL with multiple seasonal periods 3 Dynamic harmonic regression 4 TBATS model 5 Lab session 21 10

11 Dynamic harmonic regression fit <- auto.arima(calls, seasonal=false, lambda=0, xreg=fourier(calls, K=c(10,10))) fit %>% forecast(xreg=fourier(calls, K=c(10,10), h=2*169)) %>% autoplot(include=5*169) + ylab("call volume") + xlab("weeks") 400 Forecasts from Regression with ARIMA(3,1,4) errors Call volume level Weeks 11

12 Outline 1 Time series with complex seasonality 2 STL with multiple seasonal periods 3 Dynamic harmonic regression 4 TBATS model 5 Lab session 21 12

13 TBATS model TBATS Trigonometric terms for seasonality Box-Cox transformations for heterogeneity ARMA errors for short-term dynamics Trend (possibly damped) Seasonal (including multiple and non-integer periods) 13

14 TBATS model y observation at time t (y ω t 1)/ω if ω 0; log y t if ω = 0. l t 1 + φb t 1 + M i=1 l l t 1 + φb t 1 + αd t t m i + d t b (1 φ)b + φb t 1 + βd t d p i=1 k i φ i d t i + j,t q θ j ε t j + ε t j, s(i) j,t 1 cos λ(i) j + s (i) j,t 1 sin λ(i) j + γ 1 (i) d t j, s(i) j,t 1 sin λ(i) j + s (i) j,t 1 cos λ(i) j + γ (i) 2 d t 14

15 TBATS model y observation at time t (y ω t 1)/ω if ω 0; log y t if ω = 0. l t 1 + φb t 1 + M i=1 l l t 1 + φb t 1 + αd t t m i + d t Box-Cox transformation b (1 φ)b + φb t 1 + βd t d p i=1 k i φ i d t i + j,t q θ j ε t j + ε t j, s(i) j,t 1 cos λ(i) j + s (i) j,t 1 sin λ(i) j + γ 1 (i) d t j, s(i) j,t 1 sin λ(i) j + s (i) j,t 1 cos λ(i) j + γ (i) 2 d t 14

16 TBATS model y observation at time t (y ω t 1)/ω if ω 0; log y t if ω = 0. l t 1 + φb t 1 + M i=1 l l t 1 + φb t 1 + αd t t m i + d t Box-Cox transformation M seasonal periods b (1 φ)b + φb t 1 + βd t d p i=1 k i φ i d t i + j,t q θ j ε t j + ε t j, s(i) j,t 1 cos λ(i) j + s (i) j,t 1 sin λ(i) j + γ 1 (i) d t j, s(i) j,t 1 sin λ(i) j + s (i) j,t 1 cos λ(i) j + γ (i) 2 d t 14

17 TBATS model y observation at time t (y ω t 1)/ω if ω 0; log y t if ω = 0. l t 1 + φb t 1 + M i=1 l l t 1 + φb t 1 + αd t t m i + d t b (1 φ)b + φb t 1 + βd t d p i=1 k i φ i d t i + j,t q θ j ε t j + ε t j, Box-Cox transformation M seasonal periods global and local trend s(i) j,t 1 cos λ(i) j + s (i) j,t 1 sin λ(i) j + γ 1 (i) d t j, s(i) j,t 1 sin λ(i) j + s (i) j,t 1 cos λ(i) j + γ (i) 2 d t 14

18 TBATS model y observation at time t (y ω t 1)/ω if ω 0; log y t if ω = 0. l t 1 + φb t 1 + M i=1 l l t 1 + φb t 1 + αd t t m i + d t b (1 φ)b + φb t 1 + βd t d p i=1 k i φ i d t i + j,t q θ j ε t j + ε t j, Box-Cox transformation M seasonal periods global and local trend ARMA error s(i) j,t 1 cos λ(i) j + s (i) j,t 1 sin λ(i) j + γ 1 (i) d t j, s(i) j,t 1 sin λ(i) j + s (i) j,t 1 cos λ(i) j + γ (i) 2 d t 14

19 TBATS model y observation at time t (y ω t 1)/ω if ω 0; log y t if ω = 0. l t 1 + φb t 1 + M i=1 l l t 1 + φb t 1 + αd t t m i + d t b (1 φ)b + φb t 1 + βd t d p i=1 k i φ i d t i + j,t q θ j ε t j + ε t j, s(i) j,t 1 Box-Cox transformation M seasonal periods global and local trend ARMA error Fourier-like seasonal terms + s (i) j,t 1 sin λ(i) j + γ 1 (i) d t cos λ(i) j j, s(i) j,t 1 sin λ(i) j + s (i) j,t 1 cos λ(i) j + γ (i) 2 d t 14

20 TBATS model y observation at time t (yt ω 1)/ω if ω 0; log TBATS y t if ω = 0. l t 1 + φb t 1 + M i=1 l l t 1 + φb t 1 + αd t t m i + d t b (1 φ)b + φb t 1 + βd t d p q φ i d t i + i=1seasonal k i Trigonometric Box-Cox ARMA Trend j,t θ j ε t j + ε t j, s(i) j,t 1 Box-Cox transformation M seasonal periods global and local trend ARMA error Fourier-like seasonal terms + s (i) j,t 1 sin λ(i) j + γ 1 (i) d t cos λ(i) j j, s(i) j,t 1 sin λ(i) j + s (i) j,t 1 cos λ(i) j + γ (i) 2 d t 14

21 . Complex seasonality gasoline %>% tbats() %>% forecast() %>% autoplot() Forecasts from TBATS(1, {0,0},, {<52.18,12>}) 10 9 level Time 15

22 . Complex seasonality calls %>% tbats() %>% forecast() %>% autoplot() Forecasts from TBATS(0.555, {0,0},, {<169,6>, <845,4>}) level Time 16

23 . Complex seasonality telec %>% tbats() %>% forecast() %>% autoplot() Forecasts from TBATS(0.005, {4,2},, {<7,3>, <354.37,7>, <365.25,3>}) level Time 17

24 TBATS model TBATS Trigonometric terms for seasonality Box-Cox transformations for heterogeneity ARMA errors for short-term dynamics Trend (possibly damped) Seasonal (including multiple and non-integer periods) Handles non-integer seasonality, multiple seasonal periods. Entirely automated Prediction intervals often too wide Very slow on long series 18

25 Outline 1 Time series with complex seasonality 2 STL with multiple seasonal periods 3 Dynamic harmonic regression 4 TBATS model 5 Lab session 21 19

26 Lab Session 21 20

FORECASTING USING R. Dynamic regression. Rob Hyndman Author, forecast

FORECASTING 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 information

Forecasting: Principles and Practice. Rob J Hyndman. 12. Advanced methods OTexts.com/fpp/9/2/ OTexts.com/fpp/9/3/

Forecasting: Principles and Practice. Rob J Hyndman. 12. Advanced methods OTexts.com/fpp/9/2/ OTexts.com/fpp/9/3/ Rob J Hyndman Forecasting: Principles and Practice 12. Advanced methods OTexts.com/fpp/9/2/ OTexts.com/fpp/9/3/ Forecasting: Principles and Practice 1 Outline 1 Vector autoregressions 2 Neural network

More information

Transformations for variance stabilization

Transformations 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 information

Forecasting 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 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 information

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 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 information

Forecasting using R. Rob J Hyndman. 1.3 Seasonality and trends. Forecasting using R 1

Forecasting using R. Rob J Hyndman. 1.3 Seasonality and trends. Forecasting using R 1 Forecasting using R Rob J Hyndman 1.3 Seasonality and trends Forecasting using R 1 Outline 1 Time series components 2 STL decomposition 3 Forecasting and decomposition 4 Lab session 5 Forecasting using

More information

Forecasting using R. Rob J Hyndman. 2.5 Seasonal ARIMA models. Forecasting using R 1

Forecasting using R. Rob J Hyndman. 2.5 Seasonal ARIMA models. Forecasting using R 1 Forecasting using R Rob J Hyndman 2.5 Seasonal ARIMA models Forecasting using R 1 Outline 1 Backshift notation reviewed 2 Seasonal ARIMA models 3 ARIMA vs ETS 4 Lab session 12 Forecasting using R Backshift

More information

Forecasting electricity load demand using hybrid exponential smoothing-artificial neural network model

Forecasting electricity load demand using hybrid exponential smoothing-artificial neural network model International Journal of Advances in Intelligent Informatics ISSN: 2442-6571 131 Forecasting electricity load demand using hybrid exponential smoothing-artificial neural network model Winita Sulandari

More information

Automatic forecasting with a modified exponential smoothing state space framework

Automatic forecasting with a modified exponential smoothing state space framework ISSN 1440-771X Department of Econometrics and Business Statistics http://www.buseco.monash.edu.au/depts/ebs/pubs/wpapers/ Automatic forecasting with a modified exponential smoothing state space framework

More information

Short-term electricity demand forecasting in the time domain and in the frequency domain

Short-term electricity demand forecasting in the time domain and in the frequency domain Short-term electricity demand forecasting in the time domain and in the frequency domain Abstract This paper compares the forecast accuracy of different models that explicitely accomodate seasonalities

More information

Forecasting: principles and practice 1

Forecasting: 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 information

Time Series in R: Forecasting and Visualisation. Forecast evaluation 29 May 2017

Time Series in R: Forecasting and Visualisation. Forecast evaluation 29 May 2017 Time Series in R: Forecasting and Visualisation Forecast evaluation 29 May 2017 1 Outline 1 Forecasting residuals 2 Evaluating forecast accuracy 3 Forecasting benchmark methods 4 Lab session 7 5 Time series

More information

Exponentially weighted forecasts

Exponentially 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 information

LOAD FORECASTING APPLICATIONS for THE ENERGY SECTOR

LOAD FORECASTING APPLICATIONS for THE ENERGY SECTOR LOAD FORECASTING APPLICATIONS for THE ENERGY SECTOR Boris Bizjak, Univerza v Mariboru, FERI 26.2.2016 1 A) Short-term load forecasting at industrial plant Ravne: Load forecasting using linear regression,

More information

Time Series Outlier Detection

Time 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 information

Forecasting 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 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 information

Exponential smoothing and non-negative data

Exponential smoothing and non-negative data Exponential smoothing and non-negative data 1 Exponential smoothing and non-negative data Muhammad Akram Rob J Hyndman J Keith Ord Business & Economic Forecasting Unit Exponential smoothing and non-negative

More information

Modelling and Forecasting Australian Domestic Tourism

Modelling and Forecasting Australian Domestic Tourism Modelling and Forecasting Australian Domestic Tourism George Athanasopoulos & Rob Hyndman Outline Background 1 Background 2 Data 3 Regression models 4 Exponential smoothing via innovations state space

More information

Dennis Bricker Dept of Mechanical & Industrial Engineering The University of Iowa. Forecasting demand 02/06/03 page 1 of 34

Dennis Bricker Dept of Mechanical & Industrial Engineering The University of Iowa. Forecasting demand 02/06/03 page 1 of 34 demand -5-4 -3-2 -1 0 1 2 3 Dennis Bricker Dept of Mechanical & Industrial Engineering The University of Iowa Forecasting demand 02/06/03 page 1 of 34 Forecasting is very difficult. especially about the

More information

Intermittent state-space model for demand forecasting

Intermittent state-space model for demand forecasting Intermittent state-space model for demand forecasting 19th IIF Workshop 29th June 2016 Lancaster Centre for Forecasting Motivation Croston (1972) proposes a method for intermittent demand forecasting,

More information

Forecasting: The First Step in Demand Planning

Forecasting: The First Step in Demand Planning Forecasting: The First Step in Demand Planning Jayant Rajgopal, Ph.D., P.E. University of Pittsburgh Pittsburgh, PA 15261 In a supply chain context, forecasting is the estimation of future demand General

More information

STA 6857 VAR, VARMA, VARMAX ( 5.7)

STA 6857 VAR, VARMA, VARMAX ( 5.7) STA 6857 VAR, VARMA, VARMAX ( 5.7) Outline 1 Multivariate Time Series Modeling 2 VAR 3 VARIMA/VARMAX Arthur Berg STA 6857 VAR, VARMA, VARMAX ( 5.7) 2/ 16 Outline 1 Multivariate Time Series Modeling 2 VAR

More information

Functional time series

Functional time series Rob J Hyndman Functional time series with applications in demography 4. Connections, extensions and applications Outline 1 Yield curves 2 Electricity prices 3 Dynamic updating with partially observed functions

More information

Complex exponential Smoothing. Ivan Svetunkov Nikolaos Kourentzes Robert Fildes

Complex exponential Smoothing. Ivan Svetunkov Nikolaos Kourentzes Robert Fildes Complex exponential Smoothing Ivan Svetunkov Nikolaos Kourentzes Robert Fildes 18 March 2015 Introduction Exponential Smoothing methods performed very well in many competitions: M-Competitions in 1982

More information

SHORT TERM LOAD FORECASTING

SHORT TERM LOAD FORECASTING Indian Institute of Technology Kanpur (IITK) and Indian Energy Exchange (IEX) are delighted to announce Training Program on "Power Procurement Strategy and Power Exchanges" 28-30 July, 2014 SHORT TERM

More information

Lecture 5: Estimation of time series

Lecture 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 information

Antti Salonen PPU Le 2: Forecasting 1

Antti Salonen PPU Le 2: Forecasting 1 - 2017 1 Forecasting Forecasts are critical inputs to business plans, annual plans, and budgets Finance, human resources, marketing, operations, and supply chain managers need forecasts to plan: output

More information

Lecture 7: Interpolation

Lecture 7: Interpolation Lecture 7: Interpolation ECE 401: Signal and Image Analysis University of Illinois 2/9/2017 1 Sampling Review 2 Interpolation and Upsampling 3 Spectrum of Interpolated Signals Outline 1 Sampling Review

More information

PPU411 Antti Salonen. Forecasting. Forecasting PPU Forecasts are critical inputs to business plans, annual plans, and budgets

PPU411 Antti Salonen. Forecasting. Forecasting PPU Forecasts are critical inputs to business plans, annual plans, and budgets - 2017 1 Forecasting Forecasts are critical inputs to business plans, annual plans, and budgets Finance, human resources, marketing, operations, and supply chain managers need forecasts to plan: output

More information

A Comparison of the Forecast Performance of. Double Seasonal ARIMA and Double Seasonal. ARFIMA Models of Electricity Load Demand

A 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 information

An Evaluation of Methods for Very Short-Term Load Forecasting Using Minute-by-Minute British Data

An Evaluation of Methods for Very Short-Term Load Forecasting Using Minute-by-Minute British Data An Evaluation of Methods for Very Short-Term Load Forecasting Using Minute-by-Minute British Data James W. Taylor Saïd Business School University of Oxford International Journal of Forecasting, 2008, Vol.

More information

Ross Bettinger, Analytical Consultant, Seattle, WA

Ross 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 information

Capacity Planning and Headroom Analysis for Taming Database Replication Latency - Experiences with LinkedIn Internet Traffic

Capacity Planning and Headroom Analysis for Taming Database Replication Latency - Experiences with LinkedIn Internet Traffic Capacity Planning and Headroom Analysis for Taming Database Replication Latency - Experiences with LinkedIn Internet Traffic Zhenyun Zhuang, Haricharan Ramachandra, Cuong Tran, Subbu Subramaniam, Chavdar

More information

Daily 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 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 information

INDIAN INSTITUTE OF SCIENCE STOCHASTIC HYDROLOGY. Lecture -12 Course Instructor : Prof. P. P. MUJUMDAR Department of Civil Engg., IISc.

INDIAN INSTITUTE OF SCIENCE STOCHASTIC HYDROLOGY. Lecture -12 Course Instructor : Prof. P. P. MUJUMDAR Department of Civil Engg., IISc. INDIAN INSTITUTE OF SCIENCE STOCHASTIC HYDROLOGY Lecture -12 Course Instructor : Prof. P. P. MUJUMDAR Department of Civil Engg., IISc. Summary of the previous lecture Data Extension & Forecasting Moving

More information

Regression of Time Series

Regression of Time Series Mahlerʼs Guide to Regression of Time Series CAS Exam S prepared by Howard C. Mahler, FCAS Copyright 2016 by Howard C. Mahler. Study Aid 2016F-S-9Supplement Howard Mahler hmahler@mac.com www.howardmahler.com/teaching

More information

Time Series Analysis -- An Introduction -- AMS 586

Time 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 information

2.3 Oscillation. The harmonic oscillator equation is the differential equation. d 2 y dt 2 r y (r > 0). Its solutions have the form

2.3 Oscillation. The harmonic oscillator equation is the differential equation. d 2 y dt 2 r y (r > 0). Its solutions have the form 2. Oscillation So far, we have used differential equations to describe functions that grow or decay over time. The next most common behavior for a function is to oscillate, meaning that it increases and

More information

Seasonal Models and Seasonal Adjustment

Seasonal Models and Seasonal Adjustment LECTURE 10 Seasonal Models and Seasonal Adjustment So far, we have relied upon the method of trigonometrical regression for building models which can be used for forecasting seasonal economic time series.

More information

1. Fundamental concepts

1. Fundamental concepts . Fundamental concepts A time series is a sequence of data points, measured typically at successive times spaced at uniform intervals. Time series are used in such fields as statistics, signal processing

More information

Math 54. Selected Solutions for Week 10

Math 54. Selected Solutions for Week 10 Math 54. Selected Solutions for Week 10 Section 4.1 (Page 399) 9. Find a synchronous solution of the form A cos Ωt+B sin Ωt to the given forced oscillator equation using the method of Example 4 to solve

More information

Complex Exponential Smoothing

Complex Exponential Smoothing Complex Exponential Smoothing Ivan Svetunkov Management School Lancaster University A thesis submitted for the degree of Doctor of Philosophy 2016 Acknowledgements It was one of those evenings during our

More information

FORECASTING SUGARCANE PRODUCTION IN INDIA WITH ARIMA MODEL

FORECASTING 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 information

Time series analysis of activity and temperature data of four healthy individuals

Time series analysis of activity and temperature data of four healthy individuals Time series analysis of activity and temperature data of four healthy individuals B.Hadj-Amar N.Cunningham S.Ip March 11, 2016 B.Hadj-Amar, N.Cunningham, S.Ip Time Series Analysis March 11, 2016 1 / 26

More information

Short-term management of hydro-power systems based on uncertainty model in electricity markets

Short-term management of hydro-power systems based on uncertainty model in electricity markets Open Access Journal Journal of Power Technologies 95 (4) (2015) 265 272 journal homepage:papers.itc.pw.edu.pl Short-term management of hydro-power systems based on uncertainty model in electricity markets

More information

3) sin 265 cos 25 - cos 265 sin 25 C) Find the exact value by using a sum or difference identity. 4) sin 165 C) - 627

3) sin 265 cos 25 - cos 265 sin 25 C) Find the exact value by using a sum or difference identity. 4) sin 165 C) - 627 Bonus Assignment Name MULTIPLE CHOICE. Choose the one alternative that best completes the statement or answers the question. Use the given information to find the exact value of the expression. 1) sin

More information

Time Series. Anthony Davison. c

Time Series. Anthony Davison. c Series Anthony Davison c 2008 http://stat.epfl.ch Periodogram 76 Motivation............................................................ 77 Lutenizing hormone data..................................................

More information

Lecture 4. Section 2.5 The Pinching Theorem Section 2.6 Two Basic Properties of Continuity. Jiwen He. Department of Mathematics, University of Houston

Lecture 4. Section 2.5 The Pinching Theorem Section 2.6 Two Basic Properties of Continuity. Jiwen He. Department of Mathematics, University of Houston Review Pinching Theorem Two Basic Properties Lecture 4 Section 2.5 The Pinching Theorem Section 2.6 Two Basic Properties of Continuity Jiwen He Department of Mathematics, University of Houston jiwenhe@math.uh.edu

More information

Department of Econometrics and Business Statistics

Department of Econometrics and Business Statistics ISSN 1440-771X Australia Department of Econometrics and Business Statistics http://www.buseco.monash.edu.au/depts/ebs/pubs/wpapers/ Short-term load forecasting based on a Semi-parametric additive model

More information

EXAMINATIONS OF THE ROYAL STATISTICAL SOCIETY

EXAMINATIONS 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 information

Antti Salonen KPP Le 3: Forecasting KPP227

Antti Salonen KPP Le 3: Forecasting KPP227 - 2015 1 Forecasting Forecasts are critical inputs to business plans, annual plans, and budgets Finance, human resources, marketing, operations, and supply chain managers need forecasts to plan: output

More information

Name (print, please) ID

Name (print, please) ID Name (print, please) ID Operations Management I 7- Winter 00 Odette School of Business University of Windsor Midterm Exam I Solution Wednesday, ebruary, 0:00 :0 pm Last Name A-S: Odette B0 Last Name T-Z:

More information

Chapter 12: An introduction to Time Series Analysis. Chapter 12: An introduction to Time Series Analysis

Chapter 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 information

Forecasting with R A practical workshop

Forecasting with R A practical workshop Forecasting with R A practical workshop International Symposium on Forecasting 2016 19 th June 2016 Nikolaos Kourentzes nikolaos@kourentzes.com http://nikolaos.kourentzes.com Fotios Petropoulos fotpetr@gmail.com

More information

Forecasting. Operations Analysis and Improvement Spring

Forecasting. Operations Analysis and Improvement Spring Forecasting Operations Analysis and Improvement 2015 Spring Dr. Tai-Yue Wang Industrial and Information Management Department National Cheng Kung University 1-2 Outline Introduction to Forecasting Subjective

More information

STR: A Seasonal-Trend Decomposition Procedure Based on Regression

STR: A Seasonal-Trend Decomposition Procedure Based on Regression ISSN 1440-771X Department of Econometrics and Business Statistics http://www.buseco.monash.edu.au/depts/ebs/pubs/wpapers/ STR: A Seasonal-Trend Decomposition Procedure Based on Regression Alexander Dokumentov

More information

Forecasting: principles and practice. Rob J Hyndman 1.1 Introduction to Forecasting

Forecasting: principles and practice. Rob J Hyndman 1.1 Introduction to Forecasting Forecasting: principles and practice Rob J Hyndman 1.1 Introduction to Forecasting 1 Outline 1 Background 2 Case studies 3 The statistical forecasting perspective 4 What can we forecast? 2 Resources Slides

More information

Time-Series analysis for wind speed forecasting

Time-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 information

Evaluation of Some Techniques for Forecasting of Electricity Demand in Sri Lanka

Evaluation of Some Techniques for Forecasting of Electricity Demand in Sri Lanka Appeared in Sri Lankan Journal of Applied Statistics (Volume 3) 00 Evaluation of Some echniques for Forecasting of Electricity Demand in Sri Lanka.M.J. A. Cooray and M.Indralingam Department of Mathematics

More information

TMA4285 December 2015 Time series models, solution.

TMA4285 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

Analytic Trigonometry. Copyright Cengage Learning. All rights reserved.

Analytic Trigonometry. Copyright Cengage Learning. All rights reserved. Analytic Trigonometry Copyright Cengage Learning. All rights reserved. 7.4 Basic Trigonometric Equations Copyright Cengage Learning. All rights reserved. Objectives Basic Trigonometric Equations Solving

More information

Time Series Modeling. Shouvik Mani April 5, /688: Practical Data Science Carnegie Mellon University

Time Series Modeling. Shouvik Mani April 5, /688: Practical Data Science Carnegie Mellon University Time Series Modeling Shouvik Mani April 5, 2018 15-388/688: Practical Data Science Carnegie Mellon University Goals After this lecture, you will be able to: Explain key properties of time series data Describe,

More information

Construction of an Informative Hierarchical Prior Distribution: Application to Electricity Load Forecasting

Construction of an Informative Hierarchical Prior Distribution: Application to Electricity Load Forecasting Construction of an Informative Hierarchical Prior Distribution: Application to Electricity Load Forecasting Anne Philippe Laboratoire de Mathématiques Jean Leray Université de Nantes Workshop EDF-INRIA,

More information

Forecasting: Methods and Applications

Forecasting: Methods and Applications Neapolis University HEPHAESTUS Repository School of Economic Sciences and Business http://hephaestus.nup.ac.cy Books 1998 Forecasting: Methods and Applications Makridakis, Spyros John Wiley & Sons, Inc.

More information

Chapter 2: Complex numbers

Chapter 2: Complex numbers Chapter 2: Complex numbers Complex numbers are commonplace in physics and engineering. In particular, complex numbers enable us to simplify equations and/or more easily find solutions to equations. We

More information

State-space Model. Eduardo Rossi University of Pavia. November Rossi State-space Model Financial Econometrics / 49

State-space Model. Eduardo Rossi University of Pavia. November Rossi State-space Model Financial Econometrics / 49 State-space Model Eduardo Rossi University of Pavia November 2013 Rossi State-space Model Financial Econometrics - 2013 1 / 49 Outline 1 Introduction 2 The Kalman filter 3 Forecast errors 4 State smoothing

More information

8.2 Harmonic Regression and the Periodogram

8.2 Harmonic Regression and the Periodogram Chapter 8 Spectral Methods 8.1 Introduction Spectral methods are based on thining of a time series as a superposition of sinusoidal fluctuations of various frequencies the analogue for a random process

More information

MODELING INFLATION RATES IN NIGERIA: BOX-JENKINS APPROACH. I. U. Moffat and A. E. David Department of Mathematics & Statistics, University of Uyo, Uyo

MODELING 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 information

STA 6857 Signal Extraction & Long Memory ARMA ( 4.11 & 5.2)

STA 6857 Signal Extraction & Long Memory ARMA ( 4.11 & 5.2) STA 6857 Signal Extraction & Long Memory ARMA ( 4.11 & 5.2) Outline 1 Signal Extraction and Optimal Filtering 2 Arthur Berg STA 6857 Signal Extraction & Long Memory ARMA ( 4.11 & 5.2) 2/ 17 Outline 1 Signal

More information

ESSE Mid-Term Test 2017 Tuesday 17 October :30-09:45

ESSE Mid-Term Test 2017 Tuesday 17 October :30-09:45 ESSE 4020 3.0 - Mid-Term Test 207 Tuesday 7 October 207. 08:30-09:45 Symbols have their usual meanings. All questions are worth 0 marks, although some are more difficult than others. Answer as many questions

More information

USING THE MULTIPLICATIVE DOUBLE SEASONAL HOLT-WINTERS METHOD TO FORECAST SHORT-TERM ELECTRICITY DEMAND

USING THE MULTIPLICATIVE DOUBLE SEASONAL HOLT-WINTERS METHOD TO FORECAST SHORT-TERM ELECTRICITY DEMAND Department of Econometrics Erasmus School of Economics Erasmus University Rotterdam Bachelor Thesis BSc Econometrics & Operations Research July 2, 2017 USING THE MULTIPLICATIVE DOUBLE SEASONAL HOLT-WINTERS

More information

Exercises - Time series analysis

Exercises - 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 information

Two-Stage Models for Forecasting Time Series with Multiple Seasonality

Two-Stage Models for Forecasting Time Series with Multiple Seasonality IIMB-WP N0. 540 WORKING PAPER NO: 540 Two-Stage Models for Forecasting Time Series with Multiple Seasonality Anupama Lakshmanan Doctoral Student Indian Institute of Management Bangalore Bannerghatta Road,

More information

Time Series and Forecasting Using R

Time Series and Forecasting Using R Time Series and Forecasting Using R Lab Sessions 25 May 2016 Before doing any exercises in R, load the fpp package using library(fpp) and the ggplot2 package using library(ggplot2). Lab Session: Time series

More information

Multiplicative Winter s Smoothing Method

Multiplicative Winter s Smoothing Method Multiplicative Winter s Smoothing Method LECTURE 6 TIME SERIES FORECASTING METHOD rahmaanisa@apps.ipb.ac.id Review What is the difference between additive and multiplicative seasonal pattern in time series

More information

The admissible parameter space for exponential smoothing models

The admissible parameter space for exponential smoothing models The admissible parameter space for exponential smoothing models Rob J. Hyndman Department of Econometrics and Business Statistics Monash University, VIC 3800, Australia. Telephone: +613 9905 2358 Fax:

More information

A New Approach to Modeling Multivariate Time Series on Multiple Temporal Scales

A New Approach to Modeling Multivariate Time Series on Multiple Temporal Scales University of Nebraska - Lincoln DigitalCommons@University of Nebraska - Lincoln Dissertations and Theses in Statistics Statistics, Department of 5-2015 A New Approach to Modeling Multivariate Time Series

More information

Package thief. R topics documented: January 24, Version 0.3 Title Temporal Hierarchical Forecasting

Package thief. R topics documented: January 24, Version 0.3 Title Temporal Hierarchical Forecasting Version 0.3 Title Temporal Hierarchical Forecasting Package thief January 24, 2018 Methods and tools for generating forecasts at different temporal frequencies using a hierarchical time series approach.

More information

The Analysis to Tertiary-industry with ARIMAX Model

The Analysis to Tertiary-industry with ARIMAX Model Vol. 1, No. 2 ISSN: 1916-9795 The Analysis to Tertiary-industry with ARIMAX Model Jing Fan (Corresponding author) Department of Science, Yanshan University Qin Huangdao 066004, China E-mail: fanjing850328@163.com

More information

Forecasting intraday call arrivals using the seasonal moving average method

Forecasting intraday call arrivals using the seasonal moving average method Forecasting intraday call arrivals using the seasonal moving average method Barrow, DK Author post-print (accepted) deposited by Coventry University s Repository Original citation & hyperlink: Barrow,

More information

DEPARTMENT OF ECONOMETRICS AND BUSINESS STATISTICS

DEPARTMENT OF ECONOMETRICS AND BUSINESS STATISTICS ISSN 1440-771X ISBN 0 7326 1085 0 Unmasking the Theta Method Rob J. Hyndman and Baki Billah Working Paper 5/2001 2001 DEPARTMENT OF ECONOMETRICS AND BUSINESS STATISTICS AUSTRALIA Unmasking the Theta method

More information

Monitoring and data filtering II. Dynamic Linear Models

Monitoring and data filtering II. Dynamic Linear Models Monitoring and data filtering II. Dynamic Linear Models Advanced Herd Management Cécile Cornou, IPH Dias 1 Program for monitoring and data filtering Friday 26 (morning) - Lecture for part I : use of control

More information

FE570 Financial Markets and Trading. Stevens Institute of Technology

FE570 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 information

Serie temporali: analisi e decomposizione

Serie temporali: analisi e decomposizione A. A. 2016-2017 Serie temporali: analisi e decomposizione prof. ing. Antonio Comi Department of Enterprise Engineering University of Rome Tor Vergata Bibliography Forecasting: principles and practice by

More information

21 Linear State-Space Representations

21 Linear State-Space Representations ME 132, Spring 25, UC Berkeley, A Packard 187 21 Linear State-Space Representations First, let s describe the most general type of dynamic system that we will consider/encounter in this class Systems may

More information

A time series is called strictly stationary if the joint distribution of every collection (Y t

A time series is called strictly stationary if the joint distribution of every collection (Y t 5 Time series A time series is a set of observations recorded over time. You can think for example at the GDP of a country over the years (or quarters) or the hourly measurements of temperature over a

More information

Review of Exponential Relations

Review of Exponential Relations Review of Exponential Relations Integrated Math 2 1 Concepts to Know From Video Notes/ HW & Lesson Notes Zero and Integer Exponents Exponent Laws Scientific Notation Analyzing Data Sets (M&M Lab & HW/video

More information

Modelling Wind Farm Data and the Short Term Prediction of Wind Speeds

Modelling Wind Farm Data and the Short Term Prediction of Wind Speeds Modelling Wind Farm Data and the Short Term Prediction of Wind Speeds An Investigation into Wind Speed Data Sets Erin Mitchell Lancaster University 6th April 2011 Outline 1 Data Considerations Overview

More information

CHAPTER 4 FOURIER SERIES S A B A R I N A I S M A I L

CHAPTER 4 FOURIER SERIES S A B A R I N A I S M A I L CHAPTER 4 FOURIER SERIES 1 S A B A R I N A I S M A I L Outline Introduction of the Fourier series. The properties of the Fourier series. Symmetry consideration Application of the Fourier series to circuit

More information

Trend and long-run relations in electricity prices

Trend and long-run relations in electricity prices Trend and long-run relations in electricity prices Why pre-filtering is inevitable Matteo Pelagatti 2 jointly with A. Gianfreda 1 L. Parisio 2 P. Maranzano 2 1 Free University of Bozen 2 University of

More information

Solar irradiance forecasting for Chulalongkorn University location using time series models

Solar irradiance forecasting for Chulalongkorn University location using time series models Senior Project Proposal 2102499 Year 2016 Solar irradiance forecasting for Chulalongkorn University location using time series models Vichaya Layanun ID 5630550721 Advisor: Assist. Prof. Jitkomut Songsiri

More information

Periodic functions: simple harmonic oscillator

Periodic functions: simple harmonic oscillator Periodic functions: simple harmonic oscillator Recall the simple harmonic oscillator (e.g. mass-spring system) d 2 y dt 2 + ω2 0y = 0 Solution can be written in various ways: y(t) = Ae iω 0t y(t) = A cos

More information

New York University Department of Economics. Applied Statistics and Econometrics G Spring 2013

New York University Department of Economics. Applied Statistics and Econometrics G Spring 2013 New York University Department of Economics Applied Statistics and Econometrics G31.1102 Spring 2013 Text: Econometric Analysis, 7 h Edition, by William Greene (Prentice Hall) Optional: A Guide to Modern

More information

Statistical Methods for Forecasting

Statistical Methods for Forecasting Statistical Methods for Forecasting BOVAS ABRAHAM University of Waterloo JOHANNES LEDOLTER University of Iowa John Wiley & Sons New York Chichester Brisbane Toronto Singapore Contents 1 INTRODUCTION AND

More information

Forecasting using robust exponential smoothing with damped trend and seasonal components

Forecasting using robust exponential smoothing with damped trend and seasonal components Forecasting using robust exponential smoothing with damped trend and seasonal components Crevits R, Croux C. KBI_1714 Forecasting using Robust Exponential Smoothing with Damped Trend and Seasonal Components

More information

3 Time Series Regression

3 Time Series Regression 3 Time Series Regression 3.1 Modelling Trend Using Regression Random Walk 2 0 2 4 6 8 Random Walk 0 2 4 6 8 0 10 20 30 40 50 60 (a) Time 0 10 20 30 40 50 60 (b) Time Random Walk 8 6 4 2 0 Random Walk 0

More information

Autobox. Statistical Principles Guide

Autobox. Statistical Principles Guide Autobox Statistical Principles Guide 2 Contents Chapter 1 Overview Chapter 2 Analytical Capabilities Chapter 3 Learning Box-Jenkins Page # 4 10 14 Chapter 4 Statistical Tests and Their Role in Autobox

More information

Dynamic Time Series Regression: A Panacea for Spurious Correlations

Dynamic 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 information

António Ascensão Costa

António Ascensão Costa ISSN 1645-0647 Notes on Pragmatic Procedures and Exponential Smoothing 04-1998 António Ascensão Costa Este documento não pode ser reproduzido, sob qualquer forma, sem autorização expressa. NOTES ON PRAGMATIC

More information

peak half-hourly Tasmania

peak half-hourly Tasmania Forecasting long-term peak half-hourly electricity demand for Tasmania Dr Shu Fan B.S., M.S., Ph.D. Professor Rob J Hyndman B.Sc. (Hons), Ph.D., A.Stat. Business & Economic Forecasting Unit Report for

More information