Next Article in Journal
Statistical Tests of Symbolic Dynamics
Next Article in Special Issue
On Approximation of the Tails of the Binomial Distribution with These of the Poisson Law
Previous Article in Journal
Digital Image Processing Method for Characterization of Fractures, Fragments, and Particles of Soil/Rock-Like Materials
Previous Article in Special Issue
Chebyshev–Edgeworth-Type Approximations for Statistics Based on Samples with Random Sizes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Limit Theory for Stationary Autoregression with Heavy-Tailed Augmented GARCH Innovations

Department of Applied Statistics, Gachon University, Seongnam 13120, Korea
Mathematics 2021, 9(8), 816; https://doi.org/10.3390/math9080816
Submission received: 24 February 2021 / Revised: 7 April 2021 / Accepted: 8 April 2021 / Published: 9 April 2021
(This article belongs to the Special Issue Analytical Methods and Convergence in Probability with Applications)

Abstract

:
This paper considers stationary autoregressive (AR) models with heavy-tailed, general GARCH (G-GARCH) or augmented GARCH noises. Limit theory for the least squares estimator (LSE) of autoregression coefficient ρ = ρ n is derived uniformly over stationary values in [ 0 , 1 ) , focusing on ρ n 1 as sample size n tends to infinity. For tail index α ( 0 , 4 ) of G-GARCH innovations, asymptotic distributions of the LSEs are established, which are involved with the stable distribution. The convergence rate of the LSE depends on 1 ρ n 2 , but no condition on the rate of ρ n is required. It is shown that, for the tail index α ( 0 , 2 ) , the LSE is inconsistent, for α = 2 , log n / ( 1 ρ n 2 ) -consistent, and for α ( 2 , 4 ) , n 1 2 / α / ( 1 ρ n 2 ) -consistent. Proofs are based on the point process and the asymptotic properties in AR models with G-GARCH errors. However, this present work provides a bridge between pure stationary and unit-root processes. This paper extends the existing uniform limit theory with three issues: the errors have conditional heteroscedastic variance; the errors are heavy-tailed with tail index α ( 0 , 4 ) ; and no restriction on the rate of ρ n is necessary.

1. Introduction

