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
, including local vicinities of unity with condition
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
as
.
In this work, as an extension of [
15], we establish uniform limit theory for the LSEs of autoregression coefficient
in the stationary AR(1) model with heavy-tailed G-GARCH(1,1) errors, for all
, but focusing on
as
. The asymptotic distributions of the LSEs, involved with the stable distributions depending on the tail index, are established for stationary
that includes neighborhood of the unity,
, but any condition on the rate of
is not required. The tail index
of the G-GARCH(1,1) process is considered in range
in this work. Our limit theory shows that, for the tail index
, the LSE is inconsistent, for
,
-consistent and for
,
-consistent. This work might have the novelty that the condition
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
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
defined by
with initialization
and coefficient
, where
is generated by the general GARCH(1,1) process, G-GARCH(1,1):
where
,
,
,
and
are nonnegative functions,
is a sequence of i.i.d. symmetric random variables with
,
. Formally,
is a triangular array but
n is omitted for notational simplicity. Throughout the paper, we make the following assumptions:
- (A1)
.
- (A2)
There exists a such that , , and , where .
- (A3)
The density of 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
as
.
For
, model (
1) is stationary and has the following expansion:
Given observations
, the LSE
of
in model (
1) and their difference are respectively given by
The limiting distribution of
is obtained by the asymptotic behaviors of
and
. By (
3),
and
can be approximated by
respectively, as
. We need to see the behavior of the heavy-tailed G-GARCH process
as well as the convergence of
, for
. We first give the following lemma for the tail index of
, presented by [
13].
Lemma 1 (Zhang and Ling (2015) [
13])
. Under Conditions (A1)–(A3), there exists a unique such that and where In (
5) or in the following,
means that
.
For tail index
in (
5), this work considers three cases with different asymptotic results:
,
and
. As a case of
in the simple heavy-tailed GARCH noise, but not G-GARCH, Ref. [
22] recently established a Gaussian limit theory as
converges to the unity with condition
as
. This present work excludes case
in the G-GRACH model
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
, let the
-dimensional stable random vector with stable index
be denoted by
For any given integers
l and
H, we define two
-dimensional random vectors:
where
Under the conditions in (A1)–(A3),
is a varying regular random vector sequence with index
. Let
where
is the limit of sum of all components in
divided by
, and it is independent of
l by the stationarity. Note that
exists by Equation (3.5) of [
13]. According to Theorem 2.8 and Proposition 3.3 of [
21], we have
where
denotes convergence in distribution. The limiting stable distribution in (
9) has somewhat different forms for
and for
, 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
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 as , where and is a -dimensional stable random vector with stable index . 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
for
. The limit theory is uniform over all stationary values of
, including
, on which this present work focuses. The convergence rate of
depends on
, but any condition of the rate of
is not required. As seen in the asymptotic results below, if
, the LSE is inconsistent as the limiting represents a ratio of
-stable random variables, like the estimators of autocorrelation function of the heavy-tailed process with
in [
21]. If
, the LSE is
-consistent, while, if
, it is
-consistent. We now present our main results as follows:
Theorem 1. Assume (A1)–(A3) and let . If as , then we have Theorem 2. Assume (A1)–(A3) and let . If as , then we have Theorem 3. Assume (A1)–(A3) and let . If as , then we have To derive the asymptotic results in Theorems above, we need to observe the joint limiting distribution of the numerator and the denominator of
in (
4). In the following lemma, the joint asymptotic distribution is established, which unify all the three cases considered:
and
.
Lemma 3. We consider the AR(1) model with coefficient in (1) and assume (A1)–(A3). For and for given as in (6), if as , then we havewhere 5. Concluding Remarks
This work considered an AR-G-GARCH model with tail index
and established limit theory for LSE of autoregression coefficient
with condition
as
. For the tail index
of G-GARCH errors, asymptotic distributions of the LSEs are involved with the stable distribution. The convergence rate of the LSE depends on
, but no condition on the rate of
is required. It is shown that, for the tail index
, the LSE is inconsistent, for
,
-consistent, and for
,
-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
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
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
, 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
, the asymptotic result for the LSE in model (
1) provides a Gaussian limit theory with convergence rate
. Its proof for case
is completely different from arguments in this work dealing with the cases
, and might be so lengthy as in [
22]. Thus, the Gaussian limit theory for
remains for future study at this moment.
As the convergence
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;
with assumption
as
, where
c and
d are nonzero constants, and
is a sequence of i.i.d. errors. The limit theory has been established for mildly stationary case
and for mildly explosive case
. 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
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.