arxiv: v2 [stat.me] 1 Jan 2018

Transcription

1 UNIFORM RATES OF CONVERGENCE OF SOME REPRESENTATIONS OF EXTREMES : A FIRST APPROACH arxiv: v2 [stat.me] Jan 208 M. NGOM ), T. A. KPANZOU 5), C. M. M. TRAORÉ,4), M.DIALLO,4), AND G. S. LO,2,3) Abstract. Uniform convergence rates are provided for asymptotic representations of sample extremes. These bounds which are universal in the sense that they do not depend on the extreme value index are meant to be extended to arbitrary samples extremes in coming papers. Résumé. Des vitesses de convergence uniformes, universelles en ce sens qu elles ne dpendent pas de l index extremal, sont fournies pour des representations asymptotiques des extrmes d chantillon. Ces bornes sont tendre dans le cas gnral des extrmes dans des articles venir. LERSTAD, Gaston Berger University, Saint-Louis, Sngal Gane Samb Lo, Modou Ngom). 2 Affiliated to LSTA, Universit Pierre et Marie Curie, Paris, France Gane Samb LO) 3 Associated to African University of Sciences and Technology, AUST, Abuja, Nigeria Gane Samb Lo) 4 Universit des Sciences, Techniques et des Technologies de Bamako, Mali Cherif Traoré, Moumouni Diallo) 5 Universit de Kara, Togo Tchilabalo Atozou Kpanzou).. Corresponding author. Gane Samb Lo. gane-samb.lo@edu.ugb.sn, ganesamblo@ganesamblo.net Permanent address : 78 Evanston Dr NW T3P 0J9, Calgary, Alberta, Canada. keywords and phrases : Extreme value Theory; weak convergence; rate of convergence, Scheff s Theorem. AMS 200 Classification Object : Primary 60B0; 60G70; 60G20.. Introduction This paper presents a view of the rate of convergence of the univariate extremes of samples in a simple form. Rather than trying to handle the general problem for a distribution function in the extreme domain of attraction, we focus here on the simplest representations for which we give the most precise rates. The results are expected to serve as tools in general for the univariate case and later for the multivariate frame. In a few words, the extreme value theory started with the univariate case, especially with independent data. Given a sequence of independent and identically distributed randomvariablesx n ) n 0 with commoncumulative distribution function cdf) and

2 2M. NGOM ), T. A. KPANZOU 5), C. M. M. TRAORÉ,4), M.DIALLO,4), AND G. S. LO,2,3) defined on the same probability space Ω, A, P), the max-stability problem consists of finding possible limits in distribution of the sequences of partial maxima M n = maxx,...,x n ), n, when appropriately centered and normalized. Precisely, we want to find nonrandom sequences a n > 0) n R and b n ) n R such that M n b n )/a n converges in distribution to a random variable Z as n +, denoted as M n b n Z. M) a n This problem was solved around the middle of 20-th century with the contributions of many people, from whom we can cite Gnedenko 943), Fisher and Tippet 928), Frchet 927), etc. The following result is usually quoted as the Gnedenko 943) result since he had the chance to close the characterization theorem. If Z is non-degenerate, that is : Z takes at least two different values, Formula M) can hold only if Z is one of the three types in terms of its cdf: the Frchet type with parameter γ Γ = {x R, x > 0} : H γ x) φ γ x) = exp x /γ ) x 0), the Weibull type with parameter γ Γ 2 = {x R, x < 0} : H γ x) ψ γ x) = exp x) /γ ) x<0) + x 0) and the Gumbel type with parameter γ Γ 0 = {0} : H γ x) Λx) = exp e x ), x R. By type of distribution, we mean that any non-degenerate Z in Formula M) have a cdf F Z which is of the form F Z x) = H γ Ax+B), where 0 < A R, B R, and γ R. A modern account of the theory, including statistical estimation, can be found in Galambos 958), de Haan 970), de Haan and Ferreira 2006), Embrechts et al. 997), Resnick 987), Beirlant et al. 204), etc. This theory is part of the general weak convergence which is thoroughly treated in Billingsley 968) and van der Vaart and Wellner 996). But all our needs in extreme value theory and in weak convergence are already gathered in Lo et al. 206) and Lo 206), which are currently cited in the paper. We denote by Z = Frγ), Z 2 = Weiγ) and Z 0 = Λ random variables admitting the above cdf s respectively. We remark that the cdf s H γ are continuous, and the cdf of M n b n )/a n is Fa n x+b n ) n, x R. By Theorem 3 in Chapter 2 and Point 5) in Chapter 4, Section in Lo et al. 206), M) is equivalent to R n = Fa n x+b n ) n H γ x) 0 as n +.