Autoregressive (AR) time series models are not only extensively used in econometrics and financial markets but also their theories are applicable to various time series models. Inference on autoregressive coefficient ρ in an AR model of order one, or more generally, on the sum of coefficients in an AR model of order p, is an important issue in times series analysis. One of the most common inferences in the AR models is given by the least squares estimator (LSE), and its limit theory has been developed gradually in time series studies from a simple AR model with i.i.d. errors to an advanced AR model with various types of errors, for example, Martingale difference errors or conditional heteroscedastic errors.
We consider a stationary AR(1) model with conditional heteroscedastic errors. In particular, the error model adopted in this work is a more general form than the autoregressive conditional heteroscedastic (ARCH) model of Engle [1] and the generalized ARCH (GARCH) model of Bollerslev [2]. Since the introduction of the (G)ARCH model of Engle [1] and Bollerslev [2], many time series models in the GARCH family have been proposed in modeling financial time series data with applications to financial volatility: for instance, we refer to multiplicative GARCH in [3], exponential GARCH in [4], nonlinear asymmetric GARCH in [5], GJR-GARCH in [6], threshold GARCH in [7], etc.
To unify and generalize such parametric GARCH models, Ref. [8] defined the so-called augmented GARCH process and established its diffusion limit. Independently, Ref. [9] suggested a general GARCH (G-GARCH) process with nonparametric functionals and studied its moment properties, which is a first-order simple form of the augmented GARCH model of [8]. Ref. [10] presented strong approximation for the sum of squares of augmented GARCH sequences with applications to change-point analysis, and Ref. [11] developed a functional central limit theorem for a polynomial form of the augmented GARCH model. Ref. [12] studied the probabilistic structure of the augmented GARCH(1,1) sequences and asymptotic distribution of their various functionals. Meanwhile, Ref. [13] dealt with asymptotic inference for AR models with heavy-tailed G-GARCH noises, and Ref. [14] studied weak convergence of renorming volatilities in the G-GARCH models.
In the mean time, Ref. [15] developed the Gaussian limit theory of the least squares estimator (LSE) for a stationary autoregression with Martingale difference errors having homoscedastic variance, which holds uniformly over stationary values of the autoregression coefficient ρ n , including local vicinities of unity with condition ( 1 ρ n ) n as sample size n tends to infinity. In this paper, we provide a limit theory for the stationary autoregression, like [15], but by allowing for the errors to have conditional heteroscedastic variance. For the error model of the stationary AR(1) model considered in the present work, the heavy-tailed augmented GARCH (or G-GARCH process) is adopted.
Several authors studied analysis and inference for time series models with GARCH errors. For example, Ref. [16] provided a review of theoretical results for time series models with GARCH errors. The LSEs for ARMA-ARCH models have been done by [17,18] under the finite fourth moment conditions, and their results can be easily extended to the ARMA-GARCH models. Ref. [19] dealt with the OLS estimation for AR-GARCH models under non-finite fourth moments, and Ref. [20] considered the likelihood-based inference for AR-ARCH models. For the case of heavy-tailed GARCH errors, it would be expected from the results of [21] that the LSE is inconsistent if the variance of GARCH process is infinite but is consistent with a slower convergence rate if it has finite variance and an infinite fourth moment. This has been theoretically verified by [13] for the AR models with the heavy-tailed G-GARCH errors. The notable work of [13] make us extend to limit theory combined with an interesting result of [15], who dealt with the theory on condition ( 1 ρ n ) n as n .
In this work, as an extension of [15], we establish uniform limit theory for the LSEs of autoregression coefficient ρ n in the stationary AR(1) model with heavy-tailed G-GARCH(1,1) errors, for all ρ n [ 0 , 1 ) , but focusing on ρ n 1 as n . The asymptotic distributions of the LSEs, involved with the stable distributions depending on the tail index, are established for stationary ρ n that includes neighborhood of the unity, ρ n 1 , but any condition on the rate of ρ n is not required. The tail index α of the G-GARCH(1,1) process is considered in range ( 0 , 4 ) in this work. Our limit theory shows that, for the tail index α ( 0 , 2 ) , the LSE is inconsistent, for α = 2 , log n / ( 1 ρ n 2 ) -consistent and for α ( 2 , 4 ) , n 1 2 / α / ( 1 ρ n 2 ) -consistent. This work might have the novelty that the condition ρ n 1 plays a role of bridge between pure stationary and unit-root processes. In addition, this paper is a notable extension by adding three issues to the work of [15]: first, the errors have conditional heteroscedastic variance; second, the errors have heavy-tailed feature with specific tail indices; and, finally, any condition on the rate of ρ n is not required.
The remainder of the paper is organized as follows: in Section 2, the model, assumptions, and existing theory are described. Our main results are presented in Section 3 while their proofs are given in Section 4. Concluding remarks are stated in Section 5.

2. Assumption and Existing Theory

