STAT 497LECTURE NOTES 4
MODEL INDETIFICATION AND NON-STATIONARY TIME SERIES MODELS
1
MODEL IDENTIFICATION
• We have learned a large class of linear parametric models for stationary time series processes.
• Now, the question is how we can find out the best suitable model for a given observed series. How to choose the appropriate model (on order of p and q).
2
MODEL IDENTIFICATION
• ACF and PACF show specific properties for specific models. Hence, we can use them as a criteria to identify the suitable model.
• Using the patterns of sample ACF and sample PACF, we can identify the model.
3
MODEL SELECTION THROUGH CRITERIA
• Besides sACF and sPACF plots, we have also other tools for model identification.
• With messy real data, sACF and sPACF plots become complicated and harder to interpret.
• Don’t forget to choose the best model with as few parameters as possible.
• It will be seen that many different models can fit to the same data so that we should choose the most appropriate (with less parameters) one and the information criteria will help us to decide this.
4
MODEL SELECTION THROUGH CRITERIA
• The three well-known information criteria are– Akaike’s information criterion (AIC) (Akaike, 1974)– Schwarz’s Bayesian Criterion (SBC) (Schwarz, 1978).
Also known as Bayesian Information Criterion (BIC)– Hannan-Quinn Criteria (HQIC) (Hannan&Quinn,
1979)
5
AIC• Assume that a statistical model of M parameters
is fitted to data
• For the ARMA model and n observations, the log-likelihood function
.2likelihood maximumln2 MAIC
sidualSS
qpa
a Sn
L
Re
22 ,,
21
2ln2
ln
.,0~ assuming 2..
a
dii
t Na 6
AIC
• Then, the maximized log-likelihood is
constant
2 2ln12
ˆln2
nl nnL a
MnAIC a 2ˆln 2
Choose model (or the value of M) with minimum AIC.
7
SBC• The Bayesian information criterion (BIC) or
Schwarz Criterion (also SBC, SBIC) is a criterion for model selection among a class of parametric models with different numbers of parameters.
• When estimating model parameters using maximum likelihood estimation, it is possible to increase the likelihood by adding additional parameters, which may result in overfitting. The BIC resolves this problem by introducing a penalty term for the number of parameters in the model.
8
SBC
• In SBC, the penalty for additional parameters is stronger than that of the AIC.
• It has the most superior large sample properties.
• It is consistent, unbiased and sufficient.
nMnSBC a lnˆln 2
9
HQIC
• The Hannan-Quinn information criterion (HQIC) is an alternative to AIC and SBC.
• It can be shown [see Hannan (1980)] that in the case of common roots in the AR and MA polynomials, the Hannan-Quinn and Schwarz criteria still select the correct orders p and q consistently.
nMnHQIC a lnln2ˆln 2
10
THE INVERSE AUTOCORRELATION FUNCTION
• The sample inverse autocorrelation function (SIACF) plays much the same role in ARIMA modeling as the sample partial autocorrelation function (SPACF), but it generally indicates subset and seasonal autoregressive models better than the SPACF.
11
THE INVERSE AUTOCORRELATION FUNCTION
• Additionally, the SIACF can be useful for detecting over-differencing. If the data come from a nonstationary or nearly nonstationary model, the SIACF has the characteristics of a noninvertible moving-average. Likewise, if the data come from a model with a noninvertible moving average, then the SIACF has nonstationary characteristics and therefore decays slowly. In particular, if the data have been over-differenced, the SIACF looks like a SACF from a nonstationary process
12
THE INVERSE AUTOCORRELATION FUNCTION
• Let Yt be generated by the ARMA(p, q) process
• If (B) is invertible, then the model
is also a valid ARMA(q, p) model. This model is sometimes referred to as the dual model. The autocorrelation function (ACF) of this dual model is called the inverse autocorrelation function (IACF) of the original model.
.,0~ where 2attqtp WNaaBYB
tptq aBZB
13
THE INVERSE AUTOCORRELATION FUNCTION
• Notice that if the original model is a pure autoregressive model, then the IACF is an ACF that corresponds to a pure moving-average model. Thus, it cuts off sharply when the lag is greater than p; this behavior is similar to the behavior of the partial autocorrelation function (PACF).
• Under certain conditions, the sampling distribution of the SIACF can be approximated by the sampling distribution of the SACF of the dual model (Bhansali 1980). In the plots generated by ARIMA, the confidence limit marks (.) are located at 2n1/2. These limits bound an approximate 95% confidence interval for the hypothesis that the data are from a white noise process.
14
THE EXTENDED SAMPLE AUTOCORRELATION FUNCTION_ESACF• The extended sample autocorrelation function
(ESACF) method can tentatively identify the orders of a stationary or nonstationary ARMA process based on iterated least squares estimates of the autoregressive parameters. Tsay and Tiao (1984) proposed the technique.
15
ESACF
• Consider ARMA(p, q) model
or
then
follows an MA(q) model
16
tq
qtp
p aBBYBB 101 11
qtqttptptt aaaYYY 11110
tp
pt YBBZ 11
.110 qtqttt aaaZ
ESACF• Given a stationary or nonstationary time series
Yt with mean corrected form with a true autoregressive order of p+d and with a true moving-average order of q, we can use the ESACF method to estimate the unknown orders and by analyzing the sample autocorrelation functions associated with filtered series of the form
17
tt YY
m
ijt
jmittjm
jmt YYYBZ
1
,,
, ˆˆ
process. jm,ARMAan is series that the
assumption under the estimatesparameter theare s'ˆ where i
ESACF
• It is known that OLS estimators for ARMA process are not consistent so that an iterative procedure is proposed to overcome this.
• The j-th lag of the sample autocorrelation function of the filtered series is the extended sample autocorrelation function, and it is denoted as
18
1,
1,111,
11,1,
ˆˆ
ˆˆˆ
jm
m
jmmjm
ijm
ijm
i
.ˆ mj
ESACF
ESACF TABLE
19
ARMA
0 1 2 3 …
0 …
1 …
2 …
3 …
… … … … … …
01
02
03
04
34
33
32
22
12
31
11 21
13
14 24
23
ESACF
• For an ARMA(p,q) process, we have the following convergence in probability, that is, for m=1,2,… and j=1,2,…, we have
20
otherwise 0,X
0 ,0ˆ
qjpmmj
ESACF
• Thus, the asymptotic ESACF table for ARMA(1,1) model becomes
21
ARMA
0 1 2 3 4 …
0 X X X X X …
1 X 0 0 0 0 …
2 X X 0 0 0 …
3 X X X 0 0 …
4 X X X X 0 …
… … … … … … …
ESACF• In practice, we have finite samples and may not be exactly zero.
However, we can use the Bartlett’s approximate formula for the asymptotic variance of .
• The orders are tentatively identified by finding a right (maximal) triangular pattern with vertices located at (p+d, q) and (p+d, qmax) and in which all elements are insignificant (based on asymptotic normality of the autocorrelation function). The vertex (p+d, q) identifies the order.
22
qjpmmj 0,
mj
MINIMUM INFORMATION CRITERION
MINIC TABLE
23
MA AR 0 1 2 3 …0 SBC(0,0) SBC(0,1) SBC(0,2) SBC(0,3) …1 SBC(1,0) SBC(1,1) SBC(1,2) SBC(1,3) …2 SBC(2,0) SBC(2,1) SBC(2,2) SBC(2,3) …3 SBC(3,0) SBC(3,1) SBC(3,2) SBC(3,3) …… … … … … …
NON-STATIONARY TIME SERIES MODELS
• Non-constant in mean
• Non-constant in variance
• Both
24
NON-STATIONARY TIME SERIES MODELS
• Inspection of the ACF serves as a rough indicator of whether a trend is present in a series. A slow decay in ACF is indicative of a large characteristic root; a true unit root process, or a trend stationary process.
• Formal tests can help to determine whether a system contains a trend and whether the trend is deterministic or stochastic.
25
NON-STATIONARITY IN MEAN
• Deterministic trend– Detrending
• Stochastic trend– Differencing
26
DETERMINISTIC TREND• A deterministic trend is when we say that the
series is trending because it is an explicit function of time.
• Using a simple linear trend model, the deterministic (global) trend can be estimated. This way to proceed is very simple and assumes the pattern represented by linear trend remains fixed over the observed time span of the series. A simple linear trend model:
27
tt atY
DETERMINISTIC TREND
• The parameter measure the average change in Yt from one period to the another:
• The sequence {Yt} will exhibit only temporary departures from the trend line +t. This type of model is called a trend stationary (TS) model.
28
t
ttttt
YE
aattYYY 11 1
TREND STATIONARY
• If a series has a deterministic time trend, then we simply regress Yt on an intercept and a time trend (t=1,2,…,n) and save the residuals. The residuals are detrended series. If Yt is stochastic, we do not necessarily get stationary series.
29
DETERMINISTIC TREND• Many economic series exhibit “exponential
trend/growth”. They grow over time like an exponential function over time instead of a linear function.
• For such series, we want to work with the log of the series:
30
t
tt
YE
atY
ln
: is rategrowth average theSo
ln
DETERMINISTIC TREND
• Standard regression model can be used to describe the phenomenon. If the deterministic trend can be described by a k-th order polynomial of time, the model of the process
31
.,0~ where 2
2210
at
tk
kt
WNa
atttY
DETERMINISTIC TREND
• This model has a short memory. • If a shock hits a series, it goes back to
trend level in short time. Hence, the best forecasts are not affected.
• Rarely model like this is useful in practice. A more realistic model involves stochastic (local) trend.
32
STOCHASTIC TREND
• A more modern approach is to consider trends in time series as a variable. A variable trend exists when a trend changes in an unpredictable way. Therefore, it is considered as stochastic.
33
STOCHASTIC TREND
• Recall the AR(1) model: Yt = c + Yt−1 + at.
• As long as || < 1, everything is fine (OLS is consistent, t-stats are asymptotically normal, ...).
• Now consider the extreme case where = 1, i.e.Yt = c + Yt−1 + at.
• Where is the trend? No t term.
34
STOCHASTIC TREND• Let us replace recursively the lag of Yt on the
right-hand side:
35
t
ii
ttt
ttt
aYtc
aaYcc
aYcY
10
12
1
Deterministic trend
• This is what we call a “random walk with drift”. If c = 0, it is a“random walk”.
STOCHASTIC TREND
• Each ai shock represents shift in the intercept. Since all values of {ai} have a coefficient of unity, the effect of each shock on the intercept term is permanent.
• In the time series literature, such a sequence is said to have a stochastic trend since each ai shock imparts a permanent and random change in the conditional mean of the series. To be able to define this situation, we use Autoregressive Integrated Moving Average (ARIMA) models.
36
DETERMINISTIC VS STOCHASTIC TREND• They might appear similar since they both lead to
growth over time but they are quite different.• To see why, suppose that through any policies,
you got a bigger Yt because the noise at is big. What will happen next period?
– With a deterministic trend, Yt+1 = c +(t+1)+at+1. The noise at is not affecting Yt+1. Your stupendous policy had a one period impact.
– With a stochastic trend, Yt+1 = c + Yt + at+1 = c + (c + Yt−1 + at) + at+1. The noise at is affecting Yt+1. In fact, the policy will have a permanent impact.
37
DETERMINISTIC VS STOCHASTIC TREND
Conclusions:– When dealing with trending series, we are always interested
in knowing whether the growth is a deterministic or stochastic trend.
– There are also economic time series that do not grow over time (e.g., interest rates) but we will need to check if they have a behavior ”similar” to stochastic trends ( = 1 instead of || < a, while c = 0).
– A deterministic trend refers to the long-term trend that is not affected by short term fluctuations in the economy. Some of the occurrences in the economy are random and may have a permanent effect of the trend. Therefore the trend must contain a deterministic and a stochastic component.
38
AUTOREGRESSIVE INTEGRATED MOVING AVERAGE (ARIMA) PROCESSES• Consider an ARIMA(p,d,q) process
39
.0,WN~ and
rootscommon no share 1
and 1 where
1
2a
1
1
0
t
qqq
ppp
tqtd
p
a
BBB
BBB
aBYBB
ARIMA MODELS
• When d=0, 0 is related to the mean of the process.
• When d>0, 0 is a deterministic trend term.– Non-stationary in mean:
– Non-stationary in level and slope:
40
.1 10 p
tqtp aBYBB 01
tqtp aBYBB 021
RANDOM WALK PROCESS
• A random walk is defined as a process where the current value of a variable is composed of the past value plus an error term defined as a white noise (a normal variable with zero mean and variance one).
• ARIMA(0,1,0) PROCESS
41
.,0~ where
12
1
at
tttttt
WNa
aYBYaYY
RANDOM WALK PROCESS
• Behavior of stock market.• Brownian motion.• Movement of a drunken men.• It is a limiting process of AR(1).
42
RANDOM WALK PROCESS• The implication of a process of this type is that the
best prediction of Y for next period is the current value, or in other words the process does not allow to predict the change (YtYt-1). That is, the change of Y is absolutely random.
• It can be shown that the mean of a random walk process is constant but its variance is not. Therefore a random walk process is nonstationary, and its variance increases with t.
• In practice, the presence of a random walk process makes the forecast process very simple since all the future values of Yt+s for s > 0, is simply Yt.
43
RANDOM WALK PROCESS
44
RANDOM WALK PROCESS
45
RANDOM WALK WITH DRIFT
• Change in Yt is partially deterministic and partially stochastic.
• It can also be written as
46
t
Y
tt aYY
t
01
trend
stochastic
1
trendticdeterminis
00
t
iit atYY Pure model of a trend
(no stationary component)
RANDOM WALK WITH DRIFT
47
00 tYYE t
After t periods, the cumulative change in Yt is t0.
flatnot 0 sYYYE ttst
Each ai shock has a permanent effect on the mean of Yt.
RANDOM WALK WITH DRIFT
48
ARIMA(0,1,1) OR IMA(1,1) PROCESS
• Consider a process
• Letting
49
.,0~ where
112at
tt
WNa
aBYB
tt YBW 1
stationaryaBW tt 1
ARIMA(0,1,1) OR IMA(1,1) PROCESS
• Characterized by the sample ACF of the original series failing to die out and by the sample ACF of the first differenced series shows the pattern of MA(1).
• IF:
50
.1 where11
t
jjt
jt aYY
1
121 1,,
jjt
jttt YYYYE
Exponentially decreasing. Weighted MA of its past values.
ARIMA(0,1,1) OR IMA(1,1) PROCESS
51
,,1,, 2111 ttttttt YYYEYYYYE
where is the smoothing constant in the method of exponential smoothing.
REMOVING THE TREND
• Shocks to a stationary time series are temporary over time. The series revert to its long-run mean.
• A series containing a trend will not revert to a long-run mean. The usual methods for eliminating the trend are detrending and differencing.
52
DETRENDING
• Detrending is used to remove deterministic
trend.
• Regress Yt on time and save the residuals.
• Then, check whether residuals are stationary.
53
DIFFERENCING
• Differencing is used for removing the stochastic trend.
• d-th difference of ARIMA(p,d,q) model is stationary. A series containing unit roots can be made stationary by differencing.
• ARIMA(p,d,q) d unit roots
54
Integrated of order d, I(d)
dIYt ~
DIFFERENCING
• Random Walk:
55
ttt aYY 1
tt aY
Non-stationary
Stationary
DIFFERENCING
• Differencing always makes us to loose observation.
• 1st regular difference: d=1
• 2nd regular difference: d=2
56
tttt YYYYB 11
21222 2211 tttttt YYYYBBYYB
difference 2nd not the is 2 tt YY
DIFFERENCING
57
Yt Yt 2Yt YtYt-23 * * *
8 83=5 * *
5 58=3 35=8 53=2
9 95=4 4(3)=7 98=1
KPSS TEST
• To be able to test whether we have a deterministic trend vs stochastic trend, we are using KPSS (Kwiatkowski, Phillips, Schmidt and Shin) Test (1992).
58
stationary difference 1~:
stationary or trend level0~:
1
0
IYH
IYH
t
t
KPSS TESTSTEP 1: Regress Yt on a constant and trend and
construct the OLS residuals e=(e1,e2,…,en)’.STEP 2: Obtain the partial sum of the residuals.
STEP 3: Obtain the test statistic
where is the estimate of the long-run variance of the residuals.
59
t
iit eS
1
n
t
tSnKPSS
12
2
2
KPSS TEST
• STEP 4: Reject H0 when KPSS is large, because that is the evidence that the series wander from its mean.
• Asymptotic distribution of the test statistic uses the standard Brownian bridge.
• It is the most powerful unit root test but if there is a volatility shift it cannot catch this type non-stationarity.
60
PROBLEM
• When an inappropriate method is used to eliminate the trend, we may create other problems like non-invertibility.
• E.g.
61
. and circleunit
theoutside are 0 of roots thewhere
stationary Trend10
tt
tt
aB
B
tYB
PROBLEM
• But if we misjudge the series as difference stationary, we need to take a difference. Actually, detrending should be applied. Then, the first difference:
62
tt BYB 11
Now, we create a non-invertible unit root process in the MA component.
PROBLEM
• To overcome this, look at the inverse sample autocorrelation function. If it has the same ACF pattern of non-stationary process (that is, slow decaying behavior), this means that we over-differenced the series.
• Go back and detrend the series instead of differencing.
• There are also smoothing filters to eliminate the trend (Decomposition Methods).
63
NON-STATIONARITY IN VARIANCE
• Stationarity in mean Stationarity in variance• Non-stationarity in mean Non-stationarity in
variance• If the mean function is time dependent,1. The variance, Var(Yt) is time dependent.2. Var(Yt) is unbounded as t.3. Autocovariance and autocorrelation functions are
also time dependent.4. If t is large wrt Y0, then k 1.
64
VARIANCE STABILIZING TRANSFORMATION
• The variance of a non-stationary process changes as its level changes
for some positive constant c and a function f.• Find a function T so that the transformed
series T(Yt) has a constant variance.
65
tt fcYVar .
The Delta Method
VARIANCE STABILIZING TRANSFORMATION
• Generally, we use the power function
66
1964) Cox, and(Box 1
tt
YYT
Transformation1 1/Yt
0.5 1/(Yt)0.5
0 ln Yt
0.5 (Yt)0.5
1 Yt (no transformation)
VARIANCE STABILIZING TRANSFORMATION
• Variance stabilizing transformation is only for positive series. If your series has negative values, then you need to add each value with a positive number so that all the values in the series are positive. Now, you can search for any need for transformation.
• It should be performed before any other analysis such as differencing.
• Not only stabilize the variance but also improves the approximation of the distribution by Normal distribution.
67