arxiv: v2 [stat.me] 29 Jan 2018

Transcription

1 Inference on he ail process wih applicaion o financial ime series modelling arxiv: v2 [sa.me] 29 Jan 2018 Richard A. Davis Columbia Universiy Deparmen of Saisics 1255 Amserdam Ave. New York, NY 10027, USA rdavis@sa.columbia.edu Johan Segers Holger Drees Universiy of Hamburg Deparmen of Mahemaics Bundessraße 55, Hamburg, Germany holger.drees@mah.uni-hamburg.de Micha l Warcho l Universié caholique de Louvain Insiu de Saisique, Biosaisique e Sciences Acuarielles Voie du Roman Pays 20, B-1348 Louvain-la-Neuve, Belgium johan.segers@uclouvain.be, m87warchol@gmail.com January 30, 2018 Absrac To draw inference on serial exremal dependence wihin heavy-ailed Markov chains, Drees, Segers and Warcho l [Exremes , ] proposed nonparameric esimaors of he specral ail process. The mehodology can be exended o he more general seing of a saionary, regularly varying ime series. The large-sample disribuion of he esimaors is derived via empirical process heory for cluser funcionals. The finie-sample performance of hese esimaors is evaluaed via Mone Carlo simulaions. Moreover, wo differen boosrap schemes are employed which yield confidence inervals for he pre-asympoic specral ail process: he saionary boosrap and he muliplier block boosrap. The esimaors are applied o sock price daa o sudy he persisence of posiive and negaive shocks. Keywords: Financial ime series; Heavy ails; Muliplier block boosrap; Regular variaion; Shock persisence; Saionary ime series; Tail process. 1 Inroducion The ypical modelling paradigm for a ime series ofen sars by choosing a flexible class of models ha capures salien feaures presen in he daa. Of course, feaures depends on he ype of characerisics one is looking for. For a financial ime series consising of say log-reurns of some asse, he key feaures, ofen referred o as sylized facs, include heavy-ailed marginal disribuions and serially uncorrelaed bu dependen daa. These characerisics are readily deeced using sandard diagnosics such as qq-plos of he marginal disribuion and plos of he sample auocorrelaion funcion ACF of he daa and he squares of he daa. The GARCH process and is varians as well as he sochasic volailiy SV process driven by heavy-ailed noise exhibi hese aribues and ofen serve as a saring poin for building a model. More recenly, considerable aenion has been direced owards sudying he exremal behavior of boh financial and environmenal ime series, especially as i relaes o esimaing risk facors. Exremes for such ime series can occur in clusers and geing a handle on he naure of clusers boh in erms of size and frequency of occurrence is imporan for evaluaing various risk measures. Ulimaely, one wans o choose models ha adequaely describe various exremal dependence feaures observed 1