We consider an AR(1) process { y t : t = 1 , 2 , , n } defined by
y t = ρ n y t 1 + ϵ t
with initialization y 0 and coefficient ρ = ρ n [ 0 , 1 ) , where { ϵ t : t = 1 , 2 , } is generated by the general GARCH(1,1) process, G-GARCH(1,1):
ϵ t = h t η t and h t δ = g ( η t 1 ) + c ( η t 1 ) h t 1 δ
where δ > 0 , P { h t δ > 0 } = 1 , c ( 0 ) < 1 , c ( · ) and g ( · ) are nonnegative functions, { η t } is a sequence of i.i.d. symmetric random variables with E η 1 = 0 , E η 1 2 = 1 . Formally, { y t = y t , n : t = 1 , 2 , , n } is a triangular array but n is omitted for notational simplicity. Throughout the paper, we make the following assumptions:
(A1)
E log ( c ( η t ) ) < 0 .
(A2)
There exists a k 0 > 0 such that E ( c ( η t ) ) k 0 1 , E [ ( c ( η t ) ) k 0 log + ( c ( η t ) ) ] < , and E ( g ( η t ) + | η t | δ ) k 0 < , where log + ( x ) = max { 0 , log ( x ) } .
(A3)
The density f ( x ) of η 1 is positive in the neighborhood of zero.
Conditions (A1)–(A3) are adopted from [13], who derived asymptotic theory for AR models with the heavy-tailed G-GARCH noises. The remarkable work of [13] led us to extend an interesting limit theory along with one of [15], who studied uniform limit theory with condition ( 1 ρ n ) n as n .
For ρ = ρ n [ 0 , 1 ) , model (1) is stationary and has the following expansion:
y t = l = 0 ρ n l ϵ t l .
Given observations { y 0 , y 1 , , y n } , the LSE ρ ^ n of ρ n in model (1) and their difference are respectively given by
ρ ^ n = t = 1 n y t y t 1 t = 1 n y t 1 2 and ρ ^ n ρ n = t = 1 n y t 1 ϵ t t = 1 n y t 1 2 .
The limiting distribution of ρ ^ n ρ n is obtained by the asymptotic behaviors of t = 1 n y t 1 ϵ t and t = 1 n y t 1 2 . By (3), t = 1 n y t 1 ϵ t and t = 1 n y t 1 2 can be approximated by
l = 0 H ρ n l t = 1 n ϵ t 1 l ϵ t and l = 0 H m = 0 H ρ n l + m t = 1 n ϵ t 1 l ϵ t 1 m ,
respectively, as H . We need to see the behavior of the heavy-tailed G-GARCH process { ϵ t } as well as the convergence of t = 1 n ϵ t l ϵ t m , for l , m { 1 , 2 , , H } . We first give the following lemma for the tail index of { ϵ t } , presented by [13].
Lemma 1
(Zhang and Ling (2015) [13]). Under Conditions (A1)–(A3), there exists a unique α ( 0 , δ k 0 ] such that E ( c ( η t ) ) α / δ = 1 and
P ( | ϵ 1 | > x ) c 0 E | η 1 | α x α ,
where
c 0 = c 0 ( α ) = E [ { g ( η 1 ) + c ( η 1 ) h 1 δ } α / δ { c ( η 1 ) h 1 δ } α / δ ] α E [ c ( η 1 ) α / δ log + ( c ( η 1 ) ) ] .
In (5) or in the following, a n c b n means that lim n a n / b n = c .
For tail index α in (5), this work considers three cases with different asymptotic results: 0 < α < 2 , α = 2 and 2 < α < 4 . As a case of α = 4 in the simple heavy-tailed GARCH noise, but not G-GARCH, Ref. [22] recently established a Gaussian limit theory as ρ n converges to the unity with condition ( 1 ρ n ) n as n . This present work excludes case α = 4 in the G-GRACH model { ϵ t } because the asymptotic result is not a stable limit, which is a main focus of this work, but it is a Gaussian limit with such a lengthy analysis as in [22]. This work is an extension of [15], and, moreover, our limit theory plays a role of a bridge between the inferences in [13,23], who have dealt with unit-root and stationary AR models with heavy-tailed GARCH errors, respectively.
We here summarize the discussion in [13] along with results of [21]. The same notations as those in [13] are used, but some complicated notations are referred to in their paper.
For any positive integer H > 0 , let the ( H + 1 ) -dimensional stable random vector with stable index α / 2 ( 0 , 2 ] be denoted by
S H ( α ) : = ( S 0 ( α ) , S 1 ( α ) , , S H ( α ) ) .
For any given integers l and H, we define two ( H + 1 ) -dimensional random vectors:
E t , l , H ( 1 ) = ( ϵ t l , ϵ t l 1 , , ϵ t l H ) , E t , l , H ( 2 ) = ϵ t l 2 c n , ϵ t l ϵ t l 1 , , ϵ t l ϵ t l H ,
where
c n = c n ( α ) = 0 for 0 < α < 2 E ϵ 1 2 I ( | ϵ 1 | c 0 n ) = c 0 log n for α = 2 E ϵ 1 2 for 2 < α < 4 .
Under the conditions in (A1)–(A3), E t , l , H ( 1 ) , t 1 is a varying regular random vector sequence with index α . Let
U ( α ) = lim H lim n 1 a n 2 t = 1 n ( ϵ t l 2 c n ) + h = 1 H t = 1 n ϵ t l ϵ t l h
where
a n = a n ( α ) = c 0 ( α ) E | η 1 | α n 1 / α
U ( α ) is the limit of sum of all components in t = 1 n E t , l , H ( 2 ) divided by a n 2 , and it is independent of l by the stationarity. Note that U ( α ) exists by Equation (3.5) of [13]. According to Theorem 2.8 and Proposition 3.3 of [21], we have
1 a n 2 t = 1 n E t , l , H ( 2 ) d ( S 0 ( α ) , S 1 ( α ) , , S H ( α ) ) = S H ( α ) ,
where d denotes convergence in distribution. The limiting stable distribution in (9) has somewhat different forms for 0 < α < 2 and for 2 α < 4 , whose explicit expressions are given in Equations (3.2) and (3.3) of [13], respectively. See [13,21] for its details. Note that the limiting distribution S H ( α ) in (9) does not depend on l. Ref. [13] has discussed the following result.
Lemma 2
(Zhang and Ling (2015) [13]). Under Conditions (A1)–(A3), for any positive integers K and H, we have
1 a n 2 t = 1 n E t , 0 , H ( 2 ) , t = 1 n E t , 1 , H ( 2 ) , , t = 1 n E t , K , H ( 2 ) d ( S H ( α ) , S H ( α ) , , S H ( α ) ) 1 × ( K + 1 )
as n , where 0 < α < 4 and S H is a ( H + 1 ) -dimensional stable random vector with stable index α / 2 ( 0 , 2 ) .
The next section presents our main limit theories with their proofs. The stable limit in Lemma 2 is used to derive our results.

