Extremes of -norm of Vector-valued Gaussian processes with Trend
Abstract: Let be a Gaussian vector process and be a continuous function. The asymptotics of distribution of , the norm for Gaussian finite-dimensional vector, have been investigated in numerous literatures. In this contribution we are concerned with the exact tail asymptotics of with trend over . Both scenarios that is locally stationary and non-stationary are considered. Important examples include and chi-square processes with trend, i.e., . These results are of interest in applications in engineering, insurance and statistics, etc.
Keywords: Tail asymptotics; -norm; vector-valued Gaussian process; fractional Brownian motion; Pickands constant; Piterbarg constant.
AMS Classification: Primary 60G15; secondary 60G70
In engineering sciences, extreme values of non-linear functions of multivariate Gaussian processes are of interest in dealing with the safety of structures, see [EXCRA1980] and the references therein. Probabilistic structural analysis to answer the question is: what is the probability that a certain mechanical (or other) structure will survive when it is subject to a random load. The load is then usually defined by some -dimensional vector process , and one seeks the probability that exceeds some more or less well-defined safe region, which is specific for the structure as
where the time-dependent safety region is defined by
with some continuous function and , the norm, i.e.,
in the space .
Assume that where are independent copies of a centered Gaussian process which has continuous trajectories, variance function and correlation function and
and hereafter, we call the norm process.
When , for a positive constant , as in the convention is called the chi process when and the chi-square process when .
Further, as the Gaussian processes, we can introduce the stationary, locally-stationary, and non-stationary norm processes according to the stationary, locally-stationary, and non-stationary properties of , respectively.
The investigate of
is initiated by the studies of high excursions of envelope of a Gaussian process, see e.g., [BN1969] and generalized in [Lindgren1980a, Lindgren1980b, Lindgren1989]. When is stationary with and
[Albin1990, Albin1992] develop the Berman’s approach in [Berman82] to obtain an asymptotic behavior of large deviation probabilities of the stationary chi-square processes.
Further, if there exists unique satisfies and
where and are positive constants related to , the tail asymptotic behavior of the non-stationary and are investigated in [Pitchi1994] and [FatalovA1993], respectively, under the application of the so-called ”double-sum method” in [Pit96].
Some recent contributions are focused on more general scenarios of chi process and chi-square process with , i.e.,
where the continuous function is generally considered as a trend or a drift.
When are non-stationary Gaussian processes, , the non-stationary chi processes with trend, and , the non-stationary chi-square processes with trend, are studied in [EnkelejdJi2014Chi] and [PL2015], respectively.
When are locally-stationary Gaussian processes, [LJ2017] obtains the extreme of the supremum of with trend, see, e.g., [Berman92, Hus90] for more details about locally stationary Gaussian processes.
Considering both the locally stationary and non-stationary norm processes, the contribution of this paper concerns an exact asymptotic behavior of large deviation probabilities for with , constant and a continuous function, which contains the aforementioned results.
Organisation of the rest of the paper: In Section 2, the notation and some preliminaries are given. Our main results are displayed in Section 3. Following in Section 4 are two applications related to insurance and statistics. Finally, we present the proofs in Section 5 and several lemmas in Section 6.
2. Notation and preliminaries
First we introduce some notation, starting with the well-known Pickands constant defined by
where are constants and is a standard fractional Brownian motion (fBm) with Hurst index Further, define for non-negative continuous function.
The exact values of are known for and , namely,
See [PicandsA, Pit72, debicki2002ruin, DI2005, DE2014, DiekerY, DEJ14, Pit20, Tabis, DM, SBK, GeneralPit16] for various properties of and .
Through this paper means asymptotic equivalence when the argument tends to or .
We notice that denotes the tail distribution function of an random variable and .
For the norm process in (3) and a continuous function , we shall investigate the asymptotics of
with a constant. As in [FatalovA1993, Pitchi1994], for , using the duality property of norm we find
where is a centered Gaussian field defined on cylinder with
where if , if and if .
On , attains its maximum at:
(i) for at points where ( stands at the i-th position), ( stands at the i-th position), ;
(ii) for at points on , ;
(iii) for at points , where
( we take all possible combinations of signs ”+” and ”-” ), where .
The proof can be easily carried out by method of Lagrangian multipliers or referring to [FatalovA1993] [Lemma 3.1].
Next by [randomChaos2015], we have the following lemma.
3. Extremes of norm processes with trend
In this section, recall that in (3) is the norm process and ’s are independent copies of with continuous trajectories, variance functions and correlation functions .
3.1. Extremes of non-stationary norm processes with trend
As in [Threshold2016], if is non-stationary, we introduce the following assumptions:
attains its maximum on at the unique point and
for some positive constants .
for some constants and
Further, we introduce a bounded measurable trend function which satisfies
for some constants .
3.2. Extremes of locally stationary norm processes with trend
If is locally stationary, as in [Threshold2016], we shall suppose that:
where are positive continuous function on .
Before giving the scenarios with trend, we consider the extremes of the centered locally stationary norm processes.
Assume that , i.e., unit variance and covariance function satisfies assumptions (iv) and (v). Then we have for
where is the same as in Lemma 2.1.
Assume that , i.e., unit variance and correlation function satisfies assumptions (iv) and (v). Assume that is a continuous function which attains its maximum at a unique point satisfying assumption (iii) for some constants . Further, set and and is the same as in Lemma 2.1.
If , then we have as
where and if , if .
If , then we have
If , then we have
By the proof, we notice that for the case in Theorem 3.4, the result always holds for any continuous function . When , the result holds for any bounded function .
Following example is a special case of Theorem 3.4, which is corresponded with [LJ2017] [Theorem 2.1].
In Theorem 3.4, assume that , and is a continuous function, then we have
4.1. Ruin probability of a risk model
In theoretical insurance modelling a surplus process can be defined by
see [MR1458613], where is the initial reserve, is the rate of premium and the stochastic process denotes the aggregate claims process. See [rolski2009stochastic, DHJ13a, HX2007, ParisianBrownianfinite2017, Threshold2016, ParisianInfinite2018] for more studies on related risk models. Here we investigate
where is the same as in (2) and are independent fractional Brownian motions. can be considered as the sum of independent claims or payments until time . The corresponding ruin probability over a finite-time horizon is defined as
We present next approximation of this ruin probability.
We have as
Besides in risk modelling, the norm processes, especially the chi-square processes, are also widely utilized in hypothesis testing, see [HTDa1987, SCSP2014] and the reference there. Next we give an example.
4.2. The Ornstein-Uhlenbeck chi-square process in Quantitative Trait Locus detection
A Quantitative Trait Locus (QTL) denotes a gene with quantitative effect on a trait. The method used by most of geneticists in order to detect a QTL on a chromosome, is the Interval Mapping proposed by [LanderBO1989]. Using the Haldane distance and modelling in [Haldane1919], each chromosome is represented by a segment . The distance on is called the genetic distance. At each location , using the ”genome information” brought by genetic markers, a likelihood ratio test (LRT) is performed, testing the presence of a QTL at this position. [ADR2014] prove that when the number of genetic markers and the number of progenies tends to infinity, the limiting process of the LRT process is an Ornstein-Uhlenbeck chi-square process under the null hypothesis of the absence of QTL on the interval . In order to take decision about the presence of a QTL on , we need to calculate the critical value for the supremum of an Ornstein-Uhlenbeck chi-square process, i.e.,
where the Ornstein-Uhlenbeck chi-square process is
and are independent identically stationary Gaussian processes with covariance function given by
We have as
During the following proofs, are some positive constants which can be different from line by line and for interval we denote
Proof of Theorem 3.1:
We first present the proof for the case .
Set , with and for large enough
with the same as in (5) which is a centered Gaussian field.
We have for some small and large enough
We first give the upper bounds of and .
Set and . Then by Borell inequality as in [AdlerTaylor] and Lemma 2.2 for large
By assumptions (i) and (iii), we know that for some
hold for when small enough, then
Denote with . By assumption (ii), we have that
which combined with (6) imply
Now we focus on the asymptotic of , as .
Denote for any and some