Time series analysis of activity and temperature data of four healthy individuals
|
|
- Camron Curtis
- 5 years ago
- Views:
Transcription
1 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, / 26
2 Aims Fit time series models to activity data Identify periodicity in the data Produce a 24-hr ahead forecast for activity and temperature Medical applications B.Hadj-Amar, N.Cunningham, S.Ip Time Series Analysis March 11, / 26
3 Data Healthy individuals wore devices recording physical activity and skin temperature Four individuals recorded for approximately four days Physical activity recorded minutely, skin temperature every ten minutes hourly median taken Missing data B.Hadj-Amar, N.Cunningham, S.Ip Time Series Analysis March 11, / 26
4 Data 60 Rest Activity Temperature B.Hadj-Amar, N.Cunningham, S.Ip Time Series Analysis March 11, / 26
5 Autocorrelation function Activity Temperature ACF Lag Lag B.Hadj-Amar, N.Cunningham, S.Ip Time Series Analysis March 11, / 26
6 ARMA {Y 1, Y 2,, Y T } is a time series with zero mean. ɛ t is some noise with zero mean and variance σ 2. B.Hadj-Amar, N.Cunningham, S.Ip Time Series Analysis March 11, / 26
7 ARMA {Y 1, Y 2,, Y T } is a time series with zero mean. ɛ t is some noise with zero mean and variance σ 2. ARMA ARMA(p, q) model has the form Y t = p q φ i Y t i + ɛ t θ j ɛ t j (1) i=1 j=1 B.Hadj-Amar, N.Cunningham, S.Ip Time Series Analysis March 11, / 26
8 Backshift Backshift Operator B k Y t = Y t k (2) B k ɛ t = ɛ t k (3) B.Hadj-Amar, N.Cunningham, S.Ip Time Series Analysis March 11, / 26
9 Backshift Backshift Operator B k Y t = Y t k (2) B k ɛ t = ɛ t k (3) Y t = p q φ Bi i Y t + ɛ t θ Bj j ɛ t i=1 j=1 B.Hadj-Amar, N.Cunningham, S.Ip Time Series Analysis March 11, / 26
10 ARMA Characteristic Equations Y t = p q φ Bi i Y t + ɛ t θ Bj j ɛ t i=1 j=1 B.Hadj-Amar, N.Cunningham, S.Ip Time Series Analysis March 11, / 26
11 ARMA Characteristic Equations Y t = p φ Bi i Y t + ɛ t i=1 φ ( B q θ Bj j ɛ t j=1 ) ( ) Y t = θ B ɛ t (4) B.Hadj-Amar, N.Cunningham, S.Ip Time Series Analysis March 11, / 26
12 ARMA Characteristic Equations Y t = p φ Bi i Y t + ɛ t i=1 φ ( B q θ Bj j ɛ t j=1 ) ( ) Y t = θ B ɛ t (4) ARMA Characteristic Equations φ(x) = 1 θ(x) = 1 p φ i x i (5) i=1 q θ j x j (6) j=1 B.Hadj-Amar, N.Cunningham, S.Ip Time Series Analysis March 11, / 26
13 SARMA Characteristic Equations Introduce seasonality of lag s. B.Hadj-Amar, N.Cunningham, S.Ip Time Series Analysis March 11, / 26
14 SARMA Characteristic Equations Introduce seasonality of lag s. SARMA Characteristic Equations Φ(x) = 1 Θ(x) = 1 P Φ i x is (7) i=1 Q Θ j x js (8) j=1 B.Hadj-Amar, N.Cunningham, S.Ip Time Series Analysis March 11, / 26
15 SARMA SARMA SARMA(p, q) (P, Q) s model has the form ( ) ( ) ( φ B Φ B Y t = θ B ) Θ ( B ) ɛ t (9) B.Hadj-Amar, N.Cunningham, S.Ip Time Series Analysis March 11, / 26
16 SARMA SARMA SARMA(p, q) (P, Q) s model has the form ( ) ( ) ( φ B Φ B Y t = θ B ) Θ ( B ) ɛ t (9) Assume Normal noise, parameters can be estimated using maximum log likelihood. B.Hadj-Amar, N.Cunningham, S.Ip Time Series Analysis March 11, / 26
17 Information Criterions AIC and BIC can be used to select what p, q, P, Q, s to use. B.Hadj-Amar, N.Cunningham, S.Ip Time Series Analysis March 11, / 26
18 Information Criterions AIC and BIC can be used to select what p, q, P, Q, s to use. AIC and BIC Select p, q, P, Q, s which minimizes one of these AIC = T ln σ 2 + (p + q + P + Q)2 (10) BIC = T ln σ 2 + (p + q + P + Q) ln T (11) B.Hadj-Amar, N.Cunningham, S.Ip Time Series Analysis March 11, / 26
19 Experiment 3 day fit, 1 day forecast Select s = 24 hours. Fit ARMA(p, q) and select the best p, q pair. Fix p, q, fit SARMA(p, q) (P, Q) 24 and select the best P, Q pair. Fit and then assess 24 hour forecast B.Hadj-Amar, N.Cunningham, S.Ip Time Series Analysis March 11, / 26
20 Results Figure: Fit SARMA on activity time series B.Hadj-Amar, N.Cunningham, S.Ip Time Series Analysis March 11, / 26
21 Results Figure: Fit SARMA on temperature time series B.Hadj-Amar, N.Cunningham, S.Ip Time Series Analysis March 11, / 26
22 Results AIC Data p q P Q MSE ( C 2 ) Temp ± 0.09 Temp ± 0.1 Temp ± 0.04 Temp ± 0.1 BIC Data p q P Q MSE ( C 2 ) Temp ± 0.08 Temp ± 0.1 Temp ± 0.03 Temp ± 0.06 Table: Selected SARMA models for temperature with mean squared error B.Hadj-Amar, N.Cunningham, S.Ip Time Series Analysis March 11, / 26
23 Harmonic Regression Let us consider the periodic model: X t = Acos(2πωt + φ) + Z t, where A amplitude, φ phase shift, ω fixed frequency, Z t WN(0, σ 2 ) B.Hadj-Amar, N.Cunningham, S.Ip Time Series Analysis March 11, / 26
24 Harmonic Regression Let us consider the periodic model: X t = Acos(2πωt + φ) + Z t, where A amplitude, φ phase shift, ω fixed frequency, Z t WN(0, σ 2 ) Using, trigonometric identities we re-write: X t = β 1 cos(2πωt) + β 2 sin(2πωt) + Z t where β 1 = Acos(φ) and β 2 = Asin(φ). B.Hadj-Amar, N.Cunningham, S.Ip Time Series Analysis March 11, / 26
25 Harmonic Regression Let us consider the periodic model: X t = Acos(2πωt + φ) + Z t, where A amplitude, φ phase shift, ω fixed frequency, Z t WN(0, σ 2 ) Using, trigonometric identities we re-write: X t = β 1 cos(2πωt) + β 2 sin(2πωt) + Z t where β 1 = Acos(φ) and β 2 = Asin(φ). Linear in β 1 and β 2 Linear regression B.Hadj-Amar, N.Cunningham, S.Ip Time Series Analysis March 11, / 26
26 Harmonic Regression Spectral Representation Theorem states that any (weakly) stationary time series can be approximated: X t = µ + K { βk1 cos(2πω k t) + β k2 sin(2πω k t) } k=1 B.Hadj-Amar, N.Cunningham, S.Ip Time Series Analysis March 11, / 26
27 Harmonic Regression Spectral Representation Theorem states that any (weakly) stationary time series can be approximated: X t = µ + K { βk1 cos(2πω k t) + β k2 sin(2πω k t) } k=1 Issue: Find this collection of frequencies {ω k } K k=1 that drive the data Therefore, we use the periodogram I (ω j ), estimator of the frequency spectrum f (ω j ) E[I (ω j )] = f (ω j ) B.Hadj-Amar, N.Cunningham, S.Ip Time Series Analysis March 11, / 26
28 Periodogram However, the periodogram is not a consistent estimator. Generating periodogram of AR(1), for different values of T : Figure: Showing not consistency of the periodogram: black line is the periodogram, red line is the true spectrum. B.Hadj-Amar, N.Cunningham, S.Ip Time Series Analysis March 11, / 26
29 Smoothed Periodogram The smoothed periodogram is instead a consistent estimator of the frequency spectrum: 1 M ˆf (ω j ) = h k I (ω j + k 2M + 1 T ) k= M B.Hadj-Amar, N.Cunningham, S.Ip Time Series Analysis March 11, / 26
30 Temperature: Finding the frequencies that drive the data Figure: Periodogram and smoothed periodograms (uniform and Daniell weights) for temperature of patient 8. This figure is best viewed in colours. B.Hadj-Amar, N.Cunningham, S.Ip Time Series Analysis March 11, / 26
31 Periodogram and Spectrum: AR(p) approach The spectrum of any (weakly) stationary time series can be approximated by the spectrum of an AR(p) process. (a) AIC and BIC (b) Spectrum AR(2), AR(9) Figure: AR(p) approach to obtain the correct frequency representation. Main frequencies that drive the data are around 1/24, and 1/8 B.Hadj-Amar, N.Cunningham, S.Ip Time Series Analysis March 11, / 26
32 Harmonic Regression: Temperature Figure: Harmonic regression for temperature of patient 8, using 5 different harmonics. Dotted lines represent the fitting for single harmonics; thick red line is the superposition of the dotted harmonics, which is the final fit B.Hadj-Amar, N.Cunningham, S.Ip Time Series Analysis March 11, / 26
33 Forecasting: Temperature Figure: 24 hour forecast for temperature of patient 8, given by averaging the fitted model of 4 days. B.Hadj-Amar, N.Cunningham, S.Ip Time Series Analysis March 11, / 26
34 Harmonic Regression: Rest Activity Figure: Harmonic regression for rest activity of patient 8, using 5 harmonics. On top it is shown the classic harmonic fitting, whereas on the bottom the alternative model, where negative values are taken to be zero. B.Hadj-Amar, N.Cunningham, S.Ip Time Series Analysis March 11, / 26
35 Conclusions and future work Harmonic regression provided a better fit for the data In both cases skin temperature was more accurately modeled Further work Examine larger datasets Consider weekday vs weekend effects Model dependence structure between skin temperature and activity levels B.Hadj-Amar, N.Cunningham, S.Ip Time Series Analysis March 11, / 26
36 References Mark Fiecas - Spectral Analysis of Time Series Data University of Warwick - ST Jenkins, GM and Reinsel, GC - Time series analysis: forecasting and control Holden-Day, San Francisco 1976 Shumway, Robert H and Stoffer, David S - Time series analysis and its applications Springer Science & Business Media 2013 B.Hadj-Amar, N.Cunningham, S.Ip Time Series Analysis March 11, / 26
Time 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 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 informationA 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 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 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 informationMCMC analysis of classical time series algorithms.
MCMC analysis of classical time series algorithms. mbalawata@yahoo.com Lappeenranta University of Technology Lappeenranta, 19.03.2009 Outline Introduction 1 Introduction 2 3 Series generation Box-Jenkins
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 informationChapter 3: Regression Methods for Trends
Chapter 3: Regression Methods for Trends Time series exhibiting trends over time have a mean function that is some simple function (not necessarily constant) of time. The example random walk graph from
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 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 informationForecasting Area, Production and Yield of Cotton in India using ARIMA Model
Forecasting Area, Production and Yield of Cotton in India using ARIMA Model M. K. Debnath 1, Kartic Bera 2 *, P. Mishra 1 1 Department of Agricultural Statistics, Bidhan Chanda Krishi Vishwavidyalaya,
More informationElements of Multivariate Time Series Analysis
Gregory C. Reinsel Elements of Multivariate Time Series Analysis Second Edition With 14 Figures Springer Contents Preface to the Second Edition Preface to the First Edition vii ix 1. Vector Time Series
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 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 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 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 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 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 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 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 informationTime 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 informationof seasonal data demonstrating the usefulness of the devised tests. We conclude in "Conclusion" section with a discussion.
DOI 10.1186/s40064-016-3167-4 RESEARCH Open Access Portmanteau test statistics for seasonal serial correlation in time series models Esam Mahdi * *Correspondence: emahdi@iugaza.edu.ps Department of Mathematics,
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 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 informationReview Session: Econometrics - CLEFIN (20192)
Review Session: Econometrics - CLEFIN (20192) Part II: Univariate time series analysis Daniele Bianchi March 20, 2013 Fundamentals Stationarity A time series is a sequence of random variables x t, t =
More informationTime Series Examples Sheet
Lent Term 2001 Richard Weber Time Series Examples Sheet This is the examples sheet for the M. Phil. course in Time Series. A copy can be found at: http://www.statslab.cam.ac.uk/~rrw1/timeseries/ Throughout,
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 informationStatistics 349(02) Review Questions
Statistics 349(0) Review Questions I. Suppose that for N = 80 observations on the time series { : t T} the following statistics were calculated: _ x = 10.54 C(0) = 4.99 In addition the sample autocorrelation
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 informationForecasting Egyptian GDP Using ARIMA Models
Reports on Economics and Finance, Vol. 5, 2019, no. 1, 35-47 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/ref.2019.81023 Forecasting Egyptian GDP Using ARIMA Models Mohamed Reda Abonazel * and
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 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 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 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 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 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 informationChapter 4: Models for Stationary Time Series
Chapter 4: Models for Stationary Time Series Now we will introduce some useful parametric models for time series that are stationary processes. We begin by defining the General Linear Process. Let {Y t
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 informationTime Series: Theory and Methods
Peter J. Brockwell Richard A. Davis Time Series: Theory and Methods Second Edition With 124 Illustrations Springer Contents Preface to the Second Edition Preface to the First Edition vn ix CHAPTER 1 Stationary
More informationTIME SERIES DATA PREDICTION OF NATURAL GAS CONSUMPTION USING ARIMA MODEL
International Journal of Information Technology & Management Information System (IJITMIS) Volume 7, Issue 3, Sep-Dec-2016, pp. 01 07, Article ID: IJITMIS_07_03_001 Available online at http://www.iaeme.com/ijitmis/issues.asp?jtype=ijitmis&vtype=7&itype=3
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 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 informationINTRODUCTION TO TIME SERIES ANALYSIS. The Simple Moving Average Model
INTRODUCTION TO TIME SERIES ANALYSIS The Simple Moving Average Model The Simple Moving Average Model The simple moving average (MA) model: More formally: where t is mean zero white noise (WN). Three parameters:
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 PREFACE xiii 1 Difference Equations 1.1. First-Order Difference Equations 1 1.2. pth-order Difference Equations 7
More information3. ARMA Modeling. Now: Important class of stationary processes
3. ARMA Modeling Now: Important class of stationary processes Definition 3.1: (ARMA(p, q) process) Let {ɛ t } t Z WN(0, σ 2 ) be a white noise process. The process {X t } t Z is called AutoRegressive-Moving-Average
More informationApplied Time. Series Analysis. Wayne A. Woodward. Henry L. Gray. Alan C. Elliott. Dallas, Texas, USA
Applied Time Series Analysis Wayne A. Woodward Southern Methodist University Dallas, Texas, USA Henry L. Gray Southern Methodist University Dallas, Texas, USA Alan C. Elliott University of Texas Southwestern
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 informationLesson 2: Analysis of time series
Lesson 2: Analysis of time series Time series Main aims of time series analysis choosing right model statistical testing forecast driving and optimalisation Problems in analysis of time series time problems
More informationProblem Set 2 Solution Sketches Time Series Analysis Spring 2010
Problem Set 2 Solution Sketches Time Series Analysis Spring 2010 Forecasting 1. Let X and Y be two random variables such that E(X 2 ) < and E(Y 2 )
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 information4. MA(2) +drift: y t = µ + ɛ t + θ 1 ɛ t 1 + θ 2 ɛ t 2. Mean: where θ(l) = 1 + θ 1 L + θ 2 L 2. Therefore,
61 4. MA(2) +drift: y t = µ + ɛ t + θ 1 ɛ t 1 + θ 2 ɛ t 2 Mean: y t = µ + θ(l)ɛ t, where θ(l) = 1 + θ 1 L + θ 2 L 2. Therefore, E(y t ) = µ + θ(l)e(ɛ t ) = µ 62 Example: MA(q) Model: y t = ɛ t + θ 1 ɛ
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 informationKernel-based portmanteau diagnostic test for ARMA time series models
STATISTICS RESEARCH ARTICLE Kernel-based portmanteau diagnostic test for ARMA time series models Esam Mahdi 1 * Received: 01 October 2016 Accepted: 07 February 2017 First Published: 21 February 2017 *Corresponding
More informationLesson 15: Building ARMA models. Examples
Lesson 15: Building ARMA models. Examples Dipartimento di Ingegneria e Scienze dell Informazione e Matematica Università dell Aquila, umberto.triacca@ec.univaq.it Examples In this lesson, in order to illustrate
More informationEstimating AR/MA models
September 17, 2009 Goals The likelihood estimation of AR/MA models AR(1) MA(1) Inference Model specification for a given dataset Why MLE? Traditional linear statistics is one methodology of estimating
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 informationIntroduction to Signal Processing
to Signal Processing Davide Bacciu Dipartimento di Informatica Università di Pisa bacciu@di.unipi.it Intelligent Systems for Pattern Recognition Signals = Time series Definitions Motivations A sequence
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 informationX random; interested in impact of X on Y. Time series analogue of regression.
Multiple time series Given: two series Y and X. Relationship between series? Possible approaches: X deterministic: regress Y on X via generalized least squares: arima.mle in SPlus or arima in R. We have
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 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 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 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 informationForecasting 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 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 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 informationIDENTIFICATION OF ARMA MODELS
IDENTIFICATION OF ARMA MODELS A stationary stochastic process can be characterised, equivalently, by its autocovariance function or its partial autocovariance function. It can also be characterised by
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 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 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 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 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 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 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 III. Stationary models 1 Purely random process 2 Random walk (non-stationary)
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 informationLecture 1: Fundamental concepts in Time Series Analysis (part 2)
Lecture 1: Fundamental concepts in Time Series Analysis (part 2) Florian Pelgrin University of Lausanne, École des HEC Department of mathematics (IMEA-Nice) Sept. 2011 - Jan. 2012 Florian Pelgrin (HEC)
More informationForecasting. This optimal forecast is referred to as the Minimum Mean Square Error Forecast. This optimal forecast is unbiased because
Forecasting 1. Optimal Forecast Criterion - Minimum Mean Square Error Forecast We have now considered how to determine which ARIMA model we should fit to our data, we have also examined how to estimate
More informationFORECASTING OF COTTON PRODUCTION IN INDIA USING ARIMA MODEL
FORECASTING OF COTTON PRODUCTION IN INDIA USING ARIMA MODEL S.Poyyamozhi 1, Dr. A. Kachi Mohideen 2. 1 Assistant Professor and Head, Department of Statistics, Government Arts College (Autonomous), Kumbakonam
More informationUNIVARIATE TIME SERIES ANALYSIS BRIEFING 1970
UNIVARIATE TIME SERIES ANALYSIS BRIEFING 1970 Joseph George Caldwell, PhD (Statistics) 1432 N Camino Mateo, Tucson, AZ 85745-3311 USA Tel. (001)(520)222-3446, E-mail jcaldwell9@yahoo.com (File converted
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 information8.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 informationEcon 623 Econometrics II Topic 2: Stationary Time Series
1 Introduction Econ 623 Econometrics II Topic 2: Stationary Time Series In the regression model we can model the error term as an autoregression AR(1) process. That is, we can use the past value of the
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 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 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 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 informationTime Series Examples Sheet
Lent Term 2001 Richard Weber Time Series Examples Sheet This is the examples sheet for the M. Phil. course in Time Series. A copy can be found at: http://www.statslab.cam.ac.uk/~rrw1/timeseries/ Throughout,
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 informationScenario 5: Internet Usage Solution. θ j
Scenario : Internet Usage Solution Some more information would be interesting about the study in order to know if we can generalize possible findings. For example: Does each data point consist of the total
More informationChapter 9: Forecasting
Chapter 9: Forecasting One of the critical goals of time series analysis is to forecast (predict) the values of the time series at times in the future. When forecasting, we ideally should evaluate the
More informationEstadística Oficial. Transfer Function Model Identication
Boletín de Estadística e Investigación Operativa Vol 25, No 2, Junio 2009, pp 109-115 Estadística Oficial Transfer Function Model Identication Víctor Gómez Ministerio de Economía y Hacienda B vgomez@sgpgmehes
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 informationMoving Average (MA) representations
Moving Average (MA) representations The moving average representation of order M has the following form v[k] = MX c n e[k n]+e[k] (16) n=1 whose transfer function operator form is MX v[k] =H(q 1 )e[k],
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 information2 Time Series Regression and Exploratory Data Analysis
2 Time Series Regression and Exploratory Data Analysis 2.1 Introduction The linear model and its applications are at least as dominant in the time series context as in classical statistics. Regression
More informationShort-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 informationMultiplicative Sarima Modelling Of Nigerian Monthly Crude Oil Domestic Production
Journal of Applied Mathematics & Bioinformatics, vol.3, no.3, 2013, 103-112 ISSN: 1792-6602 (print), 1792-6939 (online) Scienpress Ltd, 2013 Multiplicative Sarima Modelling Of Nigerian Monthly Crude Oil
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 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 information