3. Main Results

In this section, we provide the limit theory, involved with the stable distribution, of the LSE ρ ^ n for 0 < α < 4 . The limit theory is uniform over all stationary values of ρ = ρ n [ 0 , 1 ) , including ρ n 1 , on which this present work focuses. The convergence rate of ρ ^ n depends on 1 ρ n 2 , but any condition of the rate of ρ n 1 is not required. As seen in the asymptotic results below, if 0 < α < 2 , the LSE is inconsistent as the limiting represents a ratio of α / 2 -stable random variables, like the estimators of autocorrelation function of the heavy-tailed process with 0 < α < 2 in [21]. If α = 2 , the LSE is ( log n ) / ( 1 ρ n 2 ) -consistent, while, if 2 < α < 4 , it is n 1 2 / α / ( 1 ρ n 2 ) -consistent. We now present our main results as follows:
Theorem 1.
Assume (A1)–(A3) and let 0 < α < 2 . If ρ n 1 as n , then we have
1 1 ρ n 2 ( ρ ^ n ρ n ) d U ( α ) S 0 ( α ) + 2 U ( α ) .
Theorem 2.
Assume (A1)–(A3) and let α = 2 . If ρ n 1 as n , then we have
log n 1 ρ n 2 ( ρ ^ n ρ n ) d U ( α ) .
Theorem 3.
Assume (A1)–(A3) and let 2 < α < 4 . If ρ n 1 as n , then we have
n 1 2 / α 1 ρ n 2 ( ρ ^ n ρ n ) d ( c 0 E | η 1 | α ) 2 / α E ϵ 1 2 U ( α ) .
To derive the asymptotic results in Theorems above, we need to observe the joint limiting distribution of the numerator and the denominator of ρ ^ n ρ n in (4). In the following lemma, the joint asymptotic distribution is established, which unify all the three cases considered: 0 < α < 2 , α = 2 and 2 < α < 4 .
Lemma 3.
We consider the AR(1) model with coefficient ρ n in (1) and assume (A1)–(A3). For 0 < α < 4 and for c n given as in (6), if ρ n 1 as n , then we have
1 ρ n 2 n 2 / α t = 1 n y t 1 2 c n 1 ρ n 2 , 1 n 2 / α t = 1 n y t 1 ϵ t d Z α / 2 ( 1 ) , Z α / 2 ( 2 )
where
Z α / 2 ( 1 ) = ( c 0 ( α ) E | η 1 | α ) 2 / α S 0 ( α ) + 2 U ( α ) and Z α / 2 ( 2 ) = ( c 0 ( α ) E | η 1 | α ) 2 / α U ( α ) .

4. Proofs

