Asymptotic theory for quadratic variation of harmonizable fractional stable processes
Andreas Basse-O'Connor, Mark Podolskij

TL;DR
This paper develops the asymptotic theory for the quadratic variation of harmonizable fractional stable processes, revealing non-ergodic limits and convergence to Lévy-driven Rosenblatt variables under specific conditions.
Contribution
It provides new asymptotic results for quadratic variation of harmonizable fractional stable processes, including non-ergodic limits and weak convergence to Rosenblatt distributions.
Findings
Law of large numbers with non-ergodic limit
Weak convergence to Lévy-driven Rosenblatt variable
Conditions on Hurst parameter and stability index
Abstract
In this paper we study the asymptotic theory for quadratic variation of a harmonizable fractional -stable process. We show a law of large numbers with a non-ergodic limit and obtain weak convergence towards a L\'evy-driven Rosenblatt random variable when the Hurst parameter satisfies and . This result complements the asymptotic theory for fractional stable processes investigated in e.g. \cite{BHP19,BLP17,BP17,BPT20,LP18,MOP20}.
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Taxonomy
TopicsStochastic processes and statistical mechanics · Nonlinear Partial Differential Equations · Stochastic processes and financial applications
Asymptotic theory for quadratic variation of
harmonizable fractional stable processes
Andreas Basse-O’Connor Department of Mathematics, University of Aarhus, E-mail: [email protected].
Mark Podolskij Department of Mathematics, University of Luxembourg, E-mail: [email protected] Podolskij gratefully acknowledges financial support of ERC Consolidator Grant 815703 “STAMFORD: Statistical Methods for High Dimensional Diffusions”.
Abstract
In this paper we study the asymptotic theory for quadratic variation of a harmonizable fractional -stable process. We show a law of large numbers with a non-ergodic limit and obtain weak convergence towards a Lévy-driven Rosenblatt random variable when the Hurst parameter satisfies and . This result complements the asymptotic theory for fractional stable processes investigated in e.g. [4, 5, 6, 7, 17, 21].
Key words: fractional processes, harmonizable processes, limit theorems, quadratic variation, stable Lévy motion.
AMS 2010 subject classifications: 60F05, 60F15, 60G22, 60G48, 60H05.
1 Introduction
Since the seminal work of Mandelbrot and Van Ness [20] fractional stochastic processes and fields have gained a lot of attention in the probabilistic and statistical literature. One of the key results of [20] was the introduction of the fractional Brownian motion, the most famous fractional stochastic process. The authors show that, up to scaling, fractional Brownian motion with Hurst parameter is the unique self-similar centered Gaussian process with stationary increments. When the Gaussianity assumption is replaced by -stable distributions the aforementioned class becomes much richer: Rosinski [24] shows that it consists of (mixed) moving average -stable processes, -stable harmonizable processes and processes of a third kind (which lack an explicit description). The most well-known representative of the first class is the linear fractional stable motion (lfsm), whose probabilistic and statistical properties have been studied in numerous articles. We refer to [4, 5, 6, 7, 17, 22, 23] for a detailed exposition of limit theorems for functionals of lfsm, and to [1, 10, 12, 18, 19, 21] for statistical estimation of lfsm and related processes.
In the paper [9] the authors introduced the -stable harmonizable fractional motion, defined as
[TABLE]
where , is the Hurst parameter and , , denotes an isotropic complex-valued -stable Lévy motion with . The processes is self-similar with index , i.e. , and possesses stationary increments. Probabilistic properties of model (1.1) and their extensions have been studied in numerous articles, see e.g. [2, 8, 9, 11, 15] among others. However, to the best of our knowledge, asymptotic theory for statistics of harmonizable fractional -stable processes has not yet being investigated in the literature. In contrast to lfsm, harmonizable fractional -stable processes are not ergodic, which indicates potential difficulties when studying limit theorems.
This paper aims to derive asymptotic results for the quadratic statistic
[TABLE]
and our main statement is the following theorem:
Theorem 1.1**.**
Let be the -stable harmonizable fractional motion introduced in (1.1) with parameters and .
- (i)
As ,
[TABLE]
where denotes the quadratic variation of for . 2. (ii)
For and we obtain the weak convergence
[TABLE]
Remark 1.2**.**
(A) The result of Theorem 1.1(i) is in sharp contrast to the setting of short memory processes. To make this comparison concrete, let us consider a sequence of independent -stable isotropic random variables with . According to the Stable Central Limit Theorem, cf. [13, Theorem 5.25],
[TABLE]
where is an -stable random variable. That is, for independent -stable random variables, the quadratic variation is of the order , whereas is of the order , cf. Theorem 1.1(i). This difference may be viewed as an indication of the strong dependence between the random variables .
(B) We believe that the randomness of the limit in (1.3) is connected to the non-ergodicity of the process .
(C) The weak convergence result of (1.4) is similar in spirit to the Gaussian case. Indeed, the quadratic variation of the fractional Brownian motion with Hurst parameter is asymptotically distributed according to the Rosenblatt random variable when with convergence rate ; the classical Rosenblatt random variable possesses the same representation as the right hand side of (1.4) where is replaced by a complex-valued Gaussian measure. We also recall that there is no convergence regime of the form (1.4) in the setting of lfsm (cf. [4, 5]).
(D) The proof of Theorem 1.1 heavily relies on the quadratic form of the statistic. Extensions to more general functional forms and statistical estimation of the process are therefore far from obvious.
The rest of the paper is structured as follows. In Chapter 2, we discuss the integration theory and a dominated convergence result we will use in our proof. Chapter 3 is devoted to the proof of Theorem 1.1.
Notation
Throughout the paper all random variables are defined on a probability space . By we denote convergence in probability, and denotes convergence in distribution. For a complex number , denotes its complex conjugate and its norm. Furthermore, we denote by (resp. ) the real (resp. imaginary) part of .
2 Some background
In this section we will collect some important results on stochastic integration theory. A complex-valued random variable is called isotropic if it holds that for any . Similarly, a stochastic process is isotropic when any linear combination is isotropic. It is easy to see that a complex-valued Lévy process with characteristic triplet , where denotes the Lévy measure, is isotropic if and only if
[TABLE]
In particular, and both have the same symmetric -stable distribution. For any measurable function we define
[TABLE]
whenever the four real-valued integrals on the right-hand side exist. We note that the integral is well defined if and only if is well defined, which holds if . The next result, which follows by [25, Theorem 2], gives a criterion for the existence of the double integral:
Proposition 2.1**.**
Let be a measure function satisfying , and suppose that is a two-dimensional symmetric -stable Lévy process with . The double integral of a measurable function vanishing outside exists when the following condition holds:
[TABLE]
where for and for .
We remark that [25, Theorem 2] is only stated in the case where ; however, the extension to follows directly using a polarization techniques similar to the polarization equality for real numbers . Proposition 2.4 extends directly to complex-valued functions and processes as well.
Remark 2.2**.**
Let is measurable function with and set for all . Then we have the following statements.
(A) The function satisfies condition (2.4), and hence the double integral is well defined, cf. Proposition 2.1. Indeed, (2.4) follows by choosing .
(B) The following identity holds
[TABLE]
To show (2.5) we note that exists for all . For all in we have by integration by parts that
[TABLE]
where denotes the quadratic variation of two complex-valued semimartingales and . If we let and in (2.7) we obtain
[TABLE]
By reconizing the integral on the right-hand side of (2.8) as a double integral and calculating the quadratic variation we derive the identity
[TABLE]
where we have used that . Equation (2.9) implies (2.5).
We will also use the dominated convergence result [14, Theorem 6.2(c)], which reads as follows.
Proposition 2.3**.**
Assume that , and are deterministic measurable functions supported on all satisfying condition (2.4). Suppose moreover that is a two-dimensional symmetric -stable Lévy process with . If for Lebesgue almost all points and , then
[TABLE]
3 Proof of Theorem 1.1
Throughout this section all positive constants are denoted by , although they may change from line to line. We start with the identity
[TABLE]
In particular, we have that
[TABLE]
which implies the statement
[TABLE]
By using the representations (3.1) and (2.5) on the function given by we deduce the identity
[TABLE]
where
[TABLE]
3.1 Proof of convergence in (1.3)
By the decomposition (3.4) we have
[TABLE]
and thus we only need to prove that . We first note that
[TABLE]
for any , using the geometric series. On the other hand, we have
[TABLE]
In view of (3.3) and Remark 2.2(A), we can apply Proposition 2.3 to conclude that
[TABLE]
which implies (1.3).
3.2 Proof of weak convergence
We will divide the proof into several steps. We write and define the following complex-valued functions
[TABLE]
By the representation (3.4) we have
[TABLE]
where the last equality follows from the identity in (3.7) and the definition of . Thus, to show the limit theorem (1.4) it is enough to show the convergence
[TABLE]
The following subsections will be denoted to the proof of (3.13), and to this aim we define
[TABLE]
3.2.1 Existence of the limit
In this subsection we will show that the limiting integral from (3.13) is well-defined by proving that satisfies condition (2.4). Set
[TABLE]
for all , where is any real number satisfying and is chosen such that . To show condition (2.4) with respect to the function we will use the following lemma.
Lemma 3.1**.**
Let be a complex-valued function satisfying the inequality
[TABLE]
where the numbers satisfy the conditions and . Then it holds that
[TABLE]
Proof.
According to condition (3.15) we have that and hence we need to show the integrability of the functions and . We have that
[TABLE]
since . Hence, we conclude
[TABLE]
and the latter is finite as .
Now we turn our attention to the function . We use the substitution to obtain the bound
[TABLE]
Indeed, the latter integral is finite since . Thus, we deduce
[TABLE]
due to assumptions and . This completes the proof. ∎
Now, we proceed with the proof of existence. First of all, we observe that the function satisfies the condition (3.15) with
[TABLE]
It is easy to see that assumptions and are satisfied due to and . On the other hand, for any we readily deduce that
[TABLE]
When we choose small enough, we conclude that the function
[TABLE]
also satisfies the conditions of Lemma 3.1. This proves the existence of the double integral , cf. Proposition 2.1.
3.2.2 Main decomposition
Now, we show the weak limit theorem in (1.4). For this purpose, we introduce the sets
[TABLE]
for some fixed . For we also use the notation where is the closest number to . We obtain the following estimate for the function :
[TABLE]
for any . We decomposition the double integral on the right-hand side of (3.12) as follows
[TABLE]
with
[TABLE]
In the following, we will prove that is the dominating term and show its convergence, while we will prove that . The following identity in distribution
[TABLE]
where , will be used in our proof. The identity (3.25) follows from the self-similarity of ; more precisely, from the fact that equals in finite dimensional distributions.
3.2.3 Term
We apply the distributional identity (3.25) to deduce that
[TABLE]
For any fixed with we obtain the convergence
[TABLE]
Hence, we conclude that
[TABLE]
On the other hand, we have that
[TABLE]
As in the proof of existence we deduce that the function satisfies the condition (2.4). Thus,
[TABLE]
due to Proposition 2.3.
3.2.4 Term
We will use the approximation (3.20). First of all, we have that
[TABLE]
Note that since . On the other hand, the function satisfies the condition
[TABLE]
where we use an estimate of the type (3.19) for the log term. For all we have
[TABLE]
It is easy to show that the function satisfies the condition (2.4) when . Choosing we deduce that
[TABLE]
if we choose small enough since . We, therefore, conclude that there exists a measurable function satisfying (2.4) such that for all . Since Lebesgue almost surely, cf. (3.31)–(3.32), Proposition 2.3 yields .
3.2.5 Term
The proof is similar to the one in Subsection 3.2.4. For simplicity we only consider the integration region ; the case is analogue. For all :
[TABLE]
We have and the function satisfies the condition (2.4). On the other hand, we have
[TABLE]
The bounds (3.33) and (3.34) imply the existence of a measurable function satisfying (2.4) such that for all . Since a similar bound holds on we conclude, cf. Proposition 2.3, that since (cf. (3.33)–(3.34)).
By the decomposition (3.21) and the convergence results of Subsections 3.2.3–3.2.5 we obtain the limit theorem (3.13). This completes the proof of (1.4) and of Theorem 1.1. ∎
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1[1] A. Ayachea and J. Hamoniera (2012): Linear fractional stable motion: A wavelet estimator of the α 𝛼 \alpha parameter. Statistics and Probability Letters 82, 1569–1575.
- 2[2] A. Ayachea and Y. Xiao (2016): Harmonizable fractional stable fields: Local nondeterminism and joint continuity of the local times. Stochastic Processes and Their Applications 126(1), 171–185.
- 3[3] A. Basse-O’Connor, T. Grønbæk and M. Podolskij (2021): Local asymptotic self-similarity for heavy-tailed harmonizable fractional Lévy motions. ESAIM: Probability and Statistics 25, 286–297.
- 4[4] A. Basse-O’Connor, C. Heinrich and M. Podolskij (2019): On limit theory for functionals of stationary increments Lévy driven moving averages. Electronic Journal of Probability 24(79), 1–42.
- 5[5] A. Basse-O’Connor, R. Lachièze-Rey and M. Podolskij (2017): Power variation for a class of stationary increments Lévy driven moving averages. Annals of Probability , 45(6B), 4477–4528.
- 6[6] A. Basse-O’Connor and M. Podolskij (2017): On critical cases in limit theory for stationary increments Lévy driven moving averages. Festschrift for Bernt Øksendal, Stochastics 81(1), 360–383.
- 7[7] A. Basse-O’Connor, M. Podolskij and C. Thäle (2020): A Berry-Esseén theorem for partial sums of functionals of heavy-tailed moving averages. Electronic Journal of Probability 25(31), 1–31.
- 8[8] A. Benassi, S. Cohen and J. Istas (2002): Identification and properties of real harmonizable fractional Lévy motions. Bernoulli 8(1), 97–115.