3 UNIFORM RATES OF CONVERGENCE OF SIMPLE EXTREMES 3 The general problem of the rate of convergence here is to find R n or to have universal bounds for it. As we already stressed, here we propose results to be used in the general case. To introduce this, pose that we have a sequence of independent and uniformly distributed random variables U n ) n 0 on 0,) and let us consider for each the U,n = minu,...,u n ) for each n and define and Z n,γ) = nu,n ) γ, γ Γ ; Z n 2,γ) = nu,n ) γ, γ Γ 2 Z n 0,γ) = lognu,n ), γ Γ 0. It is straightforward to show that for γ Γ i, Z n i,γ) weakly converges to Frγ), to Weiγ) and to Λ accordingly to i = 0, i = and i = 2. But the most important is that, based on Section 3.2 in Chapter in Lo 206) for example, if M) holds with Z non-degenerate, we may always get the following asymptotic representation in distribution M n b n a n = Z n i,γ)+r n, with r n 0 as n +, MG) for some i {0,,2} with suitable sequences a n > 0) n R and b n ) n R. Because of the general formula MG), we think it is important to have a special study of the rate of convergence of the Z n i,γ) s, which paves the way, further, for the handling of the general rate r n. This is what we are going to do exactly in Section Rate of convergence of Z n i,γ) We are going to see that a universal bound will be guided by the following rate of convergence e n) n ) which is estimated as follows. Lemma. For all n 2 0 e n) n ) ) C 0 2n + n 2. logn where C 0 = exp )=

4 4M. NGOM ), T. A. KPANZOU 5), C. M. M. TRAORÉ,4), M.DIALLO,4), AND G. S. LO,2,3) Proof. It is given in the Appendix. The main result of the paper is the following. Theorem. We have for all n 2 i {0,,2} γ Γ i F Zni,γ)x) F Zi x) 2+C 0 4n + C 0 2n 2 logn. Proof. We first deal with the Z n,γ),γ > 0. Set α = /γ. We have and x R), F n x) = F Zn,γ)x) = { ) n nx α if nx α > 0 otherwise x R), f n x) = F Z n,γ) x) = { αx α nx α ) n if nx α > 0 otherwise. Since, as n +, Z n,γ) Fr/γ), we have, by Scheffé s Theorem see for example Theorem 4, Chapter 3 in Lo et al. 206)) : FZn,γ)x) φ γ x) f n x) fx) dx B measurable B = f n x) fx) dx, 2 with { x R), fx) = F αx Fr/γ) x) = α e x α if x > 0 0 otherwise. For short, we also write F Fr/γ) = F. Now, we have to compute 2.) 2a n α) = f n x) fx) dx. We have f n x) fx) dx = 2.2) [0,n /α ] = Fn /α )+ f n x) fx) dx+ ]n /α,+ [ ]n /α,+ [ f n x) fx) dx. f n x) fx) dx

5 UNIFORM RATES OF CONVERGENCE OF SIMPLE EXTREMES 5 We also have f n x) fx) dx = ]n /α,+ [ ]n /α,+ [ αx α e x α e x α ) n nx α dx. Now, let us study g n,α x) = e x α ) n nx α on ]n /α,+ [. Differentiating g n,α x) gives g n,αx) = α n ex α x α ) n 2 nx α x α ) on ]n /α,+ [. For n 2, g n,αx) is positive on ]n /α, ] and negative on [,+ [ so that g n,α x) is increasing ]n /α,] and decreasing on [,+ [. Since for n 2, and, as n +, g n,α ) = e /n) n > g n,α ) + i.e., by from above ), we have a unique number x n,α ]n /α,] such that g n,α x n,α ) = and and for all x n /α, 0 g n,α x) on ]n /α,x n,α ] g n,α x) g n,α ) = e /n) n. By exploiting the sign of g n,α x) on ]n /α,] and ], + ] respectively, we have ]n /α,+ [ xn,α αx α e x n /α = 2.3) a n α,)+a n α,2). Next, we have 2.4) αx α e xα e x α nx α a n α,) = α Combining all, this leads to g n,α x))dx+ xn,α n /α fx)dx xn,α ) n ]x n,α,+ [ n /α = Fx n,α ) Fn /α ) dx αx α e x α g n,α x) )dx ) n ) nx α dx ) n nx α. n,α