2 in he daa. The heme of his paper is o provide addiional ools ha no only give measures of exremal dependence, bu can be used as a basis for assessing he qualiy of a model s fi o exremal properies presen in he daa. The exremal index θ 0, 1] Leadbeer, 1983 is one such measure of exremal dependence for a saionary ime series. I is a measure of exremal clusering 1/θ is he mean cluser size of exremes wih θ < 1 indicaing clusering and θ = 1 signifying no clusering in he limi. Unforunaely, θ is a raher crude measure and does no provide fine deail abou exremal dependence. The exremogram, developed in Davis and Mikosch 2009, is an aemp o provide a measure of serial dependence among he exremes in a saionary ime series. I was conceived o be used in much he same way as an ACF in radiional ime series modelling, bu only applied o exreme values. In his paper, we will use he specral ail process, as formulaed by Basrak and Segers 2009 for heavy-ailed ime series, o assess and measure exremal dependence. The specral ail process provides a more in-deph descripion of he srucure of exremal dependence han he exremogram. The firs objecive of his paper will be o esablish limi heory for nonparameric esimaes of he disribuion of he specral ail process for a class of heavy-ailed saionary ime series. This builds on earlier work of Drees e al for heavy-ailed Markov chains. The nonparameric esimaes provide quaniaive informaion abou exremal dependence wihin a ime series and as such can be used in boh exploraory and confirmaory phases of modelling. As an example, i provides esimaes of he probabiliy ha an exreme observaion will occur a ime, given one has occurred a ime 0, and ha is absolue value will be even larger. These esimaes can also be used for model confirmaion, in much he same way ha he ACF is used for assessing qualiy of fi for second-order models of ime series. For example, one can compue a pre-asympoic version o be defined laer of he disribuion of he specral ail process from a GARCH process, which in mos cases can be easily calculaed via simulaion. Then he esimaed disribuion of he specral ail process can be compared wih he pre-asympoic version corresponding o a model for compaibiliy. A good fi would indicae he plausibiliy of using a GARCH model for capuring serial exremal dependence. The second main objecive is hen o provide a useful way of measuring compaibiliy, which we propose using resampling mehods. Recenly, here has been increasing ineres in he economeric lieraure for esimaing quaniies relaed o exremal dependence. For sochasic processes in coninuous ime, Bollerslev e al define a χ-coefficien, derived from he exremogram, for assessing ail dependencies applied o financial ime series. In a follow-up paper ha explores ail risk premia, Bollerslev e al make a connecion beween heir esimaes of he ime-varying ail shape parameers and he exremogram. Linon and Whang 2007 see also Han e al inroduced he quanilogram, a diagnosic ool for measuring direcional predicabiliy in a ime series. In some respecs, our developmen can be viewed as he quanilogram for exreme quaniles. The heory, however, is differen in ha our quaniles are going o infiniy. Neverheless, our work does focus on a ype of direcional predicabiliy, bu only concenraed in he exremes. Tjøsheim and Hufhammer 2013 consider local Gaussian correlaion and relae i o ail index dependence and he exremogram in a ime series conex. Their mehodology is applied o financial ime series. The key objec of sudy in his paper is he ail process and in paricular, is normalized version he specral ail process. A sricly saionary univariae ime series X Z is said o have a ail process Y Z if, for all inegers s, we have L u 1 X s,..., u 1 X X 0 > u d L Ys,..., Y, u, 1.1 wih he implici undersanding ha he law of Y 0 is non-degenerae. The law of Y 0 is hen necessarily Pareoα for some α > 0 and he funcion u P[ X 0 > u] is regularly varying a infiniy wih index α: P[ X 0 > uy] lim = P[ Y 0 > y] = y α, y [1,. 1.2 u P[ X 0 > u] The exisence of a ail process is equivalen o mulivariae regular variaion of he finie-dimensional disribuions of X Z Basrak and Segers, 2009, Theorem 2.1. In many respecs, his condiion 2

3 can be viewed as he heavy-ailed analogue of he condiion ha a process is Gaussian in he sense ha all he finie-dimensional disribuions are specified o be of a cerain ype. The specral ail process is defined by Θ = Y / Y 0, for Z. By 1.1, i follows ha for all inegers s, we have L X 0 /u, X s / X 0,..., X / X 0 X 0 > u d L Y 0, Θ s,..., Θ, u. 1.3 The difference beween 1.1 and 1.3 is ha in he laer equaion, he variables X have been normalized by X 0 raher han by he hreshold u. Such auo-normalizaion allows he ail process o be decomposed ino wo sochasically independen componens, i.e., Y = Y 0 Θ, Z. Independence of Y 0 and Θ Z is saed in Basrak and Segers 2009, Theorem 3.1. The random variable Y 0 characerizes he magniudes of exremes, whereas Θ Z capures serial dependence. The specral ail process a ime = 0 yields informaion on he relaive weighs of he upper and lower ails of X 0 : since Θ 0 = Y 0 / Y 0 = signy 0, we have P[X 0 > u] p = P[Θ 0 = +1] = lim u P[ X 0 > u], 1 p = P[Θ 0 = 1]. 1.4 The disribuions of he forward ail process Y 0 and he backward ail process Y 0 muually deermine each oher Basrak and Segers, 2009, Theorem 3.1. For all i, s, Z wih s 0 and for all measurable funcions f : R s+1 R saisfying fy s,..., y = 0 whenever y 0 = 0, we have, provided he expecaions exis, [ E[f Θ s i,..., Θ i ] = E f Θs Θ i,..., Θ Θ i ] Θ i α 1{Θ i 0}. 1.5 The indicaor variable 1{Θ i 0} can be omied because of he presence of Θ i α, bu someimes, i is useful o menion i explicily in order o avoid errors arising from division by zero. By exploiing he ime-change formula 1.5, we will be able o improve upon he efficiency of esimaors of he specral ail process. Main ineres in his paper is in he cumulaive disribuion funcion cdf, F Θ, of Θ. If F Θ is coninuous a a poin x, hen lim P[X / X 0 x X 0 > u] = P[Θ x] = F Θ x. 1.6 u We consider wo esimaes of F Θ x based on forward and backward represenaions for he ail process. While hese esimaes are asympoically normal, he expressions for he asympoic variances are oo complicaed o be useful for consrucing confidence regions. To overcome his limiaion, inference procedures can be carried ou using resampling mehods. Two resampling mehods for consrucing confidence inervals, based on he saionary boosrap as used in Davis e al. 2012, and he muliplier block boosrap as described in Drees 2015, are applied o our esimaes of F Θ y. In erms of coverage probabiliies, he muliplier block boosrap performed beer han he saionary boosrap procedure in all he cases we considered. However, boh procedures require care when applied for very high hresholds. We apply he mehodology o sudy serial exremal dependence of daily log-reurns on he S&P500 index and he P&G sock price. We disinguish beween wo sources of such dependence posiive and negaive shocks poining ou an asymmeric behavior. Specifically, we consider cases when exreme values posiive or negaive follow posiive/negaive shocks ime lags laer. In erms of he specral ail process, his corresponds o he probabiliies P[±Θ > 1 Θ 0 = ±1]. We illusrae how well he GARCH1,1 model and an exension of i allowing for a leverage effec, he APARCH1,1 model, can capure serial exremal dependence, as measured by hese direcional probabiliies. These examples demonsrae how our mehodology can provide useful 3

4 informaion on he behavior of exremes ha follow boh posiive and negaive shocks, which, in urn, can be used in a model-building conex. The remainder of he paper is organized as follows: The wo esimaes of he ail process are described in Secion 2, while he companion limi heory for hese esimaors is formulaed in Secion 3. The saionary boosrap and muliplier boosrap procedures are presened in Secion 2. The validiy of he proposed boosrap mehodology is esablished in Secion 3 oo. The finie-sample performance is invesigaed hrough Mone Carlo simulaions in Secion 4. The applicaion of our mehodology is provided in Secion 5. The proofs of he main resuls are colleced in Secion 6. 2 Mehodology 2.1 Esimaors The daa consis of a srech X,..., X 1 n+, where is fixed and corresponds o he maximal lag of ineres, drawn from a regularly varying, saionary univariae ime series wih specral ail process Θ Z and index α > 0. In order o esimae p = P[Θ 0 = 1], we simply ake he empirical version of 1.4, yielding ˆp n = 1 > u n 1 > u n. For ˆp n o be consisen and asympoically normal, he hreshold sequence u n should end o infiniy a a cerain rae described in he nex secion. To esimae he cdf, F Θ, of Θ, we propose he forward esimaor ˆF f,θ n x := 1 +/ x, > u n 1 X. 2.1 i > u n This is jus he empirical version of he lef-hand side of 1.6. In equaions 1.1 and 1.3, he condiioning even is { X 0 > u}, making no disincion beween posiive exremes, X 0 > u, and negaive exremes, X 0 < u. However, hese wo cases can be disinguished by condiioning on he sign of Θ 0. In paricular, we define ˆF f,θ Θ0=±1 n x := 1 +/ ± x, ± > u n 1 ±X. 2.2 i > u n The numeraor in he esimaor is a sum of indicaor funcions, mos of which are zero. This ofen leads o a large variance. The ime-change formula 1.5 yields a differen represenaion of he law of Θ, moivaing a differen esimaor han he one above. Depending on he value of x, he new esimaor will involve more non-zero indicaors, which receive weighs insead. The simulaion sudy repored in Secion 4.1 will show ha he resuling esimaor may have a smaller variance han he one in 2.2, in paricular if x is large. Lemma 2.1. Le X Z be a saionary univariae ime series, regularly varying wih index α and specral ail process Θ Z. Then, for all ineger 0, 1 E[ Θ α 1 Θ 0 / Θ > x] if x 0, P[Θ x] = 2.3 E[ Θ α 1 Θ 0 / Θ x] if x < 0. Moreover 1 1 p E[ Θ α 1 1/Θ > x, Θ 0 = 1 ] if x 0, P[Θ x Θ 0 = 1] = 1 p E[ Θ α 1 1/Θ x, Θ 0 = 1 ] if x < 0, 2.4 4

5 and p E[ Θ α 1 1/Θ > x, Θ 0 = 1] if x 0, P[Θ x Θ 0 = 1] = 1 1 p E[ Θ α 1 1/Θ x, Θ 0 = 1] if x < If populaion quaniies are replaced by heir sample counerpars, Lemma 2.1 suggess he following backward esimaor of he cdf of Θ : 1 Xi ˆαn 1 / > x, > u n 1 f x 0, i > u n ˆF b,θ n x := Xi ˆαn 1 / x, > u n 1 > u n if x < 0. Here, ˆα n is an esimaor of he ail index, for which we will ake he Hill-ype esimaor ˆα n = 1 > u n log /u n 1 > u n Condiioning on an exreme value of a specific sign, we ge ˆF b,θ Θ0=±1 n x := 1 ±Xi ˆαn 1 ±Xi / > x, > u n 1 ± > u n Xi ˆαn 1 ±Xi / x, < u n 1 ± > u n if x 0, if x < 0. The asympoic and finie-sample disribuions of hese esimaors will be invesigaed in he following secions. 2.2 Resampling We explore wo differen boosrap schemes ha yield confidence inervals for F Θ x, or raher, for he pre-asympoic version P[X / X 0 x X 0 > u]: he saionary boosrap and he muliplier block boosrap. We apply each of he wo resampling schemes o boh he forward and backward esimaors a various levels x and a differen lags. The saionary boosrap goes back o Poliis and Romano 1994 and is an adapaion of he block boosrap by allowing for random block sizes. The resampling scheme was applied o he exremogram in Davis e al I consiss of generaing pseudo-samples X 1,..., X n, drawn from he sample X 1,..., X n by aking he firs n values in he sequence X K1,..., X K1+L 1 1, X K2,..., X K2+L 2 1,..., where K 1, K 2... is an iid sequence of random variables uniformly disribued on {1,..., n} and L 1, L 2,... is an iid sequence of geomerically disribued random variables independen of K j j N wih disribuion P[L 1 = l] = p 1 p l 1, l = 1, 2,... for some p = p n 0, 1 such ha p n 0 and np n. If he index hus obained exceeds he sample size n, we replace by 1 mod n + 1, i.e., we coninue from he beginning of he sample. The esimaors are hen applied o X 1,..., X n+. The muliplier block boosrap mehod was applied o cluser funcionals in Drees I consiss of spliing he daa se ino m n = n/r n blocks of lengh r n and muliplying he cluser funcionals of each block by a random facor. Here x denoes he ineger par of x. 5

6 Specifically, for iid random variables ξ j, independen of X Z, wih E[ξ j ] = 0 and var [ξ j ] = 1, he boosrapped forward esimaor can be wrien as mn j=1 1 + ξ j 1 Xi+ x, > u n ˆF f,θ n x := mn j=1 1 + ξ j 1 > u n where I j = {j 1r n +1,..., jr n } denoes he se of indices belonging o he jh block. Similarly, he boosrapped backward esimaor for x > 0 wih esimaed index of regular variaion is mn ˆF n b,θ j=1 1 + ξ j Xi ˆα n 1 X > x, i > u n x := 1 mn j=1 1 + ξ j > u n wih mn ˆα n j=1 := 1 + ξ j 1 > u n mn j=1 1 + ξ j log /u n 1 > u n. 2.9 If he hreshold u n is high, i may be advisable o consruc boosrap confidence inervals based on lower hresholds and hen scale accordingly; see he explanaion afer Theorem Tesing for dependence of exreme observaions For iid random variables, he specral ail process simplifies o Θ 0 a.s. for all nonzero. If his occurs for a saionary, regularly varying ime series, hen we say ha he series exhibis serial exremal independence. The opposie case is referred o as serial exremal dependence, i.e., a leas one of he variables Θ for 0 is no degenerae a 0. Since he convergence of he pre-asympoic disribuion can be arbirarily slow, one canno formally es for exremal dependence wihin he presen framework. However, if one wans o es wheher he exceedances over a given high hreshold u are independen hen one may check wheher he lower bound of a confidence inerval for, say, P[ X X 0 X 0 > u] consruced by one of he boosrap mehodologies is larger han his probabiliy under he assumpion of exac independence of X 0 1 X 0 > u and X 1 X > u, which is easily shown o equal P[ X 0 > u] /2. If one prefers o work wih exceedances of he original ime series insead of is absolue values, hen he probabiliy under he assumpion of independence depends on he relaive weighs of he upper and lower ails, and can hus no be calculaed analyically. In ha case, i seems naural o calculae his probabiliy by Mone Carlo simulaion. To his end, one generaes condiionally iid samples according o he empirical disribuion of he original ime series by drawing wih replacemen from he observaions, which corresponds o he classical boosrap procedure for iid daa. The considered probabiliy under independence can be approximaed by he peraining relaive frequency, which is hen compared wih he lower confidence bound for he probabiliy esimaed from he original ime series. If he laer is larger his indicaes ha he exceedances in he ime series exhibi a non-negligible serial dependence see Figure 4 in Secion 5. Alernaively, one may compare he pre-asympoic probabiliy esimaed from he observed ime series using eiher he forward or he backward esimaor wih quaniles of he disribuion of his esimaor under independence. In a similar way as described above, he laer can be approximaed by an empirical quanile obained in Mone Carlo simulaions wih condiionally iid samples cf. Figure 5. 3 Large-sample heory Under cerain condiions, he sandardized esimaion errors of he forward and he backward esimaors converge joinly o a cenered Gaussian process Secion 3.1. Convergence of he muliplier block boosrap follows under he same condiions Secion 3.2. In order no o overload he presenaion, we focus on nonnegaive ime series. We briefly indicae how he condiions and resuls mus be modified in he real-valued case., 6

7 3.1 Asympoic normaliy of he esimaors All esimaors under consideraion can hen be expressed in erms of generalized ail array sums. These are saisics of he form φx n,i, wih X n,i := u 1 n Xi,..., X i,..., + 1Xi > u n. 3.1 Drees and Roozén 2010 give condiions under which, afer sandardizaion, such saisics converge o a cenered Gaussian process, uniformly over appropriae families of funcions φ. From hese resuls we will deduce a funcional cenral limi heorem for he processes of forward and backward esimaors defined in 2.1 and 2.6 wih ˆα n according o 2.7, respecively. To ensure consisency, he hreshold u n mus end o infiniy in such a way ha v n := P[X 0 > u n ] ends o 0, bu he expeced number, nv n, of exceedances ends o infiniy. Moreover, we have o ensure ha observaions which are sufficienly separaed in ime are almos independen. The srengh of dependence will be assessed by he β-mixing coefficiens β n,k := sup E 1 l n k 1 [ sup B B n n,l+k+1 P [ ] B Bn,1 l ] P[B]. Here B j n,i is he σ-field generaed by X n,l i l j. We assume ha here exis sequences l n, r n and some x 0 0 such ha he following condiions hold: Ax 0 The cdf of Θ, F Θ, is coninuous on [x 0,, for {1,..., }. B As n, we have l n, l n = or n, r n = onv n 1/2, r n v n 0, and β n,ln n/r n 0. C For all k {0,..., r n }, here exiss [ X0 { Xk } ] s n k E log max log, 1X k > u n X0 > u n u n u n 3.2 rn such ha s k = lim n s n k exiss, lim n k=1 s nk = k=1 s k holds and he las sum is finie. Moreover, here exiss δ > 0 such ha r n k=1 [ E log + X 0 log + X 1+δ ] 1/1+δ k X0 > u n = O1, n. 3.3 u n u n Wihou Condiion Ax 0 one canno expec uniform convergence of he esimaed cdf of Θ o he rue cdf on [x 0,. Indeed, in his case even P[X /X 0 x X 0 > u] need no converge o F Θ x for a poin of disconinuiy x. Condiion B imposes resricions on he rae a which v n ends o 0 and hus on he rae a which u n ends o. Ofen, he β-mixing coefficiens decay geomerically, i.e., β n,k = Oη k for some η 0, 1. Then one may choose l n = Olog n, and Condiion B is fulfilled for a suiably chosen r n if log n 2 /n = ov n and v n = o1/log n. The echnical Condiion C rules ou oo large a cluser of exreme observaions. Using inegraion by pars, he righ-hand side of 3.2 can be bounded by v 1 n 1 P[X 0 > u n s, X k > u n ] + 1 P[X 0 > u n s, X k > u n ] 1 d s 1 ds. Now one can use echniques employed in Drees 2000 and Drees 2003 o verify 3.2 for specific ime series models like soluions o sochasic recurrence equaions or suiable heavy ailed linear 7

8 ime series. Typically, he upper bounds s n k are of he form ρ k + ξ n for a summable sequence ρ k and ξ n = o1/r n. The lef-hand side of 3.3 can be rewrien in he form r n k=1 1 + δ 2 v 1 n 1 1 P[X 0 > u n s, X k > u n ] log s log δ s 1 ds d 1/1+δ, which can hen be bounded by similar echniques. Under hese condiions, one can prove he asympoic normaliy of relevan generalized ail array sums see Proposiion 6.1 below and hus he join uniform asympoic normaliy of he appropriaely cenered forward and he backward esimaor of F Θ. Theorem 3.1. Le X Z be a saionary, regularly varying process. If Ax 0, B and C are fulfilled for some x 0 0 and y 0 [x 0, 0,, hen f,θ ˆF nv n 1/2 n x P[X /X 0 x X 0 > u n ] x [x 0, b,θ ˆF n y 1 E[X /X 0 α 1X 0 /X > y X 0 > u n ] y [y 0, Zφ 2,x F Θ x Zφ 1 x [x 0, Zφ 3,x F Θ y Zφ 1 + α 2 Zφ 0 αzφ 1 E[logΘ 1Θ > y ] y [y 0, {1,..., } d {1,..., } 3.4 where Z is a cenered Gaussian process, indexed by funcions defined in 6.2, whose covariance funcion is given in 6.3, and F Θ := 1 F Θ denoes he survival funcion of Θ. Asserion 3.4 means ha for suiable versions of he processes he convergence holds uniformly for all x x 0, y y 0 and {1,..., } almos surely. Addiional condiions are needed o ensure ha he biases of he forward and he backward esimaor of F Θ are asympoically negligible: [ sup P X ] x X 0 > u n F Θ x x [x 0, X 0 = o nv n 1/2, 3.5 [ X α1x0 ] sup E /X > y X 0 > u n y [y 0, X F Θ y 0 = o nv n 1/2, 3.6 E[logX0 /u n X 0 > u n ] 1/α = o nv n 1/2, 3.7 for {,..., } \ {0} as n. These condiions are fulfilled if nv n ends o sufficienly slowly, because by definiion of he specral ail process, he regular variaion of X 0 and by 2.3, he lef-hand sides in end o 0 if F Θ is coninuous on [x 0,. Corollary 3.2. Le X Z be a saionary, regularly varying process. If Ax 0, B, C, and are fulfilled for some x 0 0 and y 0 [x 0, 0,, hen nv n 1/2 ˆF f,θ n ˆF b,θ n x F Θ x x [x 0, y F Θ y y [y 0, {1,..., } Zφ 2,x F Θ x Zφ 1 x [x 0, Zφ 3,x F Θ y Zφ 1 + α 2 Zφ 0 αzφ 1 E[logΘ 1Θ > y ] y [y 0, d where Z is he cenered Gaussian process defined in Theorem 3.1. {1,..., } In general, i is difficul o compare he asympoic variances of he backward and he forward esimaor. 8

9 3.2 Consisency of he muliplier block boosrap Here we discuss he asympoic behavior of he muliplier block boosrap version of he forward and backward esimaors. For he sake of breviy, we focus on esimaors of F Θ x for a fixed x. Drees 2015 has shown convergence of boosrap versions of empirical processes of ail array sums under he same condiions needed for convergence of he original empirical processes. Le P ξ denoe he probabiliy w.r.. ξ = ξ j j N, i.e., he condiional probabiliy given X n,i 1 i n. Theorem 3.3. Le ξ j, j N, be iid random variables independen of X Z wih E[ξ j ] = 0 and var [ξ j ] = 1. Then, under he condiions of Theorem 3.1, for all x x 0, y y 0, [ sup P ξ nv n 1/2 f,θ ˆF f,θ n x ˆF n x r, r,s R 2 ] nv n 1/2 b,θ ˆF b,θ n ˆF n y s, {1,..., } [ P nv n 1/2 ˆF f,θ n y x F Θ x r, nv n 1/2 b,θ ˆF n y F Θ y in probabiliy. [ ] b,θ In paricular, if a and b are such ha P ξ ˆF n y [a, b] [ 2 b,θ ˆF n y b, 2 ] b,θ ˆF n y a s, {1,..., }] 0 = β, hen is a confidence inerval for F Θ y wih approximaive coverage level β. However, if he number of exceedances over a given hreshold is oo small, one may prefer o consruc confidence inervals based on boosrap esimaors corresponding o lower hresholds. Le ũ n denoe anoher hreshold sequence, le ṽ n = P[X 0 > ũ n ] denoe he corresponding exceedance probabiliies, and le ˆ F b,θ n y and ˆ F b,θ n y denoe he backward esimaor and he boosrap version hereof, respecively, based on he exceedances over ũ n. The condiional disribuion of nṽ n 1/2 b,θ ˆ F n y ˆ F b,θ n y given he daa is approximaely he same as he uncondiional disribuion of nv n 1/2 b,θ ˆF n y F Θ y. [ ] b,θ Hence, if a and b are such ha P ξ ˆ F n y [a, b] = β and if ˆṽ n /ˆv n is a suiable esimaor of ṽ n /v n, hen 1/2 ˆṽ n ˆ F b,θ n y b b,θ + ˆF n y, ˆv n ˆṽ n ˆv n 1/2 ˆ F b,θ n y a + ˆF b,θ n y 3.8 is a confidence inerval for F Θ y wih approximaive coverage probabiliy β. In pracice, one will ofen use large order saisics as hresholds, say he k n -h and k n -h larges observaions, respecively. In ha case, ṽ n /v n can be replaced by k n /k n. A similar approach, namely o use a variance esimaor which is based on a lower hreshold, has successfully been employed in Drees 2003, Secion 5. Of course, confidence inervals based on he boosrap version of he forward esimaor can be consruced analogously. 9

10 Remark 3.4. I is possible o generalize Theorem 3.3 o cover he join limi disribuion of he boosrap esimaors for all x x 0 and y y 0. Technically, his requires o endow he space of probabiliy measures on spaces of bounded funcions from [x 0, resp. [y 0, o R 2 wih a meric ha induces weak convergence, e.g., he bounded Lipschiz meric. This is he approach in Drees 2015 o esablish he consisency of a boosrap mehod for esimaing he exremogram, a close cousin of he ail process. Based on such a resul, one may consruc uniform confidence bands for he funcion F Θ on [x 0, or [y 0,, respecively, which in general will be considerably wider han he poinwise confidence inervals discussed above and will hus ofen be raher uninformaive. For breviy, we omi he deails. Remark 3.5. For ime series which may ake on negaive values oo, he forward and backward esimaors of F Θ can be represened in erms of generalized ail array sums consruced from X n,i = u 1 n Xi,..., X i,..., + 1 Xi > u n. When x < 0, for example, he backward esimaor ail array sums peraining o he funcions b,θ F n x is equal o he raio of he generalized y,..., y 0,..., y y /y 0 α 1y 0 / y x, y 0 > 1, y,..., y 0,..., y 1 y 0 > 1. Limi heorems can be obained by he same mehods as in he case of non-negaive observaions under obvious analogues o he condiions Ax 0, B and C wih v n := P[ X 0 > u n ]. 4 Finie-sample performance In Secion 4.1, we show resuls from a numerical simulaion sudy designed o es he performance of he forward 2.1 and he backward 2.6 esimaors. We coninue in Secion 4.2 by evaluaing he performance of wo boosrap schemes, he muliplier block boosrap and he saionary boosrap, described in Secion 2.2. The simulaions are based on pseudo-random samples from wo widely used models for financial ime series. Boh models are of he form X = σ Z where σ and Z are independen. Firs, we consider he GARCH1,1 model wih σ 2 = X σ 2 1, he innovaions Z being independen 4 random variables, sandardized o have uni variance. The second model is he sochasic volailiy SV process wih log σ = 0.9 log σ 1 +ɛ, wih independen sandard normal innovaions ɛ and independen innovaions Z wih common disribuion 2.6. The parameers have been chosen o ensure ha boh ime series are regularly varying wih index α = 2.6 Davis and Mikosch, 2001; Mikosch and Sărică, Forward and backward esimaors We esimae P[Θ x] for boh he GARCH1,1 and he SV model, for various argumens x and lags, via he forward and he backward esimaor, wih esimaed ail index α. The hreshold is se a he empirical 95% quanile of he absolue values of a ime series of lengh n = We do Mone Carlo repeiions and calculae bias, sandard [ deviaion, and roo mean squared ] error RMSE wih respec o he pre-asympoic values P X / X 0 x X 0 > F X 0.95 in he forward represenaion. The rue quanile F X 0.95 of X 0 and he rue pre-asympoic values were calculaed numerically via Mone Carlo simulaions based on ime series of lengh I was already repored in he conex of Markovian ime series ha for = 1 and x large, he backward esimaors usually have a smaller variance han he forward esimaors Drees e al., Here, numerical simulaions sugges ha his is rue for non-markovian ime series and for higher lags as well. The resuls are presened in Figure 1. 10

11 P^Θ x, GARCH P^Θ x, GARCH P^Θ x, GARCH bias =1 forward =1 backward =10 forward =10 backward sandard deviaion =1 forward =1 backward =10 forward =10 backward RMSE raio =1 = x x x P^Θ x, SV P^Θ x, SV P^Θ x, SV bias =1 forward =1 backward =10 forward =10 backward sandard deviaion =1 forward =1 backward =10 forward =10 backward RMSE raio =1 = x x x Figure 1: Performance of he forward and he backward esimaors: bias lef, sandard deviaion middle, raio of roo mean square errors righ wih respec o he pre-asympoic values in he forward esimaor, for GARCH1,1 model op and SV model boom. The righ column, which shows he RMSE of he backward esimaor divided by he RMSE of he forward esimaors boh wih respec o he pre-asympoic values of he specral ail process in he forward represenaion, shows ha he backward esimaor ouperforms he forward esimaor if x is sufficienly large in absolue value. This phenomenon was also observed a oher lags no shown. For some oher models, however, such as cerain sochasic recurrence equaion or copula Markov models, he advanage of he backward esimaor was observed only for smaller lags = 1,..., Boosrapped specral ail process We asses he performance of he wo boosrap schemes, he saionary boosrap and he muliplier block boosrap. To do so, we esimae he coverage probabiliy of he boosrapped confidence inervals wih respec o he rue pre-asympoic specral ail process in he forward represenaion. We focus on probabiliies of he form P[ Θ > 1]. This paricular value can be of ineres due o is inerpreaion as he probabiliy of a shock being followed by an even larger afershock, i.e., X being larger han X 0 condiionally on X 0 exceeding some hreshold already. The rue pre-asympoic values in he forward represenaion were calculaed numerically via Mone Carlo simulaions wih ime series of lengh In Figure 2, we plo he resuls for he GARCH1, 1 model and for he SV model, for he forward and he backward esimaors. The expeced block size for he saionary boosrap represened by gray lines was chosen as 100. For he muliplier block boosrap black lines, he block size was fixed a 100 and he muliplier variables ξ j were drawn independenly from he sandard normal disribuion. Esimaes of he coverage probabiliies are based on simulaions. In each such sample, we use boosrap samples for calculaing he confidence inervals wih nominal coverage probabiliy 95%. We use wo differen hresholds, i.e., he 95% and 98% empirical quaniles of he absolue values of a ime series of lengh n = For he higher 11

12 Coverage forward GARCH Coverage backward GARCH muliplier saionary muliplier saionary Coverage forward GARCH Coverage backward GARCH muliplier saionary muliplier rescaled saionary rescaled muliplier saionary muliplier rescaled saionary rescaled Coverage forward SV Coverage backward SV muliplier saionary muliplier rescaled saionary rescaled muliplier saionary muliplier rescaled saionary rescaled Figure 2: Coverage probabiliies of confidence inervals for P[ X /X 0 > 1 X 0 > u n] lef: forward esimaor; righ: backward esimaor based on he saionary boosrap gray and he muliplier block boosrap black. The op and he middle plos correspond o he GARCH1, 1 model wih hresholds se a he 95% and he 98% empirical quaniles, respecively. The boom plos correspond o he SV simulaion sudy wih hreshold se a he 98% empirical quanile. In he laer wo cases, he dashed lines correspond o he coverage probabiliies of he rescaled confidence inervals 3.8. The horizonal black line is he 0.95 reference line. 12

13 Coverage saionary Coverage muliplier block 5 block 10 block 100 block block 5 block 10 block 100 block Median CI widh saionary Median CI widh muliplier block 5 block 10 block 100 block block 5 block 10 block 100 block Figure 3: Coverage probabiliies op and median widhs boom of confidence inervals for P[ X /X 0 > 1 X 0 > u n] based on he saionary boosrap lef column and he muliplier block boosrap righ column for differen block lenghs. The dash doed, dashed, solid, and doed lines represen mean block lenghs 5, 10, 100, and 250 respecively. The plos correspond o he backward esimaor and he GARCH1, 1 model wih hresholds se a he 95% empirical quanile. The horizonal black lines in he op plos are he 0.95 reference lines. hreshold, he confidence inervals were calculaed eiher direcly indicaed by he solid lines or using a rescaled boosrap esimaor ha was based on he exceedances over he 95% empirical quanile as in 3.8 dashed lines. In all cases, he muliplier block boosrap produces a beer coverage probabiliy han he saionary boosrap. Moreover, he backward esimaor is more sable han he forward one, a leas for x = 1, and his ranslaes ino higher sabiliy of he boosrapped confidence inervals. The effec is especially visible for higher hresholds, e.g., a he 98% quanile, leaving insufficienly many pairs of exceedances for accurae inference. Finally, rescaled confidence inervals 3.8 based on lower hresholds can have a much beer coverage han confidence inervals based on higher hresholds. In addiion, in Figure 3 we show coverage probabiliies and median confidence inervals widhs for differen block sizes. The muliplier block boosrap is more robus o he choice of block lengh han he saionary boosrap. In conras o he saionary boosrap, he muliplier block boosrap produces confidence inervals whose coverage probabiliies are fairly sable across differen lags for a given block lengh. I is imporan no o se he block lengh oo low since i can lead o poor coverage probabiliies, especially for higher lags. On he oher hand, oo large a block lengh can resul in confidence 13

14 ω α 1 β 1 δ γ 1 S&P500 GARCH APARCH P&G GARCH APARCH Table 1: Parameers of he models fied o daily log-reurns of he S&P500 index op and he P&G sock price boom. Sandard errors in parenheses. inervals ha are oo wide. 5 Applicaion We firs consider daily log-reurns on he S&P500 sock marke index beween and aken from Yahoo Finance 1. In Figure 4, we plo he sample specral ail process probabiliies P[ Θ > 1 Θ 0 = ±1] and P[±Θ > 1 Θ 0 = ±1] based on he backward esimaor wih 98% empirical quanile aken as a hreshold and he 80% poinwise confidence inervals from he muliplier boosrap scheme rescaled via he 95% quanile as hreshold as in 3.8. The esimaed index of regular variaion is ˆα = The lef-hand plos correspond o condiioning on a posiive exreme a he curren ime insan, whereas he righ-hand plos correspond o condiioning on a negaive shock. The former plos indicae much weaker serial exremal dependence han he laer ones: negaive shocks are more persisen han posiive ones. This is indicaed by he lower bounds of he confidence inervals being above he horizonal lines which correspond o probabiliies under independence. In paricular, he paern of negaive exremes followed by posiive ones is clearly visible; see he righ-hand plo on he second row. The above menioned characerisics are shared by oher sock s daily reurns which were esed bu no repored here. Consider wo widely used financial models of he ype X = σ Z : firs, he GARCH1, 1 process, where σ 2 = ω + α 1 X β 1 σ 2 1, and second, he APARCH1,1 process Ding e al., 1993 wih σ δ = ω + α 1 X 1 γ 1 X 1 δ + β 1 σ δ 1. Boh models allow for volailiy clusering in he limi. Addiionally, he APARCH model capures asymmery in he volailiy of reurns. Tha is, volailiy ends o increase more when reurns are negaive, as compared o posiive reurns of he same magniude if γ 1 > 0. The asymmeric response of volailiy o posiive and negaive shocks is well known in he finance lieraure as he leverage effec of he sock marke reurns Black, We fi hose wo models o daily log-reurns of he S&P500 index. We use he garchfi funcion from he fgarch library available in R, he funcion being based on maximum likelihood esimaion Wuerz e al., The innovaions, Z, are assumed o be sandard normally disribued. The fied parameers are given in he op par of Table 1. In Figure 4 we plo he pre-asympoic specral ail process probabiliies based on he forward esimaor for he fied GARCH and APARCH models, ogeher wih he sample specral ail 1 hp://finance.yahoo.com/ 14

15 P^ Θ > 1 Θ 0 = 1 P^ Θ > 1 Θ 0 = S&P independence GARCH APARCH S&P independence GARCH APARCH P^Θ > 1 Θ 0 = 1 P^Θ > 1 Θ 0 = S&P independence GARCH APARCH S&P independence GARCH APARCH P^Θ < 1 Θ 0 = 1 P^Θ < 1 Θ 0 = S&P independence GARCH APARCH S&P independence GARCH APARCH Figure 4: Sample specral ail process probabiliies solid black bold line for he S&P500 daily log-reurns based on he backward esimaor and he pre-asympoic specral ail process probabiliies of he fied GARCH1,1 doed line and APARCH1,1 dashed line models. The gray area corresponds o he 80% poinwise confidence inervals for he pre-asympoic specral ail probabiliies based on he muliplier boosrap wih replicaions. The op, middle and boom rows concern he condiional probabiliies ha Θ > 1, Θ > 1 and Θ < 1, respecively, given ha Θ 0 = 1 lef column and Θ 0 = 1 righ column. The horizonal line corresponds o he pre-asympoic specral ail probabiliies under independence. 15

16 process probabiliies for S&P500 daily log-reurns esimaed by he backward esimaor. The pre-asympoic values corresponding o he fied models were calculaed numerically via Mone Carlo simulaions wih ime series of lengh Clearly, he APARCH model capures he asymmery which he GARCH model canno. As a second example, we sudy daily log-reurns on he P&G sock price beween and The ail index is esimaed a ˆα = 3.3. We fi he GARCH1,1 and APARCH1,1 models o he ime series and show he esimaed parameers in he boom par of Table 1. In Figure 5 we plo he sample specral ail process probabiliies based on he daily log-reurns hemselves and on he residuals of he fied GARCH1,1 and APARCH1,1 models obained by he backward esimaor. The op-righ plo indicaes ha here is significan serial exremal dependence in he P&G daily log-reurns riggered by he negaive shocks. Due o high asymmery in volailiy, his feaure is sill presen in he residuals of he fied GARCH model whereas i is beer removed by he APARCH filer. 6 Proofs Proof of Lemma 2.1. To prove 2.3, apply he ime-change formula 1.5 wih s = = 0, i = h, and fy 0 = 1y 0 x 10 x o see ha P[Θ h x] 10 x = E[ Θ h α 1Θ 0 / Θ h x] 10 x E[ Θ h α ]. For x 0 in 2.4, apply he ime-change formula 1.5 wih s = h, = 0, i = h and fy h,..., y 0 = 1y 0 > x, y h = 1 o ge P[Θ h > x, Θ 0 = 1] = E[ Θ h α 1Θ 0 / Θ h > x, Θ h > 0] = E [ Θ α h 11/Θ h > x, Θ 0 = 1 ], whereas for x < 0, ake fy h,..., y 0 = 1y 0 x, y h = 1 o obain P[Θ h x, Θ 0 = 1] = E [ Θ α h 1 1/Θ h x, Θ h > 0, Θ 0 = 1 ]. Similarly, in 2.5 choose fy h,..., y 0 = 1y 0 > x, y h = 1 and fy h,..., y 0 = 1y 0 x, y h = 1 for x 0 and x < 0, respecively. Nex we urn o he asympoic normaliy of he forward and backward esimaors. Recall he definiion of X n,i in 3.1. Consider he empirical process Z n ψ := nv n 1/2 n ψxn,i E[ψX n,i ], 6.1 where ψ is one of he following funcions: φ 0 y,..., y 0,..., y = log + y 0, φ 1 y,..., y 0,..., y = 1y0 > 1, φ 2,x y,..., y 0,..., y = 1y /y 0 > x, y 0 > 1, φ 3,x y,..., y 0,..., y = y /y 0 α 1y 0 /y > x, y 0 > for {1,..., } and x 0. The asympoic behavior of Z n can be derived from more general resuls by Drees and Roozén Proposiion 6.1. Le X Z be a non-negaive, saionary, regularly varying ime series wih ail process Y Z. Assume ha condiions Ax 0, B and C are fulfilled for some x 0 0. Then, for all y 0 [x 0, 0,, he sequence of processes Z n φ 0, Z [ n φ 1, Z n φ 2,x x [x0,, Z ] n φ 3,y y [y0, {1,..., } 16

17 P^ Θ > 1 Θ 0 = 1 P^ Θ > 1 Θ 0 = P&G upper CL for independence P&G upper CL for independence P^ Θ > 1 Θ 0 = 1 P^ Θ > 1 Θ 0 = GARCH residuals upper CL for independence GARCH residuals upper CL for independence P^ Θ > 1 Θ 0 = 1 P^ Θ > 1 Θ 0 = APARCH residuals upper CL for independence APARCH residuals upper CL for independence Figure 5: Sample specral ail process black line for P&G daily log-reurns op, GARCH1,1 residuals middle, and APARCH1,1 residuals boom based on he backward esimaor. Plos in he firs column represen condiioning on a posiive shock whereas in he second column one condiions on a negaive shock. The horizonal gray lines correspond o he empirical 80% quanile of he backward esimaor under independence obained from simulaions. 17

18 converges weakly o a cenered Gaussian process Z wih covariance funcion given by cov Zψ 1, Zψ 2 = j= for all ψ 1, ψ 2 { φ 0, φ 1, φ 2,x, φ 3,y x x 0, y y 0, {1,..., } }. E [ ψ 1 Y,..., Y ψ 2Y j,..., Y j+ ] =: cψ 1, ψ The weak convergence saemens in Proposiion 6.1 hold in he space of bounded funcions on { φ 0, φ 1, φ 2,x, φ 3,y x x 0, y y 0, {1,..., } } equipped wih he supremum norm; see van der Vaar and Wellner 1996, Secion 1.5 for deails. Proof of Proposiion 6.1. One can argue similarly as in he proof of Proposiion B.1 of Drees e al. 2015, because he asympoic equiconinuiy of he process can be esablished for each separaely. Noe ha he discussion in Drees and Roozén 2016 shows ha par ii of condiion B of Drees e al is no needed. By saionariy, he covariance of Zψ 1 and Zψ 2 is obained as he limi of 1 E r n v n r n ψ 1 X n,i r n j=1 ψ 2 X n,j = 1 r n 1 v n k= r n+1 1 k r n E[ψ 1 X n,0 ψ 2 X n,k ]. This sum can be shown o converge o cψ 1, ψ 2 using Pra s lemma and Condiion C, as in Drees e al Remark 6.2. The covariances can be expressed in erms of he specral ail process. For example, cφ 3,x, φ 0 = = E [ Θ α 11/Θ > x log + Y 0 Θ j ] j= j= E [ Θ α 11/Θ > x Θ α j 1 log + Θ j + α 1]. Here we have used ha Y 0 is independen of Θ s s Z wih disribuion P[Y 0 > y] = y α for y 1. Theorem 3.1 and Corollary 3.2 can now be proved in he same way as Theorem 4.5 in Drees e al We omi he deails, which can also be inferred from he more involved discussion of he boosrap esimaor below. Drees 2015 has shown ha under roughly he same condiions as used by Drees and Roozén 2010, condiionally on he daa, he following boosrap version of he empirical process Z n has he same asympoic behavior as Z n : m n Z n,ξ ψ := nv n 1/2 ξ j ψxn,i E[ψX n,i ], 6.4 j=1 wih I j := {j 1r n + 1,..., jr n } and m n := n/r n. In wha follows, he symbol E ξ denoes he expecaion w.r.. ξ = ξ j j N, i.e., he expecaion condiionally on X n,i 1 i n. Moreover, le BL 1 denoe he se of all funcions g : R 4 +2 R such ha sup z R 4 +2 gz 1 and gz 1 gz 2 z 1 z 2 for all z 1, z 2 R Proposiion 6.3. Suppose ha X Z is a non-negaive, saionary, regularly varying ime series and ha he condiions Ax 0, B and C are fulfilled for some x 0 0. Then, for all x x 0 and all y 0 [x 0, 0,, one has Z n,ξ φ 0, Z n,ξ φ 1, [ Z n,ξ φ 2,x, Z n,ξ φ 3,y ] {1,..., } d Zφ 0, Zφ 1, [ Zφ 2,x, Zφ 3,y ] {1,..., } 18

19 wih Z as defined in Theorem 3.1. Moreover, sup E ξg Z n,ξ φ 0, Z n,ξ φ 1, [ Z n,ξ φ 2,x, Z n,ξ φ 3,y ] {1,..., } g BL 1 E g Zφ 0, Zφ 1, [ Zφ 2,x, Zφ 3,y ] {1,..., } in probabiliy. Proposiion 6.3 follows immediaely from Drees 2015, Theorem 2.1, because in he proof of Proposiion 6.1 cf. he proof of Proposiion B.1 of Drees e al i is shown ha he assumpions of Drees 2015, Theorem 2.1 follow from he condiions of Proposiion 6.3. Now we are ready o prove he consisency of he muliplier block boosrap procedure. Proof of Theorem 3.3. We only prove consisency of he boosrap version of he backward esimaor, as he proof for he forward esimaor is considerably simpler. For simpliciy, we assume ha n = m n r n. Le α n := 1 E [ v log + ] n = X 0 /u n X 0 > u n E[φ 0 X n,1 ]. Recall Z n and Z n,ξ in 6.1 and 6.4 respecively, recall I j = {j 1r n + 1,..., jr n }, and recall ˆα n and ˆα n in 2.7 and 2.9, respecively. Then mn nv n 1/2 ˆα n ˆα n = nv n 1/2 j=1 ξ mn j 1 > u n ˆα n mn j=1 ξ j j=1 1 + ξ j log + /u n mn log + /u n = Z n,ξφ 1 ˆα n Z n,ξ φ 0 + r n v n 1/2 m 1/2 n j=1 ξ j1 ˆα n /α n αn m 1 mn n j=1 ξ j + nv n 1/2 { Z. n φ 0 + Z n,ξ φ 0 } Since mn 1/2 mn j=1 ξ j and Z n are sochasically bounded and ˆα n α in probabiliy, he assumpions nv n, r n v n 0, and α n α, m 1/2 n and Proposiion 6.3 ensure ha nv n 1/2 ˆα n ˆα n = αz n,ξ φ 1 α 2 Z n,ξ φ 0 + o P 1, 6.6 which converges weakly o αzφ 1 α 2 Zφ 0. Moreover, condiionally on he daa, i converges o he same limi weakly in probabiliy in he sense of 6.5. b,θ Nex, recall ˆF b,θ n y and ˆF n y in 2.6 and 2.8, respecively. For y > 0, we have 1 I follows ha ˆF n b,θ = y [ m n b,θ ˆF n y b,θ ˆF n Xi + ξ j j=1 1 > u n = y ˆα n Xi {1 ˆαn Xi n /X ˆαn i 1 / > y, > u n. 1 / > y, > u n ˆα n 1 / > y, > u n ˆF b,θ n ] m n / m n y} ξ j 1 > u n 1 + ξ j 1 > u n. j=1 j=1 19

20 For any pair α, α such ha 0 < α < α < α, here exiss a consan 0 < C < such ha for all α [α, α] and, for suiable consans λ = λ α 0, 1, we have, on he even { / > y 0 }, α α α Xi Xi Xi Xi log α α = 1 α+λ α α Xi log 2 Xi α α 2 C α α 2. 2 Hence wih ˆF n b,θ = y [ m n b,θ ˆF n Xi + ξ j j=1 1 ˆF b,θ n y α log { Xi Xi α + ˆα n ˆα n 1 / > y, > u n Xi α log Xi } ˆα n α 1 / > y, > u n ] m n / m n y ξ j 1 > u n + R n y 1 + ξ j 1 > u n 6.7 j=1 j=1 R n y Cˆα n ˆα n 2 1 / > y, > u n Consider he funcion m n + Cˆα n α 2 ξ j 1 / > y, > u n j=1 = O P nvn 1 nv n + nv n 1 m n r n v n = OP 1, n. φ 4,x y,..., y 0,..., y = y /y 0 α log y /y 0 1y 0 /y > x, y > 0, y 0 > 1. One may show as in he proof of Proposiion 6.1 ha Z n φ 4,y and Z n,ξ φ 4,y boh converge weakly o Zφ 4,y. In paricular, as n, n α nv n 1 Xi Xi log 1 / > y, > u n X i [ X α X ] = E log 1X 0 /X > y X 0 > u n + O P nvn 1/2 X 0 X 0 E [ Θ α logθ 11/Θ > y ] = E[logΘ 1Θ > y], where he las sep follows from he ime-change formula 1.5 applied wih fy 0 = logy 0 1y 0 > y and, 0, insead of s,, i. Therefore n α Xi Xi log ˆα n ˆα n 1 / > y, > u n = nv n 1/2 E[logΘ 1Θ > y] + o P 1 nv n 1/2 ˆα n ˆα n. 6.8 Likewise, one can conclude ha m n α nv n 1/2 Xi ξ j log j=1 Xi 1 / > y, > u n 20

21 m n [ X α = Z n,ξ φ 4,y + rv n 1/2 m 1/2 n ξ j E log X 0 = O P 1. As a consequence, m n ξ j j=1 { Xi α + Xi j=1 α log Xi X X 0 ] 1X 0 /X > y X 0 > u n } ˆα n α 1 / > y, > u n m n = nv n 1/2 Z n,ξ φ [ 3,y + ξ j r n v n E Θ α 11/Θ > y ] + o1 + O P 1 j=1 = nv n 1/2 Z n,ξ φ 3,y + O P rn v n 1/2 + O P nvn 1/ Moreover, we find, as r n v n 0 and m n j=1 ξ j = O P m 1/2 n, ha m n ξ j j=1 1 > u n = nv n 1/2 Z n,ξ φ 1 + r n v n m n ξ j j=1 = nv n 1/2 Z n,ξ φ 1 + o P 1, n The denominaor of 6.7 equals nv n + O P nv n 1/2. Combining and 6.6 yields nv n 1/2 b,θ ˆF b,θ n y ˆF n y = E[logΘ 1Θ > y] nv n 1/2 ˆα n ˆα n + Z n,ξ φ b,θ 3,y 1 ˆF n y Z n,ξ φ 1 + o P 1 = Z n,ξ φ 3,y F Θ y Z n,ξ φ 1 E[logΘ 1Θ > y] αz n,ξ φ 1 α 2 Z n,ξ φ 0 + o P 1. Now he asserion is a direc consequence of Proposiion 6.3 and Theorem 3.1. Acknowledgemens The auhors wish o hank he ediors and he referees for heir careful reading and for various consrucive commens and useful suggesions. J. Segers graefully acknowledges funding by conrac Proje d Acions de Recherche Concerées No. 12/ of he Communaué française de Belgique and by IAP research nework Gran P7/06 of he Belgian governmen Belgian Science Policy. The research of M. Warcho l was funded by a PhD gran of he Fonds de la Recherche Scienifique F.R.S.-FNRS. H. Drees was parially suppored by DFG research gran JA 2160/1. R. Davis was suppored in par by ARO MURI gran W911NF References Basrak, B. and J. Segers Regularly varying mulivariae ime series. Sochasic Processes and Their Applicaions 119 4, Black, F Sudies of sock price volailiy changes. In Proceedings of he 1976 Meeings of he American Saisical Associaion, Business and Economics Secion, pp Bollerslev, T., V. Todorov, and S. Z. Li Jump ails, exreme dependencies, and he disribuion of sock reurns. Journal of Economerics 172 2, Bollerslev, T., V. Todorov, and L. Xu Tail risk premia and reurn predicabiliy. Journal of Financial Economics 118 1, Davis, R. A. and T. Mikosch Poin process convergence of sochasic volailiy processes wih applicaion o sample auocorrelaion. Journal of Applied Probabiliy 38A,

22 Davis, R. A. and T. Mikosch The exremogram: a correlogram for exreme evens. Bernoulli 15 4, Davis, R. A., T. Mikosch, and I. Cribben Towards esimaing exremal serial dependence via he boosrapped exremogram. Journal of Economerics 170 1, Ding, Z., C. W. Granger, and R. F. Engle A long memory propery of sock marke reurns and a new model. Journal of Empirical Finance 1 1, Drees, H Weighed approximaions of ail processes for β-mixing random variables. The Annals of Applied Probabiliy 10 4, Drees, H Exreme quanile esimaion for dependen daa wih applicaions o finance. Bernoulli 9 4, Drees, H Boosrapping empirical processes of cluser funcionals wih applicaion o exremograms. arxiv preprin arxiv: Drees, H. and H. Roozén Limi heorems for empirical processes of cluser funcionals. The Annals of Saisics 38 4, Drees, H. and H. Roozén Correcion noe o limi heorems for empirical processes of cluser funcionals. The Annals of Saisics o appear. Drees, H., J. Segers, and M. Warcho l Saisics for ail processes of Markov chains. Exremes 18 3, Han, H., O. Linon, T. Oka, and Y.-J. Whang The cross-quanilogram: Measuring quanile dependence and esing direcional predicabiliy beween ime series. Journal of Economerics 193 1, Leadbeer, M Exremes and local dependence in saionary sequences. Zeischrif für Wahrscheinlichkeisheorie und Verwande Gebiee 65 2, Linon, O. and Y.-J. Whang The quanilogram: Wih an applicaion o evaluaing direcional predicabiliy. Journal of Economerics 141 1, Mikosch, T. and C. Sărică Limi heory for he sample auocorrelaions and exremes of a GARCH 1,1 process. The Annals of Saisics 28 5, Poliis, D. N. and J. P. Romano The saionary boosrap. Journal of he American Saisical associaion , Tjøsheim, D. and K. O. Hufhammer Local gaussian correlaion: a new measure of dependence. Journal of Economerics 172 1, van der Vaar, A. W. and J. A. Wellner Weak Convergence of Empirical Processes. New York: Springer. Wuerz, D., Y. C. wih conribuion from Michal Miklovic, C. Boud, P. Chausse, and ohers fgarch: Rmerics - Auoregressive Condiional Heeroskedasic Modelling. R package version