Cointegration, Error Correction CH 3 - 8 PDF
Cointegration, Error Correction CH 3 - 8 PDF
Cointegration, Error Correction CH 3 - 8 PDF
General Editors
C. W . J . GRANGE R G
. E . MIZO N
CO-INTEGRATION, ERROR
CORRECTION, AND
THE ECONOMETRI C
ANALYSIS O F
NON-STATIONARY DAT A
Anindya Banerjee, Juan J. Dolado,
John W. "Galbraith, and Davi d F . Hendry
OXFORD
UNIVERSITY PRES S
Preface
This boo k i s intended a s a guid e t o th e literatur e o n co-integratio n an d
modelling o f integrate d processes . Time-serie s econometric s ha s devel oped rapidl y durin g th e pas t decade , bu t especiall y s o in th e analysi s of
non-stationarity. I n particular , th e stud y o f integrate d processe s ha s
grown i n importance fro m th e statu s of a n exoti c topic, discusse d onl y in
technical journals , t o bein g a n essentia l par t o f th e econometrician' s
collection o f techniques . I t ha s thereb y develope d int o a n are a o f
interest fo r econometri c theorist s an d applie d econometrician s alike .
This boo k i s aime d a t graduat e student s i n economics , applie d econo metricians, econometri c theorists , an d th e genera l audienc e o f econo mists who use empirica l methods t o analys e tim e series.
Despite th e growin g importanc e o f th e literatur e o n integratio n an d
co-integration, mos t account s o f thi s literatur e remai n confine d t o
journals, edite d collection s o f papers , o r surve y papers. Whil e som e o f
the survey s ar e quit e detailed , spac e restriction s usuall y d o no t allo w a
full expositio n o f man y o f th e theoretica l points . Thi s boo k attempt s t o
bridge th e ga p betwee n account s suc h a s surveys , whic h ar e mainl y
descriptive, an d account s tha t ar e mainl y theoretical . I t explain s th e
important concept s informall y an d als o present s the m formally . Th e
asymptotic theor y o f integrate d processe s i s describe d an d th e tool s
provided b y thi s theor y ar e use d t o derive , i n som e detail , th e
distributions o f estimators. B y taking reader s ste p b y ste p throug h som e
of th e mai n derivations , ou r hop e i s t o mak e th e theor y readil y
accessible t o a wide audience .
We hav e trie d t o mak e th e boo k a s self-containe d a s possible . A
knowledge o f econometrics , statistics , an d matri x algebr a a t th e leve l of
a final-yea r undergraduat e o r first-yea r graduat e cours e i n econometric s
is assumed , bu t otherwis e al l o f th e importan t statistica l concept s an d
techniques ar e described .
A boo k suc h a s thi s one , whic h discusse s a n are a tha t i s developin g
rapidly, i s inevitabl y incomplet e an d run s th e ris k o f no t bein g quit e
up-to-date. T o limi t th e tim e take n i n writin g an d revising , w e di d no t
seek t o chas e a frontie r tha t wa s expanding in man y directions . Rather ,
the topic s covere d reflec t ou r view s of issues, models , an d method s tha t
are likel y t o remai n importan t fo r som e tim e t o come , man y o f whic h
will continue to provid e th e platfor m for futur e research .
Acknowledgements
Our boo k wa s writte n i n tw o continents , thre e years , an d fou r univer sities, s o th e lis t o f people , acros s time , space , an d departments , t o
whom w e ow e extensiv e debt s o f gratitud e ha s grow n formidably large.
A majo r par t o f thi s deb t i s owe d t o th e Department s o f Economic s a t
the Universitie s o f Californi a a t Sa n Diego , Florid a i n Gainesville ,
McGill, an d Oxford , an d th e Ban k o f Spain , wher e th e author s eithe r
worked o r visite d for substantia l periods. Thei r generou s suppor t o f ou r
work i s much appreciated .
The boo k ha s als o benefite d greatl y fro m th e patien t scrutin y o f
several o f ou r colleagues , wh o rea d th e entir e typescript an d mad e
detailed comments . W e hav e pleasur e i n thankin g Michae l Clements ,
Rob Engle , Neil Ericsson, Ton y Hall (an d severa l o f his students), Colin
Hargreaves, S0re n Johansen , Katarin a Juselius , Teu n Kloek , Jame s
MacKinnon, G . S . Maddala , Grayha m Mizon , Jean-Fran9oi s Richard ,
Mark Rush , Nei l Shephard , Tim o Terasvirta , an d fou r anonymou s
referees fo r thei r help . The y hav e mad e a grea t contributio n t o thi s
book, an d foun d man y infelicitie s i n earlie r versions , bu t o f cours e ar e
not responsibl e for an y that remain.
Early version s o f th e boo k wer e inflicte d b y u s upo n ou r graduat e
students. Amon g thos e wh o suffere d fro m th e confusio n cause d
by obscur e notatio n an d prose , bu t continue d unflinchingly , Hughe s
Dauphin, Caro l Dole , Jesu s Gonzalo , Catherin e Liston , Claudi o Lupi ,
Neil Rickman , an d Geet a Sing h deserve specia l thanks.
We ar e als o indebte d t o Juli a Campos , Michae l Clements , Steve n
Cook, Nei l Ericsson an d Claudi o Lup i fo r proof reading.
The financia l suppor t o f th e Economi c an d Socia l Researc h Counci l
(UK) unde r grant s B0125002 4 an d R23118 4 an d th e Fond s pou r l a
Formation de s Chercheur s e t 1'Aid e a l a Recherch e (Quebec ) i s grate fully acknowledged . Finally, w e than k Andre w Schulle r an d th e editor s
of thi s series , wh o remaine d encouragin g abou t th e projec t despit e it s
many difficulties .
Oxford A
Madrid J
Montreal J
Oxford D
.B.
. J. D .
. W. G.
. F. H.
Contents
Notational Conventions, Symbols , an d Abbreviations x
1. Introductio n and Overview 1
1.1. Equilibrium relationships and the long run 2
1.2. Stationarity and equilibrium relationships 4
1.3. Equilibrium and the specification of dynamic models 5
1.4. Estimation of long-run relationships and testing for orders
of integration and co-integration 8
1.5. Preliminary concepts an d definitions 1
1.6. Data representation an d transformations 2
1.7. Examples: typical ARM A processes 3
1.8. Empirical time series: money, prices, output, and interest
rates 4
1.9. Outline o f later chapters 4
Appendix 4
Linear Transformations , Erro r Correction , and the Lon g
Run i n Dynami c Regressio n 4
2.1. Transformations o f a simple model 4
2.2. Th e error-correction model 5
2.3. A n example 5
2.4. Bdrdsen an d Bewley transformations 5
2.5. Equivalence o f estimates from different transformations 5
2.6. Homogeneity and the ECM as a linear transformation
oftheADL 6
2.7. Variances o f estimates o f long-run multipliers 6
2.8. Expectational variables and the interpretation of
long-run solutions 6
0
8
2
0
2
3
6
8
0
2
3
5
0
1
4
9
0
1
4
6
1
5
viii Content
9
4
6
8
3
9
3
5. Co-integratio n 13
5.1. A n example 13
5.2. Polynomial matrices 14
5.3. Integration and co-integration: formal definitions and
theorems 14
5.4. Significance o f alternative representations 15
5.5. Alternative representations o f co-integrated variables:
two examples 15
5.6. Engle- Granger two-step procedure 15
6
7
0
3
7
2
4
8
2
0
4
5
6
1
4
0
6
9
0
2
4
2
5
7
1
6
8
1
5
3
Contents i
9. Conclusio n 29
9.1. Summary 29
9.2. Th e invariance o f co-integrating vectors 30
9.3. Invariance o f co-integration under seasonal adjustment 30
9.4. Structured time-series models an d co-integration 30
9.5. Recent research on integration and co-integration 30
9.6. Reinterpreting econometrics time-series problems 30
9
9
0
1
3
4
7
References 31
Author Index 32
Subject Index 32
Notational Conventions,
Symbols, and Abbreviations
The following notationa l convention s will be used throughou t the text:
Y, y endogenou
X, Z , x , z exogenou
s variables
s variables, o r vectors
containing both y an d z
Greek letters populatio
n values (parameters)
Greek letters with ~ o r ~ sampl
e values (estimates )
Bold lowe r case (Roma n o r Greek) vector s
Bold upper cas e (Roman or Greek ) matrice s
Equation numbers
Equations ar e numbere d consecutivel y i n eac h chapte r an d referre d t o
within tha t chapte r b y this number alone . Equation s fro m othe r chapter s
are referre d t o b y th e chapte r numbe r an d equatio n numbe r withi n
chapter; e.g . th e fift h equatio n i n Chapte r 2 is (5) within Chapter 2 , an d
(2.5) elsewhere .
Symbols
la
first-differenc
Kronecke
fo
g operator:
e operator :
r produc t
r al l
modulus or absolut e value of x, where x i s a scalar
determinan
t o f A, wher e A is a matrix
x
conditiona l on y
wea
k convergence
convergenc
e i n distribution
convergenc
e i n probability
Abbreviations
ADF augmente
d Dickey-Fuller
ADL autoregressive-distribute
d lag
xii Notationa
AR autoregressio
n
ARIMA autoregressiv
e integrate d movin g average
ARMA autoregressive-movin
g averag e
ARMAX ARM
A + additiona l exogenou s processe s
ASE Asymptoti
c standard erro r
BM Brownia
n motio n
Cl(d, b) co-integrate
d o f order d , b
CLT centra
l limi t theore m
COMFAC commo
n facto r error representatio n
CRDW co-integratin
g regression D W statistic
diag diagona
l matrix
d.f. degree
s o f freedom
DF Dickey-Fulle
r
DGP data-generatio
n proces s
DW Durbin-Watso
n statisti c
ECM error-correctio
n model/mechanis m
ESE (average
) estimate d standar d erro r
FCLT functiona
l centra l limi t theorem/ s
FIML full-informatio
n maximu m likelihood
GLS generalize
d least square s
GNP gros
s national produc t
\(d) integrate
d of orde r d
ID independentl
y distribute d
IID independentl
y an d identically distributed
IMA integrate
d movin g average
IN(/i, a 2 ) independentl
y and normall y distributed with mean fi an d
variance a 2
IV instrumenta
l variables
LIML limited-informatio
n maximum likelihood
MA movin
g averag e
MDS martingal
e difference sequence
MLE maximu
m likelihood estimato r
N(ju, a 2 ) normall
y distribute d wit h mean p, and variance a 2
NI near-integrate
d
OLS ordinar
y least square s
SC Schwar
z information criterion
SD standar
d deviatio n
SE standar
d erro r
SI seasonall
y integrated
SSD sampl
e standar d deviatio n
T sampl
e siz e or las t observatio n i n a time-series
TFE tota
l fina l expenditur e
VAR vecto
r autoregressio n
var varianc
e
g operator
r (Brownia n motion) process wit h increments of
variance r
Introduction an d Overvie w
This boo k consider s th e econometri c analysi s o f bot h stationar y
and non-stationar y processe s whic h ma y b e linke d b y equilibriu m
relationships. I t exposit s th e mai n tools , techniques , models , con cepts, an d distribution s involve d i n econometri c modellin g o f
possibly non-stationar y time-serie s data . Sinc e th e focu s i s o n
equilibrium concepts , includin g co-integration an d erro r correction ,
the analysi s begin s wit h a discussio n o f th e applicatio n o f thes e
concepts t o stationar y empirica l models . Late r w e wil l sho w tha t
integrated processe s ca n b e reduce d t o thi s cas e b y suitabl e
transformations tha t tak e advantag e o f co-integrating (equilibrium )
relationships. I n thi s chapte r w e wil l introduc e som e importan t
concepts fro m time-serie s analysi s an d th e theor y o f stochasti c
processes, an d i n particula r th e theor y o f Brownia n motio n pro cesses. W e als o offe r severa l empirica l example s whic h us e thes e
concepts.
A significan t re-evaluatio n o f th e statistica l basis o f econometri c model ling too k plac e durin g th e 1980s . It s analytica l basis expande d fro m th e
assumption o f stationarit y t o includ e integrate d processes . Th e effec t o f
this shif t i s fa r fro m complete , bu t i s alread y radical , influencin g th e
choice o f mode l forms , modellin g practices , statistica l inference , dis tribution theory , an d th e interpretatio n o f man y traditiona l concept s
such a s simultaneity , measurement errors , collinearity , forecasting , an d
exogeneity. Thi s boo k attempt s t o analys e thes e issues , describ e th e
tools necessar y t o investigat e integrate d processes , an d relat e th e ne w
methods t o thos e mor e familia r t o econometricians . Researc h i s con tinuing a t a rapi d pace , an d sinc e thi s boo k canno t cove r al l o f th e
techniques tha t hav e bee n explored , w e wil l concentrat e o n thos e tha t
we believe wil l remain useful .
Time-series econometric s i s concerned wit h th e estimatio n o f relation ships amon g group s of variables , eac h o f whic h is observed a t a numbe r
of consecutiv e point s i n time . Th e relationship s amon g thes e variable s
may b e complicated ; i n particular , th e valu e o f eac h variabl e ma y
depend o n th e value s take n b y man y other s i n severa l previou s tim e
periods. I n consequence , th e effec t tha t a chang e in one variabl e ha s on
another depend s upo n th e tim e horizo n tha t w e consider . I t i s eas y t o
2 Introductio
n an d Overvie w
4 Introductio
n and Overvie w
6 Introductio
n and Overvie w
8 Introductio
n and Overvie w
10 Introductio
n and Overvie w
12 Introductio
n and Overvie w
Stationarity implie s tha t y(/z ) = y ( h ) , sinc e th e autocovarianc e be tween tw o value s depend s onl y o n th e distanc e betwee n them . Th e
autocorrelation function i s defined similarly, as
with
14 Introductio
n and Overvie w
16 Introductio
n and Overvie w
as /z^>o ; tha t is , th e join t distributio n functio n o f th e tw o sub sequences o f {y,} approache s th e produc t o f th e distribution s o f eac h o f
1.5.4. Exogeneity
While ou r primar y focu s i s o n integrate d serie s an d th e problem s the y
imply fo r standar d econometri c analyses , rathe r tha n o n th e problem s
created b y a failur e o f exogeneit y (i n th e appropriat e sense) , i t wil l b e
important t o conside r exogeneit y a t several points .
Econometric analysi s often proceed s o n th e basi s o f a single-equatio n
model o f a proces s o f interest . Implicitly , w e assum e tha t knowledg e of
the processe s generatin g th e explanator y variable s woul d carr y n o
information relevan t t o th e parameter s o f interest . A s Engle , Hendry ,
18 Introductio
n and Overvie w
with
and E[y, z
tHE cORRESPONDING CONDITIONAL VARIANCE VAR
20 Introductio
n an d Overvie w
Introduction an d Overvie w 2
and
22
FlGl
23
24
Introduction an d Overvie w
25
26
Introduction an d Overview
The las t expressio n i s yi/VT, wher e y\ i s the lagge d mean . Thi s resul t
uses th e fac t that , fo r any constant c,
and hence
Thus, som e functiona l o f Wiene r processe s ar e familia r rando m vari ables i n disguis e and w e will develo p thi s aspect a s we proceed. A proo f
of (11 ) i s given in the Appendix .
7.5.7. Monte Carlo Simulation
The purpos e o f Mont e Carl o simulatio n i s t o evaluat e b y experimen t
quantities tha t woul d be ver y difficult o r impossibl e t o evaluat e analytically. Suc h experiment s typicall y begi n b y creatin g a se t o f dat a wit h
known statistica l properties . Thi s i s achieve d b y specifyin g ever y aspec t
of a data-generatin g process , o r clas s o f suc h processes , an d replacin g
the rando m error s o f th e DG P b y pseudo-rando m numbers . Pseudo random number s ar e number s generate d deterministicall y t o mimi c a
random proces s wit h a particula r distribution . A n investigato r typically
generates a large numbe r o f suc h artificial data set s (calle d replications )
to investigat e statistica l technique s whic h analys e thes e dat a a s i f th e
process generating them were no t known. Th e performanc e o f th e
statistical techniqu e i n revealin g som e characteristi c o f th e dat a se t ma y
4
28 Introductio
n and Overvie w
model wit h a constan t erro r variance . I n particular , a normall y dis tributed rando m wal k wit h drif t i s ofte n postulate d s o tha t
Axt ~ IN(jW , cr 2). Man y economi c tim e serie s (suc h a s consumption ,
national income an d expenditure , o r th e pric e level ) d o gro w over time ,
but th e amoun t b y whic h the y gro w i n eac h perio d als o tend s t o rise .
However, A.x t = x t xt-i wil l b e stationar y onl y if the absolut e amoun t
of growt h is stationary , i n whic h cas e fo r n > 0, a/x t wil l ten d t o zero .
Percentage growth , b y contrast , ofte n display s n o obviou s tendenc y t o
rise o r fall , makin g it a more likel y candidate fo r stationarity . Since th e
levels o f man y economi c variable s ar e initiall y positive , an d recallin g
that
30 Introductio
n and Overvie w
31
32 Introductio
n and Overvie w
= t [whit
e noise ]
Fig. 1.1 1 u,
= e, + 0.8e,_i [MA(1)
, stationary]
Fig. 1.1 2 u,
= e, - 0.8,_ ! [MA(1)
, stationary ]
Fig. 1.1 3 u,
= 0. 5 ,_! + e t [AR(1)
, stationary ]
Fig. 1.1 4 u,
, stationary ]
Fig. 1.1 5 u,
- 0. 5 Mr _! + e, - 0.8e t _i [ARMA(1,1)
, stationary]
Fig. 1.1 6 u,
Fig. 1.1 7 u
Fig. 1.1 8 u,
Fig. 1.1 9 u
Fig. 1.2 0 u,
, stationary ]
, stationary ]
, stationary ]
, non-stationary ]
, non-stationary ]
Observation
Observation
33
34
Observation
Observation
Introduction an d Overvie w
Observation
Observation
35
36
Introduction an d Overvie w
Observation
Observation
Introduction an d Overvie w
Observation
Observation
37
38 Introductio
n and Overvie w
Observation
For th e stationar y process , \a\ < 1, s o th e firs t ter m an d th e contribu tions o f mor e distan t error s disappea r a s n > oo , an d u, ma y b e
approximated b y a n MA(rc ) proces s wit h increasing accurac y a s oo.
If a = 1, however, the firs t ter m doe s no t disappear , an d the approxima tion fails ; thi s follow s fro m th e failur e o f th e stationarjt y conditio n
stated above . Whe n a = 1,
40
FIG 1.23. Tim e serie s o f real mone y (log M,/Pt) an d rea l TF E (lo g Yr)
42
Appendix
Equation (11)
To prove (11) , w e need t o construct a random variable X t, wher e
44 Introductio
n an d Overvie w
If
and
Introduction an d Overvie w 4
Since y/VT^> \\W(r)&r fro m result s above , w e hav e tha t y/\/T
converges t o both \\W(r}Ar an d to N(0, 1/3) . Therefor e
48 Linea
r Transformations an d ECM s
or
50 Linea
r Transformation s an d ECM s
51
with
and
52 Linea
r Transformation s an d ECM s
2.3. A n Exampl e
An exampl e o f th e us e o f th e error-correctio n mechanis m ca n b e foun d
in Davidso n e t al. (1978) , wh o us e a homogeneou s (6 = 1 ) error-correc tion mechanis m i n th e modellin g o f consumers ' expenditure . Th e 'error '
to whic h adjustment i s made i n th e mode l i s the differenc e between th e
logarithms o f consumptio n an d income , eac h lagge d fou r quarters . Th e
error-correction ter m i s significant i n a wide variety of specifications . I n
particular, usin g quarterl y seasonall y unadjuste d dat a fro m th e Unite d
Kingdom, expresse d a t constan t price s over th e sampl e perio d o f
1958(1) -1970(IV), th e author s favou r th e model 4 (standar d error s i n
parentheses):
54 Linea
r Transformation s an d ECM s
where A = (1 2;=
i ^i)"1 = Ml)]" 1- I t i s importan t t o not e tha t thi s
formula i s applicabl e onl y wher e |2)S= i ari i s strictl y les s tha n 1 .
Otherwise, n o long-ru n equilibrium ca n b e sai d t o exis t betwee n y an d
x, as thes e quantitie s ma y diverg e increasingl y a s t -* . I n particular ,
the unconditiona l expectations are no t wel l defined.
Corresponding t o th e Bewle y transform (4), we also have
5
(7)
56
so that
For example , tak e m = n = 2 and p = 1 in (2 ) so that th e matri x of the
transformation t o th e Bardse n for m (5 ) is
-1
0
0
0
0
0
0
0
0
0
0
0
-1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
1
-1
_1
0
0
0
0
0
0
(10a)
since x' t = [yt, 1 , y,^, y,_2, xt, jc,_1; *,_2] map s ont o \' t = [&yt, 1, Ayf _i,
A*<5 Ax,-i, yt-2,xt-2] i n (5) . Fo r th e Bewle y transformatio n (4 ) an d
the sam e cas e ( m = n = 2, p = 1) , the transformatio n matrix is
0
0
0
0
0
0
0
0
0
0
0
_1
0
0
0
0
0
0
^
0
0
0
0
0
0
0
0
0
0
0
0
0
1
-1
0
(f
0
0
0
(106)
_1
since
58
0 ..
.0
0
1
0
0
0
22
(15)
since A 22 i s of full rank . From (14 ) and (15) , we can deduc e that 5
(16)
5
Th e normalizin g constan t a doe s no t appea r her e becaus e i t happen s tha t u = e in (8)
(a12 = 0) fo r th e cas e o f th e Bardse n transform , s o tha t a = 1 . Fo r th e Bewle y transform ,
there is a non-zero normalizing constant .
Th e I V estimator take s this form becaus e the origina l Xs are bein g used as instruments
in th e transforme d regression model involving y and X .
60 Linea
r Transformation s an d ECM s
Comparing (19 ) with (16), it i s clear tha t onc e agai n th e estimate s fro m
the transforme d mode l ca n b e relate d bac k t o thos e fro m OL S o n th e
ADL model , o r t o thos e fro m th e othe r transformatio n A a , throug h th e
known transformatio n matrices . Moreover , comparin g (19 ) with (12),
the sam e relatio n hold s i n estimate d parameter s a s i n th e tru e para meters, s o tha t estimate s o f function s o f thes e parameter s (suc h a s th e
long-run multipliers ) wil l b e th e sam e regardles s o f th e mode l fro m
which the y ar e calculated . Here , usin g th e Bewle y transformation , th e
long-run multiplier s 9 j appea r directl y in d b; t o calculat e the m fro m th e
ADL parameters , w e would use
and
61
-1
0
0
0
0
0
0
0
0
0
0
Ah =
-1
0
0
0
0
0
0
0
0
0
0
0
0
0
0
-1
0
0
1
-1
0
0
0
0
0
0
0
0
0
0
0
0
0
o"
0
0
0
0
0
0
0
-e
_1
0
0
-8
0
0
0
0
0
0
0
0
0
1
-1
0
0
0
0
0
0
1_
0
0
0
0
0
62 Linea
r Transformation s an d ECM s
(23)
Wickens an d Breusc h (1988 : 198 ) sho w tha t thi s i s equa l t o th e
covariance matri x o f th e sam e paramete r vecto r d b<[V, calculate d (in directly) b y applying the transformation \b to the original parameter s y
and usin g th e Jacobia n o f thi s transformatio n t o approximat e th e
estimated covarianc e matri x Vg ftiv . Tha t is ,
(24)
and Vg b ca n b e reduce d fro m thi s t o th e sam e expressio n a s tha t give n
for V^ iv i n (23) .
Both th e origina l AD L mode l an d Bardsen' s transformatio n involv e a
calculation o f th e 6, a s nonlinear function s of coefficient s in th e origina l
regression. Followin g Bardsen , a standar d formul a for a n approximatio n
to th e varianc e o f a nonlinea r functio n o f element s wit h know n
variances ca n b e use d t o comput e var( ; -). Le t / = /(a 1 ,a 2 5 > //) ;
then / = /(!, 2 2, . . ., a H) and
(25a)
In th e cas e o f th e ADL , w e hav e / = , = ^^f= o^ji
wher
7
I t migh t see m impossible tha t th e Bewle y transformatio n coul d involv e a ratio, bein g a
case o f a linea r transformatio n matri x A b applie d t o th e origina l linea r regression . Recall ,
however, tha t ther e i s als o a normalizatio n facto r applied , throug h whic h divisio n b y
another linea r functio n o f the origina l coefficient s is accomplished .
we have
Wickens an d Breusch sho w that , substitutin g dj fo r / an d comparin g
the result s o f (23 ) and (24 ) with thos e o f (25a) an d (256) , th e estimate d
variances of long-run multipliers calculate d fro m th e AD L mode l ar e th e
same a s thos e provide d i n th e I V estimato r o f th e Bewle y transform .
We sho w no w tha t th e sam e i s als o tru e o f th e error-correctio n
transformation, o r an y othe r linea r transformation , usin g th e for m
(256). Tha t is , th e metho d o f proo f doe s no t us e th e feature s o f th e
ECM o r o f an y othe r particula r transformation , bu t instea d applie s t o
the estimate s fro m an y non-singula r linear transformation . Th e import ant point , a s above , i s th e equivalenc e o f result s yielde d b y differen t
linear transformations of the model .
Consider th e long-ru n multiplier vector a s calculated fro m th e ADL,
The AD L approximatio n t o th e varianc e of the multiplie r is
where J f i s th e Jacobia n o f th e transformatio n represente d b y th e
function f ( , ) . Th e long-ru n multiplie r vecto r calculate d fro m th e
ECM is
64 Linea
r Transformations an d ECM s
it follows that
var,t (0) = JfA22 var (f), % )A22 Jf , an
d var
66 Linea
r Transformation s an d ECM s
The parameter s /3, - canno t b e determine d fro m (40 ) without know ledge o f th e margina l proces s (38) ; hence, recallin g th e definitio n o f
weak exogeneit y i n Chapte r 1 , x t i s no t weakl y exogenous fo r th e f r i n
(40). If , however , w e someho w observe d x,\,-i directly , the n w e would
be abl e t o estimat e th e /3 , fro m (39) .
The proble m identifie d b y McCallu m an d Kell y i s therefor e simpl y
one aspec t o f a broader , an d well-known , one: we canno t i n genera l
count o n unbiase d estimate s fro m model s i n whic h th e explanator y
variables ar e no t (weakly ) exogenous fo r th e parameter s o f interest. Th e
solution, i n thi s circumstance , i s therefor e join t estimatio n o f (40 ) and
(38). I f (40 ) alone i s estimated , no t onl y i s th e long-ru n solutio n no t
consistently estimated , bu t th e short-ru n adjustmen t coefficient s ar e
incorrectly estimate d a s well ; wher e $,\ t-i i s omitte d fro m th e model ,
coefficients o n x t_i ar e no t /3*(L ) bu t /3 0<5(L) + )8*(L) . I f w e d o no t
have wea k exogeneity , w e canno t conduc t vali d conditiona l inference .
68 Linea
r Transformation s an d ECM s
Properties o f Integrated
Processes
A knowledg e of the fundamenta l properties o f integrated processe s
is essentia l fo r a n understandin g o f test s fo r bot h non-stationarit y
and th e existenc e o f long-ru n equilibriu m relationships . Her e w e
define an d presen t th e importan t propertie s o f integrate d pro cesses. W e dea l wit h th e issu e o f spuriou s regression s an d sho w
how a consideratio n o f th e theor y o f integrate d processe s help s u s
to understan d th e behaviou r o f standar d estimator s i n model s
involving non-stationar y data . Severa l example s illustrat e th e us e
of Wiene r distributio n theor y i n derivin g asymptoti c result s fo r
such models .
Much conventiona l asymptoti c theor y fo r least-square s estimatio n (e.g .
the standar d proof s o f consistenc y an d asymptoti c normalit y o f OL S
estimators) assume s stationarit y o f th e explanator y variables , possibl y
around a deterministi c trend . No t al l economi c tim e serie s ar e station ary, a s w e sa w i n Chapte r 1 , an d fo r man y importan t ones , includin g
aggregate consumptio n and nationa l income , stationarit y is not eve n a
sensible approximation .
Nonetheless, regressio n method s hav e ofte n appeare d t o b e effectiv e
when analysin g such series, an d i t was not clea r tha t method s develope d
for stationar y serie s woul d no t b e vali d elsewhere . T o som e extent ,
therefore, man y analyse s o f unadjuste d non-stationar y serie s hav e bee n
carried ou t o n th e assumptio n tha t th e non-stationarit y woul d no t
matter. A s som e potentia l problem s i n doing s o became clear , however ,
econometricians naturall y looked fo r method s o f transformin g their dat a
in suc h a wa y tha t th e resultin g serie s would b e stationary , an d
therefore amenabl e t o analysi s usin g 'traditional ' econometri c o r time series methods .
One illustratio n o f th e difficultie s tha t ca n aris e whe n performin g
regression wit h clearl y non-stationar y serie s i s th e proble m o f nonsense
regression, s o name d b y Yul e (1926) , o r spurious regression, i n th e
terminology o f Grange r an d Newbol d (1974) : give n tw o completel y
unrelated bu t integrate d series , regressio n o f on e o n th e othe r wil l ten d
70 Propertie
s of Integrated Processe s
1
Th e chang e i n terminology ma y be misleading , sinc e Yul e als o use d the ter m 'spuriou s
relationships', referrin g t o a correlatio n induce d betwee n tw o variable s tha t ar e causall y
unrelated bu t ar e both dependen t o n othe r commo n variables .
2
Thi s i s accomplishe d eithe r b y includin g a functio n o f tim e a s a regressor , o r b y
subtracting a functio n o f tim e fro m al l serie s used . B y th e Frisch-Waug h theorem ,
regressing al l serie s o n tim e an d usin g thei r residual s i n a furthe r regressio n i s numericall y
equivalent t o includin g tim e a s a regressor whe n usin g th e unadjuste d series.
3
Anderso n (1958 ) extend s th e standar d asymptoti c distributio n theor y t o dea l wit h th e
de-trending o f deterministi c variable s tha t ca n b e suitabl y standardized . However , her e we
are concerned wit h integrated stochasti c processes .
72 Propertie
s of Integrated Processe s
Turn no w to th e specificatio n o f an economi c hypothesis . A n econom ist ma y wish t o describ e th e relationshi p betwee n {y t} an d {z t} wit h th e
model
where fi i i s interprete d a s th e derivativ e o f y, wit h respec t t o z t.
Conventionally, equation s suc h a s (7 ) ar e estimate d b y ordinar y leas t
squares, treatin g {u,} a s a n II D proces s independen t o f z t- Sinc e y t an d
zt ar e causall y unrelate d her e b y construction , th e derivativ e fi \ i s zer o
in th e sens e tha t n o relatio n exists ; it is not tru e t o sa y that settin g fl i t o
zero i n (7 ) give s th e tru e DGP . W e wan t t o examin e th e propertie s o f
the conventiona l estimatio n an d hypothesi s testin g procedur e applie d t o
(7) when the unknow n DGP i s in fact (4)-(6) .
Standard regressio n theor y fo r model s involvin g stationar y regressor s
would sugges t tha t pli m (fi{) = fl i = 0, an d tha t th e probabilit y o f th e
absolute valu e o f th e t -statistic fo r H 0:fii = 0 exceedin g 1.9 6 i s 5 pe r
cent. Becaus e thes e regressor s ar e integrated , however , thi s i s no t so .
Reconsider (7) . Sinc e {y t} an d {z t} ar e bot h integrate d processes , (7 )
could b e a well-define d regression wit h a non-zer o j8i , i f a relationshi p
between thes e tw o variable s existed . I f howeve r /3 j = 0, a s i s tru e her e
by (4)-(6) , w e hav e y, = /3 0 + ut. Now sinc e { y t } i s 1(1), {u t} mus t b e
1(1), whic h violates the assumptio n mad e abou t {u,} above . Ther e is an
internal inconsistenc y i n conductin g hypothesi s testin g i n th e standar d
way here , becaus e i t i s no t possibl e fo r th e erro r ter m t o b e 1(0 ) whe n
/?i i s zero .
We ca n us e Mont e Carl o method s t o examin e typica l results, i n finit e
samples, o f regression s suc h a s (7 ) where {y,} an d {z t} ar e independen t
non-stationary processes . I n th e exercis e tha t follows , w e generat e {y,}
using th e DG P (4) , wit h T = 100 , a = 0, y 0 = 0, an d o = 1. Similarly,
{zt} i s generate d wit h y = 0, z 0 = 0, an d a v = 1 i n (5) . Th e rando m
errors ar e normall y distribute d an d generate d independently , consisten t
with (6) . A t eac h replication , w e recor d (i ) th e estimate d coefficients ,
(ii) th e estimate d standar d errors , (iii ) whethe r th e nul l hypothesi s
j8j = 0 i s rejecte d whe n conventiona l 5 pe r cen t critica l value s o f th e
t -distribution ar e used , (iv ) the valu e o f th e sampl e correlatio n betwee n
y an d z , an d (v ) th e valu e o f th e Durbin-Watso n statisti c fo r residua l
serial correlation. There are N = 10,00 0 replication s usin g PC-NAIVE .
In thi s experiment , th e Mont e Carl o estimat e o f the mea n valu e of / ^
for th e experimen t i s E[fii] = -0.012 , wit h Mont e Carl o standar d erro r
(that is , th e standar d erro r o f th e Mont e Carl o estimat e o f th e mea n of
fli) o f 0.006 . Becaus e w e ar e estimatin g a mea n usin g independen t
replications, a centra l limi t theore m applie s t o th e Mont e Carl o results ,
so tha t th e sampl e mea n i s asymptoticall y normall y distribute d (i.e . a s
N-* oo) . Henc e we can reject th e hypothesi s tha t [& ] = 0 at T = 100,
74
75
calculated doe s not hav e a zero mean, uni t variance distribution. I n fact ,
where I i s th e mea n f-statistic , I = -0.12 (0.07 ) an d SSD(? ) = 7.3 .
Values o f |/ | > 1.9 6 ar e ver y likely wit h suc h a large standar d deviation ,
and th e empirica l critica l value s i n th e experimen t tha t ensur e a tes t
with a siz e o f 5 pe r cen t ar e approximatel y 14.5 . However , thes e
critical values are not appropriat e a t other sample sizes .
This i s th e spurious regression problem : regressio n o f a n integrate d
series on anothe r unrelate d integrate d serie s produce s f-ratio s on the
slope paramete r whic h indicat e a relationshi p muc h mor e ofte n tha n
they shoul d a t th e nomina l test level . Th e phenomeno n i s of course no t
specific t o thi s sampl e size , an d i n particula r th e proble m wil l no t
disappear a s th e sampl e siz e i s increased. Th e distributio n o f th e f-rati o
will, however , depen d o n th e sampl e size ; Fig . 3.3 show s th e grap h o f
[/3j|r] fo r T = 20, 21, . . ., 100 , togethe r wit h 2a a t eac h T , wher e
a denote s th e Mont e Carl o standar d erro r i n th e graph . Th e bia s i s
significantly differen t fro m zer o onl y at th e large r sampl e sizes , but doe s
not chang e noticeabl y wit h T. Moreover , th e valu e o f a doe s no t fal l
greatly wit h T, whic h differ s fro m wha t on e woul d expec t i f conven tional asymptotic theory were applicable .
Figure 3. 4 record s th e mea n valu e o f th e regressio n coefficien t
together wit h th e SS D and th e mea n estimate d standar d erro r (ESE ) o f
the coefficient . Ther e i s a great differenc e betwee n th e tw o measures of
uncertainty: ES E i s th e estimate d standar d erro r o f th e coefficien t j j
that th e investigato r woul d obtai n o n average , i n a regressio n o f th e
form o f (7 ) give n th e DG P i n (4)-(6) ; th e SS D i s th e Mont e Carl o
estimate o f th e tru e standar d deviatio n o f thi s paramete r estimate . A s
Fig. 3. 4 shows , th e economis t woul d repor t a sever e underestimat e o f
the uncertaint y in the estimat e o f /J t .
The mea n valu e o f th e ^-statisti c shown i n Fig . 3.5 change s little a s T
increases fro m 2 0 t o 100 , bu t th e standar d deviatio n o f t increase s
FIG 3.3. Mea n value of the spuriou s regressio n coefficien t wit h 2c r (th e
Monte Carl o standar d error)
76
FIG 3.5. Mea n valu e o f th e 'r-test ' o f // 0:/31 = 0, wit h 2SS D (th e
Monte Carl o base d samplin g standard deviation )
rapidly. Thu s the proble m become s worse as T increases ; rejectio n of
the nul l hypothesi s of n o relatio n betwee n th e y, and z t serie s become s
more likely , despit e one' s initia l intuitio n that , i f th e serie s reall y ar e
unrelated, thi s feature shoul d eventuall y dominate a s T > oo . Figur e 3. 6
records th e rejectio n frequencie s for ever y sampl e siz e considere d i n th e
simulation exercise ; Prd^ft = 0)| ^ 2) is 0.30 at T = 20, already greate r
than th e nomina l siz e o f th e test , an d th e proble m worsen s a s T i s
increased becaus e th e rejectio n frequencie s also increase steadil y with T.
The outcome s o f th e simulation s revea l th e danger s o f usin g critica l
values justified i n on e contex t (e.g . IID processes ) t o conduc t inference s
with statistic s compute d fro m dat a generate d b y a ver y differen t
probability mechanism .
With th e DG P i n (4 ) an d (5) , the proble m o f discriminatin g betwee n
genuine interdependenc e an d spuriou s regression s i s difficul t t o solv e
77
78
79
FIG 3.9. Frequenc y distributio n o f R fo r tw o 1(2 ) processe s wit h inde pendent II D second difference s
probability o f findin g R 0 i n thi s las t case , althoug h th e populatio n
value anticipate d unde r th e nul l i s zero. Th e mos t likel y sample valu e is
R~l.
If th e degree s o f integratio n o f th e dat a serie s ar e unknown , mixtures
of case s (i)-(iii ) ar e possible . Fo r T = 100, Tabl e 3. 1 summarize s th e
outcomes.
Denote th e orde r o f integration o f y, and x t b y di an d rf 2 respectively,
and le t d = max{di, d 2}. Th e mea n o f J R i s close t o zer o i n ever y case ,
but it s standar d deviatio n increase s wit h d^ + d2. Th e estimat e o f th e
mean o f fi \ i s relativel y smal l compare d wit h th e SSD , especiall y when
TABLE 3.1. Feature s o f regression s amon g serie s wit h variou s order s o f
integration
Type3
1(0), 1(0)
1(1)
1(2), 1(2)
1(0), 1(1)
1(1), 1(0)
1(2), 1(1)
1(1), 1(2)
1(1),
R SSD(R
0
2
4
1
1
3
3
0,,0004
-0,,006
0,,004
0,,0004
0,,0008
-0,,023
-0.,013
0.101
0.490
0.818
0.099
0.101
0.613
0.610
0.101
0.102
0.103
0.031
0.384
3.84
0.0054
0.,102
0,,631
1..974
0,,033
0,.417
33,,52
0..036
0,,0493
0,,7570
0,,9406
0,,0458
0,,0486
0.,8530
0,,8444
80
81
82 Propertie
s o f Integrated Processes
The result s fro m (10 ) ar e muc h a s one woul d expec t give n those implie d
by (3 ) an d (9 ) abov e (see , again , Durlau f an d Phillip s 1988) . A s before ,
the distributio n o f c diverge s an d y tend s i n probabilit y t o zero , bu t / ?
has a non-degenerat e distributio n asymptoticall y (i.e . doe s no t converg e
to zero) . Test s fo r H 0: / ? = 0 diverg e i n distribution, tendin g t o lea d th e
investigator falsel y t o rejec t thi s nul l hypothesis . Estimatio n o f th e
regressions i n (9 ) an d (10 ) wil l produce substantia l residua l autocorrela tion. I t migh t b e though t tha t modellin g th e autoregressiv e erro r using ,
say, th e Cochrane-Orcut t algorith m shoul d remov e th e uni t roo t an d
thereby allo w vali d test s o f ft = 0 in (10) . Grange r an d Newbol d (1977 )
present Mont e Carl o evidenc e suggestin g that suc h a strategy i s ineffect ive in practice whe n based o n conventiona l critica l values.
In summary , th e proble m o f falsel y concludin g tha t a relationshi p
exists betwee n tw o unrelate d non-stationar y series , a proble m tha t
persists eve n a s th e sampl e siz e grow s without bound, i s no t alleviate d
by a n attempt t o remove a trend fro m th e underlying series .
In workin g with non-stationar y data , th e investigato r mus t b e particu larly careful . Whil e on e solutio n i s t o transfor m th e serie s t o achiev e
stationarity (a t th e cos t o f losin g som e informatio n abou t long-ru n
behaviour, a s we shal l se e below) , i t i s essential tha t th e investigato r b e
aware o f th e non-stationarit y i n th e dat a i f procedure s fo r modellin g
data o f thi s typ e ar e t o b e applie d appropriately . A s i t happens , testin g
84 Propertie
s o f Integrated Processes
6
Thi s definitio n i s simila r t o tha t o f Engl e an d Grange r (1987) , bu t rule s ou t som e
anomalies. Conside r th e stationary , I(1) , serie s z , = et ,_1; wher e e, is 1(0) . Integrat ing {z,} give s a serie s tha t i s 1(0) ; bu t i f we cal l {z, } itsel f a n 1(0 ) series , the n w e woul d
expect its integral {ej t o be 1(1).
and
In (12) , t o ensur e stationarity , le t u s assum e tha t y 0 i s draw n fro m th e
unconditional distribution o f y; that is, y0 ~ IID[0, a\/(l - p 2)].
It i s interestin g t o compar e severa l propertie s o f thes e series , viewe d
as possibl e DGPs . Tabl e 3. 2 summarize s som e o f th e difference s
between autoregressiv e serie s tha t ar e stationary , an d thos e containin g
one (o r more ) uni t root s (whic h requir e differencin g t o b e mad e
stationary). Th e propertie s i n th e right-han d colum n o f th e tabl e hol d
for integrate d serie s generally . Nonetheless , th e specificatio n (13 ) i s a
special one , an d i n a genera l treatmen t w e wan t a les s restrictiv e
TABLE 3. 2. Som e propertie s o f stationary an d integrate d processe s
Variance
Conditional variance
Autocorrelation
function a t lag i
Expected time between
crossings of y = 0
Memory3
a
DGP (12 )
(1(0))
DGP (13 )
(1(1))
Finite
(a\(l - p2)-i )
Unbounded
(grows as ta^)
Pi = P
Finite
Temporary
Pi = Vl - (i/f) - 1 V i as t - o o
Infinite
Permanent13
86 Propertie
s o f Integrated Processe s
88 Propertie
s of Integrated Processe s
90 Propertie
s of Integrated Processe s
Similarly:
Since
Thus:
, 1/3 ) T~
,1 T
2: J rdW(r) N(0
, 1/3 )
3: W(l) N(0
, 1)
1T
4: W(r)dW(r)
r f i 1-1/
5: J o W(r) 2 dH J
,1 T
e moment 5
y Sampl
(l/2)(x
2fi
^ W(r)dV(r) N(0
6: J o ( r - a)W(r)dr N(0
(l) - 1 ) T~
l 2
ly
^ y,-iu,
t =l
, 1) , T), T~
52
/ ^ ty, (a = 0)
where Y = (1/60 )
(8 - 25 + 202)
a
92 Propertie
s o f Integrated Processe s
However, b y (21) ;
by (24) ;
also;
Further,
94
From (21),
From (22) ,
Also,
and
96 Propertie
s of Integrated Processe s
98 Propertie
s of Integrated Processes
100 Testin
102 Testin
(b) Mode l
25
50
100
250
500
00
(c) Mode l
25
50
100
250
500
00
(3fl)/(8;1
0 0.9
5 0.97
5 0.9
-9.3
-9.9
-10.2
-10.3
-10.4
-10.5
-7.3
-7.7
-7.9
-8.0
-8.0
-8.1
ca
J, J
-5.5
-5.6
-5.7
-5.7
-5.7
1.01
0.97
0.95
0.93
0.93
0.93
1.40
1.35
1.31
1.28
1.28
1.28
1.79
1.70
1.65
1.62
1.61
1.60
2.28
2.16
2.09
2.04
2.04
2.03
(3b)/(8b )
-17.2
-14.6
-18.9
-15.7
-19.8
-16.3
-20.3 -16.6
-20.5
-16.8
-20.7
-16.9
-12.5
-13.3
-13.7
-14.0
-14.0
-14.1
-10.2
-10.7
-11.0
-11.2
-11.2
-11.3
-0.76
-0.81
-0.83
-0.84
-0.84
-0.85
0.01
-0.07
-0.10
-0.12
-0.13
-0.13
0.65
0.53
0.47
0.43
0.42
0.41
1.40
1.22
1.14
1.09
1.06
1.04
(3c)/(8c)
-22.5
-25.7
-27.4
-28.4
-28.9
-29.5
-17.9
-19.8
-20.7
-21.3
-21.5
-21.8
-15.6
-16.8
-17.5
-18.0
-18.1
-18.3
-3.66
-3.71
-3.74
-3.75
-3.76
-3.77
-2.51
-2.60
-2.62
-2.64
-2.65
-2.66
-1.53
-1.66
-1.73
-1.78
-1.78
-1.79
-0.43
-0.65
-0.75
-0.82
-0.84
-0.87
-11.9
-12.9
-13.3
-13.6
-13.7
-13.8
-19.9
-22.4
-23.6
-24.4
-24.8
-25.1
(3 )/(8:1
(a) Model fl
25
-2.66
50
-2.62
100
-2.60
250
-2.58
500
-2.58
00
-2.58
0 0.9
5 0.97
5 0.9
-2.26
-2.25
-2.24
-2.23
-2.23
-2.23
-1.95
-1.95
-1.95
-1.95
-1.95
-1.95
-1.60
-1.61
-1.61
-1.62
-1.62
-1.62
0.92
0.91
0.90
0.89
0.89
0.89
1.33
1.31
1.29
1.29
1.28
1.28
1.70
1.66
1.64
1.63
1.62
1.62
2.16
2.08
2.03
2.01
2.00
2.00
-3.00
-2.93
-2.89
-2.88
-2.87
-2.86
-2.63
-2.60
-2.58
-2.57
-2.57
-2.57
-0.37
-0.40
-0.42
-0.42
-0.43
-0.44
0.00
-0.03
-0.05
-0.06
-0.07
-0.07
0.34
0.29
0.26
0.24
0.24
0.23
0.72
0.66
0.63
0.62
0.61
0.60
(3c)/(8c)
-3.95
-4.38
-4.15
-3.80
-4.04
-3.73
-3.69
-3.99
-3.98
-3.68
-3.66
-3.96
-3.60
-3.50
-3.45
-3.43
-3.42
-3.41
-3.24
-3.18
-3.15
-3.13
-3.13
-3.12
-1.14
-1.19
-1.22
-1.23
-1.24
-1.25
-0.80
-0.87
-0.90
-0.92
-0.93
-0.94
-0.50
-0.58
-0.62
-0.64
-0.65
-0.66
-0.15
-0.24
-0.28
-0.31
-0.32
-0.33
-2.33
-1.65
-1.28
1.28
1.65
1.96
(c) Model
25
50
100
250
500
00
N(0, 1 )
00
-1.96
2.33
104 Testin
DGP Model
s yieldin g similar tests 1
(i) y t = Pyt-i + ut, y0 = 0 (3c)
, (36) , (3c )
(ii) y t = py,-! + ut, arbitrary y 0 (36)
, (3c )
(iii) y t = [i+ py t-i + ut, arbitrary y 0 (3c
)
(iv) y t = [a + yt + pyt-i + u f> arbitrar y y 0 Extensio n o f (3c ) necessar y
Thus, fo r example , i n cas e (i) , i f th e mode l i s give n b y (3c) , th e
appropriate critica l value s ar e give n b y Table s 4.1(c ) an d 4.2(c) . Th e
same table s ca n b e use d t o conduc t inferenc e i n (iii) , despit e a non-zer o
value o f n i n th e DGP , becaus e (3c ) yield s a simila r test . Similarit y
implies tha t th e distribution s o f p an d it s associate d ^-statisti c ar e no t
affected b y th e value , unde r th e null , o f th e nuisanc e parameter , an d
the critical value s ar e th e sam e a s the one s tha t woul d appl y fo r n = 0,
namely, those i n Tables 4.1(c ) an d 4.2(c).
There ar e a numbe r o f noteworth y additiona l points . I n cas e (i ) ther e
are n o nuisanc e parameters , s o tha t similarit y i s a trivia l property . I n
general, a s this summar y suggests , a simila r tes t havin g a Dickey-Fuller
distribution require s tha t th e mode l use d contai n more parameter s tha n
the DGP . I n order to hav e a similar test fo r (iv) , one woul d the n nee d a
model wit h a ter m suc h a s t 2, necessitatin g anothe r bloc k o f critica l
values i n eac h o f Table s 4. 1 and 4.2 . I n cas e (ii) , fo r example , w e nee d
at leas t mode l (36 ) (wit h a constant ) t o allo w fo r th e unknow n startin g
value. I n cas e (iii ) w e hav e a n unknow n constan t an d nee d th e tren d
term i n model (3c ) t o allo w for it s effect .
Each o f thes e simila r test s i s als o exac t i n finit e samples , provide d
appropriate critica l value s ar e available . I n general , however , i t wil l b e
necessary t o abando n exac t test s i n orde r t o us e variant s o f th e
Dickey-Fuller tes t wher e ther e ar e mor e unknow n parameters . Thes e
parameters ca n typicall y be estimated , s o that asymptoticall y they can b e
accounted fo r an d a tes t provided . Again , Kivie t an d Phillip s offe r
general exac t an d simila r test s fo r DGP s wher e th e dynamic s ar e
restricted t o first-order , a s wel l a s demonstratin g th e similarit y o f th e
tests just mentioned .
In th e cas e o f exac t parameterizations , suc h a s cas e (iii ) wit h mode l
(3>), w e d o no t hav e simila r test s wit h th e Dickey-Fulle r distributions .
However, a s West (1988 ) showed , the f-statistic s i n th e exactl y paramet erized cas e (wit h exogenou s item s suc h a s a constan t i n th e DGP ) ar e
asymptotically normal , jus t a s ar e f-statistic s use d fo r standar d prob lems. I n finit e samples , however , th e Dickey-Fulle r distribution s ma y
be a better approximatio n tha n th e norma l distribution . We will explor e
this asymptoti c normalit y further i n Chapte r 6 below.
Critica l value s ar e those corresponding t o the mode l use d i n Table 4.1 or 4.2 .
106 Testin
terms ar e al l 1(0 ) an d appropriat e scalin g ensure s tha t th e variance covariance matri x i s asymptoticall y block-diagonal . (Tha t is , al l cross product term s o f 1(0 ) an d 1(1 ) variable s i n th e matri x ar e asymptoticall y
negligible.) I t i s thi s asymptoti c orthogonality tha t drive s th e result ,
much as , i n a standar d regressio n model , on e use s th e orthogonalit y of
the informatio n matri x t o prov e th e statistica l independenc e o f th e
estimated coefficien t vecto r fro m th e estimat e o f the standar d error . Th e
asymptotic theor y an d th e issu e o f 'appropriate ' scalin g ar e discusse d
later i n this chapter an d i n Chapter 6 .
By allowin g the DG P t o tak e th e for m (4 ) rather tha n th e muc h mor e
restrictive AR(1 ) for m (3) , w e hav e expande d th e clas s o f model s t o
which we can validl y appl y unit-roo t test s of thi s type . Not e that , as it
will generall y b e th e cas e tha t p i s unknown even wher e y t i s strictly an
AR(p) process , i t i s generall y safe r t o tak e p t o b e a fairl y generou s
number; i f too man y lags ar e presen t i n (5) , th e regressio n i s free t o se t
them t o zer o a t th e cos t o f som e los s i n efficiency , wherea s to o fe w lags
implies som e remainin g autocorrelatio n i n (5 ) an d henc e th e inapplicab ility o f even th e asymptoti c distributions i n Tables 4. 1 an d 4.2 . On e can ,
of course , perfor m test s fo r autocorrelatio n o n th e estimate d residual s
from (5 ) i n orde r t o chec k th e acceptabilit y o f th e premis e tha t thes e
residuals ar e whit e noise . Alternatively , mode l selectio n procedure s ca n
be used t o choose p, and test fo r a unit root, jointly (see Hal l 1990) .
We have , therefore , a class o f tests fo r th e uni t root whic h can validly
be applie d t o serie s tha t follo w AR(p ) processe s containin g n o mor e
than on e uni t root . Th e nex t natura l ste p i s to attemp t t o exten d furthe r
the clas s of series t o which we can appl y such tests , ideall y in such a way
as t o allo w exogenou s variable s t o ente r th e proces s a s well . Sai d an d
Dickey (1984 ) provid e a tes t procedur e vali d fo r a genera l ARM A
process i n th e errors ; Phillip s (1987a ) an d Perro n an d Phillip s (1988 )
offer a still more genera l procedure .
While th e Said-Dicke y approac h doe s represen t a generalizatio n o f
the Dickey-Fulle r procedure , i t agai n yield s test statistic s wit h th e sam e
asymptotic critica l value s a s thos e tabulate d b y Dicke y an d Fuller . Th e
particular advantag e o f thi s tes t i s tha t w e ca n appl y i t no t onl y t o
models wit h M A part s i n th e errors , bu t als o t o model s fo r whic h (as is
typically th e case ) th e order s o f th e A R an d M A polynomial s i n th e
error proces s ar e unknown . Th e method involve s approximating the tru e
process b y a n autoregressio n i n whic h the numbe r o f lag s increases wit h
sample size .
Begin b y assuming that th e data-generatio n proces s follows :
108 Testin
or equivalently
110 Testin
and
It i s eas y t o calculat e fro m thes e regression s th e coefficien t estimate s
and th e '^-statistics ' fo r each . Fo r test s o f th e significanc e o f p,- , th e
statistics ar e the n adjuste d t o reflec t autocorrelatio n i n th e corresponding Uit series . (W e wil l omi t subscript s a , b , o r c o n u t t o simplif y
notation.) I f we defin e
and
W e trea t th e initia l observatio n a s fixe d a t zero ; not al l statistics here are invarian t t o
the initia l value. Se e Phillips (1987a) an d Perron (1988).
3
Thes e statistic s ar e vali d fo r eithe r choic e o f S 2Tt give n abov e (i.e . the Phillip s o r
Newey-West forms) .
112 Testin
or, alternatively ,
and
and
having th e limitin g distribution s tabulate d i n Table s 4.1(c ) an d 4.2(c )
respectively. Th e quantit y D x i s defined a s the determinan t o f th e inne r
product o f the dat a matri x with itself: for (8c),
114 Testin
The Phillips-Perro n correction s t o th e standar d Dickey-Fulle r statist ics mus t howeve r b e use d cautiously . Again , th e accumulate d evidenc e
of severa l Mont e Carl o simulatio n studie s suggest s tha t th e non-para metrically correcte d tes t statistic s d o no t alway s hav e th e correc t size s
even in fairl y larg e samples .
Schwert (1989 ) make s thi s poin t forcefully . Hi s results , amplifyin g
those i n th e Phillips-Perro n simulation s reporte d earlier , sho w tha t th e
critical value s o f th e augmente d Dickey-Fulle r tes t statistics , give n b y
the standar d Dickey-Fulle r tables , ar e muc h mor e robus t t o th e
presence o f movin g averag e term s i n th e error s o f th e random-wal k
process tha n ar e th e correspondin g non-parametricall y adjuste d Dickey Fuller statistics . A n example , take n fro m Schwert , i s sufficien t t o
illustrate th e point .
The data-generatio n proces s i s give n by 4 y, = yt-i + ut + du t~i,
4
TABLE4.3(a). Tes t statistics for simple hypotheses in models with drif t an d trend 3
Statistic typ e Tes
t Statistic
Critica l values for Z(TI) , Z(t2) , an d Z(T^) ar e th e sam e as those fo r TI , TI, an d 7 3 respectively and ar e tabulate d i n Table 4.4.
Note als o tha t S 2U an d S\ e ar e define d wit h respect t o th e residual s o f a particula r model , an d s o diffe r acros s models (8a), (8b),
and (8c) . c ti(j) i s the it h diagonal element of the invers e second-moment matrix of the regressors i n model j .
Sources: Dickey and Fuller (1981 ) and Perro n (1988) .
Critical values for Z(<>i), Z(<J> 2 )> and Z(<t>3) are the same as those for <!>!, <I>2, and <53 respectively and are tabulated in Table
4.5. Note also that S2U and S2T( are defined with respect to the residuals of a particular model, and so differ across models (8a),
(8b), and (8c).
Sources: Dickey and Fuller (1981) and Perron (1988).
y o f a smaller value 3
0.90 0.9
5 0.97
5 0.9
2.61
2.56
2.54
2.53
2.52
2.52
2.97
2.89
2.86
2.84
2.83
2.83
3.41
3.28
3.22
3.19
3.18
3.18
3.20
3.14
3.11
3.09
3.08
3.08
3.59
3.47
3.42
3.39
3.38
3.38
4.05
3.87
3.78
3.74
3.72
3.71
2.85
2.81
2.79
2.79
2.78
2.78
3.25
3.18
3.14
3.12
3.11
3.11
3.74
3.60
3.53
3.49
3.48
3.46
118
y of a smaller value a
0.025 0.0
0.10
0.90
0.95
0.975 0.9
r\.
4.12
3.94
3.86
3.81
3.79
3.78
; mode l (8b)
6.30
5.18
4.86
5.80
4.71
5.57
4.63
5.45
5.41
4.61
4.59
5.38
Me = 0 , yc = 0; model (8c )
6.75
4.67
5.68
5.13
5.94
4.31
5.59
4.16
4.88
4.07
4.75
5.40
4.05
4.71
5.35
4.03
4.68
5.31
Xc
=
0
;
model ( 8c)
0>3; DGP : (8c) wit h PC = 1 , '
1.33
7.24
8.65
0.90
1.08
5.91
1.37
5.61
6.73
7.81
0.93
1.11
7.44
0.94
1.12
1.38
5.47
6.49
5.39
6.34
7.25
0.94
1.13
1.39
0.94
1.39
5.36
6.30
7.20
1.13
5.34
6.25
0.94
1.39
7.16
1.13
Pc = l ,
1.10
1.12
1.12
1.13
1.13
1.13
7.88
7.06
6.70
6.52
6.47
6.43
8.21
7.02
6.50
6.22
6.15
6.09
10.61
9.31
8.73
8.43
8.34
8.27
120 Testin
Not e tha t th e firs t sequenc e too k th e smallest numbe r (i.e . 1 ) of uni t root s a s it s firs t
maintained hypothesis .
122 Testin
(21)
124 Testin
and
Hence
where
so that
where
126 Testin
The analytica l densitie s o f Vr,i,2 > ^7,2, 3 > 07,i > 07, 2 > an d 0 r> 3 ca n b e
found fro m Tabl e 3.3 . I n th e cas e o f <j) T^ w e use th e fac t tha t th e squar e
of W(l) i s distribute d a s ^ 2(1), recallin g that W(l) i s standar d normal.
The closed-for m densit y fo r th e functiona l t o whic h ^7,2, 2 converge s is
more difficul t t o derive , but a n asymptoti c expansion i s given by Abadi r
(1992).
If, a s i n thi s Dickey-Fulle r test , w e ar e particularl y interested i n th e
estimator o f p c an d it s ?-ratio , t(p c), choosin g th e appropriat e element s
from abov e gives
and
the secon d elemen t o f the 3 x1 matrix Vj1 ^. Fro m (27 ) w e note tha t
(pc 1) converge s a t rat e O p(T~l) instea d o f th e conventiona l
Op(T~^2). Similarly , fro m (28) , th e correspondin g ?-rati o ha s a non degenerate distributio n differin g fro m th e standardize d normal distribu-
128 Testin
130 Testin
4.6.4. Example: Instrumental Variables Test for Unit Roots (Hall 1989)
The non-parametri c statistic s describe d i n exampl e 4.6. 3 ar e know n no t
to perfor m wel l i n finit e sample s i n th e presenc e o f negativ e moving average error s (se e Schwer t 1989) . Hal l (1989 ) propose d estimatio n b y
instrumental variable s a s a n alternativ e t o th e us e o f non-parametri c
corrections. H e showe d tha t i n th e regressio n mode l y, = pyt~\ + ut,
where u t i s a moving-averag e proces s o f som e specifie d orde r an d p i s
equal t o 1 under H 0, the n p iv ha s the standar d Dickey-Fulle r distribu tion.
The intuitio n for thi s result ma y b e easil y described: p OLS i n th e abov e
model doe s no t hav e th e standar d Dickey-Fulle r distributio n because o f
the bia s induce d b y th e correlatio n betwee n y r _i an d u, (whe n u t i s an
ARMA(p,q) process) . I t i s therefor e necessar y t o us e a correctio n
factor t o remov e thi s bias . Thi s bia s doe s no t appea r when , say , y,_ 2 is
used a s a n instrumen t fo r y,_ i an d u t i s a n MA(1 ) process . Th e
8
A s note d above , th e finite-sampl e behaviou r o f thes e tw o estimator s ma y b e quit e
different (se e Schwer t 1989) .
Note tha t
132 Testin
g fo r a Unit Root
But a 2 = (1 + 0 2 )a 2 . Henc e
The las t equalit y follows from th e expressio n fo r a 2 give n previously, (i )
now follows routinely from (40) .
Proof of (ii).
134 Testin
Co-integration
We defin e th e concep t o f co-integratio n o f integrate d time-serie s
and giv e severa l examples . A n importan t theore m du e t o Grange r
on alternativ e representations o f a system of co-integrated variables
is state d an d it s proo f i s sketched . W e the n discus s th e Engle Granger two-ste p procedur e fo r estimatin g th e parameter s
characterizing the co-integratin g relationship.
In Chapte r 1 we discusse d ou r us e o f th e wor d 'equilibrium' . Th e ide a
that variable s hypothesize d t o b e linke d b y som e theoretica l economi c
relationship shoul d no t diverg e fro m eac h othe r i n th e lon g ru n i s a
fundamental one. 1 Suc h variable s ma y drif t apar t i n th e shor t ru n o r
because o f seasona l effects , bu t i f the y wer e t o diverg e without bound ,
an equilibriu m relationshi p amon g suc h variable s coul d no t b e sai d t o
exist. Th e divergenc e fro m a stabl e equilibriu m state must be stochastic ally bounde d and , a t som e point, diminishing over time . 'Co-integration '
may b e viewe d a s th e statistica l expressio n o f th e natur e o f suc h
equilibrium relationships.
The concep t o f co-integration is a powerful on e becaus e i t allow s us t o
describe th e existenc e o f a n equilibrium , o r stationary , relationshi p
among tw o o r mor e time-series , eac h o f whic h i s individuall y non stationary.2 Tha t is, while the componen t time-serie s ma y have moment s
such a s means , variances , an d covariance s varyin g wit h time , som e
linear combinatio n o f thes e series , whic h define s th e equilibriu m rela tionship, ha s time-invariant linear properties .
The wor d 'co-integration ' clearl y demand s a forma l definitio n o f
'integration', an d thi s wa s provided i n Chapte r 3 . Informally , a serie s is
said t o b e integrate d i f it accumulate s some pas t effects ; suc h a serie s is
non-stationary becaus e it s futur e pat h depend s upo n al l suc h pas t
influences, an d i s no t tie d t o som e mea n t o whic h i t mus t eventuall y
1
Familia r example s o f hypothesize d long-ru n relationship s includ e th e quantit y theor y
of money , th e Fishe r effect , th e permanent-incom e hypothesi s o f consumption , an d
purchasing-power parity .
2
Typically , i n economi c application s on e look s fo r th e existenc e o f co-integratin g
relationships amon g variable s individuall y integrate d o f orde r one . Th e deviatio n fro m th e
equilibrium relationshi p i s thu s integrate d o f orde r zer o (i.e . i s stationary ) whe n th e
variables ar e co-integrated .
Co-integration 13
5.1. A n example
In orde r t o illustrat e th e precedin g discussion , conside r a simpl e
example. Tw o serie s {x t} an d {y t} ar e eac h integrate d o f orde r 1 and
evolve accordin g t o th e followin g data-generation process: 4
3
Whe n regardin g a co-integratin g combination a s a n 'equilibrium ' relationship , i t i s
natural t o expec t this combination to b e integrate d o f order zero. However, definitionally,
any reductio n i n th e orde r o f integrationsay , fro m d t o d b (wher e b > 0)is
sufficient fo r th e variables to b e calle d 'co-integrated' .
4
Th e exampl e is taken fro m Engl e an d Grange r (1987).
138 Co-integratio
Co-integration 13
140 Co-integratio
(10)
The determinan t |A(L) | o f a polynomia l matri x A(L ) i s a scala r
polynomial.
A familia r exampl e o f a polynomial matri x i s A (A) = (A 0 - AI) , which
occurs i n the characteristi c equatio n
which ma y b e solve d fo r eigenvalue s o f th e matri x AQ . Ever y matri x
satisfies it s ow n characteristi c equatio n (th e Cayley-Hamilto n theorem )
in that , i f we le t /(A ) = |A(A)| , the n /(A ) = 0 (wher e thi s i s interprete d
as a matri x expression) . I n general , i f A(L ) = 2f= oA;L' , the n w e wil l
also us e the notatio n A(B ) = 2f=oA,-B', fo r a matrix argument B.
Co-integration 14
such tha t
142 Co-integratio
and
Co-integration
143
144 Co-integratio
(23)
Comparing (21 ) wit h (23) , th e analysi s can b e see n t o repeat , leadin g t o
| A O/) | after k - 1 steps. Thus ,
the laten t root s ca n b e foun d b y equatin g either expressio n t o zer o an d
solving. Sinc e A ( ) i s n x n , O i s n k x n k an d s o ha s n k eigenvalues ,
as required.
From (13) , whe n B = I, i f A(l ) ha s ran k r < n, the n |A(1) | = 0 an d
hence A(L ) ha s n r uni t roots . Conversely , i f A(l) ha s ran k n , A(L )
has none o f its eigenvalues equal to unity.
Next, derivative s o f polynomia l matrice s wit h respec t t o thei r argu ments will b e needed , an d w e have
Co-integration 14
In term s o f (16) ,
146 Co-integratio
(25)
where th e e t satisf y assumption s (3.16a)-(3.16d ) an d th e firs t k dat a
points Xj_fc , Xj-fc+i , . . ., x 0 ar e fixed . Th e mode l ca n the n b e rewritte n
in error-correction for m as
Co-integration 14
148 Co-integratio
(R4) E(a'x t) = -(
Co-integration 14
Thus,
and
0,
150 Co-integratio
Next,
Co-integration 15
152 Co-integratio
Co-integration 15
154
Co-integration
VAR representation
Co-integration 15
156 Co-integratio
n
5.5.2. Example 2
Co-integration 15
(48)
158 Co-integratio
whereas
Thus,
Co-integration
159
y, an d x t ar e co-integrate d 1(1 ) variable s wit h th e co-integratin g para meter give n b y a . I n th e contex t o f th e discussio n i n thi s chapter , th e
error-correction mechanis m i s estimate d i n (53 ) usin g th e tru e valu e
of th e co-integratin g parameter , whil e i n (54 ) a i s substitute d fo r a ,
where a i s derive d fro m th e stati c regressio n o f y t o n x t. Also ,
e
* = e + y( - oc)x t-]_. Le t zt = yt- x tWe nee d t o sho w that th e asymptoti c distributions of the estimator s f t
and y , o f / 3 an d y respectively , ar e th e sam e regardles s o f whethe r on e
uses a o r a (tha t is , whether one estimates (53 ) o r (54)). .
In standar d fashion , w e hav e fro m (53 ) (assumin g adequat e initia l
values)
The estimator s derive d fro m (54 ) ar e als o give n by (55 ) bu t wit h z t-\
and e f replacin g z t an d s t. From this , i t is easy to deduc e tha t th e resul t
will be demonstrate d if the followin g condition s are show n to be true :
160 Co-integratio
Secondly,
Using (59) ,
Co-integration 16
164 Regressio
c , implying y, = c + vt, E
t^(vt)
0, (1
with
corr(e ( , v t) = p an
d corr(e
t+; -,
v t) = 0 V; + 0.
e\P
Model (4)
1.0
0.9
0.8
0.5
0.0
1.0
0.9
0.8
0.5
0.0
30
0.99
26
22
0.98
0.95
17
12
0.90
0.00
5
(b) T = 200
0.999
29
0.99
18
13
0.98
0.95
9
0.90
7
0.00
5
24
20
17
12
9
6
20
15
15
10
8
6
11
10
8
7
6
5
7
7
7
6
6
5
60
54
50
38
28
6
45
40
37
30
22
7
36
33
30
25
19
7
16
15
14
12
10
5
6
6
5
6
6
6
23
15
10
7
6
4
20
13
9
7
6
4
10
8
7
6
6
5
5
4
5
5
6
5
61
41
29
17
10
5
48
32
24
14
9
5
38
27
20
12
8
4
18
13
11
7
6
5
5
5
6
6
7
5
(a) T = 50
0.999
166
022
1.0
0.999
1.0
0.999
0.999
0.999
0.999
1.0
0.999
1.0
0.999
0.999
0.0
0.0
=0.0
=0.0
0.0
=0.0
1(0)
NI(1)
NI(1)
NI(1)
NI(1)
NI(2)
NI(1)
1(0)
NI(1)
1(0)
NI(1)
NI(1)
3
A regressio n i s defined t o b e balance d i f and onl y i f the regressan d an d th e regressor s
(either individuall y o r collectively , a s a co-integrate d set ) ar e o f th e sam e orde r o f
integration. Th e mer e fac t tha t a regressio n i s unbalance d ma y no t b e a matte r fo r
concern; fo r example , AD F statistic s ar e compute d fro m model s that , i n thi s terminology ,
are unbalanced . The y ar e nevertheles s vali d tool s fo r inferenc e a s lon g a s th e correc t
critical value s ar e used.
regression i s balance d (ther e exist s a linea r combinatio n o f th e regres sors tha t ha s th e sam e orde r o f integratio n a s th e regressand ) an d siz e
distortions d o not appear . Cas e F resemble s C except tha t 9 U i s close t o
1, indicatin g that co-integratio n i s broken betwee n th e regressors ; non e
the less , siz e distortion s ar e no t detectabl e a s lon g a s # 12 remain s
non-zero. Thi s las t findin g demonstrate s th e difficult y o f distinguishing,
at modes t sampl e sizes , th e result s o f regression s wit h co-integrate d
regressors fro m thos e wit h regressor s o f differing , bu t bot h strictl y
positive, orders .
We see , i n summary , tha t fo r integrate d serie s (or , i n finit e samples ,
for th e borderline-stationar y serie s examine d i n thes e papers) , wit h
p=l, siz e distortion s ma y emerge whe n ther e i s no linea r combinatio n
of regressor s tha t has the sam e order of integration as the regressand. 4
For a n intuitive view of these results , le t u s return t o th e consumptio n
example an d conside r th e order s o f integratio n o f th e variable s o n th e
two side s o f th e regression . Consumptio n an d incom e ar e bot h typically
variables integrate d o f orde r one . Thus , th e regressio n (2 ) ha s a n 1(0 )
variable (difference d consumption ) regressed o n a n 1(1 ) variabl e (lagge d
income i n level ) an d th e regressio n i s unbalanced ; th e investigato r i s
attempting t o explai n a n 1(0 ) variable b y a variabl e integrate d o f highe r
order. Thi s strateg y wil l eventuall y fail , a s th e tw o variable s mus t
diverge b y ever-large r amounts . Therefore , a requiremen t o f estimatio n
with integrate d variable s mus t be balanc e i n th e order s o f integratio n of
the variable s o n th e left-han d an d right-han d side s o f th e regressio n
equation. However , ther e ma y be circumstance s i n whic h a test wil l b e
designed t o involv e regressan d an d regressor s havin g differen t order s of
integrationfor example , efficienc y test s suc h a s thos e mentione d
above. W e mus t bea r i n mind , o f course , tha t tes t statistic s fro m suc h
regressions will have non-standard distributions.
The importanc e o f th e latte r poin t follow s fro m th e observatio n that ,
even whe n th e regressan d (e.g . y t) an d th e regresso r (x t) ar e bot h
integrated o f orde r 1 an d ar e co-integrate^ , th e ^-statisti c o n th e
coefficient o f x t stil l ha s a non-standar d distributio n whic h make s
ordinary t an d norma l tables unusabl e fo r purposes o f inference. O n th e
other hand , i f th e orde r o f integratio n o f bot h side s i s zero (whic h may
be ensure d b y lookin g fo r a co-integrate d se t o f regressor s an d usin g a
sufficiently difference d ter m a s th e regressand) , th e t -statistics ca n b e
shown t o hav e asymptoticall y norma l distributions . Thi s implie s som e
advantage t o th e us e o f dynami c rathe r tha n static regressions , sinc e
lagging variable s an d includin g the m a s regressor s ofte n ha s th e sam e
effect a s providin g a co-integrate d se t o f regresso r variables . Th e
4
Campbel l an d Dufou r (1991 ) offer , a s a wa y o f overcomin g th e Mankiw-Shapir o
problem, a n alternativ e non-parametric test o f orthogonality which i s independent of some
nuisance parameters i n the DGP .
168 Regressio
170 Regressio
DFb
0 0.9
0 0.9
5 0.97 5 0.9 9
(a)n = 0
50
-3,.57
100
-3,,50
^ .47
200
400
-3,.45
(b) ni = 0.001
50
-3,.57
100
-3,,49
200
-3,,47
400
-3,,46
(c) (i = 0.010
!
,57
50
100
-3,,50
200
-3,.47
400
-3,,46
(d) /j: = 0.10
50
-3,,53
100
-3.,41
200
-3,,34
400
-3,,17
0.60
-3.22
-3.17
-3.15
-3.13
-2.92
-2.89
-2.88
-2.87
-2.60
-2.59
-2.58
-2.57
-1.55
-1.56
-1.56
-1.56
-0.40
-0.42
-0.42
-0.43
-0.03
-0.06
-0.06
-0.07
0.30
0.26
0.25
0.25
0.66
0.62
0.62
0.62
-3.22
-3.17
-3.14
-3.14
-2.93
-2.89
-2.88
-2.87
-2.60
-2.58
-2.57
-2.57
-1.55
-1.56
-1.56
-1.56
-0.40
-0.42
-0.42
-0.43
-0.03
-0.06
-0.06
-0.07
0.29
0.26
0.26
0.24
0.65
0.63
0.63
0.62
-3.22
-3.16
-3.14
-3.13
-2.93
-2.90
-2.88
-2.87
-2.60
-2.58
-2.57
-2.57
-1.55
-1.56
-1.56
-1.56
-0.40
-0.41
-0.41
-0.42
-0.03
-0.05
-0.06
-0.06
0.29
0.27
0.26
0.26
0.67
0.64
0.64
0.64
-3.17
-3.08
-3.00
-2.83
-2.87
-2.80
-2.71
-2.54
-2.54
-2.48
-2.37
-2.19
-1.46 -0.22
-1.37 -0.08
-1.20 0.11
-0.94 0.34
0.16
0.30
0.48
0.70
0.49
0.61
0.81
1.02
0.89
0.98
1.18
1.39
(e) n = 0.25
50
100
200
400
(/) J* = 0.5
50
100
200
400
(g) J* = 1
50
100
200
400
(h) n = 10
50
100
200
400
N(0, 1)
.35
-3 ,10
-2.86
-2.68
-2.95
-2.72
-2.48
-2.31
-2.64
-2.41
-2.16
-2.00
-2.28
-2.05
-1.80
-1.63
-1.,02
-0,.74
-0,.52
-0,.36
0.29
0.53
0.74
0.90
0.67
0.90
1.11
1.26
1.00
1.21
1.42
1.58
1.35
1.58
1.77
1.95
-2.94
2.70
2.60
-2.50
-2.53
-2.33
-2.21
-2.14
-2.20
-2.01
-1.89
-1.82
-1.81
-1.64
-1.53
-1.46
-0.51
-0,.35
-0..25
-0,.18
0.78
0.93
1.02
1.09
1.15
1.29
1.40
1.45
1.47
1.61
1.71
1.79
1.85
1.98
2.06
2.16
-2.65
_2.52
-2.48
-2.42
-2.26
-2.15
-2.10
-2.06
-1.93
-1.84
-1.77
-1.74
-1.55
-1.47
-1.41
-1.37
-0.24
-0.17
-0.12
-0.09
1.06
1.12
1.15
1.18
1.44
1.49
1.53
1.54
1.76
1.81
1.84
1.87
2.16
2.19
2.20
2.25
-2.44
_2.38
-2.35
-2.35
-2.32
-2.03
-2.00
-1.99
-1.98
-1.96
-1.71
-1.68
-1.67
-1.66
-1.65
-1.32
-1.31
-1.30
-1.30
-1.28
-0.02
-0.02
-0.01
-0.01
0.00
1.28
1.27
1.27
1.27
1.28
1.65
1.65
1.65
1.64
1.65
1.99
1.97
1.97
1.96
1.96
2.39
2.35
2.33
2.33
2.32
172
(7)
Now, notin g tha t y, = \n bt + ^'s=\us = fj, bt + St, i t i s possibl e t o sho w
that
173
since
The scalin g factors are suc h that an y term wit h th e stochasti c componen t
of y t, namel y S t, ha s a degenerat e asymptoti c distributio n an d ma y b e
ignored asymptotically .
From (7) , w e have tha t
where
174
3/2
(13)
Thus,
(14)
From (3.22) ,
(16)
Using (3.23 ) (als o se e Phillips 1987a : 282),
Finally,
176 Regressio
Then
1 78 Regressio
or
180
or, equivalently
where
181
182 Regressio
(36)
1
All the convergenc e result s use the scale d Z' Z matri x T^Z'ZTy ; le t us
call this scaled matri x Q .
The firs t ste p in th e proo f i s to deriv e th e limitin g matrix for Q . SSW
show that , unde r certai n regularit y conditions , Q = $ > V wher e th e elements of V may b e describe d a s follows :
(a) V u an d V 12 ar e non-rando m matrice s give n b y S7= o Fn/Fii/ an d
2F=oFii/F2iy respectively . Additionally, V ]2 = V 21.
(b) V l p = V ^ = 0, p = 3, ...,2g + l.
(c) V 22 is also non-random, give n by F22F22 + S 7=0^21/^21;
(d) V mp , wher e m, p = 3, 5 , 7 , . . ., 2 g + 1, ar e rando m matrice s
involving functionals of multivariate Wiener processes .
(e) V mp, where m = 2, 4 , 6 , . . ., 2g , p = 3, 5 , 7 , . . ., 2 g + 1, are als o
random matrice s involvin g functional s o f multivariat e Wiene r pro cesses.
(f) V mp = [2/(p + m-2)] mm'pp, p = 4, 6, . . ., 2g, m = 2, 4, 6, . . .,
2g.
This i s the firs t tim e w e have used multivariat e Wiener processes . Th e
mathematical detail s involve d i n goin g fro m univariat e t o multivariat e
Wiener processe s ar e comple x an d wil l no t b e deal t wit h her e (fo r a
good account , se e Phillip s an d Durlau f 1986) . Howeve r th e generaliza tions fro m ou r analysi s in Chapte r 3 can b e understoo d intuitivel y fairl y
easily an d the appendix sketche s th e bivariate case .
Thus, eac h elemen t o f a standardize d n x 1 multivariat e Wiene r
process W(r ) i s a univariat e Wiene r proces s an d th e element s o f W(r )
are independent . I n particular , W(l ) i s the multivariat e standar d norma l
183
density, tha t is , N(0, !). Further, W(r ) e C[0,1]", wher e C[0,1 ] is the
space of continuous function s defined on [0,1] .
Convergence result s analogou s t o (3.17) , fo r a sequence o f mean zero
random vector s {u (}, ca n b e prove d b y definin g standardize d sum s such
as
From (35) ,
184 Regressio
by substitutin g fo r s
Thus,
(38)
As note d abov e th e matri x V is the limitin g matrix of Q .
The asymptoti c distribution of
where
(a) (j) m fo r al l m ^ 3 are functional s of multivariate Wiener processes ;
(b) 0 2 = 02 i + 022 , wher e ft, 2 = vec[F 22W(l)'S1/2], W(l ) is th e multi variate standar d norma l densit y function, and
Finally,
g+i).
Consolidatin g
186
(40)
where 61 = fa, 6 2 = fa , an d 6 3 = f a + fa , Z l>t
Z 3; f = x t. It ma y also be shown that
Z 2 ( = 1 , an d
(41)
where 0(L ) = (1 + faL)' 1 an d 0*(L) = (1 - L)" 1 [0(L) Note fro m (41 ) tha t F 2 i(L) = 0. Thi s implies , b y referrin g t o th e
description o f th e V matri x above , tha t V i s block-diagonal . Th e
estimate d j o f the coefficien t on th e (differenced ) stationary ter m ha s an
asymptotically norma l distributio n wit h mea n 0 an d varianc e give n b y
Vf]1. Th e margina l distribution o f o 2) however , i s no t normal ; becaus e
F23 i s no t equa l t o zero , Z 2 ,t an d Z 3 j r ar e asymptoticall y correlated ,
and sinc e Z^ t ha s a Wiener distribution , so does the coefficien t o n Z 2:t .
If a n intercep t i s no t include d i n th e regression , w e hav e a 2 x 2
block-diagonal V matrix . Th e estimate d coefficien t o j stil l ha s a n
asymptotically norma l distribution , wit h d^ convergin g to it s probability
limit a t rat e T 1/2, whil e S 3 has a Wiene r distributio n wit h convergence
at rat e T . An y join t tes t involvin g di an d 6 3 wil l als o hav e a
non-standard distribution.
The analog y with the Stock-Wes t exampl e is direct. I n (27 ) we ha d a
series o f term s integrate d o f orde r zero . Th e coefficien t estimate s o n al l
these stationar y term s were jointly and individuall y asymptotically normally distributed . Th e join t distributio n o f 0 i n (27) , wit h an y o f th e 77, ,
was o f cours e non-standard . Thi s observatio n applie s equall y well here .
There is , however , a n importan t differenc e betwee n th e Stock-Wes t
12
Thi s transformatio n i s no t unique , an d on e coul d imagin e choosin g others ; however ,
(39) ca n be rewritte n a s x, = (f) l + /3 2)*,_i - /3 2(*,-.i ~ x t-2> + 1t> because j8 0 = 0 under th e
null, an d thi s suggest s th e decompositio n give n b y (40) . I t ha s th e advantag e o f makin g 6 l
(= /32) th e coefficien t o f a non-integrate d rando m variable , sinc e x , i s a n integrate d
series.
187
(42)
(43)
188 Regressio
regressors i n combinatio n wit h an y o f th e origina l stationar y regres sors wil l also have a standard distributio n asymptotically .
(c) Th e F-statistic s o f join t significanc e o f an y subse t o f th e se t o f
non-stationary regressor s hav e non-standar d distributions .
Moreover, a n F-statisti c tha t use s an y stationar y regressors i n
combination wit h an y non-stationar y regressor s wil l hav e a non standard distribution .
Point (a ) i s obtaine d fro m th e propert y o f th e non-stationar y regres sors formin g a co-integrate d set ; a s in Sectio n 6.2. 3 above, bot h 6 and /3
can b e writte n a s coefficient s o n mean-zer o stationar y variable s (wit h
(46) givin g on e suc h re-parameterizatio n fo r /?) . Th e nex t exampl e
reconsiders thi s poin t i n th e contex t o f modellin g practice . Poin t (b) i s
not surprisin g becaus e th e F-statistic s considere d us e onl y stationar y
regressors. Th e fac t tha t som e o f thes e stationar y regressor s ma y b e
re-parameterizations o f som e o r al l of the origina l non-stationary regres sors i s an interesting feature .
Point (c ) i s surprising in two respects. Conside r (44 ) and (46) ; the firs t
surprising featur e i s th e non-standar d behaviou r o f th e F-statisti c an d
the secon d i s that , whil e th e f-rati o o f th e coefficien t o f c t-\ ha s a
standard distributio n unde r parameterizatio n (45) , unde r th e linea r
re-parameterization give n b y (46 ) th e t -ratio ha s a Wiene r distribution .
Both result s follo w fro m th e asymptoti c singularit y o f a particula r
variance-covariance matrix. 15
Consider y i i n (46) , whic h tend s t o a non-degenerat e distributio n a t
rate T ; T l/2<fi therefor e ha s a degenerat e distribution , an d T 1/2}>2 i s
asymptotically normally distributed. Thus ,
and s o
190 Regressio
where
y* = the permanen t componen t o f disposabl e incom e whic h i s as sumed t o follo w a random wal k
ct = consumption
yst = transitory componen t o f disposabl e incom e whic h is a stationary
innovation proces s
p, = price leve l in period t.
The innovation processes u, and v t ar e uncorrelated .
Stock an d Watso n relat e th e tal e o f two econometricians tryin g to tes t
versions o f Friedman' s permanen t incom e hypothesis . Th e misguide d
econometrician, unawar e o f o r choosin g t o ignor e th e order s o f integration o f the series , estimate s the followin g regressions :
c, = <x\ + Pipt (t
ct = a 2 + $2* (t
Ac, = a 3 + !3 3Ay, (t
The firs t regressio n i s a spuriou s regressio n o f th e classica l Granger Newbold kind ; c, an d p, ar e unrelate d rando m walks , an d th e eco nometrician's findin g o f a larg e ^-statisti c fo r j8 l5 thereb y leadin g hi m t o
conclude i n favour of money illusion, 16 i s a spurious one .
The secon d regressio n i s als o spuriou s sinc e i t attempt s t o explai n a
random wal k (or, i n other words , a stochastically trending variable) b y a
deterministic trend . Nelso n an d Kan g (1981 ) pointe d ou t th e danger s of
running regression s whic h attemp t t o de-tren d stochasticall y trendin g
data i n th e vai n hop e o f achievin g stationarity aroun d a trend . I n bot h
cases th e problem s wit h th e inference s aris e becaus e th e regression s
involve variables tha t ar e no t co-integrate d (se e Chapte r 3) .
The thir d equatio n appear s t o b e correctl y specifie d bu t nevertheles s
leads t o downwardl y biased estimate s o f th e coefficien t for th e margina l
propensity t o consum e becaus e disposabl e incom e measure s th e chang e
in permanen t incom e wit h error , sinc e i t include s th e chang e i n
transitory incom e a s well . Th e fina l regressio n i s wha t w e calle d a n
'unbalanced regression ' a s i t trie s t o explai n a variabl e integrate d o f
order zer o b y a variabl e integrate d o f orde r 1 . Th e serie s o f paper s
noted abov e (Manki w an d Shapir o 1985 , 1986 ; Banerje e an d Dolad o
1988; Galbrait h e t al. 1987 ) conside r th e exten t t o whic h th e f -statistics
in suc h case s ar e biase d awa y fro m zero , leadin g t o misleadin g infer ences abou t th e significanc e of coefficients.
Stock an d Watso n compar e th e predicamen t o f thi s econometricia n
with econometricia n B , say , wh o look s a t th e result s o f th e followin g
alternative regressions :
192 Regressio
194 Regressio
where th e standar d error s o f coefficien t estimate s ar e show n i n paren theses, o i s th e equatio n standar d error , an d S C i s th e Schwar z
criterion. (Smalle r value s on balanc e produc e preferabl e models. ) Sinc e
the mea n o f Y i s 4701.0 , th e a a s a percentag e o f Y i s 3. 1 pe r cent .
However, th e coefficient s ar e no t constan t ove r th e sampl e period , a s
shown i n Fig . 6.1 fo r th e intercept , an d Fig . 6.2 fo r th e one-ste p
residuals an d o . (Se e Hendr y (1989 ) fo r details.) 21 Th e intercep t trend s
upwards, an d o increase s ove r time , eve n ignorin g the larg e shoc k i n
1919-20. O n an y constancy test, th e mode l i s rejected a t fa r beyon d th e
1 per cen t leve l (e.g. tha t of Hansen 1992) .
Next w e mode l growt h i n logs . A s before , on e lagge d differenc e
removed residua l seria l correlation, giving
21
Recursiv e estimatio n involve s estimatin g a n equatio n ove r successivel y large r sub samples, startin g fro m a minimu m sub-sampl e an d extendin g t o th e ful l sample . Paramete r
instability ma y b e tracke d b y lookin g a t th e behaviou r o f th e estimate d coefficients , a s
sample siz e i s increased , t o se e whethe r the y fluctuat e significantl y o r remai n stable .
Recursive Cho w (1960 ) test s ma y b e compute d i n a t leas t tw o ways . Th e firs t involve s
estimating th e equatio n from , say , t = 1 to ( = 7\ , wher e T l i s greater tha n th e minimu m
sample size , an d the n fro m t = I t o t = T t + 1. The one-step-ahea d Cho w tes t is based on
a compariso n o f th e residua l varianc e o f th e tw o estimate d equation s an d i s a n F-tes t
under th e nul l o f paramete r constancy . A secon d tes t i s give n b y estimatin g th e equatio n
from, say , t = 1 to ( = T } an d comparin g th e residua l varianc e o f this regressio n wit h tha t
of th e equatio n estimate d ove r th e ful l sample . A sequenc e o f thes e Cho w test s i s built u p
by augmentin g th e sub-sampl e siz e b y on e a t eac h step , e.g . T 1 + 1 t o 7 \ + 2, an d
195
comparing th e residua l varianc e o f eac h o f thes e equation s wit h th e ful l sampl e residua l
variance. Alternatively , th e sequenc e o f one-ste p residual s (o r forecas t errors ) ca n b e
examined relative to the residua l variance a t eac h sampl e size.
196
0 = exp (2v0 ){l - 2ex p [-(A + a2)] + exp [-(2A + a2)]}; and Y0 = exp (y0 )
is the startin g observatio n postulate d fo r the mode l i n levels.
Thus, i f th e logarithmi c mode l wer e true , th e mode l i n level s woul d
have bot h a drift , (5exp(Af) , an d variance , 0exp[(2A + cr 2 )/], exponen tially increasin g with time. Further , i n the regressio n
with A =j u + & 2/2, wher e f t an d a 2 ar e obtaine d fro m estimatin g (55) ,
1980
198 Regressio
22
T o se e this , simpl y substitut e A r,_1 fo r Z, . Th e instabilit y o f th e rando m wal k mode l
in level s mad e a forma l tes t i n th e level s > logarithms directio n unnecessar y i n th e
Ermini-Hendry discussion , althoug h i n principle suc h a test coul d be carrie d out .
199
or
200 Regressio
(A6)
Next, usin g a componen t b y componen t analysi s simila r t o tha t i n
Chapter 3 , fro m (A5) :
201
202
(A10)
By direct calculatio n fro m (Al ) however ,
(A13 )
Co-integration in Individua l
Equations
We firs t examin e method s o f testin g fo r co-integratio n vi a stati c
regressions, an d provid e simulatio n estimate s o f th e uppe r percen tage point s o f th e distribution s o f statistic s use d i n th e tests . Next ,
we look a t th e propertie s o f the estimator s derive d fro m suc h stati c
regressions. I n particular , w e focu s o n th e finite-sampl e biase s i n
the estimate s o f co-integratin g vector s an d th e power s o f test s t o
detect co-integration . Finally , w e conside r modifie d estimator s an d
dynamic models . I n Chapte r 8 , system s method s o f estimatin g
co-integrating relation s wil l be considered .
The previou s chapte r focuse d o n th e propertie s o f co-integrate d pro cesses an d th e implication s o f modellin g wit h co-integrate d variables .
We hav e discusse d th e 'super-consistency ' o f th e coefficien t estimate s i n
the static o r co-integratin g regression , balance d an d unbalance d regres sions, an d th e distribution s o f th e statistic s commonl y use d t o tes t fo r
the significanc e of regression coefficients .
The tw o issues o f being abl e t o tes t fo r th e existenc e o f an equilibriu m
relationship amon g variable s an d t o accuratel y estimat e suc h a relation ship ar e complementary . Indeed , a s demonstrate d i n discussin g spuriou s
regressions i n Chapte r 3 , stati c regression s amon g integrate d serie s ar e
meaningful i f an d onl y if they involve co-integrate d variables . Thus , i t i s
of interes t t o discover , first , ho w wel l th e mos t frequentl y use d test s of
co-integration perform , an d second , ho w accuratel y th e correspondin g
equilibrium relationship i s estimated.
The objectiv e o f thi s chapte r i s t o develo p test s applicabl e t o singl e
equations whic h ma y b e use d t o detec t a long-ter m relationshi p o f th e
form discusse d an d exploite d i n earlie r chapters . W e als o attemp t t o
formulate som e recommendation s fo r efficien t estimatio n o f co-integrat ing parameter s an d testin g fo r co-integratio n i n finit e samples . I t wil l
become clea r fro m th e discussio n that th e asymptotic propertie s o f static
regression estimator s ar e ofte n rathe r differen t fro m thei r behaviou r i n
empirically relevan t sampl e sizes . Further , lac k o f wea k exogeneit y du e
to co-integratin g vector s enterin g severa l equation s als o alter s finit e
sample behaviour . I t therefor e become s important , i n th e fac e o f dat a
206 Co-integratio
n i n Individua l Equation s
*t = (yt,z ty.
The modification s are necessar y because, whil e the test s for uni t root s
discussed i n Chapte r 4 us e th e origina l series , sa y {w t}, th e co-integra tion test s ar e base d o n th e estimated, o r derived, residual series ,
208 Co-integratio
n i n Individua l Equations
and
where y, an d z t satisf y (multivariate ) unit-roo t processes . Th e asymp totic distribution s o f a numbe r o f residual-base d test s ar e discussed ,
from whic h we wil l conside r fiv e (thi s analysi s is of cours e relate d t o th e
CRDW
ADF(l)
ADF(4)
50
100
200
0.72
0.38
0.20
-3.43
-3.38
-3.37
-3.29
-3.17
-3.25
50
100
200
0.89
0.48
0.25
-3.82
-3.76
-3.74
-3.75
-3.62
-3.78
50
100
200
1.05
0.58
0.30
-4.18
-4.12
-4.11
-3.98
-4.02
-4.13
50
100
200
1.19
0.68
0.35
-4.51
-4.48
-4.42
-4.15
-4.36
-4.43
where
210 Co-integratio
n i n Individual Equation s
212 Co-integratio
n i n Individua l Equations
Model
Point (% )
000
SE
0i
02
No constant ,
no tren d
Constant,
no trend
Constant
+ tren d
Constant,
no tren d
Constant
+ tren d
Constant,
no trend
Constant
+ tren d
Constant,
no tren d
Constant
+ tren d
Constant,
no trend
Constant
+ tren d
Constant,
no tren d
Constant
+ tren d
-2.5658
-1.9393
-1.6156
-3.4336
-2.8621
-2.5671
-3.9638
-3.4126
-3.1279
-3.9001
-3.3377
-3.0462
-4.3266
-3.7809
-3.4959
-4.2981
-3.7429
-3.4518
-4.6676
-4.1193
-3.8344
-4.6493
-4.1000
-3.8110
-4.9695
-4.4294
-4.1474
-4.9587
-4.4185
-4.1327
-5.2497
-4.7154
-4.4345
-5.2400
-4.7048
-4.4242
-5.5127
-4.9767
-4.6999
(0.0023)
(0.0008)
(0.0007)
(0.0024)
(0.0011)
(0.0009)
(0.0019)
(0.0012)
(0.0009)
(0.0022)
(0.0012)
(0.0009)
(0.0022)
(0.0013)
(0.0009)
(0.0023)
(0.0012)
(0.0010)
(0.0022)
(0.0011)
(0.0009)
(0.0023)
(0.0012)
(0.0009)
(0.0021)
(0.0012)
(0.0010)
(0.0026)
(0.0013)
(0.0009)
(0.0024)
(0.0013)
(0.0010)
(0.0029)
(0.0018)
(0.0010)
(0.0033)
(0.0017)
(0.0011)
-1.960
-0.398
-0.181
-5.999
-2.738
-1.438
-8.353
-4.039
-2.418
-10.534
-5.967
-4.069
-15.531
-9.421
-7.203
-13.790
-8.352
-6.241
-18.492
-12.024
-9.188
-17.188
-10.745
-8.317
-22.504
-14.501
-11.165
-22.140
-13.641
-10.638
-26.606
-17.432
-13.654
-26.278
-17.120
-13.347
-30.735
-20.883
-16.445
-10.04
1
5
10
1
5
10
1
5
10
1
5
10
1
5
10
1
5
10
1
5
10
1
5
10
1
5
10
1
5
10
1
5
10
1
5
10
1
5
10
0.0
0.0
-29.25
-8.36
-4.48
-47.44
-17.83
-7.58
-30.03
-8.98
-5.73
-34.03
-15.06
-4.01
-46.37
-13.41
-2.79
-49.35
-13.13
-4.85
-59.20
-21.57
-5.19
-50.22
-19.54
-9.88
-37.29
-21.16
-5.48
-49.56
-16.50
-5.77
-41.65
-11.17
0.0
-52.50
-9.05
0.0
Source: MacKinno n (1991) . W e ar e gratefu l t o Jame s MacKinno n fo r permis sion t o reproduce thi s table.
214 Co-integratio
n i n Individua l Equations
216
Sample size
Sample size
Fio7.1(&). Constan t i n model, estimate d bias v . sampl e size, s = 16
3.5. A s anothe r example , w e se e i n Fig . 7.2(a ) (01/02 = 4), wit h
p2 = 0.6, th e biase s a t th e sam e se t o f sampl e size s ar e 0.017 , 0.010 ,
0.005, 0.0026. 6 Her e a n eightfol d increas e i n sampl e siz e reduce s th e
6
Thes e number s ar e take n fro m th e experimenta l outpu t rathe r tha n rea d fro m th e
figures. Th e standar d erro r o f th e smalles t o f these number s i s roughly 5 x 10~ 5.
217
Sample size
Sample siz e
Sample si/.e
Sample size
219
Sample size
Sample size
220 Co-integratio
n i n Individua l Equations
that b y equatio n (10) , u t i s th e discrepanc y fro m thi s long-ru n relation ship. Smalle r value s o f p 2 an d smalle r value s o f o 2 (large r value s o f
Oi/o2) mak e the varianc e of u t relativel y small , and so we obtai n smalle r
biases a s p 2 fall s o r a s 0i/o 2 rises .
The fac t tha t thes e biase s d o disappea r les s quickl y tha n T, an d ma y
remain substantia l fo r sampl e size s larg e relativ e t o man y foun d i n
economics, suggest s tha t th e result s fro m pur e stati c model s mus t b e
treated wit h caution . W e wil l late r examin e way s i n whic h w e ca n
improve upo n simpl e stati c estimatio n eithe r b y includin g dynami c
elements, adjustin g th e result s o f th e stati c model , o r estimatin g a
system o f equations .
Finally, th e biase s ar e strongl y positivel y correlate d wit h ( 1 R2),
which indicate s tha t co-integratin g regression s wit h value s o f R 2 wel l
below unit y should b e viewe d with caution. 7 However , i n th e contex t of
multivariate regressions , a high value of R 2 i s not sufficien t t o guarante e
that th e biase s ar e small . Thi s i s because th e R 2 o f a n equatio n canno t
fall whe n a n additiona l variabl e i s adde d t o it . Thus , th e inferenc e tha t
high value s o f the R 2 impl y low biases, especially wher e th e forme r may
have bee n achieve d b y a n a d ho c additio n o f regressors , i s no t valid .
Banerjee e t al. (1986 ) explor e th e relationshi p betwee n bia s an d
(1 - R 2} i n more detail .
It i s usefu l t o conside r a n informa l explanatio n fo r th e existenc e o f
biases i n stati c regressions . Th e effec t o f usin g stati c regression s t o
estimate th e co-integratin g slope / ? is to allo w th e residua l u t t o captur e
all th e dynami c adjustmen t terms . Accordin g t o th e super-consistenc y
theorem, thi s i s certainl y permissibl e asymptotically. I t i s importan t t o
emphasize tha t th e proble m w e ar e discussin g her e i s strictl y a finite sample one ; the omissio n o f th e dynamic s ma y b e justifie d asymptotic ally b y observin g that , a s the y ar e o f a lowe r orde r o f magnitud e tha n
the non-stationar y term s i n th e regression , the y ma y b e ignore d i n th e
limit. However , th e omitte d dynamics , despit e bein g o f a lower orde r o f
magnitude, ca n matte r considerabl y i n determinin g biase s eve n i n fairl y
large bu t finit e samples. 8 Henc e i t seem s appropriat e t o pa y attentio n t o
modelling th e omitte d terms .
The dynami c terms ca n al l be parameterize d i n term s o f 1(0 ) serie s o f
the for m A >>,,, Az,_ ; , an d ( y yz)t-k wher e th e value s o f i , j, and k
1
W e ar e gratefu l t o To m Rothenber g fo r pointin g out tha t R 2 i s a rando m variabl e in
the presen t context . However, i t remain s a usefu l descriptiv e statistic.
8
Th e proble m o f finit e sampl e biase s wa s als o demonstrate d b y Hendr y an d Neal e
(1987). Usin g recursiv e procedure s fo r OL S estimation , the y estimate d a bivariat e stati c
regression fo r sampl e size s rangin g fro m 4 0 t o 200 , considering th e bia s o f th e coefficien t
estimate fo r eac h sampl e size . Th e result s indicate d that, eve n fo r sampl e size s o f 200, the
long-run coefficien t fro m th e stati c regressio n wa s approximatel y 0. 7 whil e th e tru e
long-run coefficien t wa s 1.0 . Convergenc e t o th e tru e valu e wa s no t nearl y a s fas t i n
practice a s T~ ! whic h dominate s for sufficientl y larg e T: se e (18 ) below.
222 Co-integratio
n i n Individua l Equations
The extr a lagge d variable , z t-\, i s include d t o avoi d imposin g homo geneity (se e Chapte r 2) , a s i t woul d b e unrealisti c t o assum e tha t th e
investigator know s th e precis e for m o f th e data-generatio n process . Th e
co-integrating coefficien t i s estimate d b y computin g th e expressio n
1 - d/c: se e Sect . 2.4 . Th e stati c regressio n give n b y (16 ) i s als o
estimated.
The stron g exogeneit y propert y require d o f z t i s guaranteed , i n th e
design o f th e experiment , b y drawin g e lt an d e 2t fro m uncorrelate d
pseudo-normal distributions . Th e value s o f y , ( i = 1, . . ., 3 ) ar e varie d
as i n Tabl e 7.3 , while ensurin g tha t long-ru n homogeneit y i s preserved .
The sampl e size s an d th e rati o o f the standar d deviation s o f e lr an d e 2t
are als o varied , t o giv e a se t o f 9 0 experiments . Th e simulation s ar e al l
conducted with 5000 replications .
The purpos e o f th e firs t par t o f thi s exercise i s to compar e th e biase s
in th e estimate s o f th e co-integratin g paramete r obtaine d fro m dynami c
regression wit h thos e obtaine d fro m th e stati c regression . (Th e tru e
value o f th e co-integratin g paramete r i s 1. ) Som e o f th e result s fo r
different configuration s o f th e y , parameter s an d standard-deviatio n
ratios ar e give n i n Tabl e 7.3 . We repor t th e estimate d biases , fo r fou r
different sampl e sizes , i n th e stati c model . Th e correspondin g estimate d
biases fro m th e dynami c regressio n (wher e th e co-integratin g paramete r
is calculated a s (1 d/c)) ar e i n almost al l cases so small a s to b e withi n
2 Monte Carl o standar d error s o f zero an d s o ar e no t reported . W e wil l
return t o th e compariso n o f these estimator s (stati c an d dynamic ) below ;
for th e tim e being , th e noteworth y point i s simply that substantia l biases
remain i n stati c estimate s fo r paramete r combination s a t whic h th e
biases i n dynami c estimate s ar e zero , o r ver y clos e t o zero , sinc e th e
dynamic model ha s been specifie d s o a s to mak e A close t o zero .
While th e dynami c estimate s contai n negligibl e biase s i n thes e ex amples, Z t is strongly exogenou s fo r th e paramete r o f interest . Whil e i t
is fairl y straightforwar d t o exten d thi s specificatio n t o includ e weakl y
exogenous z t , th e usefulnes s o f estimate s fro m dynami c single equation s
is reduce d substantiall y i f th e regressor s ar e no t weakl y exogenous . I t
also become s difficul t t o mak e unambiguou s comparison s betwee n
= 0 ,-5,
= 0 ,,5,
= 0 , 1 ,
= 0 .1,
0 10
0 20
0 40
-0.,39
-0.25
-0.15
-0.07
-0.,04
-0,.32
-0.22
-0.14
-0.08
-0..04
-0,,23
-0.13
-0.07
-0.03
-0,,02
-0.,21
-0.12
-0.06
-0.03
-0,,02
224 Co-integratio
n i n Individua l Equation s
A
B
C
D
Propert
1
1
1.0
1.0
1(0)
1(0)
1(1)
1(1)
y o f {z t} Propert
y o f {y,}
1(0)
nearly 1(1)
nearly 1(2)
1(1)
226
^ ^ ^^^ ^
____^
<M<QU<;m<:mQU<;<Q
^
.&1
Vj
1
ON
^
o f) rT-H ,* ON ON CM o *
l ^- - 0 0 , . O 4 O O , - N O O C O 0 0 ^ f I O T- H
ro^qooo^Hi/-)Oi-<mo<^^HOO
~I
O r- H
^ i 1' : -? <:: ^
c cr c
O^ <N O, 00 O, O^ <d, K O^ K O, O^ O^ O^ O_
^
pj
H
_, *g
^ ^
l f
'
^ ^^~, ^^ ^^ ^- v
^^^ ^^ ^^
S'SSS-S'SS'SSS'S-S-S-SS
^o-*^tNtNONvcmcocoror-~oom
o)O<-n(siinomO'-H(NrMOooom
ocNooodoo'OO^oOv-idoo
i g i
?''1 0 9
2
^
O
*-
S
~'
, O O i - H O C N C N O O ^ O f N C N O i - H
'
II
c/)
^ ^~~
\.
I!
~^.
ty)
CH
^O
Z
^
[/I "
CD <<S5
-<-*
Bj I
a S
O
|-2 S
e s
I*
>,
-a o
"O
ci'~ H
S
ra
03 c/r
"
-
S S 52
oo O
^
_I
QJ
|u- i
C\3 fl
.^4
_) I
1 7 1 1
0 ) 0 , ? 'Q - ^ o o r - O N O ( N ^ H a \ o v o o ^ o r s i v o u - i - *
Sr n O ^ ^ CT\or~a\moooNON^oos(Noor^(N i
5^o
I
$..-*-
~~~"
n. ^ ^ II
O O O O O O O O O O O O ^ t ^ ^ J ^ " ..
. C-IC3INO<NlO^HO^OvH
O
q - i C S I m ^<
NI
N(
N^
H rl ii
S,
-S-
o.
Q.
^^nin
oocooooooooqooooooqo
dddoT-H^dddd^H^ddT-H
ddddddddddddddd
ON ON ON ON ON ON ON ON ON CTN ON ON ON ON ON
o o p o o o o o o o o o *_**
ddddddddddddddd
ddddddddddddddd
^inw-jinir)'n>o'Oioi^i^i^;
to +
<O
-i- ^ N ~O i o u ~ i i r ^ i n t / ~ ) i o i n i r ) i r ^ i o i o i r ) i n r i i r
)
c N
vo
cJddddcJddddcJdddc)
S ^-^ 'C
o.
I I I I I I
^- .. d>
S ^ u
o-
CO
o O S
r>
^S
rj
-a
g2
"3
G
o
on
'g
<u ^ ?*%
y o
< "
O
^
o
in
W1
?*^
oxS
fi
t^ 0
r^
"
ui
a
o
<l=a
I
'Q.
T-H
o^
O
's O^
O^
^O
rH r-- Tj^-v CN
O*l CO O
V^
OO O
OO ON OO O> ON O> in
'v ^O
^O
(N rO
ONON^OfN
rH f* rH ^H
^H 00 C-4
<<<<!<<pqQQQU<!<<ffl<;<<:mQQQU
CN U"j ^
CO
oo
f^^^l
CN * C4
o o "^
f^^1^
OS
^ ^ '-J
^^^
CN (N i-H
r*]<Noooc>o^Hv.oooorr><-Hr-HTj-inoO'Hioooc>r-:
S'S'^i^ SS-S^ S-S2^2-S-SS^ SSS-^ SS-S-S-
I
I
I ^ I
I ^ I
II
I l l ^ l l l
^HO<N-*OOO^t-OOO^HOOO(NOOOir)OOOT-(
cN^)-^-^oost--^-moooN<^r--c^o^ooX:5sar'i/Tsi/-r^sCHs\osc7~oo"
S2'2'2'2'2'2'2'S-2'2'S'S-2'2'SSS-2'S2'2^S-2-
^ H O O O O O ^ H O O O O C ^ O O O O O O v H O O O C N f N
COOOCOOOOONrHCNlrOV^<NT-HQ\r--CNO^OO>^HCNCO^O^-^H
T-H
i i i
iI
in
i i i i i v
'
<I
rH
"n
i r~--o s i i
i i i
T-H
i i
C Noo
T-H
iifNOO-^HCNlOOTic<)OOT<c^)ioOTic^iOOTirO^nO
iI O
T-HO
'^O
if*/~)
'I *O
Tiv/^
vn
| i<Noocnvoa\oicN
cOcNi-HON<N'-H^OOO^toO'-H>OOO^-OONI--'^-a\CNUlON^H
,,
S~*^
55.
^H
Sj"
<CQ.
i i i
' ima\a\i<sit--a-i
<S2.
II
-
11
II
cT I
r- ir
^ S
II ^
^Q.
<
vc
II
iO
N ^
oooooooooooov")inini/^v")v/-)ir)i/^i/"}ir>i/^i/"i
O~
ovoONO\ov~io^ONOw^oaNOino%ooi/~jaN^o>oONO
ooooooocJoooooocsooooooocJo
N ^
^
^-
O O O O O O O O O O O O O O O C 5 O O O O O O O O O O O O O O O C 5
OOOOOOOOT-H^Hi-H^-5ooOOOOOC5T-HT-H^H^-H
\O
. O C 5 O O O C 5 O O O O O O O O O O O O O O O O O
,0
Q,
>
T3
(0
.s
'^
OJ
o>
6O
PH
CQ
t-l
CO
>
M
4)
"H3
T3
O
"Q-
" CT
o
'
ca
a
>>
Q
^
&
228
u?
C/3
in
'-H O
OO O-^ ^*~~s
^~
MD CN) T1 ON
^) f^ i^j- CO
Tf ii |
1
m s*-'
CO CN
O O
00 -<f
O O CO in oo
r- m co
o o o
1/1 CN rH
in rH in in NO in
1
NO
f? US'
rH
m CN o o o o
Ti CN ON ,,CO ON ^^oo"
C"-~ iI Tf rH "si" O O 00s O CN CO
CO <
<! <!
,-^G
NO
&o
NO
rH
'
^^
ON ON ON ON O
CN
r~ ON ON oo o
1
1
1 1
xV
yV
t- ON
in ON <o
ON
CN
G"^
oo ^t
^
CN CO
ON CN
m co
O O
in o ON in o
^ co ^|- Tf CO
88 8 8 8 88
ooo oo oo
ooo oo
CO ^
ON
in
in
in
o o o CN O O CN OO CN
10 in in
m rH m in rH
CN CN
rH
O O o o o
in in in ON i/}O
"^ CN
rH O
CO O ,, /^~s W)
^^ CO "!:T
in in o
CN OO ON rH O ^O fO T4
1 1
in 1/1
rH^f^
1/1 ON
rH
CN
O
1
in ONON m m m m in in 1/1 m m
o o o o o o o o oo o o
O O o o O O O
ON
ON ON
O O
rH
^^ xv
1 1
o
o
<o
CD
O
.~
, I
1
rH
&
O
Q,
<^
cd
.-i
'"'
IB
rH
c/l 'O
/* !
>H
3 g> en
^ aj
2
(^
0 -^
4H
u .an
M
C3
^ 'fe
Q^ &
CD
TO 43
a CD
M) S1
II
.g
3
CD
^O
CS
00
PH
CD
O
Q
^
X
^H
CD
>H
CTJ
Q. CD
cu 1)
^
CTj TJ
CD "
CD
-^
-g
^
Q
SH
CQ
.Si
cd
CD
^ "
S
^J
+j CD
SP i3
^
'||
<D
. ^
-D
HH
F-3
1 1
15
o "
C3 "^ S --^
o-5 o *
C -3 C c3
'** m .52 D
23
& "-rj CQ
CD
K
.3 uJ O
cu 73
o 8 T)CD CDrt
03 3w 'CN T30
CJ J^H
18
SH
S3 g>
^
0
-1 S
*" "CD 0 T3
?%
9 2
CD
J2 '55
x 2
* .
CD "*'
g CD
<D
04
rH
rH
O *
O O rH CO
O ON
m o
^_^^
^__^
CO ^-v
1 1
00 ON
1
1
co r-~ co o
rH
O ON
rH O
IT)
T)
rH ON rH O rH in
rH ^f
1 *> CN 1/1
i 1
1
~ ^
CO
in r-~
oo co
7
CN
rH
O O
ON
o o o oo 00 00 o o o o o O 00 oo oo o o o
rH O O o o O O rH T1 1I
o o o o o o
in ONON in
o o o o
o o
in m
o o o o
1 1 1 1
in m
m m o in in o
CO
rH
^vG
0s!--
QQ.
<CQ.
j
QQ.
T1
"*r1
["T']
**~s
--
^i.
o'
o
c^
II
o.
a
II
<
ii
n
^
"2
.s
CD
*5
CD
O.
CO
tn
ta "O
"c
3
^^ "~^I
c o
8 ~
i-^
^'ea3
^ a
w >->
<j .^
H -S
230 Co-integratio
n i n Individual Equations
)
)
232 Co-integratio
n i n Individua l Equations
(29)
where r\ = (|J#e (r)dW(r)) (^K E(r)2dr)~l. Whe n e = 0, w e reproduc e
the distributio n unde r H 0. Otherwise , fo r e < 0 , th e distributio n i s
shifted t o th e lef t b y e(\\K (r)2 dr) 1/2 . Whe n T = 100, e = - 1 implie s
that p= 0.99, an d e = -5 implie s tha t p = 0.95; a s e , th e powe r
tends t o 1 .
Kremers e t al. (1992 ) argu e tha t simila r consideration s sho w tha t th e
non-centrality paramete r o f th e ECM-base d tes t fo r co-integratio n i s
larger tha n tha t o f th e non-parametri c statistic s discusse d i n Chapte r 4 .
Their Mont e Carl o result s support thes e asymptoti c results .
Return no w t o th e Mont e Carl o experimen t give n b y equation s
(14fl)-(146). On e appealin g tes t fo r co-integratio n tha t w e hav e men tioned consist s i n usin g th e mode l (15a) , where , unde r th e nul l o f n o
co-integration, j l = 1 so tha t th e secon d coefficien t i s equa l t o zero . A
f-test fo r thi s conditio n i s therefor e a tes t fo r co-integration . Whil e w e
would expec t th e distributio n o f thi s tes t statisti c t o b e non-standard , i t
is a straightforwar d tes t an d woul d therefor e b e especially usefu l i f it s
power wer e high . I n particular , fo r strongl y exogenou s regressor s i t i s
similar (se e Kivie t an d Phillip s 1992) .
We examin e th e tes t wit h a small Monte Carl o experiment , comparin g
its powe r wit h tha t o f th e AD F tes t base d o n a static mode l t o estimat e
the co-integratin g parameter , i n th e DG P give n b y (I4a)-(l4b). Th e
first tes t i s th e AD F tes t wit h on e lag , compute d fro m th e residual s of
the static regressio n (16) . Th e secon d tes t i s base d o n th e ^-statisti c fo r
c i n (19) . A s note d earlier , i f the nul l of no co-integratio n i s true, c = 0 .
Under th e nul l (i.e . y^ = 1 , y 2 = y 3 = 0, o l = o 2 = 1 in (14a)-(14&)) ,
c=0 ha s a Wiener distribution . The critica l values of this distribution an d
the AD F wer e compute d b y simulatin g th e nul l mode l fo r 500 0
replications usin g PC-NAIV E (Hendry , Neale , an d Ericsso n 1990) . Thes e
critical value s wer e the n use d fo r computin g th e tes t power , an d ar e
shown i n Tabl e 7. 6 fo r regression s lik e (19 ) wit h a n intercept . (Th e
population constan t i s zero. ) Th e sam e critica l value s resul t fo r
72 + 7 s = 0 when thes e parameter s ar e individuall y non-zero, s o the tes t
25
50
100
0.10
-2.99
-2.95
-2.93
0.05
-3.42
-3.33
-3.28
s of ADF(l)
0.01
-4.22
-4.06
-3.95
0.10
-3.15
-3.10
-3.09
0.05
-3.51
-3.41
-3.39
0.01
-4.30
-4.08
-4.00
234 Co-integratio
n i n Individua l Equations
WADF
();l/i =
T=
25
50
100
0.13/0.13
0.21/0.17
0.44/0.31
0.05
WADF
0.01
0.06/0.06
0.10/0.10
0.26/0.20
0.01/0.01
0.02/0.02
0.07/0.05
0.06/0.05
0.10/0.09
0.30/0.19
0.01/0.01
0.02/0.02
0.08/0.04
0.07/0.05
0.12/0.07
0.40/0.14
0.02/0.01
0.03/0.01
0.13/0.03
0.45/0.20
0.97/0.72
1.00/1.00
0.16/0.05
0.78/0.34
1.00/0.97
0.66/0.18
1.00/0.67
1.00/1.00
0.29/0.04
0.94/0.28
1.00/0.96
0.87/0.12
1.00/0.60
1.00/1.00
0.64/0.03
1.00/0.22
1.00/0.94
S = CTi/0-2.
(a)
-1.15
-1.19
(*)
-1.15
-1.28
(c)
-1.15
-1.52
(d)
-2.89
-3.25
()
-2.89
-3.88
(/)
-2.89
-5.32
236 Co-integratio
n i n Individua l Equations
237
price homogeneit y i s imposed an d Ap adde d a s a regressor, correspond ing t o allowin g m an d p t o be 1(2), wit h ( m - p ) an d Ap bein g 1(1). In
that last case , R 2 fo r real mone y is equal to onl y 0.68.
We assum e no w that Ap, x S5, an d R n ar e weakl y exogenou s fo r th e
parameters i n th e conditiona l mone y deman d model . Th e outcom e o f
estimating a dynami c equatio n i n th e level s o f th e variable s wit h fiv e
lags o n eac h o f m p, Ap , * 85, an d R n (plu s a constant ) b y leas t
squares i s shown in Table 7.8.
TABLE 7.8. Empirica l result s
Variable
Lag
1
0
m p
-1.000
xss
-0.041
0.115
-0.411
0.117
-0.757
0.210
-0.124
0.169
SE
SE
Rn
SE
Ap
SE
CONSTANT
SE
0.
Sum o f
lags
A
0.164
.147
0.549
0.,240
0,,251 0 .152
0,,132
0.,135 0 ,131
0.109
0,,028
0.118
0.087
0,.162
0.293 -0,,067 -0.,240
0,,130
0 .139
0.119
0 .026
0.135
0..139
0..139
-0.361 -0,,122 -0.,046 -0 .084 -0.045 -1.070
0.130
0 .187
0.178
0.,185
0.,176 0 .175
0.069 -1,.102
0.020
0,,307 -0.,412 -0 .329
0 .222
0.255
0,.253
0,,246
0 .246
0.203
- -0.12
4
0 .169
R2 = 0.9966 a
= 0.0130 F(23
, 76) = 975.3 8 D
W = 1.976
SC = -7.85 3 Mea
n = 10.89613 1 S
D = 0.19617 3
Normality % 2(2) = 4.29
AR 1- 5 F[5, 71]
=
0.2 0 ARC H 4 F[4
, 68] = 0.22
Xj F[37,38]
= 0.6
6 RESE
T F[l,75
] = 0.98
COMFACF[15,76] = 3.14
Tests on the significance of each variable
Variable
Ffnum., denom. ]
Value
Probability
Unit-root
Mest
m p
F[5,76]
F[6, 76]
F[6, 76]
F [6, 76]
F[l,76]
340.201
7.801
12.127
6.846
0.536
0.000
0.000
0.000
0.000
0.466
-5.168
6.171
-5.719
-4.963
-0.732
*85
Rn
Ap
CONSTANT
238 Co-integratio
n i n Individual Equations
5
4
3
2
1
0.691
1.615
1.654
1.416
12.967
F [4, 76]
F [4, 76]
F [4, 76]
F [4, 76]
F[4, 76]
y
0.600
0.179
0.170
0.237
0.000
2
5
d.f . Valu
e Incrementa
0.086
0.196
4.176
8.101
47.128
3
3
3
3
3
l x 2 d.f . Valu
0.086
0.110
3.980
3.925
39.028
When th e {u it} ar e autocorrelate d an d intercorrelated , a stati c regres sion o f yit o n y 2(, b y no t usin g an y informatio n abou t th e proces s
generating y 2t, provide s a n estimat e o f y 3 whic h ca n b e quit e severel y
biased eve n i n fairl y larg e samples . Phillip s e t al. therefor e recommen d
full-system maximu m likelihood estimatio n o f co-integrate d systems . A s
an alternativ e t o estimatio n o f th e ful l system , the y propos e correctin g
the single-equatio n estimate s non-parametricall y i n orde r t o obtai n
median-unbiased an d asymptoticall y norma l estimates . Thes e re commended corrections , fo r simultaneit y bia s an d residua l autocorrela tion, us e expression s derive d fro m th e asymptoti c distribution s o f th e
estimators althoug h th e correction s ar e mad e t o estimator s fro m finit e
samples. Phillip s an d Hanse n (1990 ) sho w tha t thes e correction s wor k
effectively i n sampl e size s a s smal l a s 50. 15 Thei r exampl e i s presente d
in Sectio n 7.10. 4 below.
The estimate s obtaine d fro m full y modifie d an d full-informatio n
methods ar e asymptoticall y equivalent . Thi s equivalenc e i s o f interes t
because i t link s th e discussio n wit h a thir d possibl e metho d o f reducin g
finite-sample biases , namely , estimatin g single-equatio n dynamic regres sions. Th e ai m o f th e analysi s i n thi s sectio n i s t o compar e th e
non-parametrically corrected estimate s (whic h ar e als o asymptoticall y
efficient an d median-unbiased ) wit h estimate s obtaine d fro m dynami c
regressions i n eithe r thei r AD L o r EC M forms . Th e for m o f th e
autocorrelation i n th e erro r proces s i n (31 ) an d (32 ) i s crucia l t o thi s
comparison. Fo r som e specification s o f th e erro r process , a dynami c
15
Whil e i t i s possible t o deriv e exac t expression s fo r th e biase s i n finit e sample s t o an y
desired leve l o f accuracy , usin g Edgeworth-typ e expansions , thi s i s a complicate d pro cedure .
240 Co-integratio
n i n Individua l Equation s
where o)
result:
112
242 Co-integratio
n in Individua l Equations
where
The varianc e o f v} t i s give n b y cr n 2 = a\\ O2io22, and r] t i s orthogona l
to 2, as well a s t o th e entir e histor y o f e, given b y (f,_i , r _ 2 > )
18
Estimating th e regressio n (39 ) is asymptotically equivalen t t o maxi mizing th e conditiona l likelihoo d functio n o f (s u, e 12, . . ., CIT), give n
(s2t, t = 1, 2, . . ., T) . Assumin g invertibility of A(L) in (41), we have
which involves
}
(45)
where / ? i s th e estimat e o f th e coefficien t o f y 2t i n (44) . Bv(r) an d
B2(r) compris e a bivariate Brownian motion process with a well-defined
variance-covariance matrix .
The questio n pose d a t th e beginnin g o f thi s sub-sectio n ca n no w b e
answered. Comparin g (37 ) and (45) , the full y modifie d estimato r fi +
and th e dynami c single-equation least-squares estimator ar e equivalen t if
and onl y i f B v(r) = BI ,2(r). Thes e tw o Brownia n motion processe s ar e
not necessaril y equa l t o eac h other . Thi s i s becaus e B v(r) ca n b e
correlated wit h B 2(r), despit e it s constructio n i n (40) . The generatin g
mechanism fo r u 2t ma y therefor e b e informative , and optima l inference
then require s join t estimatio n wit h th e error-correctio n model . Phillip s
(1988) describe s thi s a s a failur e o f wea k exogeneit y or vali d conditioning. If , o n th e othe r hand , B v(r) an d B 2(r) ar e uncorrelate d a t al l
frequencies, th e conditiona l proces s i s completel y informativ e fo r th e
purposes o f estimation o f f t an d th e margina l process generating u2t ma y
be ignored . In suc h a case, B v(r) B\ 2(r).
244 Co-integratio
n i n Individual Equations
Hence
where,
246 Co-integratio
n i n Individua l Equations
Then
where CTH 2= au - o\ 2a22. The expression fo r Sln.2 follow s from appli cation o f th e conditional-expectation s formul a an d fro m inspectio n o f
(53). t], an d u 2, ar e agai n incoherent , an d th e limi t Brownia n motion s
are
247
The Brownia n motion s B^ an d B 2 ar e correlate d an d th e single equation dynami c estimato r an d th e full y modifie d estimato r ar e n o
longer equivalent , unles s $ 21 =0. Fo r th e structur e o f th e correlatio n
between B n an d B 2 (se e Phillips 1988a):
where B^ 2(r) i s a univariat e Brownia n motio n proces s wit h varianc e
given by crn 2 - oli^d^H' 1 an d is independent o f B2(r). Further ,
248 Co-integratio
n i n Individua l Equation s
modelling. I t i s eviden t tha t th e conditio n 621 ^ 0 violate s wea k exo geneity20 a s ma y b e verifie d fro m (57) ; an d onc e again , i t ma y b e see n
that th e issue s o f a full y modifie d estimation an d dynami c specification
are closel y related . Thi s exampl e form s th e basi s fo r th e simulatio n
exercise discusse d i n the fina l sub-section .
7.10.4. Simulation Example (Phillips and Hansen 1990: 116)
The data-generatio n proces s fo r thei r simulation study is given by
021 = -0. 8
i = 0. 0
OLS
D
FM
0-21 = -0. 4
OLS
D
FM
-0.137 (0.125 )
-0.062 (0.106 )
-0.025 (0.127 )
-0.090 (0.089 )
-0.021 (0.066)
-0.028 (0.079 )
-0.055 (0.061 )
-0.003 (0.041 )
-0.025 (0.052 )
-0.067 (0.081)
-0.051 (0.086 )
-0.042 (0.094 )
-0.057 (0.079 )
-0.030 (0.077 )
-0.027 (0.081 )
-0.040 (0.061 )
-0.007 (0.060 )
-0.015 (0.063 )
OLS
D
FM
0-21 = 0. 8
OLS
D
FM
-0.024 (0.040 )
-0.023 (0.046 )
-0.023 (0.048 )
-0.020 (0.046 )
-0.019 (0.053 )
-0.012 (0.052 )
-0.011 (0.050)
-0.009 (0.060 )
0.004 (0.060 )
-0.015 (0.025 )
-0.009 (0.024 )
-0.016 (0.028 )
-0.010 (0.028 )
-0.008 (0.030 )
-0.005 (0.030 )
-0.004 (0.036 )
-0.005 (0.039 )
0.015 (0.043 )
a21 = 0.4
250 Co-integratio
n i n Individua l Equations
251
02i = -0. 8
OLS
D
FM
021 = -0. 4
OLS
D
FM
02! = 0. 4
OLS
D
FM
CT21 = 0. 8
OLS
D
FM
02i = 0. 4
92i = 0. 0
-1.616 (1.268)
-1.259 (2.040 )
-0.388 (1.432 )
-1.240 (1.105 )
-0.563 (1.701 )
-0.449 (1.092 )
-0.930 (1.00 )
-0.003 (1.40)
-0.025 (0.896 )
-1.156 (1.32)
-1.058 (1.69)
-0.729 (1.49 )
-0.986 (1.25)
-0.636 (1.57 )
-0.516 (1.35 )
-0.754 (1.149)
-0.163 (1.388)
-0.335 (1.193)
-0.711 (1.19)
-0.664 (1.29 )
-0.606 (1.26 )
-0.520 (1.21)
-0.478 (1.34 )
-0.267 (1.30 )
-0.267 (1.24 )
-0.213 (1.37)
0.096 (1.36 )
-0.575 (0.955 )
-0.445 (1.15)
-0.519 (0.922 )
-0.302 (0.979 )
-0.339 (1.25 )
-0.102 (0.962 )
-0.098 (1.04 )
-0.184 (1.36 )
0.418 (1.12 )
where
Then
or
252 Co-integratio
n i n Individua l Equation s
253
(A5)
1
254 Co-integratio
n i n Individua l Equation s
and:
(A10)
Following Phillip s an d Durlau f (1986) , consider a genera l 1(1 ) vecto r
process:
and v t i s a weakl y stationar y stochasti c proces s wit h unconditiona l
covariance E(v tv't) = G an d long-ru n covarianc e Q = G + A + A'. Fro m
(A4), A ca n be writte n as:
Then:
Co-integration i n System s of
Equations
We hav e s o fa r considere d onl y single-equatio n estimatio n an d
testing. Whil e th e estimatio n o f singl e equation s i s convenien t an d
often efficient , fo r som e purpose s onl y estimatio n o f a syste m
provides sufficien t information . Thi s i s true, fo r example , whe n we
consider th e estimatio n o f multipl e co-integratin g vectors , an d
inference abou t th e numbe r o f suc h vectors. Traditionally , system s
have bee n estimate d whe n ther e i s a failur e o f weak exogeneit y i n
a singl e equation , an d thes e consideration s als o appl y here . Thi s
chapter examine s method s o f findin g th e co-integratin g rank ,
considers eircumstance s whe n dynami c single-equatio n method s
will be asymptoticall y equivalen t t o system s methods , an d provide s
examples t o illustrat e thes e issues . Asymptoti c distribution s ar e
also derived .
In earlie r chapters , w e investigate d dat a serie s containin g uni t root s i n
their scala r autoregressiv e representation s (i.e . thei r margina l distribu tions), an d denote d suc h serie s a s 1(1). I n thi s chapter w e will consider a
vector tim e serie s of dimensio n n, a, = (*u,*2o > x nt)' (generalizin g
the analysi s t o an y numbe r o f variables) , wher e x , i s 1(1 ) s o tha t Ax r i s
1(0). Generally , an y arbitrar y linea r combinatio n o f th e element s o f x f ,
say w ( = a'x t, wil l als o b e 1(1) , an d suc h linea r combination s impl y o r
give ris e t o spurious regressions. However , ther e ma y exis t vector s a ,
such tha t
in whic h case th e relevant component s o f \t are co-integrated .
In th e simples t bivariat e case , a s w e hav e seen , w e ma y tak e
xf = (y t, z ty, wher e y t an d z t ar e individuall y 1(1). Th e arbitrar y linea r
combination (y, - Kz t) wil l als o b e 1(1) , bu t i f there exist s a value i q of
K suc h tha t (y, - jqz, ) ~ 1(0) , the n y t an d z t ar e co-integrated . Lettin g
a{ = (1, iq) b e th e co-integratin g vecto r i n thi s case , a^ mus t b e
unique, sinc e fo r an y othe r valu e K*, the n y t K*zt = yt ~ *q r +
(jq - K*)z t = w t + (KI n*)zt, whic h i s the su m of a n 1(0 ) proces s an d
an 1(1) process, an d therefor e 1(1 ) unless j q = ie* .
256 Co-integratio
n i n System s o f Equations
where fi = C(l)m .
The ke y assumption s needed t o deriv e th e autoregressiv e representa tion o f th e proces s ar e give n below . A s i n Chapte r 5 , th e proo f follow s
Johansen (1991a) .
ASSUMPTION Bl. Th e characteristi c polynomial,
258 Co-integratio
n i n System s o f Equations
We therefor e hav e
with
From th e definitio n of Ae ?,
260 Co-integratio
n i n System s of Equations
(11')
The for m i n (11 ) i s th e 'canonical ' representatio n i n 1(0 ) space , an d
Phillips (1991 ) focuse s o n estimatio n o f thi s system . Whe n
E(uitu2t) + 0, a 'simultaneit y problem ' i s present, bu t thi s ca n be deal t
with b y th e inclusio n of A.x 2t a s a regressor i n th e firs t equatio n o f (11).
The functiona l central-limi t theorem s fo r Wiene r processe s note d i n
earlier chapter s appl y despit e th e seria l dependenc e i n u ( = [ lf , u 2t]',
and direc t estimatio n o f K in th e firs t equatio n o f (11) ca n b e see n a s th e
method originall y propose d b y Engl e an d Grange r (1987) . Inferenc e
must, however, allo w for the seria l dependenc e i n ut.
The latte r system , (11') , highlight s the 'structural ' form . At leas t one
of Yi o r 7 2 mus t b e non-zero , sinc e otherwis e th e syste m ca n b e
expressed i n term s o f difference d variable s alone . Wea k exogeneit y i s
violated b y (among other possibilities ) YiY2 ^ 0. Since we are unlikely to
know a priori whic h other equations ar e influence d by an y give n ECM,
we tur n no w t o a metho d o f estimatin g th e co-integratin g rank r o f a
system, which will also allo w tests o f this aspect o f weak exogeneity.
262 Co-integratio
n i n System s of Equations
3
Phillip s an d Durlau f (1986 ) deriv e th e limitin g distributio n o f th e least-square s
estimator o f (the equivalen t of) n , allowin g fo r more genera l error processes .
4
Not e that we use th e upper-cas e n fo r th e rati o of the circumferenc e of a circle to it s
diameter, a s opposed to the lower-case n define d earlie r a s the matrix product yo'.
264 Co-integratio
n i n System s o f Equation s
266 Co-integratio
n i n System s o f Equation s
268 Co-integratio
n i n System s of Equations
Th e table s i n Osterwald-Lenu m (1992 ) giv e critical values for value s o f n runnin g fro m
1 to 1 1 and are therefor e mor e extensiv e than those i n Johansen (1988) . W e ar e gratefu l t o
Michael Osterwald-Lenu m fo r permission t o reproduc e thi s table.
% 95
2
3
4
5
6
7
8
9
10
11
2.86
9.52
15.59
21.58
27.62
33.62
38.98
44.99
50.65
56.09
61.96
1
2
3
4
5
6
7
8
9
10
11
2.86
10.47
21.63
36.58
55.44
78.36
104.77
135.24
169.45
206.05
248.45
% 97.5
t,r (A-max )
3.84
11.44
17.89
23.80
30.04
36.36
41.51
47.99
53.69
59.06
65.30
t\r (trace )
3.84
12.53
24.31
39.89
59.46
82.49
109.99
141.20
175.77
212.67
255.27
% 99
4.93
13.27
20.02
26.14
32.51
38.59
44.28
50.78
56.55
61.57
68.35
6.51
15.69
22.99
28.82
35.17
41.00
47.15
53.90
59.78
65.21
72.36
4.93
14.43
26.64
42.30
62.91
86.09
114.22
146.78
181.44
219.88
261.71
6.51
16.31
29.75
45.58
66.52
90.45
119.80
152.32
187.31
226.40
269.81
270 Co-integratio
n i n System s of Equations
n 2 = r =0
n i_ ~ r = j
16.1
0.41
14..1
3..76
16 .5
0 .41
15,.4
3,.76
ph
pn 1.00
ph -1.06
0 -1.07
3 1.00
7
0
3 -0.00
2 -0.01
7
9
8.5. Extension s
The precedin g result s hol d fo r a simpl e model . Severa l possibl e exten sions an d othe r consideration s aris e i n thi s mode l an d w e shal l briefl y
consider eigh t o f these:
1. dumm y variables (suc h a s constants an d trends) ;
2. linea r restriction s o n co-integrating vectors ;
3. power s o f tests ;
4. forecastin g in co-integrated processes ;
5. finite-sampl e properties;
6. selectin g la g length;
7. 1(2 ) variables;
8. wea k exogeneity an d conditional models .
8.5.1. Dummy Variables
The firs t issu e o f practica l importanc e i s th e potentia l presenc e o f
intercepts i n th e equations . Th e inclusio n o f intercept s i n th e estimate d
system alter s th e critica l values of th e test s fro m thos e tha t obtai n whe n
no intercept s ar e presen t (a s a compariso n o f Tabl e 8. 1 (n o constant )
with Tabl e 8. 5 belo w shows) . Unde r th e nul l o f n o co-integratin g
vectors, non-zer o intercept s woul d generat e trends . However , eve n i n
equations wit h ECMs , tw o possibilitie s arise : tha t th e intercep t enter s
only i n th e ECM , o r tha t i t als o enter s a s an autonomou s growt h facto r
in th e equation . Bot h case s ar e considere d b y Osterwald-Lenu m (1992 )
and Johanse n an d Juselius (1990) . I n term s o f (12), th e mode l become s
272 Co-integratio
n i n System s o f Equations
hi s
so that
which determine s th e growt h i n th e system . Sinc e n*y= ( I + y')y =
y(I + a'y ) = y^ > where , matchin g th e structur e o f C , ip = (I + 'y) , i t
follows tha t jr* s y= yi/^ . Bu t sinc e C define s th e 1(0 ) representation ,
tys 0 a s s o o , s o tha t JT * ha s som e root s equa l t o unit y an d a
convergent componen t ip. I n a bivariate case , ty would b e th e stationar y
root o f JT* .
The matri x K i s non-symmetri c an d idempoten t wit h a' K = 0' an d
K y = 0 s o that,?r* K = K. Also , whe n y = 0 the n K = I. Sinc e th e
condition tha t fi fall s i n th e co-integratin g spac e i s fi = y 0 wher e 0 i s
r x 1 , then
confirming th e absenc e o f any linear tren d i n x, when fi = y 0.
Further, th e unconditiona l varianc e matri x o f w t , var[w, ] = G , i s
G = CGC' + , or
274 Co-integratio
n i n System s o f Equation s
% 95
2
3
4
5
6
7
8
9
10
11
2.69
12.07
18.60
24.73
30.90
36.76
42.32
48.33
53.98
59.62
65.38
1
2
3
4
5
6
7
8
9
10
11
2.69
13.33
26.79
43.95
64.84
89.48
118.50
150.53
186.39
225.85
269.96
% 97.5
t,r (A-max )
3.76
14.07
20.97
27.07
33.46
39.37
45.28
51.42
57.12
62.81
68.83
r\r (trace )
3.76
15.41
29.68
47.21
68.52
94.15
124.24
156.00
192.89
233.13
277.71
% 99
4.95
16.05
23.09
28.98
35.71
41.86
47.96
54.29
59.33
65.44
72.11
6.65
18.63
25.52
32.24
38.77
45.10
51.57
57.69
62.80
69.09
75.95
4.95
17.52
32.56
50.35
71.80
98.33
128.45
161.32
198.82
239.46
284.87
6.65
20.04
35.65
54.46
76.07
103.18
133.57
168.36
204.95
247.18
293.44
90%
1
2
3
4
5
6
7
8
9
10
11
6.50
12.91
18.90
24.78
30.84
36.35
42.06
48.43
54.01
59.19
65.07
1
2
3
4
5
6
7
8
9
10
11
6.50
15.66
28.71
45.23
66.49
90.39
118.99
151.38
186.54
226.34
269.53
0 ^ 0 e
. ^_
i.
f~IN(0,
95%
r (A-max)
8.18
14.90
21.07
27.14
33.32
39.43
44.91
51.07
57.00
62.42
68.27
r]r (trace )
8.18
17.95
31.52
48.28
70.60
95.18
124.25
157.11
192.84
232.49
277.39
fi);
97.5%
ft + e,
99%
9.72
17.07
22.89
29.16
35.80
41.86
47.59
53.85
59.80
64.98
70.69
11.65
19.19
25.75
32.14
38.78
44.59
51.30
57.07
63.37
68.61
74.36
9.72
20.08
34.48
51.54
74.04
99.32
129.75
162.75
198.06
238.26
283.84
11.65
23.52
37.22
55.43
78.87
104.20
136.06
168.92
204.79
246.27
292.65
276 Co-integratio
n i n System s o f Equations
% 95
2
3
4
5
6
7
8
9
10
11
7.52
13.75
19.77
25.56
31.66
37.45
43.25
48.91
54.35
60.25
66.02
1
2
3
4
5
6
7
8
9
10
11
7.52
17.85
32.00
49.65
71.86
97.18
126.58
159.48
196.37
236.54
282.45
% 97.5
r (A-max )
9.24
15.67
22.00
28.14
34.40
40.30
46.45
52.00
57.42
63.57
69.74
77, (trace )
9.24
19.96
34.91
53.12
76.07
102.14
131.70
165.58
202.92
244.15
291.40
% 99
10.80
17.63
24.07
30.32
36.90
43.22
48.99
54.71
60.50
66.24
72.64
12.97
20.20
26.81
33.24
39.79
46.82
51.91
57.95
63.71
69.94
76.63
10.80
22.05
37.61
56.06
80.06
106.74
136.49
171.28
208.81
251.30
298.31
12.97
24.60
41.07
60.16
84.45
111.01
143.09
177.20
215.74
257.68
307.64
(39)
In plac e o f (27) , w e mus t solv e fo r th e eigenvalue s A f s= A f & . . . ^ Af
from th e equatio n
using th e principle s applie d above . A likelihood-rati o tes t agains t th e
unrestricted valu e o f a ca n b e calculate d an d amount s t o testin g H }
within H r, an d is therefore based o n
278 Co-integratio
n i n System s o f Equations
280 Co-integratio
n in System s of Equations
= 0.02, an d
Note no w tha t consumptio n an d incom e ar e bot h 1(1 ) variables , con sumption an d incom e ar e co-integrated , an d savin g i s a stationar y
variable. Th e equation s
281
1i
283
284 Co-integratio
n i n System s o f Equation s
where
The conditiona l expectatio n E[x T+h X T] a t T i s
when m = min(/, h) .
When th e syste m i s expresse d i n difference s t o forecas t Ax
outcomes ar e give n by
T+h,
286 Co-integratio
n i n System s o f Equation s
Reimers (1991 ) compare s th e power s o f variou s test s fo r co-integra tion fo r bivariat e an d trivariat e processes . H e find s tha t th e Johanse n
procedure over-reject s whe n th e nul l i s true , i n smal l samples ,
and suggest s correctin g thi s usin g ( T - p)log( l - A,- ) instea d o f
T log (1 A,-) fo r th e tes t statistic s wher e p = nk take s accoun t o f th e
number o f estimate d parameters . Whil e nk/T i s asymptoticall y negli gible, i t ca n b e larg e i n smal l samples . Th e powe r o f th e test s i s
dependent o n th e specificatio n of the DGP , bu t Reimer s doe s no t relat e
his simulation finding s t o th e typ e of analysis in Section 8.5.3 .
8.5.6. Selecting Lag Length
Both Gonzalo' s (1990 ) an d Reimers' s (1991 ) studie s conside r th e effect s
on th e ML E o f usin g incorrec t la g length s fo r th e short-ru n dynamics .
Gonzalo find s tha t th e los s o f efficienc y fro m choosin g to o lon g a la g is
small, an d tha t th e ML E perform s best eve n i f a la g o f fou r period s i s
used fo r th e short-ru n dynamic s instea d o f th e correc t valu e o f 0 .
However, i f to o shor t a la g lengt h i s use d (fo r example , zer o lag s
instead o f one ) the n th e ML E i s n o longe r th e bes t method . Mor e
practical experienc e is required befor e a fina l judgemen t can b e reache d
on th e relativ e cost s o f under-specifyin g versu s over-specifyin g th e
lag-length, bu t Gonzalo' s simulatio n evidenc e seem s intuitivel y reason able sinc e under-specificatio n wil l induc e residua l autocorrelation .
Reimers find s tha t th e Schwar z criterio n doe s wel l i n a data-base d
lag-length selectio n exercise . However , sinc e th e rol e o f th e Ax ( _, i s t o
whiten th e error , i t i s no t clea r tha t th e us e o f th e Schwar z criterion ,
which penalize s th e additio n o f lag s strongly , will prov e optima l i n thi s
context.
8.5.7. The Analysis of 1(2) Variables
Reconsider th e basi c autoregressive system with lag length k, written as
where A 0 = I, s o that
we see tha t
288 Co-integratio
n i n System s o f Equations
(56)
Since w ( _! = (xJ-jtrAx^-j)', al l th e informatio n o n th e co-integratin g
vectors i s retained i n Wj_j . Consequently , Ax 2f i s weakly exogenous fo r
<j> i f (jt depend s o n A t alone , an d A : an d A 2 ar e variatio n free , s o tha t
A = A j x A 2. Wea k exogeneit y o f Ax 2( fo r (j> canno t occu r whe n A ! an d
A2 bot h depen d o n commo n component s o f a .
As a consequenc e o f th e normalit y assumption , an d usin g the expres sion in (55) for ya'x^, conditionin g Ax lf o n Ax 2, lead s t o th e mea n of
the conditiona l density:
290
292 Co-integratio
n i n System s o f Equation s
0.013817
Statistics
-riog(i-ft.;) ,(0.05
n 4= r =
n 3= r = 1
n - 2 = r =2
n - 1 =r = 3
72.82
28.73
6.22
1.39
0.060350
30.33
23.78
16.87
3.74
0.249694
0.517240
-riog(l - M, ;) ?
109.17
36.34
7.62
1.39
n - r (0.05)
54.64
34.55
18.17
3.74
m p
1.0000
0.0311
-0.2633
0.9838
'l
l>2
R,,
6.3966
1.0000
0.9435
4.5659
-0.8938
-0.3334
1.0000
-0.7701
7.6838
-0.1377
-1.2117
1.0000
293
7i
72
m- p
Ap
-0.0952
0.0048
-0.0210
-0.0001
0.4268
-0.5147
0.2578
-0.2253
*85
Rn
-0.0300
-0.0013
-0.0318
0.0796
-0.0076
0.0024
0.0116
0.0069
294
Further,
where
296 Co-integratio
n in System s o f Equations
Hence
using th e decomposition BI(S) = a)i 2a)22B2(s) + (<ui i where eo, y i s th e (i , /)th elemen t o f Q an d W(s ) i s a Brownia n motio n
process independen t o f B 2 sinc e aj 12co2~21 = 0cr 1/a2(l - p ) an d
The importan t propertie s o f stati c estimatio n o f f l ar e apparen t fro m
(69) an d (70) . The estimato r i s super-consistent , bu t contain s second order biase s reflecte d i n (706 ) an d (70c) ; thes e term s als o mak e
standard distribution s inappropriat e fo r hypothesi s testing . Estimatio n
will b e improve d t o th e exten t tha t (706 ) and (70c ) can b e reduce d o r
eliminated, an d the y will vanish if the short-ru n an d long-ru n multipliers
are equal , sinc e i/ > = / 3 implie s 6 = 0, an d s o A 2 = A% = 0. Whil e th e
limiting distributio n abov e i s specifi c t o th e DG P (63) , i/ > = / ? wil l
typically onl y aris e becaus e o f a n absenc e o f lagge d value s o f z t an d y t
from th e DGP ; if fo r exampl e y, = yzt + Y\yt-\ + Y2Zt-i + error, the n
the long-ru n multiplie r i s / ? = ( V + 72)/( l ~ 7i) > m whic h cas e
7i 72 = 0 i s sufficien t fo r fi = ty . A commo n facto r (y 2 = VYi) i s
necessary an d sufficient .
The term s A 2 an d A 3 abov e ca n b e eliminate d whe n if> = / ? by th e us e
of othe r estimatio n methods, a s will be see n below .
(ii) Non-linear least squares (Stock 1987). Thi s method , whic h elimin ates th e bia s containe d i n (70c) , consist s i n minimizin g th e su m o f
squared residual s defined as
298 Co-integratio
n i n System s of Equations
Conclusion
We briefl y summariz e th e mai n theme s o f th e book , an d the n
consider th e invarianc e o f th e matri x o f co-integrating vectors i n a
linear syste m unde r bot h linea r transformation s an d seasona l
adjustment. Next , co-integratio n i s related t o structure d time-serie s
models, whic h offe r a n alternativ e approac h t o modellin g inte grated data . Recen t researc h o n integratio n an d co-integratio n i s
described, an d th e boo k conclude s b y re-interpretin g som e ol d
econometric problem s i n the ligh t of co-integration theory .
9.1. Summar y
Many economi c tim e serie s appea r t o b e non-stationar y and to drif t ove r
time. Efficien t inferenc e i n time-serie s econometric s require s takin g
account o f thi s phenomenon . Thi s boo k describe d th e modellin g o f
economic variable s a s integrate d processes , allowin g fo r th e possibilit y
that variable s ma y b e linke d i n th e lon g run , implyin g tha t linea r
combinations of them ar e co-integrated .
We firs t presente d th e backgroun d t o th e theor y o f integrate d series ,
building o n concept s fro m time-serie s analysi s an d th e theor y o f sto chastic processes . Th e resultin g distribution s o f estimator s an d test s
applied t o integrate d dat a wer e functional s o f Wiene r processes , whic h
when combine d wit h a functional central-limi t theorem le d to a powerfu l
and genera l metho d fo r derivin g their limitin g distributions. These wer e
different fro m th e limitin g distribution s conventionall y applie d t o sta tionary processes , bot h becaus e th e normalizatio n facto r was the sampl e
size rathe r tha n it s squar e root , an d becaus e th e for m o f the asymptoti c
distribution wa s non-normal . A n importan t implicatio n wa s tha t th e
critical value s o f tes t statistic s differe d betwee n 1(0 ) an d 1(1 ) data .
Although th e asymptoti c distributio n theor y involve d ne w type s o f
derivations, i t wa s feasibl e t o maste r th e logi c o f Wiene r processe s
without excessiv e effort ; th e pay-of f wa s tha t th e approac h simplifie d
other derivation s (suc h a s constanc y tests , a s i n Hanse n 1992) , and , i n
addition, wa s very general.
The Wiene r proces s tool s the n allowe d u s t o analys e suc h divers e
problems a s spuriou s (o r nonsense ) regressions , spuriou s detrending ,
300 Conclusio
Conclusion 30
(1
)
where e ( ~IN(0,i;). Th e syste m i n (1 ) ha s parameter s (T , y, a, fi, E).
Then, \, is 1(1 ) i f an d onl y i f rank (yl^aj j = n r wher e * P i s th e
mean la g matrix defined i n Chapter 8 . Here (y : y ) has rank n, with y
being n X (n r) suc h tha t y i y = 0 an d (a:a ) ha s ran k n wit h
^ = 0 fo r _ L o f siz e nx(n r). Pre-multiplyin g (1 ) b y a know n
n x n non-singula r matri x B (s o | B = 0),
t
302 Conclusio
When Sf(-L ) i s a scala r time s th e uni t matri x (th e sam e filte r fo r al l x it),
vat = ef. I n (6) , i t look s a s i f y i s als o a n invariant , bu t a s o at involve s
lagged, current , an d futur e difference s of x, o f dth o r highe r order , a s
well a s e", the n on e o f v at o r e t i s likel y t o b e autocorrelated . Sinc e
'x?_i i s a n 1(0 ) variable , conventiona l seria l correlatio n biase s appl y t o
it, an d henc e y will usuall y b e affecte d b y whethe r o r 'not th e dat a ar e
seasonally adjusted . Th e short-ru n dynamic s wil l be change d whe n e t i s
an innovation , becaus e v" i s correlate d wit h Ax?_i , an d additiona l lag s
are neede d t o remov e it s autocorrelation .
Conclusion 30
304 Conclusio
by (14) .
Thus, a s long as cr 2 + 0,
Hence fro m (7) ,
When cr ^ = 0, we have t = t, t_v = 0, say, so that
and C o i s th e mea n growt h rat e [Ay r ] = g y = 0- Whe n a 2 =0, (18 )
entails changes in [Ay r ] = g y (f) over tim e an d generate s y , a s 1(2).
The alternativ e possibilit y to evolvin g growt h rate s i s tha t o f change s
in mean s ove r time , s o tha t g y(t) take s differen t value s i n differen t
epochs. Suc h behaviour coul d b e approximate d b y a mode l i n which th e
distribution D n(r]t) wa s non-normal, wit h a large mass a t zer o an d smal l
probabilities o f larg e values . The n r woul d usuall y b e constant , bu t
would occasionall y jum p t o a ne w level . Thus , i t i s unsurprisin g tha t
discrimination betwee n integrate d an d regime-chang e model s i s difficul t
(see Perro n 1989) . Conversely , ther e ar e clos e affinitie s betwee n struc tured time-serie s an d econometri c model s fo r integrate d data . Indeed ,
several researcher s hav e suggested switchin g from a unit-root nul l to on e
of 1(0 ) o r co-integration . Fo r example , on e migh t see k t o tes t a 2, = 0
when a ^ = 0 (s o r = W) a s a tes t fo r a uni t roo t (se e e.g . Kwiatkowski, Phillips , an d Schmid t (1991) an d Leybourn e an d McCab e (1992)) .
Conclusion 30
306 Conclusio
Conclusion 30
9.6.1. Multi-collinearity
When x , ~ 1(1 ) an d a'x , ~ 1(0) , the n includin g all the element s o f x ( o r
\t-i a s regressors i n a singl e equatio n wil l induc e a n apparentl y seriou s
collinearity problem . Th e secon d momen t matri x (X'X ) will b e O(T 2),
whereas th e linea r combinatio n (a'X'Xa ) wil l b e O(T). Consequently ,
(T~ 2 X'X) will converge on a singular matrix . Generally, it is inadvisable
to 'solve ' thi s proble m b y deletin g variables ; fo r 1(1 ) data , doin g s o
jeopardizes th e possibilit y of co-integration . I f th e dependen t variabl e i s
1(0), the n th e solutio n i s to fin d th e co-integratin g combination a'x t o r
'x,-i an d us e tha t a s a n explanator y variable . Thi s strateg y cor responds t o th e usua l recommendatio n o f transformin g t o near-ortho gonal an d interpretabl e variables . I n othe r cases , wher e th e dependen t
variable i s 1(1) bu t i s co-integrated wit h a subset o f \t, say, elimination
may b e sensible , bu t Wiener-base d critica l value s shoul d b e use d fo r
variables tha t canno t b e writte n implicitl y a s a n 1(0 ) functio n (se e
Chapter 7) . Thes e idea s ar e relate d t o th e earlie r techniqu e o f con fluence analysi s in Hendry an d Morga n (1989) .
308 Conclusio
Conclusion 30
310 Conclusio
References
ABADIR, K . M . (1992) , 'Th e Limitin g Distributio n o f th e Autocorrelatio n
Coefficient Unde r a Unit Root' , Annals o f Statistics, forthcoming .
AHN, S . K. , an d REINSEL , G . C . (1988) , 'Neste d Reduced-Ran k Autoregressiv e
Models fo r Multipl e Tim e Series' , Journal o f th e American Statistical Association, 83: 849-56.
ANDERSON, T . W . (1958) , A n Introduction t o Multivariate Statistical Analysis,
John Wiley , New York.
(1976), 'Estimatio n o f Linea r Functiona l Relationships : Approximat e Distributions an d Connection s wit h Simultaneou s Equation s i n Econometric s
(with discussion)' , Journal of th e Royal Statistical Society B,38 : 1-36 .
ANDREWS, D . W . K. , an d PLOBERGER , W . (1991) , 'Optima l Test s o f Paramete r
Constancy', mimeo. , Yale University Press.
BANERJEE, A. , an d DOLADO , J . (1987) , 'D o W e Rejec t Rationa l Expectation s
Models Too Often ? Interpretin g Evidence using Nagar Expansions', Economics Letters, 24: 27-32.
(1988), 'Test s o f th e Lif e Cycle-Permanen t Incom e Hypothesi s i n th e
Presence o f Rando m Walks : Asymptoti c Theor y an d Smal l Sampl e Interpre tations', Oxford Economic Papers, 40: 610-33.
-and GALBRAITH , J . W . (1990a) , 'Orthogonalit y Test s wit h De-trende d
Data: Interpretin g Mont e Carl o Result s using Nagar Expansions' , Economics
Letters, 32: 19-24.
-HENDRY, D . F. , an d SMITH , G . W . (1986) , 'Explorin g Equilibriu m
Relationships i n Econometric s throug h Stati c Models : Som e Mont e Carl o
Evidence', Oxford Bulletin of Economics an d Statistics, 48: 253-77.
-GALBRAITH, J . W. , an d DOLADO , J . (19906) , 'Dynami c Specificatio n with
the Genera l Error-Correctio n Form' , Oxford Bulletin o f Economics an d
Statistics, 52: 95-104.
-and HENDRY , D . F . (eds. ) (1992) , Testing Integration an d Cointegration,
special issue of th e Oxford Bulletin of Economics and Statistics, 54, 225-55.
BARDSEN, G . (1989) , 'Th e Estimatio n o f Long-Ru n Coefficient s fro m Error Correction Models' , Oxford Bulletin of Economics and Statistics, 51: 345-50.
BEWLEY, R . A . (1979) , 'Th e Direct Estimatio n of the Equilibriu m Response i n a
Linear Model' , Economics Letters, 3 : 357-61.
BEWLEY, R . A. , ORDEN , D. , an d FISHER , L . (1991) , 'Box-Tia o an d Johanse n
Canonical Estimator s o f Cointegratin g Vectors' , Universit y o f Ne w Sout h
Wales, Economics Discussion Paper, 91/5 .
BHARGAVA, A . (1986) , 'O n th e Theor y o f Testin g fo r Uni t Root s i n Observe d
Time Series' , Review of Economic Studies, 53 : 369-84.
BILLINGSLEY, P . (1968) , Convergence of Probability Measures, John Wiley , New
York.
BOSSAERTS, P . (1988) , 'Commo n Non-Stationar y Components o f Asse t Prices' ,
Journal o f Economic Dynamics an d Control, 12 : 347-64.
312 Reference
BOSWIJK, H . P . (1991) , 'Testin g fo r Cointegratio n i n Structura l Models', Univer sity o f Amsterdam, Econometric s Discussio n Pape r AE7/91 .
(1992), 'Efficien t Inferenc e on Cointegratio n Parameter s i n Structural Erro r
Correction Models' , Universit y o f Amsterdam , Econometric s Discussio n
Paper,
-and FRANSES , P . H . (1992) , 'Dynami c Specificatio n an d Cointegration' ,
Oxford Bulletin o f Economics an d Statistics, 54: 369-81.
Box, G . E . P. , an d JENKINS , G. M . (1970) , Time Series Analysis Forecasting and
Control, Holden-Day , Sa n Francisco.
and TIAO , G . C . (1977) , ' A Canonica l Analysi s o f Multipl e Tim e Series' ,
Biometrika, 64: 355-65.
BRANDNER, P. , an d KUNST , R . (1990) , 'Forecastin g Vecto r Autoregressions : Th e
Influence o f Cointegration', Memorandu m 265 , IAS , Vienna .
CAMPBELL, B. , an d DUFOUR , J.-M . (1991) , 'Over-Rejection s i n Rationa l Expec tations Models : A Non-Parametri c Approac h t o th e Mankiw-Shapir o Prob lem', Economics Letters, 35 : 285-90.
CAMPBELL, J . Y. , an d PERRON , P . (1991) , 'Pitfall s an d Opportunities : Wha t
Macroeconomists Shoul d Kno w Abou t Uni t Roots' , i n Blanchard , O . J . an d
Fischer, S . (eds) , NBER Economics Annual 1991, MIT Press .
and SHILLER , R . J . (1991) , 'Cointegratio n an d Test s o f Presen t Valu e
Models', Journal o f Political Economy, 95 : 1062-88.
CHAMBERS, M . J . (1991) , ' A Not e o n Forecastin g i n Co-Integrate d Systems' ,
Department o f Economics, Universit y of Essex .
CHAN, N . H. , an d WEI , C. Z . (1988) , 'Limitin g Distribution s o f Least-Square s
Estimates o f Unstabl e Autoregressiv e Processes' , Annals o f Statistics, 16 :
367-401.
CHOI, I . (1992) , 'Durbin-Hausma n Test s fo r Uni t Roots' , Oxford Bulletin o f
Economics an d Statistics, 54: 289-304.
CHONG, Y . Y. , an d HENDRY , D . F . (1986) , 'Econometri c Evaluatio n o f Linea r
Macroeconomic Models' , Review o f Economic Studies, 53 : 671-90.
CHOW, G . C . (1960) , 'Test s o f Equalit y Betwee n Set s o f Coefficient s i n Tw o
Linear Regressions' , Econometrica, 52: 211-22.
CHU, C.-S . J. , an d WHITE , H . (1991) , 'Testin g fo r Structura l Chang e i n som e
Simple Tim e Serie s Models' , Discussio n Pape r 91-6 , Universit y of California,
San Diego, Dept . o f Economics .
(1992) ' A Direc t Tes t fo r Changin g Trend' , Journal o f Business an d
Economic Statistics, 10: 289-99.
CLEMENTS, M . P. , an d HENDRY , D . F . (1991) , 'O n th e Limitation s o f Mea n
Square Erro r Forecas t Comparisons' , Discussio n pape r 138 , Oxfor d Institut e
of Economic s an d Statistics . Forthcoming, Journal o f Forecasting.
(1992), 'Forecastin g i n Cointegrate d Systems' , Discussio n pape r 139 ,
Oxford Institut e o f Economics an d Statistics .
DAVIDSON, J . E . H. , HENDRY , D . F. , SRBA , F. , an d YEO , S. (1978) , 'Economet ric Modellin g of th e Aggregat e Time-Serie s Relationshi p Between Consumers '
Expenditure an d Incom e i n th e Unite d Kingdom' , Economic Journal, 88 :
661-92.
DAVIDSON, R. , an d MACKINNON , J . G . (1992) , Estimation an d Inference i n
Econometrics, Oxfor d University Press.
DEATON, A . S. , an d MUELLBAUER , J . N . J . (1980) , Economics an d Consumer
References 31
314 Reference
References 31
316 Reference
References 31
318 Reference
References 31
320 Reference
Author Inde x
Abadir, K . M . 126 , 128
Abrahamse, A . P . J . 10 4
Ahn, S . K . 30 5
Anderson, G . J . 5 , 50, 140
Anderson, H . 27 2
Anderson, T. W . 70n. , 26 5 n., 285
Andrews, D . W . K . 31 0
Banerjee, A . 55 , 95, 97, 163 , 166, 177n.,
187, 191 , 192, 214, 215, 220, 222, 230 ,
233, 306 , 307
Bardsen, G . 47 , 53, 56, 62, 235
Bewley, R. 47 , 49, 53, 152 , 305
Bhargava, A. 101 , 104, 155, 176, 207, 209
Billingsley, P . 24 , 89
Birchenhall, C . R . 12 2
Bossaerts, P . 29 8
Boswijk, H . P . 235 , 305, 307, 310
Box, G . E . P . 10 , 13, 121, 305
Brandner, P . 28 2
Breusch T . S . 47 , 55 , 56 , 59 , 62 , 63 , 64
Campbell, B . 167n .
Campbell, J . Y . 30 6
Campos, J . 23 6
Chan, N . H . 91 , 96 n.
Chiu, A . P . L . 12 2
Choi, I . 30 6
Chong, Y . Y . 28 2
Chow, G . C . 194n .
Chu, C.-S . J. 31 0
Clements, M . P . 282 , 283, 285
Davidson, J . E . H . 5 , 50, 52, 140, 300
Davidson, R. 16 , 28
Deaton, A. S . 5 3
Dickey, D . A . 8 , 24, 82, 100 , 103, 107,
108, 112-23 , 169
Dolado, J. J . 55 , 97, 163, 166, 177n., 187,
191, 192 , 230
Dufour, J.-M . 167n.
Durlauf, S . N . 82 , 92 , 93 , 182 , 203, 238 ,
254, 262n .
Engle, R . F . 6 , 7, 17 , 18, 19, 43, 67, 84n.,
121, 122 , 137 n., 145 , 146, 152, 157-9,
163, 205n. , 208, 209, 211, 215, 231, 242 ,
256, 261, 278, 279, 282, 283, 287, 288 ,
305, 30 9
324
Author Index
Yeo, S . 5
Yoo, B . S . 121 , 152 , 208 , 209 , 278 , 279 ,
282, 283 , 287 , 305
Yule, G . U . 69 , 70n., 71, 77, 138
Subject Inde x
absolute summabilit y 15 8
adjustment:
coefficient 15 5
disequilibrium 51 , 52, 55, 61
speed of 26 8
approximation theore m 12 3
asymptotic:
convergence 15 8
independence 16 , 17
normality 105 , 126, 134, 163, 177, 178,
180, 185 ; and drif t ter m 169-7 4
asymptotic standar d erro r (ASE ) 235
Augmented Dickey-Fulle r tes t (ADF ) 106 ,
108, 109 , 207-12, 232-4 , 238, 239 n.
asymptotic distributio n 127 , 128
comparison wit h non-parametrically ad justed D F 114- 9
use o f IV i n 11 9
autocorrelation 13 , 71-2, 83 , 129, 163, 191,
206, 207, 212, 221 n., 238-42, 244,
286, 29 2
function 12 , 1 3
autocovariance functio n 12 , 13
autoregressive:
-distributed lag (ADL) model 47-55 ,
60-4, 224 , 239, 242
error 83 , 114 , 191, 291
process 12 , 72, 251, 257-60; see also
autoregressive moving-average
(ARMA) proces s
representation (VAR) , see co-integrat ing: representations o f co-integrate d
systems
autoregressive integrate d moving-averag e
(ARIMA) process 13 , 38, 39, 221
autoregressive moving-averag e (ARMA )
process 12 , 13, 39, 84 , 85 , 88 , 107,
108
examples o f 32- 8
Bardsen transformation , se e transformation: Bardse n
Bartlett windo w 24 8
Bewley:
representation 152 , 153
transformation, se e transformation :
Bewley
bias 67 , 68, 191 , 244, 246-8, 249, 250, 290 ,
309
in AR(1 ) parameter 100 , 101
correction ter m 241 , 246
326
Subject Inde x
Subject Index
generalized co-integratin g vector 17 9
general-to-specific modellin g 168 , 192
Granger causalit y 18 , 291
Granger Representatio n Theore m 48 ,
146-53, 300
homogeneity 47 , 51, 52, 60, 61, 221, 222 ,
231, 23 6
impact matri x 151 , 260
inconsistent regressio n 164-8 , 190 , 191,
229, 230
innovation sequenc e 12 , 85-7, 183
instrumental variable s (IV) 55 , 59, 62, 63,
119, 130- 3
integrated process 1 , 6, 7, 11 , 12, 21, 39 ,
69-71, 73, 136-8, 162-9 9
asymptotic theory o f 86-9 1
near-, see near-integrated process
properties o f 84- 6
see also non-stationar y proces s
integration:
order of , se e ordej r o f integration
seasonal, see seasonal integratio n
intercept 72 , 151 , 210, 232, 234, 271, 272 ,
273, 274
interim multiplie r representation 15 3
invariance 20 , 282, 283
principle 22 ; see also functiona l central
limit theore m
invertibility 13 , 84, 108 , 242
invertible system 148 , 149, 258, 259, 266
Jacobian 62 , 63
Johansen maximum-likelihoo d procedure 211 , 262-9, 285, 286, 300
power o f 277 , 278
Kronecker product 18 1
lag 9 , 11 , 47, 50, 52, 66, 106-8 , 123 , 225 ,
248, 250, 251, 286, 303
length 248 , 286
mean 28 7
polynomial 22 9
structure 208 , 222, 229
truncation paramete r 110 , 111, 113
latent roo t 13 , 104 , 142, 144, 158, 224
law o f large numbers 86 , 90
life-cycle hypothesi s 164 , 188
likelihood rati o test s 153 , 277, 278, 294,
295
limited-information maximum-likelihoo d
(LIML) 264 , 28 5
linear system 30 0
logarithms v. level s 29-32 , 193- 7
327
long-run:
covariance matrix 240 , 241, 245-7, 252,
290
multiplier 8 , 47-9, 51 , 54, 57, 59-64,
188, 230 , 235, 293, 295, 296; variance
of estimate s o f 61- 4
relationship 2 , 7, 8 , 140 , 220; see also
co-integrating: vecto r
response 15 3
solution 50 , 64-8
marginal:
distribution 18 , 19 , 290, 295
process 240 , 243-5, 248n .
marginalization 30 4
market clearing 3
martingale difference sequence (MDS ) 11,
12, 21, 163 , 179n., 185, 242, 244, 245 ,
247
maximal-eigenvalue statisti c 267 , 273
maximum-likelihood 159 , 241-5, 256 , 262 ,
264, 265, 266, 267, 269, 277, 283, 285,
286, 288
full-information, se e full-information
maximum-likelihood
limited-information, se e limited-information maximum-likelihood
mean la g 144 , 287, 301
memory 8 5
mixing:
coefficient 8 7
strong 16 , 17 , 87
uniform 16 , 17
mixingale 17 9 n.
Monte Carlo :
method 9 , 27, 28
response surface s 28 , 211, 213, 214
results 73-83 , 101 , 106 , 108, 114,
117-19, 133, 165, 214, 215, 222-3,
225-9, 232-5, 248-51, 279, 282, 283 ,
285, 291, 298
standard erro r 7 5
moving-average 12 , 88; see also auto regressive moving-average (ARMA)
process
component o f errors 10 7
negative components 113 , 119, 250, 304
parameter 24 8 n.
representation 133 , 153, 155, 156
seasonal filte r 12 1
multiple roots 119-2 2
multiplier, long-run, se e long-run: multiplier
near-integrated process 95-7 , 99 , 164, 166,
225, 231, 277
nearly-inconsistent regressio n 229 , 230
non-centrality parameter 97 , 98
328
Subject Inde x
non-parametric:
correction/test 9 , 108-10 , 114-9 , 130 ,
208, 210 , 211 , 238-40 , 25 1 asymptotic
theory o f 129-3 0
estimation 244 , 248 , 249
nonsense regressio n 69 , 80, 138
see also spuriou s regressio n
non-stationarity 4 , 8, 9, 65, 67, 72, 81-4 ,
134, 150 , 21 5
transformation t o stationarit y 69 , 70, 82,
83, 99 , 134 , 14 7
non-stationary process 5 , 6 , 9 , 38 , 39, 70,
71, 81 , 163 , 24 4
v. integrate d proces s 1 2
normality 180 , 28 9
asymptotic, se e asymptotic : normalit y
normalization 57-9 , 265 , 285
nuisance parameter s 100 , 104-6 , 172 , 176 ,
207, 21 0
order:
of magnitud e 14 , 15 , 21 , 9 0
in probabilit y 14 , 1 5
order of integration 6-9 , 48 , 79-80, 84 , 85,
147, 151 , 190-2 , 258
defined 8 4
first 137 , 17 7
higher 138 , 157 , 16 3
zero 13 7
Ornstein-Uhlenbeck proces s 9 6
orthogonal complemen t 14 7
orthogonality 86 , 149 , 151 , 242 , 244 , 245 ,
258n., 259,260 , 273
asymptotic 10 7
testing 164- 8
over-identification tes t 278 , 30 0
over-rejection 206 , 210 , 28 6
parameterization 48 , 207 , 208 , 250 , 274 ,
275
of dynamic s 22 1
exact 105 , 224
of nearly-integrate d processe s 9 5
over-/under- 224-9 , 262
permanent incom e hypothesi s 164 , 177 ,
178, 188 , 19 0
Perron-Phillips/Phillips test, se e non-para metric: correction/tes t
polynomial matrice s 140-5 , 152 , 257
isomorphism wit h companion mat rices 142- 4
power serie s expansio n 9 7
power o f tests 8 , 15 , 96, 101 , 108 , 113 , 198 ,
208, 214 , 223-4 , 230-5, 277 , 278 , 28 6
pre-determinedness 1 9
random wal k 11 , 21, 22, 24-9 , 38 , 71, 72,
82, 87 , 93 , 100 , 101 , 114 , 191 , 220 , 272
in logarithm s o r level s 19 3 n.
see also unit root
rank:
co-integrating, se e co-integrating: ran k
full 56 , 58 , 59 , 144 , 147 , 151 , 181 , 258 ,
260, 28 7
reduced 144 , 147 , 151 , 256 , 257 , 264 ,
285, 287 , 288 , 30 1
recursive estimatio n 194n. , 221 n.
re-parameterization 67 , 157 , 168 , 189 , 191 ,
222
see also transformatio n
representation theorem, see Granger Rep resentation Theore m
Said-Dickey tes t 107 , 108
compared wit h Perron-Phillip s tes t 11 3
Sargan-Bhargava test , se e Durbin-Watson
test (CRD W test)
Schwarz Criterion 194 , 28 6
seasonal adjustmen t filte r 301 , 303
seasonal integratio n 121- 3
sequential cu t 18 , 19
similar test s 100 , 104 , 105 , 16 9 n.
size distortion s 113 , 133 , 166 , 16 7
Slutsky's theore m 89 , 173
spurious:
correlation 70 , 71; in de-trended rando m
walks 82 , 8 3
regression 69-81 , 83 , 92-5, 134 , 138-9 ,
158, 159 , 162 , 191 , 230 , 25 5
stacked form , se e companion for m
static regression 162 , 163 , 167 , 205 , 214 ,
220-3, 231 , 238 , 246 , 251 , 29 6
comparison wit h dynami c 167 , 168 ,
224-30
example o f 23 6
see also Engle-Granger: two-ste p pro cedure
stationarity 1 , 4, 12 , 13 , 17 , 69, 212 , 26 2
stationary proces s 4 , 5, 6 , 7, 9, 11 , 29, 38,
39, 47 , 85 , 86 , 134 , 138 , 256 , 257 , 267 ,
279
strictly 11 , 1 2
weakly/second-order/covariance 11 , 1 2
stochastic:
differential equatio n 9 6
trend, se e trend, stochasti c
structural representatio n 261 , 30 3
super-consistency 158 , 176 , 191 , 214 , 220 ,
230, 251 , 294 , 296
total effect 142 , 25 7
trace 267 , 273
transformation 6 , 28-32, 88, 111 , 125 ,
178-80, 185
ADL 51 , 59
ADL t o EC M 60 , 61 300, 301
Subject Inde x
transformation (cont.):
Bardsen 51 , 54-9, 62 , 63
Bewley 51 , 53-6, 58n. , 59, 60, 62, 63
equivalence of , 54-60 , 62 , 64
linear 47 , 51 , 60, 61, 63, 64, 145 , 152 ,
178, 224 ; in dynamic regression 167-8 ,
177, 178 ; o f polynomial matrice s 144 ,
145
logarithmic 99 , 192- 9
trend (inclusio n of) 5 , 9, 82, 100, 101 , 106 ,
125, 185 , 211 , 212 , 213 , 214 , 236
non-stochastic (deterministic ) 6 , 20, 21,
69-72, 82, 84, 125 , 146 , 151 , 172 , 173 ,
185, 187 , 27 5
stochastic 153 , 169 , 172 , 174 , 179 , 180 ,
185, 187 , 191 ; se e also commo n trend ;
unit roo t
sums of powers o f 2 0
unit circl e 13 , 104 , 123 , 141 , 149 , 15 8
unit root 8 , 9 , 13 , 38, 72, 83-6, 95 , 96,
133, 144 , 147 , 163 , 177 , 185 , 215 , 236 ,
255, 258-60, 267, 270 , 287 , 289
multiple 12 2
329