This section gives proofs of the main results. Thanks to Lemma 3, which includes all the three cases of tail index α considered, the results in Theorems 1–3 are derived neatly and concisely, and thus we state the proof of Lemma 3 prior to those of Theorems 1–3. Indeed, the proof of Lemma 3 is a main contribution to the limit theory in this study.
Proof of Lemma 3.
To prove the desired limiting distribution, by Cramer’s Wold Device, we may show that, for any real numbers ( u , v ) ,
u 1 ρ n 2 n 2 / α t = 1 n y t 1 2 c n 1 ρ n 2 + v 1 n 2 / α t = 1 n y t 1 ϵ t d u Z α / 2 ( 1 ) + v Z α / 2 ( 2 ) .
Recalling (8), we will verify that A n ( 1 ) ( u ) + A n ( 2 ) ( v ) : =
u 1 ρ n 2 a n 2 t = 1 n y t 1 2 c n 1 ρ n 2 + v 1 a n 2 t = 1 n y t 1 ϵ t d u S 0 ( α ) + 2 U ( α ) + v U ( α ) .
First, for A n ( 1 ) ( u ) , we write
t = 1 n y t 1 2 c n 1 ρ n 2 = t = 1 n l = 0 m = 0 ρ n l + m ϵ t 1 l ϵ t 1 m c n 1 ρ n 2
= t = 1 n l = 0 H m = 0 H ρ n l + m ϵ t 1 l ϵ t 1 m c n 1 ρ n 2 + I t ( 1 ) ( H )
where I t ( 1 ) ( H ) is given as three double-summations as follows:
I t ( 1 ) ( H ) = l = H + 1 m = 0 H ρ n l + m ϵ t 1 l ϵ t 1 m +
l = 0 H m = H + 1 ρ n l + m ϵ t 1 l ϵ t 1 m + l = H + 1 m = H + 1 ρ n l + m ϵ t 1 l ϵ t 1 m .
Note that, for any δ > 0 , P 1 a n 2 t = 1 n I t ( 1 ) ( H ) > δ
P 1 a n 2 l = H + 1 ρ n l m = 0 H ρ n m t = 1 n ϵ t 1 l ϵ t 1 m > δ / 3
+ P 1 a n 2 l = 0 H ρ n l m = H + 1 ρ n m t = 1 n ϵ t 1 l ϵ t 1 m > δ / 3
+ P 1 a n 2 l = H + 1 ρ n l m = H + 1 ρ n m t = 1 n ϵ t 1 l ϵ t 1 m > δ / 3 0 as n , ( H )
and thus
1 a n 2 t = 1 n I t ( 1 ) ( H ) p 0 as n .
In addition, (11) is written as
t = 1 n l = 0 H ρ n 2 l ϵ t 1 l 2 c n 1 ρ n 2 + B t ( H ) + I t ( 1 ) ( H )
where B t ( H ) = l m ρ n l + m ϵ t 1 l ϵ t 1 m
= l > m ρ n l + m ϵ t 1 l ϵ t 1 m + l < m ρ n l + m ϵ t 1 l ϵ t 1 m = : B t ( 1 ) ( H ) + B t ( 2 ) ( H ) .
In the summation of B t ( 1 ) ( H ) , set j = l m , and in the summation of B t ( 2 ) ( H ) , set i = m l , then we get B t ( H ) = B t ( 1 ) ( H ) + B t ( 2 ) ( H ) =
= j = 1 H m = 0 H ρ n 2 m + j ϵ t 1 m ϵ t 1 m j + i = 1 H l = 0 H ρ n 2 l + i ϵ t 1 l ϵ t 1 l i
= 2 j = 1 H m = 0 H ρ n 2 m + j ϵ t 1 m ϵ t 1 m j .
The summation from t = 1 to t = n of the first two expressions in (13) is given by
t = 1 n l = 0 H ρ n 2 l ( ϵ t 1 l 2 c n ) + l = 0 H ρ n 2 l c n c n 1 ρ n 2
= l = 0 H ρ n 2 l t = 1 n ϵ t 1 l 2 c n + t = 1 n c n l = 0 H ρ n 2 l 1 1 ρ n 2
= l = 0 H ρ n 2 l t = 1 n ϵ t 1 l 2 c n + n c n ρ n 2 H + 2 1 ρ n 2 .
Therefore, we see that the first term A n ( 1 ) ( u ) of the left-hand side in (10) becomes, by (11)–(15),
u 1 ρ n 2 a n 2 l = 0 H ρ n 2 l t = 1 n ϵ t 1 l 2 c n + n c n ρ n 2 H + 2 1 ρ n 2 + 2 t = 1 n j = 1 H m = 0 H ρ n 2 m + j ϵ t 1 m ϵ t 1 m j + o p ( 1 )
Second, for A n ( 2 ) ( v ) , we write
t = 1 n y t 1 ϵ t = t = 1 n ϵ t l = 0 ρ n l ϵ t 1 l = l = 0 H ρ n l t = 1 n ϵ t ϵ t 1 l + t = 1 n I t ( 2 ) ( H )
where I t ( 2 ) ( H ) = l = H + 1 ρ n l ϵ t ϵ t 1 l . It is clear that 1 a n 2 t = 1 n I t ( 2 ) ( H ) p 0 as n , ( H ).
Thus, the second term A n ( 2 ) ( v ) of the left-hand side in (10) is
v 1 a n 2 l = 0 H ρ n l t = 1 n ϵ t ϵ t 1 l + o p ( 1 ) .
Hence, by (16) and (17), in (10), A n ( 1 ) ( u ) + A n ( 2 ) ( v ) =
u ( 1 ρ n 2 ) l = 0 H ρ n 2 l 1 a n 2 t = 1 n ϵ t 1 l 2 c n + 2 j = 1 H m = 0 H ρ n 2 m + j 1 a n 2 t = 1 n ϵ t 1 m ϵ t 1 m j
+ v l = 0 H ρ n l 1 a n 2 t = 1 n ϵ t ϵ t 1 l + O n c n ρ n 2 H a n 2 + o p ( 1 ) .
Now, we choose H = H n , α , depending on n and α , such that H and O n c n ρ n 2 H / a n 2 0 as n and ρ n 1 . When 0 < α < 2 , c n = 0 and thus any sequence H may be taken. If α = 2 , then O n c n ρ n 2 H / a n 2 = O ( ρ n 2 H log n ) by (6). Hence, we choose H = log log n / log ρ n , the integer part of log log n / log ρ n . Since log ρ n < 0 and log ρ n 0 , H as n , and it follows that 2 H log ρ n = 2 log log n = log ( log n ) 2 and thus ρ n 2 H log n = 1 / log n 0 . Therefore, on behalf of the following discussion, if 0 < α 2 , then we choose
H = log log n log ρ n .
If 2 < α < 4 , then we choose
H = 1 2 α log n log ρ n ,
which tends to and satisfies O n c n ρ n 2 H / a n 2 = O ( n 1 2 / α ρ n 2 H ) = O ( 1 / n 1 2 / α ) 0 .
By (9) and Lemma 2, Ref. (18) has the same limiting distribution as
u ( 1 ρ n 2 ) l = 0 H ρ n 2 l S 0 ( α ) + 2 m = 0 H ρ n 2 m j = 1 H ρ n j S j ( α ) + v l = 0 H ρ n l S l + 1 ( α ) .
We see that ( 1 ρ n 2 ) l = 0 H ρ n 2 l = 1 ρ n H + 2 1 as ρ n 1 and H , since log ρ n H = H log ρ n for chosen H in both (19) and (20), and thus ρ n H 0 . Letting j = l + 1 in the third term of (21), we see that, as H , (21) tends to u S 0 ( α ) + 2 j = 1 S j ( α ) + v j = 1 S j ( α ) . By (7) and (9), we obtain the desired result in (10). □
As mentioned above, the main contribution of the limit theory in the present paper is the proof of Lemma 3, which is somewhat lengthy and challenging. Due to the combination of Lemma 3 and Slutsky’s theorem, we can see simple validations of Theorems 1–3. When we prove Theorems 1–3 in the following, Slutsky’s theorem is used.
Proof of Theorem 1.
For 0 < α < 2 , c n = 0 and
1 1 ρ n 2 ( ρ ^ n ρ n ) = t = 1 n y t 1 ϵ t / n 2 α ( 1 ρ n 2 ) t = 1 n y t 1 2 / n 2 α d Z α / 2 ( 2 ) Z α / 2 ( 1 ) = U ( α ) S 0 ( α ) + 2 U ( α )
by Lemma 3, the desired result holds. □
Proof of Theorem 2.
For α = 2 , c n = c 0 log n , and we write
log n 1 ρ n 2 ( ρ ^ n ρ n ) = 1 n t = 1 n y t 1 ϵ t 1 ρ n 2 n log n t = 1 n y t 1 2 .
The first component of the left side in Lemma 3 implies
1 ρ n 2 n t = 1 n y t 1 2 c n 1 ρ n 2 d Z α / 2 ( 1 )
and thus
1 ρ n 2 n log n t = 1 n y t 1 2 = Z α / 2 ( 1 ) log n + c n log n + o p ( 1 ) .
By Lemma 3, we have
1 n t = 1 n y t 1 ϵ t 1 ρ n 2 n log n t = 1 n y t 1 2 = Z α / 2 ( 2 ) Z α / 2 ( 1 ) / log n + c 0 + o p ( 1 ) + o p ( 1 ) d 1 c 0 Z α / 2 ( 2 ) = U ( α ) ,
noting E | η 1 | 2 = 1 . The desired result is obtained. □
Proof of Theorem 3.
For 2 < α < 4 , c n = E ϵ 1 2 . We see
n 1 2 / α 1 ρ n 2 ( ρ ^ n ρ n ) = n 1 2 / α t = 1 n y t 1 ϵ t / n 2 α ( 1 ρ n 2 ) t = 1 n y t 1 2 / n 2 α .
By the first component of the left side in Lemma 3, we have
1 ρ n 2 n 2 / α t = 1 n y t 1 2 = Z α / 2 ( 1 ) + c n n 1 2 / α + o p ( 1 )
and thus (22) is
n 1 2 / α Z α / 2 ( 2 ) Z α / 2 ( 1 ) + c n n 1 2 / α + o p ( 1 ) + o p ( 1 ) = Z α / 2 ( 2 ) Z α / 2 ( 1 ) / n 1 2 / α + E ϵ 1 2 + o p ( 1 ) + o p ( 1 ) d 1 E ϵ 1 2 Z α / 2 ( 2 ) .
The desired result is obtained. □

