The Convergence Rate of Fréchet Distribution under the Second-Order Regular Variation Condition

Abstract

In this article we consider the asymptotic behavior of extreme distribution with the extreme value index γ>0 . The rates of uniform convergence for Fréchet distribution are constructed under the second-order regular variation condition.

Share and Cite:

Dai, X. (2024) The Convergence Rate of Fréchet Distribution under the Second-Order Regular Variation Condition. Journal of Applied Mathematics and Physics, 12, 1597-1605. doi: 10.4236/jamp.2024.125098.

1. Introduction

The central limit theory focuses on the extreme behavior of sample partial sums, but in nature and human society, there are also a class of extreme risk events, such as floods, earthquakes, precipitation and economic crises. Although these events are rare, once they occur, they will bring significant losses to society. Therefore, studying the laws of extreme value occurrence is extremely important. Extreme value theory emerged in this context, as an important branch of probability theory, mainly focusing on the tail behavior of extreme value distributions. In recent years, the application range of extreme value theory has been very extensive. For example, predicting the probability of extreme events such as the above, estimating the percentile of extreme value distribution, and applying it to fields such as financial risk management, see de Haan and Ferreira (2006) [1] .

Let { X n , n 1 } be independent, identically distributed (iid) random variables with common distribution function

F ( x ) = P [ X 1 x ] , x R .

Denote the extreme value by

M n = i = 1 n X i ,

and suppose there exist normalizing constants a n > 0 and b n R such that

M n b n a n

has a nondegenerate limit distribution as n , i.e.

P [ M n a n x + b n ] G ( x ) . (1.1)

Fisher and Tippett (1928) [2] , Gnedenko (1943) [3] proposed the extreme value distribution G ( x ) takes the form of

G ( x ) = G γ ( x ) = exp { ( 1 + γ x ) 1 / γ } , γ R , 1 + γ x 0 , (1.2)

where the parameter γ in (1.2) is called the extreme value index. This also means that F is in the domain of attraction of extreme value distribution.

Under the special case of the extreme value index γ > 0 , the extreme value distribution can be written as