6 6M. NGOM ), T. A. KPANZOU 5), C. M. M. TRAORÉ,4), M.DIALLO,4), AND G. S. LO,2,3) 2.5) a n α,2) g n,α ) ) Fx n,α )) ) e /n) n ) α n 3) In total, by putting Formulas , we get 2α n α) exp x α n,α ) ) n +α n 3) exp x α n,α) expx α n,α) exp x α n,α) exp x α n g n,α x n,α ) n,α) x α n,αexp x α n,α) +α n 3), ) +α n 3), ) n ) )) +α n 3) L2) )) +α n 3) wherewe identified g n,α x n,α ), which is equal to one, in Line 2 in the abovegroup of formulas. Now, the term between the big parentheses is bounded by the remum ofthe function lx) = x α e x α on 0,)whosederivative, αx α e x α x α ), is non-negative on 0,). Hence our bound is l) = /e. We conclude that for n 2 2α n α) ) e /n) n + n. By using Lemma 2, we finally get α n α) ) 2n +C 0 4n + 2 n 2. logn From there, we use Lemma 2 to conclude the proof for Z n,γ). To finish for the other cases, it is enough to see that we have the following two formulas: and F Zn2,γ)x) F Z2 x) = F Zn, /γ) x ) F Z x ) which puts an end to the proof. F Zn0,γ)x) F Z0 x) = F Zn,)e x ) F Z e x ),

7 UNIFORM RATES OF CONVERGENCE OF SIMPLE EXTREMES 7 Remark. The bounds provided here are to be used in the general theory of a coming paper. Proof. Fix n 2. We have 3. Appendix n )log /n) = + k k ) k + n k, which in turn implies g n, α ) = e n) n ) = exp kk +)n k. F) By the comparison methods between series with monotonic terms, say a n = an) with the appropriate improper Riemann integral of ax), we have k k kk +)n k + 2n + dx 2 n x 2n + n 2 logn. F2) By combining F) and F2), we have, for 0 θ, 0 e n) n ) ) exp 2n + n 2 logn ) )) 2n + n 2 exp θ logn 2n + n 2 logn ) exp ) 2n + n 2. logn References Beirlant, J., Goegebeur, Y. Teugels, J. 2004). Statistics of Extremes Theory and Applications. Wiley. MR20803) Billingsley, P. 968). Convergence of Probability measures. John Wiley, New-York. de Haan, L. 970). On regular variation and its application to the weak convergence of sample extremes. Mathematical Centre Tracts, 32, Amsterdam. MR028656) de Haan, L. and Ferreira A. 2006). Extreme value theory : An introduction. Springer. MR223456)

8 8M. NGOM ), T. A. KPANZOU 5), C. M. M. TRAORÉ,4), M.DIALLO,4), AND G. S. LO,2,3) Embrechts, P., Kűppelberg C. and Mikosh T. 997). Modelling extremal events for insurance and Finance. Springer Verlag. Fisher, R. and Tippet, L. 928) Limiting Forms of the Frequency Distribution of the Largest or Smallest Member of a Sample. Proceedings of the Cambridge Philosophical Society, 24, Frchet, M. 927). Sur la loi de probabilit de l cart maximum. Annales de la Socit Polonaise de Mathmatique, 927, 6, Galambos, J. 985). The Asymptotic theory of Extreme Order Statistics. Wiley, Nex-York. MR ) Gnedenko, B. 943) Sur la distribution limite du terme maximum d une srie alatoire. Annals of Mathematics, 44, Lo, G.S., Ngom, M. and Kpanzou T.A. 206). Weak Convergence IA). Sequences of random vectors. SPAS Books Series. Saint-Louis, Senegal - Calgary, Canada. Doi : /sbs/ Arxiv : ISBN : Lo, G.S. 206). Univariate Theory of Extreme Value Theory and Statistical Estimation in the Extreme Value Domain. SPAS Books Series. Saint-Louis, Senegal - Calgary, Canada. Doi : /sbs/ Arxiv : ISBN : ISBN Resnick, S.I. 987). Extreme Values, Regular Variation and Point Processes. Springer-Verlag, New-York. MR090080) van der Vaart, A. W. and Wellner, J. A. 996). Weak Convergence and Empirical Processes With Applications to Statistics. Springer. address: ngomodou@gmail.com, kpanzout@gmail.com, cheriftraore75@gmail.com, moudiallo@gmail.com, gane-samb.lo@ugb.edu.sn