5. Concluding Remarks

This work considered an AR-G-GARCH model with tail index α ( 0 , 4 ) and established limit theory for LSE of autoregression coefficient ρ = ρ n [ 0 , 1 ) with condition ρ n 1 as n . For the tail index α ( 0 , 4 ) of G-GARCH errors, asymptotic distributions of the LSEs are involved with the stable distribution. The convergence rate of the LSE depends on 1 ρ n 2 , but no condition on the rate of ρ n is required. It is shown that, for the tail index α ( 0 , 2 ) , the LSE is inconsistent, for α = 2 , log n / ( 1 ρ n 2 ) -consistent, and for α ( 2 , 4 ) , n 1 2 / α / ( 1 ρ n 2 ) -consistent. This paper extends the uniform limit theory of [15] with three issues: the errors have conditional heteroscedastic variance; the errors are heavy-tailed; and no restriction on the rate of ρ n is necessary. The results can be applicable to financial volatility analysis. Volatility reveals heavy-tail and long-memory features as seen in [24], who discussed volatility of realized volatility. Ref. [24] claimed that adopting a GARCH process in the errors improves the goodness-of-fit, especially for realized volatility with non-Gaussianity and heavy-tails. Empirical analysis with financial market data, based on the results in this work, would be interesting and challenging.
This work focuses on convergence condition of ρ n 1 but similar derivations for fixed coefficients ρ can be done with results of uniform limit theory, which is a special case of [13]. If tail index α = 4 , Ref. [22] investigated the limit theory of the LSE of AR(1)-GARCH(1,1) models, which is a simple linear case of (2). For the case of G-GARCH process in (2) with tail index α = 4 , the asymptotic result for the LSE in model (1) provides a Gaussian limit theory with convergence rate ( n / log n ) 1 / 2 / ( 1 ρ n 2 ) . Its proof for case α = 4 is completely different from arguments in this work dealing with the cases 0 < α < 4 , and might be so lengthy as in [22]. Thus, the Gaussian limit theory for α = 4 remains for future study at this moment.
As the convergence ρ n 1 with a slower rate is accounted for, the asymptotic results provide a bridge between pure stationary or explosive process and unit root (or local to unity) process as mentioned by [25], who has developed limit theory for a mildly integrated process with intercept; y t = d + ρ n y t 1 + u t with assumption ρ n = 1 + c / k n , k n = o ( n ) as n , where c and d are nonzero constants, and { u t } is a sequence of i.i.d. errors. The limit theory has been established for mildly stationary case c < 0 and for mildly explosive case c > 0 . Earlier, Ref. [26] developed limit theory for a mildly integrated process, including both mildly stationary and mildly explosive cases. The extension of the mildly integrated process with nonzero intercept by adopting the heavy-tailed G-GARCH errors with tail index α ( 0 , 4 ) will be an interesting study. It will make a bridge between the inferences in [13,23], who have dealt with, respectively, unit-root and stationary AR models with heavy-tailed G-GARCH errors. This extension will be considered as future research.