G ( x ) = Φ 1 / γ ( x ) = { 0 , x < 0 , exp { x 1 / γ } , x 0 , (1.3)

which is also called the Fréchet distribution, and the convergence in (1.1) can be rewritten as

P [ M n a n x ] Φ 1 / γ ( x ) . (1.4)

Based on theoretical studies, many scholars focus on the first-order asymptotic analysis in extreme events. But with the widespread application of extreme value theory, several authors discovered the first-order asymptotic results obtained by using the limits of extreme value distributions are relatively rough, and often requires a more accurate approximate representation. It is necessary to know the further expansion of first-order convergence. Therefore, research on the convergence rate of first-order asymptotic result in extreme value theory has attracted the attention of many scholars. de Haan and Peng (1997) [4] considers the convergence rate of two-dimensional extreme value distribution. In the research on the convergence speed of one-dimensional extreme value distribution, de Haan and Resnick (1996) [5] established the rates of convergence of the distribution of the extreme order statistics to its limit distribution under the second-order von Mises condition with γ R . Cheng and Jiang (2001) [6] focuses on the rates of the uniform convergence for distributions of extreme values ( F n ( a n x + b n ) to G γ ( x ) ) under the condition of generalized regular variation of second-order. For the speed at which the extreme value distribution converges to its limit distribution in special distributions, Liao et al. (2014) [7] derived the asymptotic behavior of the distribution of the maxima for samples obeying skew-normal distribution. Peng et al. (2010) [8] established the limiting distribution of the extremes and the associated convergence rates for the mixed exponential distributions. Chen and Huang (2014) [9] construsted the exact uniform convergence rate of the asymmetric normal distribution of the maximum and minimum to its extreme value limit. Chen and Feng (2014) [10] considered the rates of convergence of extremes for short-tailed symmetric distribution under power normalization. Chen et al. (2012) [11] studied the rates of convergence of extremes for general error distribution under power normalization.

The second-order asymptotic result can provide a more accurate approximate expression, and it can characterize the speed of first-order convergence, which can provide a better guidance for the prediction, risk management, and control of extreme events, see Lin (2012) [12] , Mao and Hu (2013) [13] . The focus of this paper is on rates of convergence in (1.4). We set out to explain our condition. For a nondecreasing function f, define the left-continuous inverse of f is

f ( t ) = inf { x : f ( x ) t } .

Let f = ( 1 / log F ) . Necessary and sufficient condition for the convergence for (1.4) is that f is regulary varying, i.e.

l i m t f ( t x ) f ( t ) = x γ (1.5)

holds for x > 0 and γ > 0 , written as f R V ( γ ) . Regarding regular variation refer to Resnick (1987) [14] . So in order to get the convergence rate of (1.4), we need to require a rate of convergence condition in (1.5). The condition as follows.

Supppose the second-order regular variation condition

l i m t f ( t x ) / f ( t ) x γ B ( t ) = κ ( x ) (1.6)

holds for all x > 0 , where B has constant sign near infinity and satisfying lim t B ( t ) = 0 . The function κ ( x ) should not be a multiple of x γ . By Theorem B.3.1 of de Haan and Ferreira (2006) [1] , We know that | B | R V ρ and κ ( x ) = x γ ( x ρ 1 ) / ρ for ρ 0 .

For convenience, let G 0 ( x ) = exp { e x } and its derivative G 0 ( x ) = G 0 ( x ) e x , define

J n ( x ) = G 0 ( γ 1 log f ( n x ) / f ( n ) ) G 0 ( log x ) ; J ( x ) = γ 1 x γ G 0 ( log x ) κ ( x )

and a n : = f ( n ) . Moreover, for any function g on ( 0 , ) , denote g ( ) = lim t g ( t ) and g ( 0 ) = lim t 0 g ( t ) if the limits exist.

In the following, we will provide the rates of convergence in (1.4) under the second-order regular variation condition (1.6). The rest of the paper is organized as follows. In Section 2, we present the auxiliary lemmas. Theorem and its proof are given in Section 3.

2. Lemmas

Before presenting the main conclusion, we first provide the following lemmas. Recall that a measureable function f on R + is said to be generalized regular varying with prameter γ R and auxiliary a, denote f G R ( γ , a ) , if

l i m n f ( t x ) f ( t ) a ( t ) = x γ 1 γ , x R + . (2.1)

Define

a * ( t ) = { γ f ( t ) , γ > 0 ; γ ( f ( ) f ( t ) ) , γ < 0 ; f ( t ) 1 t 0 t f ( u ) d u , γ = 0

and

p n ( x ) = f ( n x ) f ( n ) x γ .

Lemma 2.1 (cf. de Haan (1970) [15] ). if f R V ( γ ) , for any ξ , δ 0 , there exists a t 0 = t 0 ( ξ , δ ) 0 such that

( 1 ξ ) min { x γ + δ , x γ δ } < f ( t x ) f ( t ) < ( 1 + ξ ) max { x γ + δ , x γ δ } , t , t x t 0 . (2.2)

Lemma 2.2 (Cheng and Jiang (2001) [6] , Proposition 1.2). If f G R ( γ , a ) , then a * ( t ) ~ a ( t ) as n and for all ξ , δ > 0 , there exists t 0 = t 0 ( ξ , δ ) 1 such that,

| f ( t x ) f ( t ) a * ( t ) x γ 1 γ | ξ max { x γ + δ , x γ δ } , t , t x t 0 (2.3)

holds for t , t x t 0 .

Lemma 2.3 (de Haan and Ferreira (2006) [1] , Theorem 2.3.9). If f satisfies the second-order condition (1.6), then for all ξ , δ 0 , there exists t 0 = t 0 ( ξ , δ ) 1 such that

| f ( t x ) f ( t ) x γ B ( t ) κ ( x ) | ξ x γ + ρ max { x δ , x δ } (2.4)

holds for t , t x t 0 , where

B ( t ) = { ρ ( 1 lim s s γ f ( s ) t γ f ( t ) ) , ρ < 0 ; 1 0 t s γ f ( s ) d s t 1 γ f ( t ) , ρ = 0.

Proof.

| f ( t x ) f ( t ) x γ B ( t ) κ ( x ) | = x γ | f ( t x ) f ( t ) x γ 1 B ( t ) x ρ 1 ρ | = x γ | ( t x ) γ f ( t x ) t γ f ( t ) t γ f ( t ) B ( t ) x ρ 1 ρ | .

Obviously, t γ f ( t ) G R ( ρ , t γ f ( t ) B ( t ) ) , then by Lemma 2.2 the lemma is complete.

Lemma 2.4. Suppose f satisfies the second-order regular variation condition (1.6), then

l i m n sup x [ α n , β n ] x ( 1 + γ ) G 0 ( log ( x ) ) ( B 1 ( n ) p n ( x ) κ ( x ) ) = 0 , (2.5)

l i m n sup x [ α n , β n ] B 1 ( n ) x 2 γ p n 2 ( x ) = 0, (2.6)

where α n = 1 / log B 2 ( n ) , β n = 1 / B 2 ( n ) .

Proof. Note that B ( t ) R V ρ and B 2 ( t ) R V 2 ρ with ρ 0 . From (2.2) in Lemma 2.1 there exists a constant C 0 > 0 and an integer n 0 > 0 such that B 2 ( n ) C 0 n 2 ρ 1 for all n n 0 . Hence n α n n [ ( 2 ρ 1 ) log n + log C 0 ] 1 as n . This implies that (2.4) holds for all x [ α n , β n ] . Therefore, for any δ ( 0,1 ) , we have

sup x [ α n , β n ] x ( 1 + γ ) G 0 ( log ( x ) ) ( B 1 ( n ) p n ( x ) κ ( x ) ) ξ sup x [ α n , β n ] x ( 1 + γ ) exp ( x 1 ) max { x γ + ρ + δ , x γ + ρ δ } = ξ sup x R + exp ( x 1 ) max { x ρ 1 + δ , x ρ 1 δ } < .

Hence (2.5) holds. For (2.6), choosing δ ( 0, 1 / 4 ) ,

B 1 ( n ) x 2 γ p n 2 ( x ) = x 2 γ B ( n ) B 2 ( n ) p n 2 ( x ) = x 2 γ B ( n ) ( B 1 ( n ) p n ( x ) κ ( x ) + κ ( x ) ) 2 2 x 2 γ B ( n ) ( B 1 ( n ) p n ( x ) κ ( x ) ) 2 + 2 x 2 γ B ( n ) κ 2 ( x ) 2 ξ B ( n ) sup x [ α n , β n ] max { x 2 ρ + 2 δ , x 2 ρ 2 δ } + 2 x 2 γ B ( n ) κ 2 ( x ) = Q n 1 + Q n 2 .

The first part Q n 1 2 ξ B ( n ) max { ( 1 / log B 2 ( n ) ) 2 ρ 2 δ , ( B ( n ) ) 4 ρ 4 δ } 0 as n . Similarly, the second part Q n 2 = sup x [ α n , β n ] 2 B ( n ) ( ( x ρ 1 ) / ρ ) 2 0 as n . The lemma is proved.

Lemma 2.5. If f satisfies the second-order regular variation condition (1.6), then

l i m n sup x [ 0, ] | B 1 ( n ) J n ( x ) J ( x ) | = 0.

Proof. We only prove this lemma for the case that B ( n ) is positive, the proof of another case is similar.

J n ( x ) = G 0 ( 1 γ log f ( n x ) f ( n ) ) G 0 ( log x ) = q n ( x ) G 0 ( log x + θ ^ 1 n q n ( x ) ) ,

where q n ( x ) = γ 1 log f ( n x ) / f ( n ) log x and θ ^ 1 n ( 0,1 ) . Note that for some θ ^ 2 n ( 0,1 ) ,

q n ( x ) = γ 1 x γ p n ( x ) ( 1 + θ ^ 2 n x γ p n ( x ) ) 1 .

Therefore,

| B 1 ( n ) J n ( x ) J ( x ) | = | B 1 ( n ) q n ( x ) G 0 ( log x + θ ^ 1 n q n ( x ) ) J ( x ) | B 1 ( n ) | q n ( x ) γ 1 x γ p n ( x ) | G 0 ( log x + θ ^ 1 n q n ( x ) ) + B 1 ( n ) γ 1 x γ p n ( x ) | G 0 ( log x + θ ^ 1 n q n ( x ) ) G 0 ( log x ) | + | B 1 ( n ) γ 1 x γ p n ( x ) G 0 ( log x ) J ( x ) | = E n 1 + E n 2 + E n 3 .

By Lemma 2.4, we know that

sup x [ α n , β n ] E n 3 ( x ) = sup x [ α n , β n ] γ 1 x γ G 0 ( log x ) ( B 1 ( n ) p n ( x ) κ ) ( x ) = sup x [ α n , β n ] γ 1 x ( γ + 1 ) G 0 ( log x ) ( B 1 ( n ) p n ( x ) κ ) ( x ) 0.

Letting M = max { sup x R G 0 ( log x ) , sup x R G 0 ( log x ) } , we have from (2.6) that

sup x [ α n , β n ] E n 1 ( x ) sup x [ α n , β n ] | M γ 1 B 1 ( n ) x γ p n ( x ) [ ( 1 + θ ^ 2 n x γ p n ( x ) ) 1 1 ] | sup x [ α n , β n ] | M γ 1 B 1 ( n ) x 2 γ p n 2 ( x ) ( 1 + θ ^ 2 n x γ p n ( x ) ) 1 | 0.

For the second part E n 2 ( x ) , by mean value theorem and also from (2.6),

sup x [ α n , β n ] E n 2 ( x ) = sup x [ α n , β n ] | γ 1 B 1 ( n ) x γ p n ( x ) G 0 ( log x + θ ^ 3 n θ ^ 1 n q n ( x ) ) θ ^ 1 n q n ( x ) | sup x [ α n , β n ] M γ 2 B 1 ( n ) x 2 γ p n 2 ( x ) ( 1 + θ ^ 2 n x γ p n ( x ) ) 1 0.

So we obtain

l i m n sup x [ α n , β n ] | B 1 ( n ) J n ( x ) J ( x ) | = 0. (2.7)

It remains to deal with the parts of the integral near ± . For x β n ,

sup x [ β n , ] | B 1 ( n ) J n ( x ) | B 1 ( n ) sup x [ β n , ] ( | 1 G 0 ( γ 1 log f ( n x ) / f ( n ) ) | + | 1 G 0 ( log x ) | ) B 1 ( n ) | 1 G 0 ( γ 1 log f ( n β n ) / f ( n ) ) | + | 1 G 0 ( log β n ) | B 1 ( n ) | G 0 ( γ 1 log f ( n β n ) / f ( n ) ) G 0 ( log β n ) | + 2 B 1 ( n ) | 1 G 0 ( log β n ) | .

Noting J ( ) = 0 , the first part goes to zero by (2.7). The second part goes to zero because of 1 G 0 ( log β n ) ~ B 2 ( n ) and B ( n ) 0 as n . So we have lim n sup x [ β n , ] | B 1 ( n ) J n ( x ) J ( x ) | = 0 . Similarly, for x α n ,

sup x [ 0, α n ] | B 1 ( n ) J n ( x ) | B 1 ( n ) sup x [ 0, α n ] ( | G 0 ( γ 1 log f ( n x ) / f ( n ) ) | + | G 0 ( log x ) | ) B 1 ( n ) | G 0 ( γ 1 log f ( n α n ) / f ( n ) ) G 0 ( log α n ) | + 2 B 1 ( n ) | G 0 ( log α n ) | 0.

Combing J ( 0 ) = 0 we have lim n sup x [ 0, α n ] | B 1 ( n ) J n ( x ) J ( x ) | = 0 . The proof of the lemma is completed.

3. Theorem and Its Proof

Theorem 3.1. If f satisfies the second-order regular variation condition (1.6), then

l i m n F n ( a n x ) Φ 1 / γ ( x ) B ( n ) = J ( x 1 / γ ) (3.1)

holds uniformly on x > 0 .

Proof. Let z n ( x ) = [ n log F ( a n x ) ] 1 , equivalently, F n ( a n x ) = G 0 ( log z n ( x ) ) . So we have

K n ( x ) : = B 1 ( n ) ( F n ( a n x ) Φ 1 / γ ( x ) ) + J ( x 1 / γ ) = B 1 ( n ) ( G 0 ( log z n ( x ) G 0 ( 1 / γ log x ) ) + J ( z n ( x ) ) + J ( x 1 / γ ) J ( z n ( x ) ) ) = K n 1 ( x ) + K n 2 ( x ) .

In order to show (3.1), we need only to prove

lim n sup 0 < F ( a n x ) < 1 | K n i ( x ) | = 0 for i = 1 , 2 ; (3.2)

l i m n sup F ( a n x ) = 0 | K n ( x ) | = 0 ; (3.3)

lim n sup F ( a n x ) = 1 | K n ( x ) | = 0. (3.4)

If 0 < F ( a n x ) < 1 , according to the definition of the f, for any δ > 0 we have

f ( n z n ( x ) ) a n x f ( n z n ( x ) + δ ) a n .

Hence

B 1 ( n ) ( G 0 ( log z n ( x ) ) G 0 ( 1 / γ log ( f ( n z n ( x ) + δ ) / a n ) ) + J ( z n ( x ) ) ) K n 1 ( x ) B 1 ( n ) ( G 0 ( log z n ( x ) ) G 0 ( 1 / γ log ( f ( n z n ( x ) ) / a n ) ) + J ( z n ( x ) ) ) .

Then by Lemma 2.5, we can obtain lim n sup 0 < F ( a n x ) < 1 | K n 1 ( x ) | = 0 . It is obvious that z n ( x ) x 1 / γ . Since J ( x ) is continuous on ( 0, ) and J ( 0 ) = J ( ) = 0 , we can also obtain lim n sup 0 < F ( a n x ) < 1 | K n 2 ( x ) | = 0 .

For the situation of F ( a n x ) = 0 , we have x f ( 0 ) / a n . The left of (3.1) is

l i m n sup F ( a n x ) = 0 B 1 ( n ) Φ 1 / γ ( x ) l i m n sup F ( a n x ) = 0 B 1 ( n ) Φ 1 / γ ( a n 1 f ( 0 ) ) = l i m n sup F ( a n x ) = 0 ( B 1 ( n ) J n ( 0 ) J ( 0 ) ) 0.

Note that lim n f ( 0 ) / a n 0 . For any δ , there exists n 0 such that x f ( 0 ) / a n δ for all n n 0 . Therefore,

l i m n sup F ( a n x ) = 0 | J ( x 1 / γ ) | l i m n sup x δ | J ( x 1 / γ ) | = 0

by letting δ 0 . Then,

l i m n sup F ( a n x ) = 0 | B 1 ( n ) ( F n ( a n x ) Φ 1 / γ ( x ) ) + J ( x 1 / γ ) | l i m n sup F ( a n x ) = 0 B 1 ( n ) Φ 1 / γ ( x ) + l i m n sup F ( a n x ) = 0 J ( x 1 / γ ) 0.

If F ( a n x ) = 1 , then x f ( ) / a n . The left of (3.1) is

l i m n sup F ( a n x ) = 1 | B 1 ( n ) ( 1 Φ 1 / γ ( x ) ) | l i m n sup F ( a n x ) = 1 | B 1 ( n ) ( 1 Φ 1 / γ ( f ( ) / a n ) ) | l i m n sup F ( a n x ) = 1 | B 1 ( n ) ( G 0 ( log ) Φ 1 / γ ( f ( ) / a n ) ) | = l i m n sup F ( a n x ) = 1 | B 1 ( n ) J n ( ) J ( ) | 0.

For sufficiently large number M 0 , there exists n 0 such that x f ( ) / a n M 0 for all n n 0 . Hence

l i m n sup F ( a n x ) = 1 | J ( x 1 / γ ) | l i m n sup x M 0 | J ( x 1 / γ ) | = 0

by letting M 0 . Furthermore,

l i m n sup F ( a n x ) = 1 | B 1 ( n ) ( F n ( a n x ) Φ 1 / γ ( x ) ) + J ( x 1 / γ ) | 0.

So we prove this theorem.

Remark 3.1. Uniform limit in Theorem 3.1 gives an Edgeworth expansion as follows:

F n ( a n x ) = Φ 1 / γ ( x ) + 1 γ B ( n ) ( log Φ 1 / γ ( x ) ) 1 + γ Φ 1 / γ ( x ) κ ( log 1 Φ 1 / γ ( x ) ) + o ( B ( n ) )

holds uniformly on x > 0 .

Conflicts of Interest

The author declares no conflicts of interest regarding the publication of this paper.

References

[1] Haan, L. and Ferreira, A. (2006) Extreme Value Theory. Springer, New York.
https://doi.org/10.1007/0-387-34471-3
[2] Fisher, R.A. and Tippett, L.H.C. (1928) Limiting Forms of the Frequency Distribution of the Largest or Smallest Member of a Sample. Mathematical Proceedings of the Cambridge Philosophical Society, 24, 180-190.
https://doi.org/10.1017/S0305004100015681
[3] Gnedenko, B.V. (1943) Sur La Distribution Limite Du Terme Maximum D’Une Serie Aleatoire. Annals of Mathematics, 44, 423-453.
https://doi.org/10.2307/1968974
[4] de Haan, L. and Peng, L. (1997) Rates of Convergence for Bivariate Extremes. Journal of Multivariate Analysis, 61, 195-230.
https://doi.org/10.1006/jmva.1997.1669
[5] de Haan, L. and Resnick, S. (1996) Second Order Regular Variation and Rates of Convergence in Extreme Value Theory. Annals of Probability, 24, 97-124.
https://doi.org/10.1214/aop/1042644709
[6] Cheng, S. and Jiang, C. (2001) The Edgeworth Expansion for Distributions of Extreme Values. Science in China Series A: Mathematics, 44, 427-437.
https://doi.org/10.1007/BF02881879
[7] Liao, X., Peng, Z., Nadarajah, S. and Wang, X. (2014) Rates of Convergence of Extremes from Skew-Normal Samples. Statistics and Probability Letters, 84, 40-47.
https://doi.org/10.1016/j.spl.2013.09.027
[8] Peng, Z., Weng, Z. and Nadarajah, S. (2010) Rates of Convergence of Extremes for Mixed Exponential Distributions. Mathematics and Computers in Simulation, 81, 92-99.
https://doi.org/10.1016/j.matcom.2010.07.003
[9] Chen, S. and Huang, J. (2014) Rates of Convergence of Extreme for Asymmetric Normal Distribution. Statistics and Probability Letters, 84, 158-168.
https://doi.org/10.1016/j.spl.2013.10.003
[10] Chen, S. and Feng, B. (2014) Rates of Convergence of Extreme for STSD under Power Normalization. Journal of the Korean Statistical Society, 43, 609-619.
https://doi.org/10.1016/j.jkss.2014.02.001
[11] Chen, S., Wang, C. and Zhang, G. (2014) Rates of Convergence of Extreme for General Error Distribution under Power Normalization. Statistics and Probability Letters, 82, 385-395.
https://doi.org/10.1016/j.spl.2011.10.019
[12] Lin, J. (2012) Second Order Asymptotics for Ruin Probabilities in a Renewal Risk Model with Heavy-Tailed Claims. Insurance: Mathematics and Economics, 51, 422-429.
https://doi.org/10.1016/j.insmatheco.2012.07.001
[13] Mao, T. and Hu, T. (2013) Second-Order Properties of Risk Concentrations without the Condition of Asymptotic Smoothness. Extremes, 16, 383-405.
https://doi.org/10.1007/s10687-012-0164-z
[14] Resnick, S.I. (1987) Extreme Values, Regular Variation, and Point Processes. Springer, New York.
https://doi.org/10.1007/978-0-387-75953-1
[15] de Haan, L. (1970) On Regular Variation and Its Application to the Weak Convergence of Sample Extremes. Mathematisch Centrum, Amsterdam.

Copyright © 2024 by authors and Scientific Research Publishing Inc.

Creative Commons License

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.