A Note On The Stationarity and The Existence of Moments of The Garch Model
A Note On The Stationarity and The Existence of Moments of The Garch Model
A Note On The Stationarity and The Existence of Moments of The Garch Model
1. Introduction
Consider the following non-linear time series model
xt = t ht 2
ht = 0 + 1 x2t1 + + p x2tp + 1 ht1 + + q htq ,
(1.1)
506
for the GARCH(p, q) model in the literature. The purpose of this paper is to give
some sucient conditions for the strict stationarity and the existence of moments
for the GARCH(p, q) model.
2. The Strict Stationarity of the GARCH Model
First, we introduce some notation. Let {xt } conform to model (1.1). Dene yt = xt 2 , t = t 2 , Y t = (yt , . . . , ytp+1 , ht , . . . , htq+1 ) , Bt = (0 t , 0, . . . ,
0, 0 , 0, . . . , 0) , where denotes the transposition of a matrix, and
1 t p1 t p t 1 t q1 t q t
1
0
0
0 0
0
. .
.
.
.
.
.
.
.
..
..
..
..
..
.. . .
.
0
At =
1
..
.
1
p1
0
..
..
.
.
0 0
0
q1
0
..
..
.
.
0 1
0
p
0
..
.
0
1
1
..
.
0
q
0
..
.
1 p1 p 1 q1 q
1 0 0 0 0 0
. .
.. .. . .
.. ..
. . . ..
.
. . . . .
.
0 1 0 0 0 0
A.
EAt =
p1 p 1 q1 q
1
0 0 0 1 0 0
.
.. .. . .
. .
.. . .
. ..
.
. . . .. ..
0 0 1
xt = t
0 +
i xti +
i=1
which implies
yt =
p
i=1
j htj
j=1
i t 2 xti 2 +
i=1
j t 2 htj + 0 t 2 ,
(2.1)
j t htj + 0 t .
(2.2)
j=1
i t yti +
q
j=1
507
Thus, {yt } is a solution of (2.2) if and only if {Yt } is a solution of the following
stochastic dierence equation
Yt = At Yt1 + Bt .
Lemma 2.1. If
i=1 i
j=1 j
(2.3)
k1
Atj Btk
j=0
k=1
k1
k=1
Atj Btk ,
(2.4)
j=0
k1
Atj Btk =
k1
j=0
EAtj EBtk = Ak B.
j=0
Ak < ,
k=1
k1
j=0
k=1
and hence
k1
k=1
a.s.
j=0
Yt = Bt + At Bt1 +
k1
k=2 j=0
Atj )Btk
508
= Bt + At Bt1 +
l1
At1l )Bt1l
l=1 j=0
= Bt + At Yt1 .
Lemma 2.2. If (2.3) admits a strictly stationary solution with finite first mo
ment, then pi=1 i + qj=1 j < 1. Moreover, the strictly stationary solution of
(2.3) is unique.
Proof. By (2.3), we have
Y0 = A0 Y1 + B0
= B0 + A0 B1 + A0 A1 Y2
..
.
= B0 +
n1
k1
k=1 j=0
(2.5)
n1
Aj )Bk + (
Aj )Yn .
j=0
n1
k=1
k1
j=0
Aj Bk =
n1
Ak B.
k=1
Ak B < .
k=1
Therefore,
lim An B = 0.
(2.6)
(2.7)
then (2.7) implies that limn An = 0, which again implies that (A) < 1. As
we showed before the later is equivalent to pi=1 i + qj=1 j < 1, which leads
to the rst part of this lemma.
In fact, since B = 0 (1 + p+1 ) and 0 < 0 < , by (2.6) and the denition
of matrix A, (2.7) holds for i = 1 and i = p + 1. Again by the denition of A
and B, Ap+q = q (1 + p+1 ). If q = 0, then Ap+q = 0, hence (2.7) holds. If
q > 0, by the above equalities,
lim An p+q = lim An1 q (1 + p+1 ) = 0.
509
where x is a random variable.We need the Kronecker product (), the direct
operations (vec operations ), the notation Am = A A A, and the
basic identity vec(ABC) = (C A)vec(B). We denote the (i, j)th element of a
matrix D by (D)ij and dene m = E(At m ).
Theorem 3.1. Let {xt } be specified to be the strictly stationary solution of model
(1.1) and pi=1 i + qj=1 j < 1.
(i) If E|t |4 < and (2 ) < 1, then |xt |2 L2 .
(ii) If E|t |4(s1) < , for some integer s > 2, and (s ) < 1, then |xt |2 Ls .
Proof. Let Yt be generated according to (2.3) with the starting value Y0 = 0.
Let Y be a random variable having the same distribution as that of (2.4) with
t = 0, which is the (marginal) distribution of the unique stationary solution of
(2.3). It is clear that Yt Y in distribution. Let (Y ) be a random variable.
From weak convergence theory (Billingsley (1968)), it is known that to show
E(Y ) < , it suces to show lim inf E(Yt ) < . Let V (t) = E(Yt ). Then,
taking expectation on both sides of (2.3), we get
V (t) = AV (t 1) + B.
510
It is well know that lim V (t) exists and is nite if the spectral radius of A is less
than 1(cf. Subba Rao (1981)). Hence if the spectral radius of A is less than 1,
{V (t)} is bounded. Next, let V1 (t) = E(vec(Yt Yt )). Then
V1 (t) = E(At At )V1 (t 1) + E(Bt At )V1 (t 1)
+E(At Bt )V1 (t 1) + vec(E(Bt Bt )).
Note that the matices E(At At ), E(Bt At ), E(At Bt ) and vec(E(Bt Bt ))
are constant and nite matrices. As {V (t)} is bounded, it is clear that lim V1 (t)
exists and is nite if the spectral radius of 2 = E(At At ) is less than 1. Note
x4t ). Because x
4t converges in distribution to
that the rst element of V1 (t) is E(
4
4
x0 and lim V1 (t) is nite, E(x0 ) < . This completes the proof of (i). The proof
of (ii) is similar and hence omitted.
Acknowledgement
We are very grateful to referees for their careful reading and valuable suggestions which have enabled us to improve the earlier version of this paper. The
work was supported by NNCF of China and Prob. Lab. Inst. Applied Math.
Academia Sinica.
References
Billingsley, P. (1968). Convergence of Probability Measure. Wiley, New York.
Bollerslev, T. (1986). Generalized autoregressive conditional heteroscedasticity. J. Econometrics 31, 307-327.
Bollerslev, T., Chou, R. Y. and Kroner, K. F. (1992). ARCH modeling in finance. J. Econometrics 52, 5-59.
Bougeral, P. and Picard, N. (1992). Stationarity of GARCH processes and of some nonnegative
time series. J. Econometrics 52, 115-127.
Engle, R. F. (1982). Autoregressive conditional heteroscedasticity with estimates of the variance
of U. K. inflation. Econometrica 50, 987-1007.
Guegan, D. and Diebolt, J. (1994). Probabilistic properties of the -ARCH model. Statist.
Sinica 4, 71-87.
Horn, R. A. and Johnson, C. R. (1985). Matrix Analysis. Cambridge Univ. Press.
Subba Rao, T. (1981). On the theory of bilinear time series models. J. Roy. Statist. Soc. Ser.
B 43, 244-255.
Institute of Applied Mathematics, Academia Sinica, Beijing, 100080.
E-mail: mchen@iss01.iss.ac.cn
E-mail: hzan@sun.ihep.ac.cn
(Received October 1995; accetped March 1997)