INDIAN INSTITUTE OF SCIENCE STOCHASTIC HYDROLOGY. Lecture -20 Course Instructor : Prof. P. P. MUJUMDAR Department of Civil Engg., IISc.
|
|
- Dominick Hancock
- 5 years ago
- Views:
Transcription
1 INDIAN INSTITUTE OF SCIENCE STOCHASTIC HYDROLOGY Lecture -20 Course Instructor : Prof. P. P. MUJUMDAR Department of Civil Engg., IISc.
2 Summary of the previous lecture Case study -3: Monthly streamflows ate KRS reservoir Plots of Time series, Correlogram, Partial Autocorrelation function and Power spectrum Standardization to remove periodicities Candidate ARMA models : (contiguous and noncontiguous) Log Likelihood Mean square error 2
3 CASE STUDIES - ARMA MODELS 3
4 Case study 3 Monthly Stream flow (cum/sec) Statistics( ) for Cauvery River at Krishna Raja Sagar Reservoir is considered in the case study. Time series of the data, auto correlation function, partial auto correlation function and the power spectrum are plotted. The series indicates presence of periodicities. The series is standardized to remove the periodicities. 4
5 Case study 3 (Contd.) Standardized series is considered for fitting the ARMA models Total length of the data set N = 480 Half the data set (240 values) is used to construct the model and other half is used for validation. Both contiguous and non-contiguous models are studied Non-contiguous models consider the most significant AR and MA terms leaving out the intermediate terms 5
6 Case study 3 (Contd.) Contiguous models: Sl. No Model Likelihood values 1 ARMA(1,0) ARMA(2,0) ARMA(3,0) ARMA(4,0) ARMA(5,0) ARMA(6,0) ARMA(1,1) ARMA(1,2) ARMA(2,1) ARMA(2,2) ARMA(3,1) L N = ln ( σ ) n 2 i i i 6
7 Case study 3 (Contd.) Non-contiguous models*: Sl. No Model Likelihood values 1 ARMA(2,0) ARMA(3,0) ARMA(4,0) ARMA(5,0) ARMA(6,0) ARMA(7,0) ARMA(2,2) ARMA(2,3) ARMA(3,2) ARMA(3,3) ARMA(4,2) *: The last AR and MA terms correspond to the 12 th lag 7
8 Case study 3 (Contd.) For this time series, the likelihood values for contiguous model = non-contiguous model = Hence contiguous ARMA(4,0) can be used. The parameters for the selected model are as follows φ 1 = φ 2 = φ 3 = φ 4 = Constant =
9 Case study 3 (Contd.) Contiguous models: Forecasting Models Sl. No Model Mean square error values 1 ARMA(1,0) ARMA(2,0) ARMA(3,0) ARMA(4,0) ARMA(5,0) ARMA(6,0) ARMA(1,1) ARMA(1,2) ARMA(2,1) ARMA(2,2) ARMA(3,1)
10 Case study 3 (Contd.) Non-contiguous models: Sl. No Model Mean square error values 1 ARMA(2,0) ARMA(3,0) ARMA(4,0) ARMA(5,0) ARMA(6,0) ARMA(7,0) ARMA(2,2) ARMA(2,3) ARMA(3,2) ARMA(3,3) ARMA(4,2)
11 Case study 3 (Contd.) The simplest model AR(1) results in the least value of the MSE For one step forecasting, quite often the simplest model is appropriate Also as the number of parameters increases, the MSE increases which is contrary to the common belief that models with large number of parameters give better forecasts. AR(1) model is recommended for forecasting the series and the parameters are as follows φ 1 = and C =
12 Case study 3 (Contd.) Validation tests on the residual series Significance of residual mean Significance of periodicities Cumulative periodogram test or Bartlett s test White noise test Whittle s test Portmanteau test Residuals, Residual m1 m 2 et = Xt φjxt j + θjet j + C j= 1 j= 1 Data Simulated from the model 12
13 Case study 3 (Contd.) Significance of residual mean: Sl. No Model η(e) t 0.95 (239 ) 1 ARMA(1,0) ARMA(2,0) ARMA(3,0) ARMA(4,0) ARMA(5,0) ARMA(6,0) ARMA(1,1) ARMA(1,2) ARMA(2,1) ARMA(2,2) η ( e) = N ρ 1/2 ˆ1/2 e η(e) < t(0.95, 240 1); All models pass the test 13
14 Case study 3 (Contd.) Significance of periodicities: η e ( ) = γ 2 k N ( ) 4 ˆ ρ 1 2 γ k2 = α k 2 + β k 2 1 N ˆ ρ ˆ 1 = ˆcos α ω βsin ω N t= 1 n 2 α = e { ( ) ( )} 2 et kt kt cos ( ω t) k t k N t= 1 14
15 Case study 3 (Contd.) n 2 β = e sin ( ω t) k t k N t= 1 2π/ω k is the periodicity for which test is being carried out. η(e) < F α (2, N 2 ) Model passes the test 15
16 Case study 3 (Contd.) Significance of periodicities: Sl. No Model η Value for the periodicity 1 st 2 nd 3 rd 4 th F 0.95 (2,238 ) 1 ARMA(1,0) ARMA(2,0) ARMA(3,0) ARMA(4,0) ARMA(5,0) ARMA(6,0) ARMA(1,1) ARMA(1,2) ARMA(2,1) ARMA(2,2)
17 Case study 3 (Contd.) Significance of periodicities by Bartlett s test : (Cumulative periodogram test) N 2 N k = etcos( kt) + etsin ( kt) N t= 1 N t= 1 γ ω ω g k = k j= 1 N /2 k = 1 γ γ 2 j 2 k k = 1,2, N/2 17
18 Case study 3 (Contd.) Cumulative periodogram for the original series without standardizing gk Lag, k
19 Case study 3 (Contd.) The confidence limits (1.35/(N/2) 1/2 = 0.087) are plotted for 95% confidence. (N = 480 in this case) Most part of the cumulative periodogram lies outside the significance bands confirming the presence of periodicity in the data. for k=40, a spike in the graph is seen indicating the significant periodicity This k corresponds to a periodicity of 12 months (480/40) k=80, corresponds to a periodicity of 6 months 19
20 Case study 3 (Contd.) Cumulative periodogram for the residual series of ARMA(4, 0) model gk Lag, k 20
21 Case study 3 (Contd.) The confidence limits (+1.35/(N/2) 1/2 = ) are plotted for 95% confidence. (N = 240, in this case) Cumulative periodogram lies within the significance bands confirming that no significant periodicity present in the residual series. The model passes the test. 21
22 Whittle s test for white noise White noise test (Whittle s test): This test is carried out to test the absence of correlation in the series. The covariance r k at lag k of the error series {e t } r N 1 k = ejej k N k j= k+ 1 k = 0, 1, 2,.k max The value of k max is normally chosen as 0.15N Ref: Kashyap R.L. and Ramachandra Rao.A, Dynamic stochastic models from empirical data, Academic press, New York,
23 Whittle s test for white noise The covariance matrix is r r r.. r r r r.. r r max k 1 2 Γ n1 =.. rk r max k r max max 1 0 k k max x k max 23
24 Whittle s test for white noise A statistic η(e) is defined as η e N ˆ ρ = 1 n1 1 ˆ ρ1 ( ) 0 n1 = k max Where ˆρ 0 is the lag zero correlation =1, and ˆ ρ = 1 det Γ det Γ n1 n11 The matrix Γ n1-1 is constructed by eliminating the last row and the last column from the Γ n1 matrix. 24
25 Whittle s test for white noise The statistic η(e) is approximately distributed as F α (n1, N n1 ), where α is the significance level at which the test is being carried out. If the value of η(e) < F α (n1, N n1 ), then the residual series is uncorrelated. 25
26 Case study 3 (Contd.) N ˆ ρ 0 η ( e) = 1 n1 1 ˆ ρ1 F 0.95 (2,239 ) n1 = 73 n1 = 49 n1 = 25 Model η η η ARMA(1,0) ARMA(2,0) ARMA(3,0) ARMA(4,0) ARMA(5,0) ARMA(6,0) ARMA(1,1) ARMA(1,2) ARMA(2,1) ARMA(2,2) Whittle s test - white noise model fails 26
27 Case study 3 (Contd.) Portmanteau test for white noise: 2 n1 r k = k= 1 r0 ( e) ( N n1) χ 2 kmax = 48 kmax = 36 kmax = 24 kmax = (k max ) Model η η η η ARMA(1,0) ARMA(2,0) ARMA(3,0) ARMA(4,0) ARMA(5,0) ARMA(6,0) ARMA(1,1) ARMA(1,2) ARMA(2,1) ARMA(2,2) model fails η 27
28 Case study 4 Monthly Stream flow ( ) of a river (37 years of monthly data); Flow in cumec Time (months) 28
29 Case study 4 (Contd.) Correlogram Correlation coefficient Lag 29
30 Case study 4 (Contd.) Partial Auto Correlation function 1 Fig 10 Partial Autocorrelation Function of flows PAC Partial Autocorrelations Lag Lag 30
31 Case study 4 (Contd.) Power Spectrum 1.2E+09 1E I(k) E+08 W(k) 31
32 Case study 4 (Contd.) Model name constant φ 1 φ 2 φ 3 φ 4 φ 5 φ 6 φ 7 φ 8 θ 1 θ 2 ARMA(1,0) ARMA(2,0) ARMA(3,0) ARMA(4,0) ARMA(5,0) ARMA(6,0) ARMA(7,0) ARMA(8,0) ARMA(1,1) ARMA(2,1) ARMA(3,1) ARMA(4,1) ARMA(1,2) ARMA(2,2) ARMA(0,1) ARMA(0,2)
33 Case study 4 (Contd.) Sl. No Model Mean Square Error Likelihood value 1 ARMA(1,0) ARMA(2,0) ARMA(3,0) ARMA(4,0) ARMA(5,0) ARMA(6,0) ARMA(7,0) ARMA(8,0) ARMA(1,1) ARMA(2,1) ARMA(3,1) ARMA(4,1) ARMA(1,2) ARMA(2,2) ARMA(0,1) ARMA(0,2)
34 Case study 4 (Contd.) Significance of residual mean: e η ( e) = N ρ 1/2 ˆ1/2 e is the estimate of the residual mean ˆρ is the estimate of the residual variance All models pass the test Model Test t η(e) t(α, N 1) ARMA(1,0) E ARMA(2,0) E ARMA(3,0) E ARMA(4,0) E ARMA(5,0) E ARMA(6,0) E ARMA(7,0) E ARMA(8,0) E ARMA(1,1) ARMA(2,1) ARMA(3,1) ARMA(4,1) ARMA(1,2) ARMA(2,2) ARMA(0,1) ARMA(0,2)
35 Case study 5 Sakleshpur Annual Rainfall Data ( ) Rainfall in mm Time (years) 35
36 Case study 5 (Contd.) Correlogram 36
37 Case study 5 (Contd.) PAC function 37
38 Case study 5 (Contd.) Model Likelihood AR(1) AR(2) AR(3) AR(4) AR(5) AR(6) ARMA(1,1) ARMA(1,2) ARMA(2,1) ARMA(2,2) ARMA(3,1) ARMA(3,2)
39 Case study 5 (Contd.) ARMA(5,0) is selected with highest likelihood value The parameters for the selected model are as follows φ 1 = φ 2 = φ 3 = φ 4 = φ 5 = Constant =
40 Case study 5 (Contd.) Significance of residual mean Model η(e) t 0.95 (N ) ARMA(5,0)
41 Case study 5 (Contd.) Significance of periodicities: Periodicity η F 0.95 (2,239 ) 1 st nd rd th th th
42 Case study 5 (Contd.) Whittle s white noise test: Model η F 0.95 (2, N-2 ) ARMA(5,0)
43 Case study 5 (Contd.) Model MSE AR(1) AR(2) AR(3) AR(4) AR(5) AR(6) ARMA(1,1) ARMA(1,2) ARMA(2,1) ARMA(2,2) ARMA(3,1) ARMA(3,2)
44 Case study 5 (Contd.) ARMA(1, 2) is selected with least MSE value for one step forecasting The parameters for the selected model are as follows φ 1 = θ 1 = θ 2 = Constant =
45 Case study 5 (Contd.) Significance of residual mean Model η(e) t 0.95 (N ) ARMA(1, 2)
46 Case study 5 (Contd.) Significance of periodicities: Periodicity η F 0.95 (2,239 ) 1 st nd rd th th th
47 Case study 5 (Contd.) Whittle s white noise test: Model η F 0.95 (2, N-2 ) ARMA(1, 2)
48 Markov Chains
49 Markov Chains A Markov chain is a stochastic process with the property that value of process X t at time t depends on its value at time t-1 and not on the sequence of other values (X t-2, X t-3,. X 0 ) that the process passed through in arriving at X t-1. [,,... ] = [ ] PX X X X PX X t t 1 t 2 0 t t 1 Single step Markov chain 49
50 Markov Chains P X = a X = a t j t 1 i This conditional probability gives the probability at time t will be in state j, given that the process was in state i at time t-1. The conditional probability is independent of the states occupied prior to t-1. For example, if X t-1 is a dry day, we would be interested in the probability that X t is a dry day or a wet day. This probability is commonly called as transition probability 50
51 Markov Chains P X = a X = a = P t t j t 1 i ij P ij t Usually written as indicating the probability of a step from a i to a j at time t. If P ij is independent of time, then the Markov chain is said to be homogeneous. P t ij = t P + τ ij i.e., v t and τ the transition probabilities remain same across time 51
52 Markov Chains Transition Probability Matrix(TPM): P t+1 t m P P P.. P P P P.. P P m m 31 =.. m.. P P P m1 m2 mm m x m 52
53 Markov Chains m j= 1 P ij = 1 V i Elements in any row of TPM sum to unity TPM can be estimated from observed data by enumerating the number of times the observed data went from state i to j P j (n) is the probability of being in state j in time step n. 53
54 Markov Chains p j (0) is the probability of being in state j in period t = 0. ( 0) ( 0) ( 0) ( 0) p = p1 p2.. p. Probability m ( n) ( n) ( n) ( n) 1 2 If p (0) is given and TPM is given 1 m 1 m vector at time 0 p = p p.. p. Probability ( 1) ( 0) p = p P m vector at time n 54
55 Markov Chains P P P.. P P P P.. P ( 1) ( 0) ( 0) ( 0) p p p p P m m = m 31 And so on ( 0) ( 0) ( 0). P P P m1 m2 mm = p1 P11+ p2 P pm Pm 1. Probability of ( 0) ( 0) ( 0) going to state 1 = p1 P12+ p2 P pm Pm2. Probability of going to state 2 55
56 Markov Chains Therefore In general, ( 1) ( 1) ( 1) ( 1) p p p p = m 1 m ( 2) ( 1) p = p P ( 0) = p P P ( ) = p P 0 2 ( n) ( 0) p = p P n 56
57 Markov Chains As the process advances in time, p j (n) becomes less dependent on p (0) The probability of being in state j after a large number of time steps becomes independent of the initial state of the process. The process reaches a steady state at large n p= p P n As the process reaches steady state, the probability vector remains constant 57
58 Example 1 Consider the TPM for a 2-state first order homogeneous Markov chain as TPM = State 1 is a non-rainy day and state 2 is a rainy day Obtain the 1. probability of day 1 is non-rainfall day / day 0 is rainfall day 2. probability of day 2 is rainfall day / day 0 is non-rainfall day 3. probability of day 100 is rainfall day / day 0 is non-rainfall day 58
59 Example 1 (contd.) 1. probability of day 1 is non-rainfall day / day 0 is rainfall day TPM No rain rain No rain = rain The probability is probability of day 2 is rainfall day / day 0 is nonrainfall day ( ) ( ) p = p P
60 Example 1 (contd.) ( 2) p = [ ] = [ ] The probability is probability of day 100 is rainfall day / day 0 is nonrainfall day ( n) ( 0) p = p P n 60
61 Example 1 (contd.) 2 P = P P = = P = P P = P = P P = P = P P =
62 Example 1 (contd.) Steady state probability p = [ ] For steady state, n p= p P = [ ] = [ ] 62
INDIAN INSTITUTE OF SCIENCE STOCHASTIC HYDROLOGY. Lecture -18 Course Instructor : Prof. P. P. MUJUMDAR Department of Civil Engg., IISc.
INDIAN INSTITUTE OF SCIENCE STOCHASTIC HYDROLOGY Lecture -18 Course Instructor : Prof. P. P. MUJUMDAR Department of Civil Engg., IISc. Summary of the previous lecture Model selection Mean square error
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 informationINDIAN INSTITUTE OF SCIENCE STOCHASTIC HYDROLOGY. Lecture -33 Course Instructor : Prof. P. P. MUJUMDAR Department of Civil Engg., IISc.
INDIAN INSTITUTE OF SCIENCE STOCHASTIC HYDROLOGY Lecture -33 Course Instructor : Prof. P. P. MUJUMDAR Department of Civil Engg., IISc. Summary of the previous lecture Regression on Principal components
More informationStochastic models of streamflow: some case studies
Hydrological Sciences Journal ISSN: 0262-6667 (Print) 2150-3435 (Online) Journal homepage: http://www.tandfonline.com/loi/thsj20 Stochastic models of streamflow: some case studies P. P. MUJUMDAR & D. NAGESH
More informationINDIAN INSTITUTE OF SCIENCE STOCHASTIC HYDROLOGY. Course Instructor : Prof. P. P. MUJUMDAR Department of Civil Engg., IISc.
INDIAN INSTITUTE OF SCIENCE STOCHASTIC HYDROLOGY Course Instructor : Prof. P. P. MUJUMDAR Department of Civil Engg., IISc. Course Contents Introduction to Random Variables (RVs) Probability Distributions
More informationINDIAN INSTITUTE OF SCIENCE STOCHASTIC HYDROLOGY. Lecture -35 Course Instructor : Prof. P. P. MUJUMDAR Department of Civil Engg., IISc.
INDIAN INSTITUTE OF SCIENCE STOCHASTIC HYDROLOGY Lecture -35 Course Instructor : Prof. P. P. MUJUMDAR Department of Civil Engg., IISc. Summary of the previous lecture Multivariate stochastic models Multisite
More informationINDIAN INSTITUTE OF SCIENCE STOCHASTIC HYDROLOGY. Lecture -27 Course Instructor : Prof. P. P. MUJUMDAR Department of Civil Engg., IISc.
INDIAN INSTITUTE OF SCIENCE STOCHASTIC HYDROLOGY Lecture -27 Course Instructor : Prof. P. P. MUJUMDAR Department of Civil Engg., IISc. Summary of the previous lecture Frequency factors Normal distribution
More informationINDIAN INSTITUTE OF SCIENCE STOCHASTIC HYDROLOGY. Lecture -30 Course Instructor : Prof. P. P. MUJUMDAR Department of Civil Engg., IISc.
INDIAN INSTITUTE OF SCIENCE STOCHASTIC HYDROLOGY Lecture -30 Course Instructor : Prof. P. P. MUJUMDAR Department of Civil Engg., IISc. Summary of the previous lecture IDF relationship Procedure for creating
More informationWater Resources Systems: Modeling Techniques and Analysis
INDIAN INSTITUTE OF SCIENCE Water Resources Systems: Modeling Techniques and Analysis Lecture - 9 Course Instructor : Prof. P. P. MUJUMDAR Department of Civil Engg., IISc. Summary of the previous lecture
More informationModule 8. Lecture 3: Markov chain
Lecture 3: Markov chain A Markov chain is a stochastic rocess having the roerty that the value of the rocess X t at time t, deends only on its value at time t-1, X t-1 and not on the sequence X t-2, X
More informationINDIAN INSTITUTE OF SCIENCE STOCHASTIC HYDROLOGY. Lecture -29 Course Instructor : Prof. P. P. MUJUMDAR Department of Civil Engg., IISc.
INDIAN INSTITUTE OF SCIENCE STOCHASTIC HYDROLOGY Lecture -29 Course Instructor : Prof. P. P. MUJUMDAR Department of Civil Engg., IISc. Summary of the previous lecture Goodness of fit Chi-square test χ
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 informationWater Resources Systems: Modeling Techniques and Analysis
INDIAN INSTITUTE OF SCIENCE Water Resources Systems: Modeling Techniques and Analysis Lecture - 20 Course Instructor : Prof. P. P. MUJUMDAR Department of Civil Engg., IISc. 1 Summary of the previous lecture
More informationForecasting of meteorological drought using ARIMA model
Indian J. Agric. Res., 51 (2) 2017 : 103-111 Print ISSN:0367-8245 / Online ISSN:0976-058X AGRICULTURAL RESEARCH COMMUNICATION CENTRE www.arccjournals.com/www.ijarjournal.com Forecasting of meteorological
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 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 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 informationSwitching Regime Estimation
Switching Regime Estimation Series de Tiempo BIrkbeck March 2013 Martin Sola (FE) Markov Switching models 01/13 1 / 52 The economy (the time series) often behaves very different in periods such as booms
More information1. 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 informationSTAT 153: Introduction to Time Series
STAT 153: Introduction to Time Series Instructor: Aditya Guntuboyina Lectures: 12:30 pm - 2 pm (Tuesdays and Thursdays) Office Hours: 10 am - 11 am (Tuesdays and Thursdays) 423 Evans Hall GSI: Brianna
More informationTime series models in the Frequency domain. The power spectrum, Spectral analysis
ime series models in the Frequency domain he power spectrum, Spectral analysis Relationship between the periodogram and the autocorrelations = + = ( ) ( ˆ α ˆ ) β I Yt cos t + Yt sin t t= t= ( ( ) ) cosλ
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 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 informationIf we want to analyze experimental or simulated data we might encounter the following tasks:
Chapter 1 Introduction If we want to analyze experimental or simulated data we might encounter the following tasks: Characterization of the source of the signal and diagnosis Studying dependencies Prediction
More informationModeling daily precipitation in Space and Time
Space and Time SWGen - Hydro Berlin 20 September 2017 temporal - dependence Outline temporal - dependence temporal - dependence Stochastic Weather Generator Stochastic Weather Generator (SWG) is a stochastic
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 informationMinitab Project Report Assignment 3
3.1.1 Simulation of Gaussian White Noise Minitab Project Report Assignment 3 Time Series Plot of zt Function zt 1 0. 0. zt 0-1 0. 0. -0. -0. - -3 1 0 30 0 50 Index 0 70 0 90 0 1 1 1 1 0 marks The series
More informationTime Series Analysis
Time Series Analysis hm@imm.dtu.dk Informatics and Mathematical Modelling Technical University of Denmark DK-2800 Kgs. Lyngby 1 Outline of the lecture Input-Output systems The z-transform important issues
More informationINDIAN INSTITUTE OF SCIENCE STOCHASTIC HYDROLOGY. Lecture -36 Course Instructor : Prof. P. P. MUJUMDAR Department of Civil Engg., IISc.
INDIAN INSTITUTE OF SCIENCE STOCHASTIC HYDROLOGY Lecture -36 Course Instructor : Prof. P. P. MUJUMDAR Department of Civil Engg., IISc. Summary of the previous lecture Multivariate stochastic models Matalas
More informationVolatility. Gerald P. Dwyer. February Clemson University
Volatility Gerald P. Dwyer Clemson University February 2016 Outline 1 Volatility Characteristics of Time Series Heteroskedasticity Simpler Estimation Strategies Exponentially Weighted Moving Average Use
More informationAkaike criterion: Kullback-Leibler discrepancy
Model choice. Akaike s criterion Akaike criterion: Kullback-Leibler discrepancy Given a family of probability densities {f ( ; ψ), ψ Ψ}, Kullback-Leibler s index of f ( ; ψ) relative to f ( ; θ) is (ψ
More informationAppendix of the paper: Are interest rate options important for the assessment of interest rate risk?
Appendix of the paper: Are interest rate options important for the assessment of interest rate risk? Caio Almeida,a, José Vicente b a Graduate School of Economics, Getulio Vargas Foundation b Research
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 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 informationChapter 12: An introduction to Time Series Analysis. Chapter 12: An introduction to Time Series Analysis
Chapter 12: An introduction to Time Series Analysis Introduction In this chapter, we will discuss forecasting with single-series (univariate) Box-Jenkins models. The common name of the models is Auto-Regressive
More informationSeasonal 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 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 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 information(x t. x t +1. TIME SERIES (Chapter 8 of Wilks)
45 TIME SERIES (Chapter 8 of Wilks) In meteorology, the order of a time series matters! We will assume stationarity of the statistics of the time series. If there is non-stationarity (e.g., there is a
More informationRegression 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 information6. The econometrics of Financial Markets: Empirical Analysis of Financial Time Series. MA6622, Ernesto Mordecki, CityU, HK, 2006.
6. The econometrics of Financial Markets: Empirical Analysis of Financial Time Series MA6622, Ernesto Mordecki, CityU, HK, 2006. References for Lecture 5: Quantitative Risk Management. A. McNeil, R. Frey,
More informationAutomatic Autocorrelation and Spectral Analysis
Piet M.T. Broersen Automatic Autocorrelation and Spectral Analysis With 104 Figures Sprin ger 1 Introduction 1 1.1 Time Series Problems 1 2 Basic Concepts 11 2.1 Random Variables 11 2.2 Normal Distribution
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 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 informationPure Random process Pure Random Process or White Noise Process: is a random process {X t, t 0} which has: { σ 2 if k = 0 0 if k 0
MODULE 9: STATIONARY PROCESSES 7 Lecture 2 Autoregressive Processes 1 Moving Average Process Pure Random process Pure Random Process or White Noise Process: is a random process X t, t 0} which has: E[X
More informationTHE ROYAL STATISTICAL SOCIETY 2009 EXAMINATIONS SOLUTIONS GRADUATE DIPLOMA MODULAR FORMAT MODULE 3 STOCHASTIC PROCESSES AND TIME SERIES
THE ROYAL STATISTICAL SOCIETY 9 EXAMINATIONS SOLUTIONS GRADUATE DIPLOMA MODULAR FORMAT MODULE 3 STOCHASTIC PROCESSES AND TIME SERIES The Society provides these solutions to assist candidates preparing
More informationChapter 5 Identifying hydrological persistence
103 Chapter 5 Identifying hydrological persistence The previous chapter demonstrated that hydrologic data from across Australia is modulated by fluctuations in global climate modes. Various climate indices
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 informationMGR-815. Notes for the MGR-815 course. 12 June School of Superior Technology. Professor Zbigniew Dziong
Modeling, Estimation and Control, for Telecommunication Networks Notes for the MGR-815 course 12 June 2010 School of Superior Technology Professor Zbigniew Dziong 1 Table of Contents Preface 5 1. Example
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 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 informationTime Series Analysis
Time Series Analysis hm@imm.dtu.dk Informatics and Mathematical Modelling Technical University of Denmark DK-2800 Kgs. Lyngby 1 Outline of the lecture Chapter 9 Multivariate time series 2 Transfer function
More 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 informationReliability and Risk Analysis. Time Series, Types of Trend Functions and Estimates of Trends
Reliability and Risk Analysis Stochastic process The sequence of random variables {Y t, t = 0, ±1, ±2 } is called the stochastic process The mean function of a stochastic process {Y t} is the function
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 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 informationTAKEHOME FINAL EXAM e iω e 2iω e iω e 2iω
ECO 513 Spring 2015 TAKEHOME FINAL EXAM (1) Suppose the univariate stochastic process y is ARMA(2,2) of the following form: y t = 1.6974y t 1.9604y t 2 + ε t 1.6628ε t 1 +.9216ε t 2, (1) where ε is i.i.d.
More informationModelling Non-linear and Non-stationary Time Series
Modelling Non-linear and Non-stationary Time Series Chapter 7(extra): (Generalized) Hidden Markov Models Henrik Madsen Lecture Notes September 2016 Henrik Madsen (02427 Adv. TS Analysis) Lecture Notes
More informationStochastic models of streamflow: some case studies
Stochastic models of streamflow: some case studies t P. P. MUWAR* & D. NAGESH KUMAR Depamnent of Civil En@eenn& Indian Im'hue of Science, Bangalore, 560 012, India. - Abstract Ten candidate models of the
More informationStationary Stochastic Time Series Models
Stationary Stochastic Time Series Models When modeling time series it is useful to regard an observed time series, (x 1,x,..., x n ), as the realisation of a stochastic process. In general a stochastic
More informationWhite Noise Processes (Section 6.2)
White Noise Processes (Section 6.) Recall that covariance stationary processes are time series, y t, such. E(y t ) = µ for all t. Var(y t ) = σ for all t, σ < 3. Cov(y t,y t-τ ) = γ(τ) for all t and τ
More information18.S096 Problem Set 4 Fall 2013 Time Series Due Date: 10/15/2013
18.S096 Problem Set 4 Fall 2013 Time Series Due Date: 10/15/2013 1. Covariance Stationary AR(2) Processes Suppose the discrete-time stochastic process {X t } follows a secondorder auto-regressive process
More informationBasic Probability space, sample space concepts and order of a Stochastic Process
The Lecture Contains: Basic Introduction Basic Probability space, sample space concepts and order of a Stochastic Process Examples Definition of Stochastic Process Marginal Distributions Moments Gaussian
More informationLecture 2: ARMA(p,q) models (part 2)
Lecture 2: ARMA(p,q) models (part 2) Florian Pelgrin University of Lausanne, École des HEC Department of mathematics (IMEA-Nice) Sept. 2011 - Jan. 2012 Florian Pelgrin (HEC) Univariate time series Sept.
More information1. How can you tell if there is serial correlation? 2. AR to model serial correlation. 3. Ignoring serial correlation. 4. GLS. 5. Projects.
1. How can you tell if there is serial correlation? 2. AR to model serial correlation. 3. Ignoring serial correlation. 4. GLS. 5. Projects. 1) Identifying serial correlation. Plot Y t versus Y t 1. See
More informationVolume 11 Issue 6 Version 1.0 November 2011 Type: Double Blind Peer Reviewed International Research Journal Publisher: Global Journals Inc.
Volume 11 Issue 6 Version 1.0 2011 Type: Double Blind Peer Reviewed International Research Journal Publisher: Global Journals Inc. (USA) Online ISSN: & Print ISSN: Abstract - Time series analysis and forecasting
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 informationChapter 3 - Temporal processes
STK4150 - Intro 1 Chapter 3 - Temporal processes Odd Kolbjørnsen and Geir Storvik January 23 2017 STK4150 - Intro 2 Temporal processes Data collected over time Past, present, future, change Temporal aspect
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 informationSTAT Financial Time Series
STAT 6104 - Financial Time Series Chapter 2 - Probability Models Chun Yip Yau (CUHK) STAT 6104:Financial Time Series 1 / 34 Agenda 1 Introduction 2 Stochastic Process Definition 1 Stochastic Definition
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 informationTime Series 2. Robert Almgren. Sept. 21, 2009
Time Series 2 Robert Almgren Sept. 21, 2009 This week we will talk about linear time series models: AR, MA, ARMA, ARIMA, etc. First we will talk about theory and after we will talk about fitting the models
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 informationAsian Economic and Financial Review. SEASONAL ARIMA MODELLING OF NIGERIAN MONTHLY CRUDE OIL PRICES Ette Harrison Etuk
Asian Economic and Financial Review journal homepage: http://aessweb.com/journal-detail.php?id=5002 SEASONAL ARIMA MODELLING OF NIGERIAN MONTHLY CRUDE OIL PRICES Ette Harrison Etuk Department of Mathematics/Computer
More informationAPPLIED ECONOMETRIC TIME SERIES 4TH EDITION
APPLIED ECONOMETRIC TIME SERIES 4TH EDITION Chapter 2: STATIONARY TIME-SERIES MODELS WALTER ENDERS, UNIVERSITY OF ALABAMA Copyright 2015 John Wiley & Sons, Inc. Section 1 STOCHASTIC DIFFERENCE EQUATION
More informationTime Series Analysis Model for Rainfall Data in Jordan: Case Study for Using Time Series Analysis
American Journal of Environmental Sciences 5 (5): 599-604, 2009 ISSN 1553-345X 2009 Science Publications Time Series Analysis Model for Rainfall Data in Jordan: Case Study for Using Time Series Analysis
More informationModeling climate variables using time series analysis in arid and semi arid regions
Vol. 9(26), pp. 2018-2027, 26 June, 2014 DOI: 10.5897/AJAR11.1128 Article Number: 285C77845733 ISSN 1991-637X Copyright 2014 Author(s) retain the copyright of this article http://www.academicjournals.org/ajar
More informationEconometric Forecasting
Robert M. Kunst robert.kunst@univie.ac.at University of Vienna and Institute for Advanced Studies Vienna October 1, 2014 Outline Introduction Model-free extrapolation Univariate time-series models Trend
More informationTesting for IID Noise/White Noise: I
Testing for IID Noise/White Noise: I want to be able to test null hypothesis that time series {x t } or set of residuals {r t } is IID(0, σ 2 ) or WN(0, σ 2 ) many such tests exist, including informal
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 informationLecture 3-4: Probability models for time series
Lecture 3-4, page 1 Lecture 3-4: Probability models for time series Outline of lesson 3-4 (chapter 3) The most heavy (and theoretical) part (and lesson) AR(p) processes MA(q) processes ARMA(p,q) processes
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 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 informationMODELLING TIME SERIES WITH CONDITIONAL HETEROSCEDASTICITY
MODELLING TIME SERIES WITH CONDITIONAL HETEROSCEDASTICITY The simple ARCH Model Eva Rubliková Ekonomická univerzita Bratislava Manuela Magalhães Hill Department of Quantitative Methods, INSTITUTO SUPERIOR
More informationCCNY. BME I5100: Biomedical Signal Processing. Stochastic Processes. Lucas C. Parra Biomedical Engineering Department City College of New York
BME I5100: Biomedical Signal Processing Stochastic Processes Lucas C. Parra Biomedical Engineering Department CCNY 1 Schedule Week 1: Introduction Linear, stationary, normal - the stuff biology is not
More informationEcon 427, Spring Problem Set 3 suggested answers (with minor corrections) Ch 6. Problems and Complements:
Econ 427, Spring 2010 Problem Set 3 suggested answers (with minor corrections) Ch 6. Problems and Complements: 1. (page 132) In each case, the idea is to write these out in general form (without the lag
More informationMidterm Suggested Solutions
CUHK Dept. of Economics Spring 2011 ECON 4120 Sung Y. Park Midterm Suggested Solutions Q1 (a) In time series, autocorrelation measures the correlation between y t and its lag y t τ. It is defined as. ρ(τ)
More informationFRACTIONAL VECTOR AUTOREGRESSIVE MODELS FOR ENVIRONMENTAL DATA
FRACTIONAL VECTOR AUTOREGRESSIVE MODELS FOR ENVIRONMENTAL DATA Mariagrazia Granturco Eugenia Nissi Dipartimento di Metodi Quantitativi e Teoria Economica Università degli Studi G. D Annunzio Viale Pindaro,
More informationRainfall Drought Simulating Using Stochastic SARIMA Models for Gadaref Region, Sudan
MPRA Munich Personal RePEc Archive Rainfall Drought Simulating Using Stochastic SARIMA Models for Gadaref Region, Sudan Hisham Moahmed Hassan and Tariq Mahgoub Mohamed University of Khartoum, Khartoum
More informationRainfall Disaggregation Methods: Theory and Applications
Workshop on Statistical and Mathematical Methods for Hydrological Analysis Rome, May 23 Universitá degli Studi di Roma La Sapienza CNR - IRPI Istituto di Ricerca per la Protezione Idrogeologica GNDCI Gruppo
More informationTIME SERIES MODELING OF MONTHLY RAINFALL IN ARID AREAS: CASE STUDY FOR SAUDI ARABIA
American Journal of Environmental Sciences 10 (3): 277-282, 2014 ISSN: 1553-345X 2014 Science Publication doi:10.3844/ajessp.2014.277.282 Published Online 10 (3) 2014 (http://www.thescipub.com/ajes.toc)
More informationTime Series Analysis of Monthly Rainfall data for the Gadaref rainfall station, Sudan, by Sarima Methods
International Journal of Scientific Research in Knowledge, 2(7), pp. 320-327, 2014 Available online at http://www.ijsrpub.com/ijsrk ISSN: 2322-4541; 2014 IJSRPUB http://dx.doi.org/10.12983/ijsrk-2014-p0320-0327
More information9. Multivariate Linear Time Series (II). MA6622, Ernesto Mordecki, CityU, HK, 2006.
9. Multivariate Linear Time Series (II). MA6622, Ernesto Mordecki, CityU, HK, 2006. References for this Lecture: Introduction to Time Series and Forecasting. P.J. Brockwell and R. A. Davis, Springer Texts
More informationRainfall Disaggregation Methods: Theory and Applications
Workshop on Statistical and Mathematical Methods for Hydrological Analysis Rome, May 23 Universitá degli Studi di Roma La Sapienza CNR - IRPI Istituto di Ricerca per la Protezione Idrogeologica GNDCI Gruppo
More informationPsychology 282 Lecture #4 Outline Inferences in SLR
Psychology 282 Lecture #4 Outline Inferences in SLR Assumptions To this point we have not had to make any distributional assumptions. Principle of least squares requires no assumptions. Can use correlations
More information1 Linear Difference Equations
ARMA Handout Jialin Yu 1 Linear Difference Equations First order systems Let {ε t } t=1 denote an input sequence and {y t} t=1 sequence generated by denote an output y t = φy t 1 + ε t t = 1, 2,... with
More informationLesson 8: Testing for IID Hypothesis with the correlogram
Lesson 8: Testing for IID Hypothesis with the correlogram Dipartimento di Ingegneria e Scienze dell Informazione e Matematica Università dell Aquila, umberto.triacca@ec.univaq.it Testing for i.i.d. Hypothesis
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 informationWeek 5 Quantitative Analysis of Financial Markets Characterizing Cycles
Week 5 Quantitative Analysis of Financial Markets Characterizing Cycles Christopher Ting http://www.mysmu.edu/faculty/christophert/ Christopher Ting : christopherting@smu.edu.sg : 6828 0364 : LKCSB 5036
More information5 Transfer function modelling
MSc Further Time Series Analysis 5 Transfer function modelling 5.1 The model Consider the construction of a model for a time series (Y t ) whose values are influenced by the earlier values of a series
More information