Classical Decomposition Model Revisited: I
|
|
- Scott Little
- 5 years ago
- Views:
Transcription
1 Classical Decomposition Model Revisited: I recall classical decomposition model for time series Y t, namely, Y t = m t + s t + W t, where m t is trend; s t is periodic with known period s (i.e., s t s = s t for all t 2 Z) satisfying P s j=1 s j = 0; and W t is a stationary process with zero mean m t & s t often taken to be deterministic as we have seen, SARIMA processes can model stochastic s t, but can also handle deterministic m t & s t through di erencing di erencing to eliminate deterministic m t and/or s t can lead to undesirable overdi erencing of W t alternative approach: regard ( ) to be a regression model with stationary errors ( ) BD 23, CC 27, 30, 32, SS 48 XVII 1
2 Classical Decomposition Model Revisited: II will now consider a linear regression approach in which m t and/or s t depend linearly on a small number of parameters note: nonparametric regression another option (expands upon ideas discussed earlier of using filtering operations to extract trends see overhead III 25 and discussion following it) BD 23, CC 27, 30, 32, SS 48 XVII 2
3 Regression with Stationary Errors: I in context of time series analysis, standard linear regression model would take form where Y = X + Z, Y = [Y 1,..., Y n ] 0 is vector containing series; X is an n k design matrix whose tth row x 0 t has values of explanatory variables for Y t ; = [ 1,..., k] 0 is vector of regression coe cients; and Z = [Z 1,..., Z n ] 0 is vector of WN(0, 2 ) RVs example: Y t = t + Z t, for which tth row of n 2 design matrix X would be x 0 t = [1, t] second example: Y t = cos (2 ft) + 3 sin (2 ft) + Z t BD 211, CC 30, SS 50 XVII 3
4 Regression with Stationary Errors: II ordinary least squares (OLS) estimator of is vector ˆOLS minimizing sum of squared errors: nx S( ) = (Y t x 0 t ) 2 = (Y X ) 0 (Y X ) t=1 ˆOLS is solution to so-called normal equations: X 0 X if X 0 X has full rank = X 0 Y and hence ˆOLS = (X 0 X) 1 X 0 Y ˆOLS is best linear unbiased estimator of (best in that, if ˆ is any other unbiased estimator, then var{c 0 ˆOLS } apple var{c 0 ˆ} for any vector c of constants) covariance matrix for ˆOLS given by 2 (X 0 X) 1 BD 211, SS 50 XVII 4
5 Regression with Stationary Errors: III for time series, uncorrelated errors Z are usually unrealistic often a more realistic model is Y = X + W, where W contains RVs from a stationary process with zero mean, an example being a causal ARMA(p, q) process: (B)W t = (B)Z t, {Z t } WN(0, under this alternative model, ˆOLS is still an unbiased estimator of, but best linear unbiased estimator is generalized least squares (GLS) estimator in the full rank case, this estimator takes the form ˆGLS = (X 0 n 1 X) 1 X 0 n 1 Y, where n is covariance matrix for W 2 ) BD 211, 212 XVII 5
6 Regression with Stationary Errors: IV ˆGLS is minimizer of weighted sum of squares: S( ) = (Y X ) 0 1 n (Y X ) to motivate this estimator, recall innovations representation for stationary process says that W = C n U, where U = W Ŵ (U are the innovations, and Ŵ are 1-step-ahead predictions) recall that, if V has covariance matrix n, then covariance matrix for AV is A n A 0 since U has diagonal covariance matrix D n, get n = C n D n Cn 0 let Dn 1/2 be diagonal matrix such that Dn 1/2 Dn 1/2 = Dn 1 covariance matrix for Dn 1/2 Cn 1 W is identity matrix I n since Dn 1/2 Cn 1 n(cn 1 ) 0 Dn 1/2 = Dn 1/2 Cn 1 C n D n Cn(C 0 n 1 ) 0 Dn 1/2 = I n BD 212, 213 XVII 6
7 Regression with Stationary Errors: V returning now to model Y = X + W, multiplication of both sides by Dn 1/2 Cn 1 yields Dn 1/2 Cn 1 Y = Dn 1/2 Cn 1 X + Dn 1/2 Cn 1 W, which can be reexpressed as Ỹ = X + Z, {Z t } WN(0, 1) for this regression model, best linear unbiased estimator is OLS = ( X 0 X) 1 X0 Ỹ = (X 0 (Cn 1 ) 0 Dn 1/2 Dn 1/2 Cn 1 X) 1 X 0 (Cn 1 ) 0 Dn 1/2 Dn 1/2 Cn 1 Y = (X 0 (Cn 1 ) 0 Dn 1 Cn 1 X) 1 X 0 (Cn 1 ) 0 Dn 1 Cn 1 Y = (X 0 n 1 X) 1 X 0 n 1 Y = ˆGLS since n = C n D n Cn 0 says 1 n = (Cn) 0 1 D n Cn 1 = (Cn 1 ) 0 D n Cn 1 BD 212, 213 XVII 7
8 Regression with Stationary Errors: VI in principle, can use ML under a Gaussian assumption to estimate all parameters in model Y = X + W (i.e., both and parameters associated with, e.g., ARMA model) in practice, following simpler (but sub-optimal) iterative scheme for parameter estimation often works well 1. compute ˆOLS and form residuals Y t x 0 t ˆOLS 2. fit ARMA(p,q) or other stationary model to residuals 3. using fitted model, compute ˆGLS and form residuals Y t x 0 t ˆGLS 4. fit same model to residuals again 5. repeat steps 3 and 4 until parameter estimates have stabilized BD 212, 213 XVII 8
9 1st Di erence of Atomic Clock Data Revisited: I let s reconsider modeling 1st di erence of atomic clock data X t, i.e., rx t 2nd di erence r 2 X t well-modeled by ARMA(1,1) process (overheads XIII 124 to 133), with support for need for 2nd di erencing coming from augmented Dickey Fuller (ADF) unit root test (overhead XIV 37) modeling rx t as a stationary fractionally di erenced (FD) process proved to be questionable (overheads XV 75 to 80) apparent linear increase in rx t suggests exploring model of linear trend + stationary noise more seriously (looked briefly at this approach in context of ADF unit root test, where we found that linear detrending obviated need for 2nd di erencing see overhead XIV 37) XVII 9
10 1st Di erence of Atomic Clock Time Series x t t XVII 10
11 Residuals from OLS Fit of Linear Regression x t ACF PACF h h XIV 34
12 1st Di erence of Atomic Clock Data Revisited: II sample ACF and PACF suggest three models for residuals: AR(p), with p somewhere around 7 or 8 ARMA(p,q), with p + q small ARFIMA(p,,q), with p + q small (= 0 gives FD( )) let s consider various possibilities, using maximum likelihood to fit each model AICC to evaluate di erent models start by considering AR(p) models of orders p = 1,..., 20 XVII 11
13 AICC ML-Based AICC for Residuals Modeled by AR(p) p (AR model order) XVII 12
14 1st Di erence of Atomic Clock Data Revisited: III AR(11) best amongst AR(p) with AICC of here are AICCs for ARMA(p,q) models with p + q apple 4 p q AICC XVII 13
15 1st Di erence of Atomic Clock Data Revisited: IV let s compare AICC for ARMA(2,1) with ones for ARMA(p,q) models with p + q = 5 p q AICC let s now see if we can find an ARFIMA(p,,q) model that has a better AICC than one for ARMA(2,1) model XVII 14
16 1st Di erence of Atomic Clock Data Revisited: V p q AICC conclusion: will go with ARMA(2,1) model XVII 15
17 1st Di erence of Atomic Clock Data Revisited: VI using R function arima to fit model rx t = a + bt + W t, {W t } ARMA(2,1), ML estimates of model parameters a, b, 1, 2 and 1 (along with standard errors (SEs)) are â ˆb ˆ1 ˆ1 ˆ 1 value SE estimate of 2 is ˆ2. = ; next plots show rx t & lines fitted by OLS (red) & GLS (black dashed) comparison of ARMA(2,1) ACF & PACF with sample ACF & PACF for rx t â ˆbt residuals ˆR t returned by arima and associated diagnostics XVII 16
18 1st Di erence of Atomic Clock Time Series x t t XVII 17
19 ACF Sample ACF for rx t â ˆbt & ARMA(2,1) ACF h (lag) XVII 18
20 PACF Sample PACF for rx t â ˆbt & ARMA(2,1) PACF h (lag) XVII 19
21 ARMA(2,1) Residuals ˆR t for Atomic Clock R^ t t XVII 20
22 ACF Sample ACF for ARMA(2,1) Residuals ˆR t h (lag) XVII 21
23 PACF Sample PACF for ARMA(2,1) Residuals ˆR t h (lag) XVII 22
24 Q LB and 0.95% quantile Portmanteau Tests of ARMA(2,1) Residuals ˆR t h XVII 23
25 Diagnostics Tests of ARMA(2,1) Residuals ˆR t expected test test value statistic p-value turning point di erence-sign rank runs AR method AICC order AIC order Yule Walker 0 0 Burg 0 0 OLS 0 4 MLE 0 0 XVII 24
26 Regression with Stationary Errors: VII cautionary note: Y 0 Y 6= Ỹ 0 Ỹ = Y 0 1 n Y, so portion of sum of squares explained by transformed model cannot be related directly to sum of squares for untransformed data assuming W in model Y = X + W has covariance matrix n, covariance matrices for ˆOLS and ˆGLS are, respectively, (X 0 X) 1 X 0 nx(x 0 X) 1 and (X 0 n 1 X) 1 assuming W for rx t clock data is ARMA(2,1) process, can assess standard errors (SEs) for estimated slopes OLS GLS slope ( 3% di erence between estimates) SE ( 3% increase for OLS over GLS) if assume W is WN, SE for OLS slope estimate would be taken as (5 smaller than when correlation is accounted for) XVII 25
27 Harmonic Regression and CO 2 Series: I atomic clock rx t illustrates handling deterministic trend m t in model Y t = m t + W t via a parametric regression approach (m t = a + bt) as a second example, reconsider CO 2 series from Mauna Loa, Hawaii (subject of Problem 11) appropriate model here is full classical decomposition model: Y t = m t + s t + W t where m t is trend; s t is seasonal component with period s = 12 (i.e., s t 12 = s t for all t 2 Z) satisfying P 12 j=1 s j = 0; and W t is a stationary process with zero mean for illustrative purposes, will take m t & s t to be deterministic XVII 26
28 2nd Example: CO 2 Series from Mauna Loa, Hawaii y t III 2 year
29 Preliminary Nonparametric Estimate of Trend y t and m^ t year XVII 27
30 Preliminary Detrended Series u t = y t m^ t year XVII 28
31 Estimated Seasonal Component yearly pattern of s^t Jan Feb Mar Apr May June July Aug Sep Oct Nov Dec month XVII 29
32 Harmonic Regression and CO 2 Series: II factoid: can represent any deterministic s t with period 12 by 6X 6X s t = A j cos (2 f j t) + B j sin (2 f j t) = D j cos (2 f j t + ' j ), j=1 j=1 where f j j/12 is a frequency with associated period 1/f j : treating t momentarily as t 2 R rather than t 2 Z, have cos (2 f j (t + 1 f j )) = cos (2 f j t + 2 ) = cos (2 f j t) with a similar result holding for sin (2 f j t). since P 12 j=1 s j = 0, any eleven s j s determine remaining twelfth six A j s and six B j s seem one too many, but in fact there are only five relevant B j s: B 6 doesn t enter in play because B 6 sin (2 f 6 t) = B 6 sin ( t) = B 6 sin ( t) = 0 for all t XVII 30
33 Harmonic Regression and CO 2 Series: III f 1 is known as the fundamental frequency, whereas f 2,..., f 6 are called first,..., fifth harmonics since f j j/12, have f j = jf 1 for j = 2,..., 6, so harmonics are integer multiples of fundamental frequency if seasonal component s t slowly varying from one month to next, can often get by with using just fundamental frequency and a small number of harmonics: JX s t A j cos (2 f j t) + B j sin (2 f j t), 1 apple J < 6 j=1 following plots show approximations of orders J = 1,... 5 to seasonal component ŝ t estimated in Problem 11 for CO 2 series XVII 31
34 J = 1 Approximation to Seasonal Component yearly pattern of s^t Jan Feb Mar Apr May June July Aug Sep Oct Nov Dec month XVII 32
35 J = 2 Approximation to Seasonal Component yearly pattern of s^t Jan Feb Mar Apr May June July Aug Sep Oct Nov Dec month XVII 33
36 J = 3 Approximation to Seasonal Component yearly pattern of s^t Jan Feb Mar Apr May June July Aug Sep Oct Nov Dec month XVII 34
37 J = 4 Approximation to Seasonal Component yearly pattern of s^t Jan Feb Mar Apr May June July Aug Sep Oct Nov Dec month XVII 35
38 J = 5 Approximation to Seasonal Component yearly pattern of s^t Jan Feb Mar Apr May June July Aug Sep Oct Nov Dec month XVII 36
39 J = 6 Perfect Fit to Seasonal Component yearly pattern of s^t Jan Feb Mar Apr May June July Aug Sep Oct Nov Dec month XVII 37
40 Harmonic Regression and CO 2 Series: IV J = 2 approximation looks reasonable, so will entertain model Y t = m t + s t + W t = a + bt + ct 2 + 2X A j cos (2 f j t) + B j sin (2 f j t) + W t j=1 (A j s & B j s are coe cients for so-called harmonic regressors) start analysis by fitting model using OLS to get estimates ˆm t = â + ˆbt 2X + ĉt 2 and ŝ t = Â j cos (2 f j t) + ˆB j sin (2 f j t) following plots look at j=1 ˆm t + ŝ t as compared to Y t cw t = Y t ˆm t ŝ t (surrogates for W t ) and also r c W t XVII 38
41 Parametric Estimate ˆm t + ŝ t and CO 2 Series y t and m^ t + s^t year XVII 39
42 Residuals c W t from Parametric Fit of m t + s t W^ t ACF PACF h h XVII 40
43 First Di erence of Residuals c W t W^ t ACF PACF h h XVII 41
44 Harmonic Regression and CO 2 Series: V sample ACFs & PACFs for f W t & r f W t do not point to obvious simple model more work needed to find reasonable model starting with f W t, consideration of AR(p), with p = 1,..., 31 (ar.mle bombs for p = 32) ARMA(p,q), with p + q small ARFIMA(p,,q), with p + q small using maximum likelihood to fit each model AICC to evaluate individual models leads to an unsatisfying AR(26) model XVII 42
45 Harmonic Regression and CO 2 Series: VI ADP unit root test on f W t suggests need to di erence, but not strongly so (p-value just a bit above 0.05) need for di erencing also hinted at by close to upper limit of 1/2 estimates, which are sample ACF & PACF for r f W t have large values at lag h = 12, suggesting that deterministic s t might be too simplistic SARIMA model with either d = 1 or D = 1 (or both) worthy of consideration bottom line: finding a suitable model for W t for the atomic clock data was easy, but finding one for the CO 2 series poses more of a problem! XVII 43
Differencing Revisited: I ARIMA(p,d,q) processes predicated on notion of dth order differencing of a time series {X t }: for d = 1 and 2, have X t
Differencing Revisited: I ARIMA(p,d,q) processes predicated on notion of dth order differencing of a time series {X t }: for d = 1 and 2, have X t 2 X t def in general = (1 B)X t = X t X t 1 def = ( X
More informationARMA Models: I VIII 1
ARMA Models: I autoregressive moving-average (ARMA) processes play a key role in time series analysis for any positive integer p & any purely nondeterministic process {X t } with ACVF { X (h)}, there is
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 informationTechnical note on seasonal adjustment for Capital goods imports
Technical note on seasonal adjustment for Capital goods imports July 1, 2013 Contents 1 Capital goods imports 2 1.1 Additive versus multiplicative seasonality..................... 2 2 Steps in the seasonal
More informationShort-run electricity demand forecasts in Maharashtra
Applied Economics, 2002, 34, 1055±1059 Short-run electricity demand forecasts in Maharashtra SAJAL GHO SH* and AN JAN A D AS Indira Gandhi Institute of Development Research, Mumbai, India This paper, has
More informationMultiple Regression Analysis
1 OUTLINE Basic Concept: Multiple Regression MULTICOLLINEARITY AUTOCORRELATION HETEROSCEDASTICITY REASEARCH IN FINANCE 2 BASIC CONCEPTS: Multiple Regression Y i = β 1 + β 2 X 1i + β 3 X 2i + β 4 X 3i +
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 informationPart II. Time Series
Part II Time Series 12 Introduction This Part is mainly a summary of the book of Brockwell and Davis (2002). Additionally the textbook Shumway and Stoffer (2010) can be recommended. 1 Our purpose is to
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 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 informationStochastic Processes: I. consider bowl of worms model for oscilloscope experiment:
Stochastic Processes: I consider bowl of worms model for oscilloscope experiment: SAPAscope 2.0 / 0 1 RESET SAPA2e 22, 23 II 1 stochastic process is: Stochastic Processes: II informally: bowl + drawing
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 informationAkaike criterion: Kullback-Leibler discrepancy
Model choice. Akaike s criterion Akaike criterion: Kullback-Leibler discrepancy Given a family of probability densities {f ( ; ), 2 }, Kullback-Leibler s index of f ( ; ) relativetof ( ; ) is Z ( ) =E
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 informationEstimating trends using filters
Estimating trends using filters... contd. 3. Exponential smoothing of data to estimate the trend m[k] ˆm[k] = v[k]+(1 )ˆm[k 1], k =2,, n ˆm[1] = v[1] The choice of has to be fine tuned according to the
More informationNew Introduction to Multiple Time Series Analysis
Helmut Lütkepohl New Introduction to Multiple Time Series Analysis With 49 Figures and 36 Tables Springer Contents 1 Introduction 1 1.1 Objectives of Analyzing Multiple Time Series 1 1.2 Some Basics 2
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 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 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 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 informationTime Series Models and Inference. James L. Powell Department of Economics University of California, Berkeley
Time Series Models and Inference James L. Powell Department of Economics University of California, Berkeley Overview In contrast to the classical linear regression model, in which the components of the
More informationMODELING MAXIMUM MONTHLY TEMPERATURE IN KATUNAYAKE REGION, SRI LANKA: A SARIMA APPROACH
MODELING MAXIMUM MONTHLY TEMPERATURE IN KATUNAYAKE REGION, SRI LANKA: A SARIMA APPROACH M.C.Alibuhtto 1 &P.A.H.R.Ariyarathna 2 1 Department of Mathematical Sciences, Faculty of Applied Sciences, South
More informationTesting for IID Noise/White Noise: I
Testing for IID Noise/White Noise: I want to be able to test null hypothesis time series {x t } or set of residuals {r t } is IID(0, 2 ) or WN(0, 2 ) there are many such tests, including informal test
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 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 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 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 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 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 informationLevinson Durbin Recursions: I
Levinson Durbin Recursions: I note: B&D and S&S say Durbin Levinson but Levinson Durbin is more commonly used (Levinson, 1947, and Durbin, 1960, are source articles sometimes just Levinson is used) recursions
More informationProf. Dr. Roland Füss Lecture Series in Applied Econometrics Summer Term Introduction to Time Series Analysis
Introduction to Time Series Analysis 1 Contents: I. Basics of Time Series Analysis... 4 I.1 Stationarity... 5 I.2 Autocorrelation Function... 9 I.3 Partial Autocorrelation Function (PACF)... 14 I.4 Transformation
More 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 informationLevinson Durbin Recursions: I
Levinson Durbin Recursions: I note: B&D and S&S say Durbin Levinson but Levinson Durbin is more commonly used (Levinson, 1947, and Durbin, 1960, are source articles sometimes just Levinson is used) recursions
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 informationSystem Identification
System Identification Arun K. Tangirala Department of Chemical Engineering IIT Madras July 26, 2013 Module 6 Lecture 1 Arun K. Tangirala System Identification July 26, 2013 1 Objectives of this Module
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 informationTIME SERIES ANALYSIS. Forecasting and Control. Wiley. Fifth Edition GWILYM M. JENKINS GEORGE E. P. BOX GREGORY C. REINSEL GRETA M.
TIME SERIES ANALYSIS Forecasting and Control Fifth Edition GEORGE E. P. BOX GWILYM M. JENKINS GREGORY C. REINSEL GRETA M. LJUNG Wiley CONTENTS PREFACE TO THE FIFTH EDITION PREFACE TO THE FOURTH EDITION
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 informationGAMINGRE 8/1/ of 7
FYE 09/30/92 JULY 92 0.00 254,550.00 0.00 0 0 0 0 0 0 0 0 0 254,550.00 0.00 0.00 0.00 0.00 254,550.00 AUG 10,616,710.31 5,299.95 845,656.83 84,565.68 61,084.86 23,480.82 339,734.73 135,893.89 67,946.95
More informationITSM-R Reference Manual
ITSM-R Reference Manual George Weigt February 11, 2018 1 Contents 1 Introduction 3 1.1 Time series analysis in a nutshell............................... 3 1.2 White Noise Variance.....................................
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 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 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 informationTrending Models in the Data
April 13, 2009 Spurious regression I Before we proceed to test for unit root and trend-stationary models, we will examine the phenomena of spurious regression. The material in this lecture can be found
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 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 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 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 informationIntroduction to Time Series Analysis. Lecture 11.
Introduction to Time Series Analysis. Lecture 11. Peter Bartlett 1. Review: Time series modelling and forecasting 2. Parameter estimation 3. Maximum likelihood estimator 4. Yule-Walker estimation 5. Yule-Walker
More information1 Augmented Dickey Fuller, ADF, Test
Applied Econometrics 1 Augmented Dickey Fuller, ADF, Test Consider a simple general AR(p) process given by Y t = ¹ + Á 1 Y t 1 + Á 2 Y t 2 + ::::Á p Y t p + ² t (1) If this is the process generating the
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 informationLecture 3: Autoregressive Moving Average (ARMA) Models and their Practical Applications
Lecture 3: Autoregressive Moving Average (ARMA) Models and their Practical Applications Prof. Massimo Guidolin 20192 Financial Econometrics Winter/Spring 2018 Overview Moving average processes Autoregressive
More 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 informationLecture 2: Univariate Time Series
Lecture 2: Univariate Time Series Analysis: Conditional and Unconditional Densities, Stationarity, ARMA Processes Prof. Massimo Guidolin 20192 Financial Econometrics Spring/Winter 2017 Overview Motivation:
More informationAnalysis. Components of a Time Series
Module 8: Time Series Analysis 8.2 Components of a Time Series, Detection of Change Points and Trends, Time Series Models Components of a Time Series There can be several things happening simultaneously
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 informationECON/FIN 250: Forecasting in Finance and Economics: Section 7: Unit Roots & Dickey-Fuller Tests
ECON/FIN 250: Forecasting in Finance and Economics: Section 7: Unit Roots & Dickey-Fuller Tests Patrick Herb Brandeis University Spring 2016 Patrick Herb (Brandeis University) Unit Root Tests ECON/FIN
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 informationParameter estimation: ACVF of AR processes
Parameter estimation: ACVF of AR processes Yule-Walker s for AR processes: a method of moments, i.e. µ = x and choose parameters so that γ(h) = ˆγ(h) (for h small ). 12 novembre 2013 1 / 8 Parameter estimation:
More informationSTOCHASTIC MODELING OF MONTHLY RAINFALL AT KOTA REGION
STOCHASTIC MODELIG OF MOTHLY RAIFALL AT KOTA REGIO S. R. Bhakar, Raj Vir Singh, eeraj Chhajed and Anil Kumar Bansal Department of Soil and Water Engineering, CTAE, Udaipur, Rajasthan, India E-mail: srbhakar@rediffmail.com
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 informationARMA models with time-varying coefficients. Periodic case.
ARMA models with time-varying coefficients. Periodic case. Agnieszka Wy lomańska Hugo Steinhaus Center Wroc law University of Technology ARMA models with time-varying coefficients. Periodic case. 1 Some
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 informationMultivariate Time Series Analysis and Its Applications [Tsay (2005), chapter 8]
1 Multivariate Time Series Analysis and Its Applications [Tsay (2005), chapter 8] Insights: Price movements in one market can spread easily and instantly to another market [economic globalization and internet
More informationEconometrics II Heij et al. Chapter 7.1
Chapter 7.1 p. 1/2 Econometrics II Heij et al. Chapter 7.1 Linear Time Series Models for Stationary data Marius Ooms Tinbergen Institute Amsterdam Chapter 7.1 p. 2/2 Program Introduction Modelling philosophy
More information1. The Multivariate Classical Linear Regression Model
Business School, Brunel University MSc. EC550/5509 Modelling Financial Decisions and Markets/Introduction to Quantitative Methods Prof. Menelaos Karanasos (Room SS69, Tel. 08956584) Lecture Notes 5. The
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 information1.4 Properties of the autocovariance for stationary time-series
1.4 Properties of the autocovariance for stationary time-series In general, for a stationary time-series, (i) The variance is given by (0) = E((X t µ) 2 ) 0. (ii) (h) apple (0) for all h 2 Z. ThisfollowsbyCauchy-Schwarzas
More informationCOMPUTER SESSION: ARMA PROCESSES
UPPSALA UNIVERSITY Department of Mathematics Jesper Rydén Stationary Stochastic Processes 1MS025 Autumn 2010 COMPUTER SESSION: ARMA PROCESSES 1 Introduction In this computer session, we work within the
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 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 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 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 informationStatistics of stochastic processes
Introduction Statistics of stochastic processes Generally statistics is performed on observations y 1,..., y n assumed to be realizations of independent random variables Y 1,..., Y n. 14 settembre 2014
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 informationMAT3379 (Winter 2016)
MAT3379 (Winter 2016) Assignment 4 - SOLUTIONS The following questions will be marked: 1a), 2, 4, 6, 7a Total number of points for Assignment 4: 20 Q1. (Theoretical Question, 2 points). Yule-Walker estimation
More informationChristopher Dougherty London School of Economics and Political Science
Introduction to Econometrics FIFTH EDITION Christopher Dougherty London School of Economics and Political Science OXFORD UNIVERSITY PRESS Contents INTRODU CTION 1 Why study econometrics? 1 Aim of this
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 informationINDIAN 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 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 informationChapter 3. ARIMA Models. 3.1 Autoregressive Moving Average Models
Chapter 3 ARIMA Models Classical regression is often insu cient for explaining all of the interesting dynamics of a time series. For example, the ACF of the residuals of the simple linear regression fit
More informationTechnical note on seasonal adjustment for M0
Technical note on seasonal adjustment for M0 July 1, 2013 Contents 1 M0 2 2 Steps in the seasonal adjustment procedure 3 2.1 Pre-adjustment analysis............................... 3 2.2 Seasonal adjustment.................................
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 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 informationFinancial Time Series Analysis: Part II
Department of Mathematics and Statistics, University of Vaasa, Finland Spring 2017 1 Unit root Deterministic trend Stochastic trend Testing for unit root ADF-test (Augmented Dickey-Fuller test) Testing
More informationSTATISTICS OF CLIMATE EXTREMES// TRENDS IN CLIMATE DATASETS
STATISTICS OF CLIMATE EXTREMES// TRENDS IN CLIMATE DATASETS Richard L Smith Departments of STOR and Biostatistics, University of North Carolina at Chapel Hill and Statistical and Applied Mathematical Sciences
More informationTime Series and Forecasting
Time Series and Forecasting Introduction to Forecasting n What is forecasting? n Primary Function is to Predict the Future using (time series related or other) data we have in hand n Why are we interested?
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 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 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 informationTime Series Econometrics 4 Vijayamohanan Pillai N
Time Series Econometrics 4 Vijayamohanan Pillai N Vijayamohan: CDS MPhil: Time Series 5 1 Autoregressive Moving Average Process: ARMA(p, q) Vijayamohan: CDS MPhil: Time Series 5 2 1 Autoregressive Moving
More informationSTAD57 Time Series Analysis. Lecture 23
STAD57 Time Series Analysis Lecture 23 1 Spectral Representation Spectral representation of stationary {X t } is: 12 i2t Xt e du 12 1/2 1/2 for U( ) a stochastic process with independent increments du(ω)=
More informationPeter Hoff Linear and multilinear models April 3, GLS for multivariate regression 5. 3 Covariance estimation for the GLM 8
Contents 1 Linear model 1 2 GLS for multivariate regression 5 3 Covariance estimation for the GLM 8 4 Testing the GLH 11 A reference for some of this material can be found somewhere. 1 Linear model Recall
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 informationNATCOR Regression Modelling for Time Series
Universität Hamburg Institut für Wirtschaftsinformatik Prof. Dr. D.B. Preßmar Professor Robert Fildes NATCOR Regression Modelling for Time Series The material presented has been developed with the substantial
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 informationParametric Signal Modeling and Linear Prediction Theory 1. Discrete-time Stochastic Processes (cont d)
Parametric Signal Modeling and Linear Prediction Theory 1. Discrete-time Stochastic Processes (cont d) Electrical & Computer Engineering North Carolina State University Acknowledgment: ECE792-41 slides
More informationarxiv: v1 [stat.co] 11 Dec 2012
Simulating the Continuation of a Time Series in R December 12, 2012 arxiv:1212.2393v1 [stat.co] 11 Dec 2012 Halis Sak 1 Department of Industrial and Systems Engineering, Yeditepe University, Kayışdağı,
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 information