Unit root problem, solution of difference equations Simple deterministic model, question of unit root
|
|
- Jesse Lawrence
- 5 years ago
- Views:
Transcription
1 Unit root problem, solution of difference equations Simple deterministic model, question of unit root (1 φ 1 L)X t = µ, Solution X t φ 1 X t 1 = µ X t = A + Bz t with unknown z and unknown A (clearly X t 1 = A + Bz t 1 ) 0 = X t φ 1 X t 1 µ = A + Bz t φ 1 (A + Bz t 1 ) µ = A + Bz t φ 1 ( A + B z zt ) µ ( = µ + A(1 φ 1 ) +B 1 φ ) 1 }{{} z =0 }{{} z t
2 z = φ 1, A = µ 1 φ 1 Elimination B by the initial condition X t0 = µ 1 φ 1 + Bφ t 0 1 B = X t = µ + 1 φ 1 ( X t µ ) 1 φ 1 X t = µ 1 φ 1 + Bφ t 1 [ X t0 µ 1 φ 1 ] [ X t0 µ = φ t φ 1 ( X t0 µ 1 φ 1 X t µ 1 φ 1 X t0 µ = φ t t φ 1 ] φ t 0 1 φ t 1 ) φ t t 0 1 Qualitative difference for φ 1 < 1 and φ 1 > 1 (unit root means z = φ 1 = 1) φ 1 = 0.8: t 1 t 0 = 0, φ 0 1 = 1; t t 0 = 1, φ 1 1 = 0.8; t t 0 = 2, φ 2 1 = 0.64; t t 0 = 3; φ 3 1 φ 1 = 1.2: t t 0 = 10, φ 10 1 = 0.52; φ101 = 6.19; φ100 = 0.107, φ1001 = ; 1 = ;
3 (1 φ 1 (1/z) φ 2 (1/z) 2 φ 3 (1/z) 3... φ p (1/z) p )z t X 0 = 0 Unit root signature of nonstationary behaviour z = φ 1 = 1 Take homogeneous equation (1 φ 1 L)X t = 0, (1 L)X t = 0 formally via the "trick" L = 1/z we get characteristic equation (1 φ 1 /z) = 0 where the unit root is achieved for z = φ 1 = 1 Generalization for AR(p) without stochastic component: p X t φ i X t i = 0 i=1 X t φ 1 X t 1 φ 2 X t 2 φ 3 X t 3... φ p X t p = 0 (1 φ 1 L φ 2 L 2 φ 3 L 3... φ p L p )X t = 0 expected form of solution X t = z t X 0, t = 0, 1, 2...
4 [ 1 φ 1 (1/z) φ 2 (1/z) 2 φ 3 (1/z) 3... φ p (1/z) p] z t X 0 = 0 characteristic equation z p φ 1 z p 1 φ 2 z p 2 φ 3 z p 3... φ p = 0 generally has exactly p roots (z λ 1 )(z λ 2 )... (z λ p ) = 0 p (z λ k ) k = 0 k=1 For non-identical (also complex) roots X t = B 1 λ t 1 + B 2 λ t B p λ t p when roots are identical for e.g. p = 2 X t = B 1 λ t 1 + B 2 tλ t 1
5 Augmented Dickey-Fuller (ADF) unit root test (Rozšírený DF test) The test is widely used by practioners. In statistics, the traditional DickeyFuller test tests (1979) whether a unit root is present in an autoregressive model. In time series models in econometrics (the application of statistical methods to economics), a unit root is a feature of processes that evolve through time that can cause problems in statistical inference if it is not adequately dealt with. Such a process is non-stationary. (If the other roots of the characteristic equation lie inside the unit circle that is, have a modulus (absolute value) less than one then the first difference of the process will be stationary). Other words: series of the first difference (previously considered integrated of order one) turned to be stationary.
6 R: test available in adf.test() tseries() [package for Time series analysis and computational finance] The testing procedure is applied to the model structure X t = α }{{} drift + δ t }{{} lin.trend + ϱ }{{} X t 1 + par. tested p ψ j X t j X t j = X t j X t j 1, X t = X t X t 1 number of lags of autoregression; default: p = (length of d.s. 1) 1/3 ; (example: length of data series = 65, p = 4); Swert criterion: 12 (length of d.s./100) 1/4 Alternative form: p X t = α + δt + (1 + ϱ) X }{{} t 1 + ψ j X t j φ 1 in former notation j=1 p j=1 ψ j X t j the autoregression terms absorbing dynamic j=1
7 H 0 : ϱ = 0 for given model has a unit root, process is I (1), i.e. the data needs to be differenced to make it stationary; H 1 : ϱ < 1 i.e. the data is trend stationary, process is I (0) and needs to be analyzed by means of using a trend in the regression model instead of differencing the data; The unit root test is carried out under H 0 : ϱ = 0 against the alternative hypothesis of ϱ < 0. Note 1: Caution: Sometimes if data is exponentially trending then you might need to take the log of the data first before differencing it. Note 2: ϱ = 0 means X t X t 1 (i.e. random walk process)
8 test statistic [one-sided left tail test]: ADF τ = ˆϱ SE(ˆϱ) (remind that t score = ˆϱ SE(ˆϱ) ) SE(ˆϱ) is the usual standard error estimate critical value: t 1 α(n) If the test statistic is less than the critical value t 1 α(n), then the H 0 ( ϱ = 0) is rejected and no unit root is present. Example: A model that includes a constant and a time trend is estimated using sample of n = 50 observations and yields the ADF τ statistic of This is more negative than the tabulated t = 3.50, so at the 95 per cent level the null hypothesis of a unit root will be rejected.
9 R implementation > x=arima.sim(n = 1000, list(ma = c(-0.2, 0.2)),sd=1) > adf.test(x) Augmented Dickey-Fuller Test data: x Dickey-Fuller = , Lag order = 9, p-value = 0.01 alternative hypothesis: stationary Warning message: In adf.test(x) : p-value smaller than printed p-value [Note: rejection of null hypothesis that data has a unit root]
10 > x= arima.sim(n = 1000, list(ma = c(-0.2, 0.2)),sd=1) + 0.1*(c(1:1000)^3) > adf.test(x) Augmented Dickey-Fuller Test data: x Dickey-Fuller = , Lag order = 9, p-value = 0.99 alternative hypothesis: stationary Warning message: In adf.test(x) : p-value greater than printed p-value
11 KPSS test Kwiatkowski, Phillips, Schmidt, Shin (1992) observed time series x t is expressed as the sum x t = T t + RW t + u t : - deterministic trend T t - pure random walk RW t (i.e. I (1)), RW t = RW t 1 + ɛ t - stationary error term u t (i.e. I (0)) H 0 : σ 2 ɛ = 0 observable time series is stationary around a deterministic trend, x t is I (0) H 1 : σ 2 ɛ > 0
12 in R: > s<-rnorm(1000) > kpss.test(s) > > KPSS Test for Level Stationarity > > data: s > KPSS Level = , Truncation lag parameter = 7, > p-value = stationarity cannot be rejected here...
13 Second order Stochastic Difference Equation and Business Cycles Samuelson (1939) applied only second order difference equations to explain the business cycles: Investment I t is dependent on the changes in income rather then the level of income I t = b 1 (Y t 1 Y t 2 ) Consumption C t depends on the past income by the first order nonhomogeneous equation C t = b 2 Y t 1 + b 3 Production/consumption balance with error term Y t = C t + I t + ɛ t
14 Consequence - accelerator principle Y t = b 1 (Y t 1 Y t 2 ) + b 2 Y t 1 + b 3 + ɛ t simulation in R b3<-1.04 b2<-0.2 b1<-1.01 Y1<-1.0 Y2<-1.04 Y0<-1 for(i in 1:200){ Y2<-Y1 Y1<-Y0 Y0<-(b1*(Y1-Y2)+b2*Y1+b3) print(paste(y0)) }
15 Y "./Simul_Samuelson0.txt" time b3<-1.04 b2<-0.2 b1<-1.01 Y1<-1.0 Y2<-1.04 Y0<-1 for(i in 1:200){ Y2<-Y1 Y1<-Y0 Y0<-(b1*(Y1-Y2)+b2*Y1+b3) + rnorm(1,mean=0,sd=0.3) }
16 8 6 "./Simul_Samuelson1.txt" "./Simul_Samuelson2.txt" 4 Y time
17 Solution of the homogeneous (non-stochastic) Samuelson s system z =..., b 3 = 0 Y t = z t Y 0 Y t = b 1 (Y t 1 Y t 2 ) + b 2 Y t 1 z t Y 0 = b 1 (z t 1 Y 0 z t 2 Y 0 ) + b 2 Y 0 z t 1 ( 1 1 = b 1 z 1 ) 1 z 2 + b 2 z lag operator formally produces 1/z ( 1 b 1 z 1 ) 1 z 2 + b 2 z 1 = 0 b 1 + b 2 b = 0
18 z 2 (b 1 + b 2 )z + b 1 = 0 z 1,2 = b 1 + b 2 ± (b 1 + b 2 ) 2 4b 1 2 b 1 = 1.01, b 2 = 0.2, b 3 = 1.04 z 1,2 = 1.21 ± = ± i z 1,2 = , arccos (0.605/ z 1,2 ) = (2π) Y t = Y 0 z t = ( ) t [cos(±(2π)0.147t) + i sin(±(2π)0.147t)] ( = ( ) t exp(i2π0.147t) = ( ) t exp i2π t ) 6.8
19 X Qualitative check of the period x x x ACF Lag
20 Comparison Samuelson s accelerator (b 3 = 0) with ARIMA(p,d,.) (1 φ 1 φ 2 L)(1 L)Y t = (1 L Lφ 1 + L 2 φ 1 L 2 φ 2 + L 3 φ 2 ) Y t = b 1 (Y t 1 Y t 2 ) + b 2 Y t 1 Try p = 1, d = 1 ( ) p 1 φ i L i (1 L) d Y t = 0 i=1 (1 φ 1 L)(1 L)Y t = (1 φ 1 L)(Y t Y t 1 ) = Y t Y t 1 φ 1 (LY t LY t 1 ) Try p = 2, d = 1 = Y t Y t 1 φ 1 (Y t 1 Y t 2 )
21 = Y t Y t 1 φ 1 Y t 1 + φ 1 Y t 2 φ 2 Y t 2 + φ 2 Y t 3 ) Try p = 2, d = 0 (1 φ 1 L φ 2 L 2 )Y t = Y t φ 1 Y t 1 φ 2 Y t 2 Y t φ 1 Y t 1 φ 2 Y t 2 = 0 Y t = φ 1 Y t 1 + φ 2 Y t 2 Y t = b 1 (Y t 1 Y t 2 ) + b 2 Y t 1 = (b 1 + b 2 ) Y }{{} t 1 + ( b 1 ) }{{} φ 1 Conclusion: Model is transformable to ARIMA(2,0,q) i.e. d = 0 model φ 2 Y t 2
22 Box Jenkins metodology 1. identification of model, mostly AR(1) making sure that model is stationary using plots of ACF, and PACF 2. estimation of the model parameters: say x t = 0.68x t 1 + ɛ t, ˆσ = using the maximum likelihood non-linear least-squares 3. diagnostics of the model: prediction abilities, Ljung-Box test Model identification by means of autocorrelation function: U-shaped after first q AR(p) MA(q) ARMA(p, q) ρ k k 0 does not exist k 0 = q k 0 does not exist ACF ρ k is U shaped ρ k becomes ρ kk k 0 = p k 0 does not exist k 0 does not exist PACF ρ kk is limited ρ kk becomes by U-shaped curve U-shaped after the first
23 Table: U denotes the curve of geometrically or sin-like decaying amplitudes. Rule 1 If the series has positive autocorrelations out to a high number of lags, then it probably needs a higher order of differencing. Rule 2 If the lag-1 autocorrelation is zero or negative, or the autocorrelations are all small a patternless, then the series does not need a higher order of differencing. If the lag-1 autocorrelation is -0.5 or more negative, the series may be overdifferenced. Beware of overdifferencing. Rule 3 The optimal order of differencing is then often the order at which the standard deviation is lowest. Rule 4 A model with no orders of differencing assumes that the original series is stationary. Rule 5 A model with no orders of differencing normally includes a constant term Rule 6 If the part of PACF of the differenciated series displays a sharp cutoff and/or the lag 1 autocorrelation is positive, i.e. the
24 series appears slightly underdifferenced then consider adding one or more AR terms to the model Rule 7 The autocorrelation function ACF of the differenced series displays a sharp cutoff and or the lag-1 autocorrelation function is negative. Rule 8 It is possible for an AR term and MA term to cancel each other s effects, so if a mixed AR-MA model seems to fit data also try a model with one fewer AR term. Rule 9 If there is a unit root in the AR part of the model the sum of AR coefficients is almost exactly 1 Rule 10 If the long-term forecasts appear erratic or unstable, there may be a unit root. Rule 11 If the series has a strong and consistent seasonal pattern, then you should use an order of seasonal differencing.
25 Autoregression of simulated time series of the price formation process Model of random process: p 0, p 1, p 2, p 3,..., p t,... if p t > p fun { p t+1 = p t 0.99} else p t+1 = p t if p t+1 < p fun {p t+1 = p t } else p t+1 = p t+1 p t+1 = p t+1 exp [ N(0, 0.01) ]
26 R project: simulation of price dynamics pfun<-1 # "fundamental" price p<-(pfun+0.01) # initial price setting vp<-p # auxilliary initial vector setting for pl Ns<-1000 # number of steps time=c(0:ns) # vector of times for(i in 1:Ns){ if(p>pfun) p<-(p * 0.99) if(p<pfun) p<-(p * 1.01) # deterministic reduction # deterministic increase p<-(p*exp(rnorm(1,mean=0,sd=0.01))) #randomization vp<-cbind(vp,p) } vp<-as.vector(vp)
27 pdf("tp_simul.pdf") plot(time,vp, type="l", pch=21, xlab="time", ylab="p") dev.off() p time
28 Testing of simulated price series for stationarity adf.test(vp) Augmented Dickey-Fuller Test data: vp Dickey-Fuller = , Lag order = 9, p-value = 0.01 alternative hypothesis: stationary Warning message: In adf.test(vp) : p-value smaller than printed p-value conclusion: reject unit root
29 Selection of the "best" arima model based on the Aikake information criterion The Akaike information criterion is a measure of the relative goodness of fit of a statistical model. It is grounded in the concept of information entropy, in effect offering a relative measure of the information lost when a given model is used to describe reality. In the general case, the AIC is AIC = 2k 2 ln(l), where k is the number of parameters in the statistical model, and L is the maximized value of the likelihood function for the estimated model. To apply AIC in practice, we start with a set of candidate models, and then find the models corresponding AIC values. There will almost always be information lost due to using one of the candidate models to represent the "true" model. We wish to select, from among n candidate models, the model that minimizes the estimated information loss. Denote the AIC values of the candidate models by AIC 1, AIC 2, AIC 3,..., AIC n. Let AIC min be the minimum of those values. Then e (AIC min AIC i )/2 can be interpreted as the relative probability that the i-th model minimizes the (estimated) information loss.
30 ########################################################################## Fit an ARIMA model to a univariate time series. Usage arima(x, order = c(0, 0, 0), seasonal = list(order = c(0, 0, 0), period = NA),... optim.control = list(), kappa = 1e6) # Example: R code for selection of the "best" model # Aikake s information criterion, select among 10 candidate models a1=arima(vp,order=c(0,0,0)) a2=arima(vp,order=c(0,0,1)) a3=arima(vp,order=c(0,1,0)) a4=arima(vp,order=c(1,0,0)) a5=arima(vp,order=c(1,1,0)) a6=arima(vp,order=c(0,1,1)) a7=arima(vp,order=c(1,0,1)) a8=arima(vp,order=c(1,1,1)) a8=arima(vp,order=c(0,0,2)) a9=arima(vp,order=c(0,2,0)) a10=arima(vp,order=c(2,0,0)) veca=as.vector(cbind(a1$aic,a2$aic,a3$aic,a4$aic, a5$aic,a6$aic,a7$aic,a8$aic,a9$aic,a10$aic)) which.min(veca)
31 Model Outputs: > which.min(veca) [1] 4 Check the best: > a4 Call: arima(x = vp, order = c(1, 0, 0)) Coefficients: ar1 intercept s.e sigma^2 estimated as : log likelihood = , aic = !!! intercept indicates fundamental price p_fun=1
32 ### Check another one = not the best > a3 Call: arima(x = vp, order = c(0, 1, 0)) sigma^2 estimated as : log likelihood = , aic =
33 Forecast by Arima predict.arima {stats} R Documentation Forecast from ARIMA fits Description: Forecast from models fitted by arima. predict(object, n.ahead = 1, newxreg = NULL, se.fit = TRUE,...) Arguments: object The result of an arima fit. n.ahead The number of steps ahead for which prediction is required.
34 Make the prediction: > a4=arima(x=vp, order=c(1,0,0)) > predict(a4,n.ahead=6) \$pred Time Series: Start = 1002 End = 1007 Frequency = 1 [1] \$se Time Series: Start = 1002 End = 1007 Frequency = 1 [1]
35 Testing of forecast vp1=vp[1:990] > a4=arima(vp1,order=c(1,0,0)) > vp2=vp[991:1000] > pre=predict(a4,n.ahead=10) Time Series: Start = 991 End = 1000 Frequency = 1 [1] [9] \$se Time Series: Start = 991 End = 1000 Frequency = 1 [1] [7] > vp2 [1] [8]
36 vpred=pre\$pred[1:10] > vp2-vpred [1] [6] compare errors with a mean standard deviation (vp2-vpred)/mean(se1) [1] [7]
Forecasting using R. Rob J Hyndman. 2.4 Non-seasonal ARIMA models. Forecasting using R 1
Forecasting using R Rob J Hyndman 2.4 Non-seasonal ARIMA models Forecasting using R 1 Outline 1 Autoregressive models 2 Moving average models 3 Non-seasonal ARIMA models 4 Partial autocorrelations 5 Estimation
More informationUnivariate, Nonstationary Processes
Univariate, Nonstationary Processes Jamie Monogan University of Georgia March 20, 2018 Jamie Monogan (UGA) Univariate, Nonstationary Processes March 20, 2018 1 / 14 Objectives By the end of this meeting,
More informationBCT Lecture 3. Lukas Vacha.
BCT Lecture 3 Lukas Vacha vachal@utia.cas.cz Stationarity and Unit Root Testing Why do we need to test for Non-Stationarity? The stationarity or otherwise of a series can strongly influence its behaviour
More informationTesting for non-stationarity
20 November, 2009 Overview The tests for investigating the non-stationary of a time series falls into four types: 1 Check the null that there is a unit root against stationarity. Within these, there are
More informationLecture 6a: Unit Root and ARIMA Models
Lecture 6a: Unit Root and ARIMA Models 1 2 Big Picture A time series is non-stationary if it contains a unit root unit root nonstationary The reverse is not true. For example, y t = cos(t) + u t has no
More informationNon-Stationary Time Series and Unit Root Testing
Econometrics II Non-Stationary Time Series and Unit Root Testing Morten Nyboe Tabor Course Outline: Non-Stationary Time Series and Unit Root Testing 1 Stationarity and Deviation from Stationarity Trend-Stationarity
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 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 information7. Integrated Processes
7. Integrated Processes Up to now: Analysis of stationary processes (stationary ARMA(p, q) processes) Problem: Many economic time series exhibit non-stationary patterns over time 226 Example: We consider
More informationPackage TimeSeries.OBeu
Type Package Package TimeSeries.OBeu Title Time Series Analysis 'OpenBudgets' Version 1.2.2 Date 2018-01-20 January 22, 2018 Estimate and return the needed parameters for visualisations designed for 'OpenBudgets'
More informationDecision 411: Class 9. HW#3 issues
Decision 411: Class 9 Presentation/discussion of HW#3 Introduction to ARIMA models Rules for fitting nonseasonal models Differencing and stationarity Reading the tea leaves : : ACF and PACF plots Unit
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 information7. Integrated Processes
7. Integrated Processes Up to now: Analysis of stationary processes (stationary ARMA(p, q) processes) Problem: Many economic time series exhibit non-stationary patterns over time 226 Example: We consider
More informationLecture 5: Unit Roots, Cointegration and Error Correction Models The Spurious Regression Problem
Lecture 5: Unit Roots, Cointegration and Error Correction Models The Spurious Regression Problem Prof. Massimo Guidolin 20192 Financial Econometrics Winter/Spring 2018 Overview Stochastic vs. deterministic
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 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 informationLecture 4a: ARMA Model
Lecture 4a: ARMA Model 1 2 Big Picture Most often our goal is to find a statistical model to describe real time series (estimation), and then predict the future (forecasting) One particularly popular model
More informationUnivariate linear models
Univariate linear models The specification process of an univariate ARIMA model is based on the theoretical properties of the different processes and it is also important the observation and interpretation
More informationNon-Stationary Time Series and Unit Root Testing
Econometrics II Non-Stationary Time Series and Unit Root Testing Morten Nyboe Tabor Course Outline: Non-Stationary Time Series and Unit Root Testing 1 Stationarity and Deviation from Stationarity Trend-Stationarity
More informationStat 153 Time Series. Problem Set 4
Stat Time Series Problem Set Problem I generated 000 observations from the MA() model with parameter 0.7 using the following R function and then fitted the ARMA(, ) model to the data. order = c(, 0, )),
More 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 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 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 informationNon-Stationary Time Series and Unit Root Testing
Econometrics II Non-Stationary Time Series and Unit Root Testing Morten Nyboe Tabor Course Outline: Non-Stationary Time Series and Unit Root Testing 1 Stationarity and Deviation from Stationarity Trend-Stationarity
More informationUnivariate ARIMA Models
Univariate ARIMA Models ARIMA Model Building Steps: Identification: Using graphs, statistics, ACFs and PACFs, transformations, etc. to achieve stationary and tentatively identify patterns and model components.
More informationFORECASTING 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 informationNonstationary Time Series:
Nonstationary Time Series: Unit Roots Egon Zakrajšek Division of Monetary Affairs Federal Reserve Board Summer School in Financial Mathematics Faculty of Mathematics & Physics University of Ljubljana September
More informationOn Consistency of Tests for Stationarity in Autoregressive and Moving Average Models of Different Orders
American Journal of Theoretical and Applied Statistics 2016; 5(3): 146-153 http://www.sciencepublishinggroup.com/j/ajtas doi: 10.11648/j.ajtas.20160503.20 ISSN: 2326-8999 (Print); ISSN: 2326-9006 (Online)
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 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 information11/18/2008. So run regression in first differences to examine association. 18 November November November 2008
Time Series Econometrics 7 Vijayamohanan Pillai N Unit Root Tests Vijayamohan: CDS M Phil: Time Series 7 1 Vijayamohan: CDS M Phil: Time Series 7 2 R 2 > DW Spurious/Nonsense Regression. Integrated but
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 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 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 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 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 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 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 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 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 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 informationARDL Cointegration Tests for Beginner
ARDL Cointegration Tests for Beginner Tuck Cheong TANG Department of Economics, Faculty of Economics & Administration University of Malaya Email: tangtuckcheong@um.edu.my DURATION: 3 HOURS On completing
More informationBox-Jenkins ARIMA Advanced Time Series
Box-Jenkins ARIMA Advanced Time Series www.realoptionsvaluation.com ROV Technical Papers Series: Volume 25 Theory In This Issue 1. Learn about Risk Simulator s ARIMA and Auto ARIMA modules. 2. Find out
More informationHomework 5, Problem 1 Andrii Baryshpolets 6 April 2017
Homework 5, Problem 1 Andrii Baryshpolets 6 April 2017 Total Private Residential Construction Spending library(quandl) Warning: package 'Quandl' was built under R version 3.3.3 Loading required package:
More 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 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 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 informationCointegrated VARIMA models: specification and. simulation
Cointegrated VARIMA models: specification and simulation José L. Gallego and Carlos Díaz Universidad de Cantabria. Abstract In this note we show how specify cointegrated vector autoregressive-moving average
More informationForecasting using R. Rob J Hyndman. 2.3 Stationarity and differencing. Forecasting using R 1
Forecasting using R Rob J Hyndman 2.3 Stationarity and differencing Forecasting using R 1 Outline 1 Stationarity 2 Differencing 3 Unit root tests 4 Lab session 10 5 Backshift notation Forecasting using
More informationEconometrics of Panel Data
Econometrics of Panel Data Jakub Mućk Meeting # 9 Jakub Mućk Econometrics of Panel Data Meeting # 9 1 / 22 Outline 1 Time series analysis Stationarity Unit Root Tests for Nonstationarity 2 Panel Unit Root
More informationSTA 6857 ARIMA and SARIMA Models ( 3.8 and 3.9)
STA 6857 ARIMA and SARIMA Models ( 3.8 and 3.9) Outline 1 Building ARIMA Models 2 SARIMA 3 Homework 4c Arthur Berg STA 6857 ARIMA and SARIMA Models ( 3.8 and 3.9) 2/ 34 Outline 1 Building ARIMA Models
More informationLecture Notes of Bus (Spring 2017) Analysis of Financial Time Series Ruey S. Tsay
Lecture Notes of Bus 41202 (Spring 2017) Analysis of Financial Time Series Ruey S. Tsay Simple AR models: (Regression with lagged variables.) Motivating example: The growth rate of U.S. quarterly real
More informationSTA 6857 ARIMA and SARIMA Models ( 3.8 and 3.9) Outline. Return Rate. US Gross National Product
STA 6857 ARIMA and SARIMA Models ( 3.8 and 3.9) Outline 1 Building ARIMA Models 2 SARIMA 3 Homework 4c Arthur Berg STA 6857 ARIMA and SARIMA Models ( 3.8 and 3.9) 2/ 34 Return Rate Suppose x t is the value
More 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 informationRomanian Economic and Business Review Vol. 3, No. 3 THE EVOLUTION OF SNP PETROM STOCK LIST - STUDY THROUGH AUTOREGRESSIVE MODELS
THE EVOLUTION OF SNP PETROM STOCK LIST - STUDY THROUGH AUTOREGRESSIVE MODELS Marian Zaharia, Ioana Zaheu, and Elena Roxana Stan Abstract Stock exchange market is one of the most dynamic and unpredictable
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 informationThe data was collected from the website and then converted to a time-series indexed from 1 to 86.
Introduction For our group project, we analyzed the S&P 500s futures from 30 November, 2015 till 28 March, 2016. The S&P 500 futures give a reasonable estimate of the changes in the market in the short
More informationChapter 2: Unit Roots
Chapter 2: Unit Roots 1 Contents: Lehrstuhl für Department Empirische of Wirtschaftsforschung Empirical Research and undeconometrics II. Unit Roots... 3 II.1 Integration Level... 3 II.2 Nonstationarity
More informationApplied time-series analysis
Robert M. Kunst robert.kunst@univie.ac.at University of Vienna and Institute for Advanced Studies Vienna October 18, 2011 Outline Introduction and overview Econometric Time-Series Analysis In principle,
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 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 informationOil price volatility in the Philippines using generalized autoregressive conditional heteroscedasticity
Oil price volatility in the Philippines using generalized autoregressive conditional heteroscedasticity Carl Ceasar F. Talungon University of Southern Mindanao, Cotabato Province, Philippines Email: carlceasar04@gmail.com
More informationPackage atsa. July 8, 2015
Type Package Title Alternative Time Series Analysis Version 3.1.2 Date 2015-06-19 Author Maintainer Package atsa July 8, 2015 Contains some tools for testing, analyzing time series data
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 informationThe Identification of ARIMA Models
APPENDIX 4 The Identification of ARIMA Models As we have established in a previous lecture, there is a one-to-one correspondence between the parameters of an ARMA(p, q) model, including the variance of
More information10) Time series econometrics
30C00200 Econometrics 10) Time series econometrics Timo Kuosmanen Professor, Ph.D. 1 Topics today Static vs. dynamic time series model Suprious regression Stationary and nonstationary time series Unit
More informationSOME BASICS OF TIME-SERIES ANALYSIS
SOME BASICS OF TIME-SERIES ANALYSIS John E. Floyd University of Toronto December 8, 26 An excellent place to learn about time series analysis is from Walter Enders textbook. For a basic understanding of
More informationReport for Fitting a Time Series Model
Report for Fitting a Time Series Model 1. Introduction In this report, I choose the stock price of Starbucks Corporation (SBUX) as my data to analyze. The time period for the data is from Jan 2009 to Dec
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 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 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 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 information1 Quantitative Techniques in Practice
1 Quantitative Techniques in Practice 1.1 Lecture 2: Stationarity, spurious regression, etc. 1.1.1 Overview In the rst part we shall look at some issues in time series economics. In the second part we
More informationConsider the trend-cycle decomposition of a time series y t
1 Unit Root Tests Consider the trend-cycle decomposition of a time series y t y t = TD t + TS t + C t = TD t + Z t The basic issue in unit root testing is to determine if TS t = 0. Two classes of tests,
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 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 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 informationE 4160 Autumn term Lecture 9: Deterministic trends vs integrated series; Spurious regression; Dickey-Fuller distribution and test
E 4160 Autumn term 2016. Lecture 9: Deterministic trends vs integrated series; Spurious regression; Dickey-Fuller distribution and test Ragnar Nymoen Department of Economics, University of Oslo 24 October
More informationFORECASTING ENERGY USAGE IN THE INDUSTRIAL SECTOR IN SWEDEN USING SARIMA AND DYNAMIC REGRESSION
FORECASTING ENERGY USAGE IN THE INDUSTRIAL SECTOR IN SWEDEN USING SARIMA AND DYNAMIC REGRESSION Submitted by Carl Anners A thesis submitted to the Department of Statistics in partial fulfillment of the
More informationTime Series I Time Domain Methods
Astrostatistics Summer School Penn State University University Park, PA 16802 May 21, 2007 Overview Filtering and the Likelihood Function Time series is the study of data consisting of a sequence of DEPENDENT
More 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 informationAutoregressive Integrated Moving Average Model to Predict Graduate Unemployment in Indonesia
DOI 10.1515/ptse-2017-0005 PTSE 12 (1): 43-50 Autoregressive Integrated Moving Average Model to Predict Graduate Unemployment in Indonesia Umi MAHMUDAH u_mudah@yahoo.com (State Islamic University of Pekalongan,
More informationAutomatic seasonal auto regressive moving average models and unit root test detection
ISSN 1750-9653, England, UK International Journal of Management Science and Engineering Management Vol. 3 (2008) No. 4, pp. 266-274 Automatic seasonal auto regressive moving average models and unit root
More informationIS THE NORTH ATLANTIC OSCILLATION A RANDOM WALK? A COMMENT WITH FURTHER RESULTS
INTERNATIONAL JOURNAL OF CLIMATOLOGY Int. J. Climatol. 24: 377 383 (24) Published online 11 February 24 in Wiley InterScience (www.interscience.wiley.com). DOI: 1.12/joc.13 IS THE NORTH ATLANTIC OSCILLATION
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 informationForecasting Gold Price. A Comparative Study
Course of Financial Econometrics FIRM Forecasting Gold Price A Comparative Study Alessio Azzutti, University of Florence Abstract: This paper seeks to evaluate the appropriateness of a variety of existing
More 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 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 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 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 informationAuthor: Yesuf M. Awel 1c. Affiliation: 1 PhD, Economist-Consultant; P.O Box , Addis Ababa, Ethiopia. c.
ISSN: 2415-0304 (Print) ISSN: 2522-2465 (Online) Indexing/Abstracting Forecasting GDP Growth: Application of Autoregressive Integrated Moving Average Model Author: Yesuf M. Awel 1c Affiliation: 1 PhD,
More informationThe ARIMA Procedure: The ARIMA Procedure
Page 1 of 120 Overview: ARIMA Procedure Getting Started: ARIMA Procedure The Three Stages of ARIMA Modeling Identification Stage Estimation and Diagnostic Checking Stage Forecasting Stage Using ARIMA Procedure
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 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 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 informationPart 1. Multiple Choice (40 questions, 1 point each) Part 2. Problems/Short Answer (10 questions, 6 points each)
GROUND RULES: This exam contains two parts: Part 1. Multiple Choice (40 questions, 1 point each) Part 2. Problems/Short Answer (10 questions, 6 points each) The maximum number of points on this exam is
More 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 informationARIMA Modelling and Forecasting
ARIMA Modelling and Forecasting Economic time series often appear nonstationary, because of trends, seasonal patterns, cycles, etc. However, the differences may appear stationary. Δx t x t x t 1 (first
More informationResponse surface models for the Elliott, Rothenberg, Stock DF-GLS unit-root test
Response surface models for the Elliott, Rothenberg, Stock DF-GLS unit-root test Christopher F Baum Jesús Otero Stata Conference, Baltimore, July 2017 Baum, Otero (BC, U. del Rosario) DF-GLS response surfaces
More information