Funding

This work was supported by the National Research Foundation of Korea (NRF-2018R1D1A1 B07048745).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Engle, R.F. Autoregressive conditional heteroscedasticity with estimates of the variance of U.K. inflation. Econometrica 1982, 50, 987–1008. [Google Scholar] [CrossRef]
  2. Bollerslev, T. Generalized autoregressive conditional heteroskedasticity. J. Econom. 1986, 31, 307–327. [Google Scholar] [CrossRef] [Green Version]
  3. Geweke, J. Comment: Modeling the persistence of conditional variances. Econom. Rev. 1986, 5, 57–61. [Google Scholar] [CrossRef]
  4. Nelson, D. Conditional heteroskedasticity in asset returns: A new approach. Econometrica 1991, 59, 347–370. [Google Scholar] [CrossRef]
  5. Engle, R.F.; Ng, V. Measuring and testing the impact of news on volatility. J. Financ. 1993, 48, 1749–1778. [Google Scholar] [CrossRef]
  6. Glosten, L.; Jagannathan, R.; Runkle, D. Relationship between the expected value and the volatility of the nomial excess return on stocks. J. Financ. 1993, 48, 1779–1801. [Google Scholar] [CrossRef]
  7. Zakoian, J.M. Threshold heteroskedastic models. J. Econ. Dyn. Control 1994, 18, 931–955. [Google Scholar] [CrossRef]
  8. Duan, J. Augmented GARCH(p.q) process and its diffusion limit. J. Econom. 1997, 79, 97–127. [Google Scholar] [CrossRef]
  9. He, C.; Teräsvirta, T. Properties of moments of a family of GARCH processes. J. Econom. 1999, 92, 173–192. [Google Scholar] [CrossRef]
  10. Aue, A.; Berkes, I.; Horvath, L. Strong approximations for the sums of squares of argmented GARCH sequences. Bernoulli 2006, 12, 583–608. [Google Scholar] [CrossRef]
  11. Berkes, I.; Hörmann, S.; Horvath, L. The functional central limit theorem for a family of GARCH observations with applications. Stat. Probablity Lett. 2008, 78, 2725–2730. [Google Scholar] [CrossRef] [Green Version]
  12. Hörmann, S. Augmented GARCH sequences: Dependence structure and asymptotics. Bernoulli 2008, 14, 543–561. [Google Scholar] [CrossRef]
  13. Zhang, R.; Ling, S. Asymptotic inference for AR models with heavy-tailed G-GARCH noises. Econom. Theory 2015, 31, 880–890. [Google Scholar] [CrossRef]
  14. Li, D.; Wu, W. Renorming volatilities in a family of GARCH models. Econom. Theory 2018, 34, 1370–1382. [Google Scholar] [CrossRef]
  15. Giraitis, L.; Philips, P.C.B. Unifrom limit theory for stationary autoregression. J. Time Ser. Anal. 2006, 27, 51–60, Corrigendum in J. Time Ser. Anal. 2006, 27, i–ii. [Google Scholar] [CrossRef] [Green Version]
  16. Li, W.K.; Ling, S.; MacAleer, M. A Survey of Recent Theoretical Results for Time Series Models with GARCH Errors; ISER Discussion Paper; Institue of Social and Economic Research, Osaka University: Osaka, Japan, 2001; p. 545. [Google Scholar]
  17. Weiss, A.A. Asymptotic theory for ARCH models: Estimation and testing. Econom. Theory 1986, 2, 107–131. [Google Scholar] [CrossRef] [Green Version]
  18. Pantula, S.G. Estimation of autoregressive models with ARCH errors. Sankhya B 1989, 50, 119–138. [Google Scholar]
  19. Lange, T. Tail behavior and OLS estimation in AR-GARCH models. Stat. Sin. 2011, 21, 1191–1200. [Google Scholar] [CrossRef] [Green Version]
  20. Lange, T.; Rahbek, A.; Jensen, S.T. Estimation and asymptotic inference in the AR-ARCH model. Econom. Rev. 2011, 30, 129–153. [Google Scholar] [CrossRef]
  21. Davis, R.A.; Mikosch, M. The sample autocorrelations of heavy-tailed process with applications to ARCH. Ann. Stat. 1998, 26, 2049–2080. [Google Scholar] [CrossRef]
  22. Hwang, E. A note on limit theory for mildly stationary autoregression with a heavy-tailed GARCH error process. Stat. Probab. Lett. 2019, 152, 59–68. [Google Scholar] [CrossRef]
  23. Chan, N.H.; Zhang, R.M. Inference for unit-root models with infinite variance GARCH errors. Stat. Sin. 2010, 20, 1363–1393. [Google Scholar]
  24. Corsi, F.; Mittnik, S.; Pigorsch, C.; Pigorsch, U. The volatility of realized volatility. Econom. Rev. 2008, 27, 46–78. [Google Scholar] [CrossRef]
  25. Fei, Y. Limit theory for mildly integrated process with intercept. Econ. Lett. 2018, 163, 98–101. [Google Scholar] [CrossRef]
  26. Phillips, P.C.B.; Magdalinos, T. Limit theory for moderate deviations from a unit root. J. Econom. 2007, 136, 115–130. [Google Scholar] [CrossRef] [Green Version]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Hwang, E. Limit Theory for Stationary Autoregression with Heavy-Tailed Augmented GARCH Innovations. Mathematics 2021, 9, 816. https://doi.org/10.3390/math9080816

AMA Style

Hwang E. Limit Theory for Stationary Autoregression with Heavy-Tailed Augmented GARCH Innovations. Mathematics. 2021; 9(8):816. https://doi.org/10.3390/math9080816

Chicago/Turabian Style

Hwang, Eunju. 2021. "Limit Theory for Stationary Autoregression with Heavy-Tailed Augmented GARCH Innovations" Mathematics 9, no. 8: 816. https://doi.org/10.3390/math9080816

APA Style

Hwang, E. (2021). Limit Theory for Stationary Autoregression with Heavy-Tailed Augmented GARCH Innovations. Mathematics, 9(8), 816. https://doi.org/10.3390/math9080816

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop