Nonparametric estimation of extreme risks from heavy-tailed distributions
|
|
- Constance Franklin
- 5 years ago
- Views:
Transcription
1 Nonparametric estimation of extreme risks from heavy-tailed distributions Laurent GARDES joint work with Jonathan EL METHNI & Stéphane GIRARD December 2013
2 1 Introduction to risk measures 2 Introduction to Extreme Value Analysis 3 Estimation of risk measures 4 Extrapolation 5 Application 2 / 41
3 Plan 1 Introduction to risk measures 2 Introduction to Extreme Value Analysis 3 Estimation of risk measures 4 Extrapolation 5 Application 3 / 41
4 The Value-At-Risk Let Y R be a random loss variable. For example, Y can represent the loss of a portfolio during some fixed interval time. A risk measure is a real-valued function R(.) quantifying the risk associated to a given loss variable. The most known risk measure is the Value-at-Risk of level α (0, 1) is the α-quantile defined by R(Y ) = VaR α(y ) := F (α) = inf{t, F (t) α}, where F (.) is the generalized inverse of the survival function of Y. 4 / 41
5 Illustration of the Value-At-Risk 5 / 41
6 Drawback of the Value-At-Risk Let Y 1 and Y 2 be two loss variables with survival functions F 1(.) and F 2(.) represented below : Full line : F 1(.), dashed line : F 2(.). For α 0.22, VaR α(y 1) = VaR α(y 2). Clearly, the loss Y 1 is more dangerous than Y 2. 6 / 41
7 Other risk measures To overcome this drawback, other risk measures have been proposed : The Conditional Tail Expectation of level α (0, 1) defined by CTE α(y ) := E(Y Y > VaR α(y )). The Conditional-Value-at-Risk of level α (0, 1) introduced by Rockafellar et Uryasev [2000] defined by with 0 λ 1. CVaR λ,α (Y ) := λvar α(y ) + (1 λ)cte α(y ), The Conditional Tail Variance of level α (0, 1) introduced by Valdez [2005] is defined by CTV α(y ) := E{(Y CTE α(y )) 2 Y > VaR α(y )}. 7 / 41
8 First contribution : definition of the Conditional Tail Moment The first goal of this work is to unify the definitions of the previous risk measures. To this end, the Conditional Tail Moment of level α (0, 1) is introduced : CTM a,α(y ) := E(Y a Y > VaR α(y )), where a 0 is such that the moment of order a of Y exists. All the previous risk measures of level α can be rewritten as CTE α(y ) = CTM 1,α(Y ), CVaR λ,α (Y ) = λvar α(y ) + (1 λ)ctm 1,α(Y ), CTV α(y ) = CTM 2,α(Y ) CTM 2 1,α(Y ). = All the previous risk measures depend on the VaR and the CTM. 8 / 41
9 Second contribution : Risk measure for extreme losses in a regression case Generalization of risk measures : 1 For extreme losses : the fixed level α (0, 1) is replaced by a sequence α n 0. n 2 Regression case : a covariate X R p is recorded simultaneously with Y. Denoting by F (. x) the conditional survival distribution function of Y given X = x, the Regression Value-at Risk is defined by : RVaR αn (Y x) := F (α n x) = inf{t, F (t x) α n}, and the Regression Conditional Tail Moment of order a is defined by : RCTM a,αn (Y x) := E(Y a Y > RVaR αn (Y x), X = x), where a > 0 is such that the moment of order a of Y exists. 9 / 41
10 Second contribution : Risk measure for extreme losses in a regression case This yields the following risk measures : RCTE αn (Y x) = RCTM 1,αn (Y x), RCVaR λ,αn (Y x) = λrvar αn (Y x) + (1 λ)rctm 1,αn (Y x), RCTV αn (Y x) = RCTM 2,αn (Y x) RCTM 2 1,α n (Y x). All the risk measures depend on the RVaR and the RCTM. To deal with extreme losses (i.e. when α n 0), we need to use Extreme Value Analysis results. 10 / 41
11 Plan 1 Introduction to risk measures 2 Introduction to Extreme Value Analysis 3 Estimation of risk measures 4 Extrapolation 5 Application 11 / 41
12 Extreme Value Analysis Let Z 1,..., Z n be n independent copies of a random value Z R with survival function F (.). We denote by Z 1,n... Z n,n the ordered statistics. The main goal of Extreme Value Analysis is to estimate extreme quantiles defined by : F (α n) = VaR αn (Z), for α n 0. The simplest way to estimate quantile consists in using the empirical survival function ˆ F n(z) := 1 n I{Z i > z}. n i=1 Clearly, this is not a good idea since when α n < 1/n, ˆ F n (α n) = Z n,n. 12 / 41
13 Extreme Value Analysis We thus need some information beyond the maximum observation. This information is given by the following result. Fisher-Tippett theorem If there exist two sequences (a n) > 0 and (b n) and a non-degenerate distribution function H(.) such that ( ) Zn,n b n P z H(z), as n goes to infinity, a n then { H(y) = where z + = max(z, 0). [ exp { 1 + γ z µ σ exp { exp ( z µ σ } 1/γ )} + ] if γ 0, if γ = 0, The distribution H(.) is called the Extreme Value Distribution (EVD). The shape parameter γ R is called the extreme value index. It controls the heaviness of the tail. 13 / 41
14 Extreme Value Analysis If γ > 0, F (z) decreases at a polynomial rate to zero as z goes to infinity. More precisely, F (z) = z 1/γ L(z), where L(.) is a slowly varying function satisfying for all t > 0 : L(tz) lim z L(z) = 1. Such distributions are known as heavy-tail distributions (Pareto, student, log-gamma, etc...) If γ < 0, F (z) decreases at a polynomial rate to zero as z goes to the right endpoint y F < (Uniform, Beta, etc...) If γ = 0, F (z) decreases at an exponential rate to zero as z goes to the right endpoint (Normal, exponential, gamma, etc...) 14 / 41
15 Plan 1 Introduction to risk measures 2 Introduction to Extreme Value Analysis 3 Estimation of risk measures 4 Extrapolation 5 Application 15 / 41
16 Recall... Recall that we are interested in the quantities : and RVaR αn (Y x) := F (α n x) = inf{t, F (t x) α n}, RCTM a,αn (Y x) := E(Y a Y > RVaR αn (Y x), X = x) = 1 E (Y a I{Y > RVaR αn (Y x)} X = x), α n where a > 0 is such that the moment of order a of Y exists. 16 / 41
17 Conditional moment We introduce the conditional moment of order a 0 of Y given X = x ϕ a(y x) = E (Y a I{Y > y} X = x), where I{.} is the indicator function. Since ϕ 0(y x) = F (y x), it follows RVaR αn (Y x) = ϕ 0 (α n x), RCTM a,αn (Y x) = 1 ϕ a(ϕ 0 (α n x) x). α n Goal : estimate ϕ a(. x) and ϕ a (. x). 17 / 41
18 Inference Estimator of ϕ a(. x) : We propose to use a classical kernel estimator given by n ( ) / n ( x Xi ϕ a,n(y x) = K Yi a x Xi I{Y i > y} K i=1 h n h n is a sequence called the window-width such that h n 0 as n, K is a bounded density on R p with support included in the unit ball of R p. Estimator of ϕ a (. x) : i=1 h n ). Since ˆϕ a,n(. x) is a non-increasing function, an estimator of ϕ a defined for α (0, 1) by (α x) can be ˆϕ a,n(α x) = inf{t, ˆϕ a,n(t x) < α}. 18 / 41
19 Inference The Regression Value-At-Risk of order α n is thus estimated by : RVaR αn (Y x) = ϕ 0 (α n x), and the Regression Conditional Tail Moment of order a is estimated by : RCTM a,αn (Y x) = 1 α n ϕ a,n ( ϕ 0 (α n x) x) We thus can now estimate all the above mentioned risk measures. To do... Find the asymptotic joint distribution of {( ) } RCTM aj,α n (Y x), j = 1,..., J, RVaR αn (Y x), where 0 a 1 < a 2 <... < a J. 19 / 41
20 Heavy-tail assumptions (F.1) The conditional survival distribution function of Y given X = x is assumed to be heavy-tailed i.e. for all λ > 0, F (y x) = y 1/γ(x) L(y x), where for x fixed, L(. x) is a slowly varying function. In this context, γ(.) is a positive function of the covariate x and is referred to as the conditional tail index. Condition (F.1) also implies that for a [0, 1/γ(x)) and all α (0, 1), RCTM a,α(y x) exists and, RCTM a,α(y x) = α aγ(x) l a(α 1 x), where for x fixed, l a(. x) is a slowly-varying function 20 / 41
21 Heavy-tail assumptions (F.2) l a(. x) is normalized for all a [0, 1/γ(x)). In such a case, the Karamata representation of the slowly-varying function can be written as ( y ) ε a(u x) l a(y x) = c a(x) exp du, u where c a(.) is a positive function and ε a(y x) 0 as y. 1 (F.3) ε a(. x) is continuous and ultimately non-increasing for all a [0, 1/γ(x)). 21 / 41
22 Regularity assumptions A Lipschitz condition on the probability density function g of X is also required : (L) There exists a constant c g > 0 such that g(x) g(x ) c g d(x, x ). where d(x, x ) is the Euclidean distance between x and x. Finally, for y > 0 and ξ > 0, the largest oscillation of the conditional moment of order a [0, 1/γ(x)) is defined by { } ϕ a(z x) ω n(y, ξ) = sup ϕ a(z x ) 1, z [(1 ξ)y, (1 + ξ)y] and d(x, x ) h. 22 / 41
23 Main result Theorem 1 : Suppose (F.1), (F.2) and (L) hold. Let Then, 0 a 1 < a 2 < < a J, x R p such that g(x) > 0 and 0 < γ(x) < 1/(2a J ), α n 0 and nh p nα n as n, ξ > 0 such that nh p nα n {h n ω n(rvar(α n x), ξ)} 0, ) RCTM aj,α n (Y x) nhp α n ( RCTM aj,α n (Y x) 1 j {1,...,J} is asymptotically Gaussian, centered, with covariance matrix K 2 2γ 2 (x)σ(x)/g(x) where 23 / 41 Σ(x) = a i a j (2 (a i +a j )γ(x)) (1 (a i +a j )γ(x)), ( ) RVaR αn (Y x) RVaR α(y x) 1 a 1. a J a 1 a J 1.
24 Conditions on the sequences α n and h n nh p nα n : Necessary and sufficient condition for the almost sure presence of at least one point in the region B(x, h n) [RVaR αn (Y x), + ) of R p R. nh p n α n {h n ω n(rvar αn (Y x), ξ)} 0 : The biais induced by the smoothing is negligible compared to the standard-deviation. 24 / 41
25 Consequences Suppose the assumptions of Theorem 1 hold. Then, if 0 < γ(x) < 1/2, ( ) ( RCTE αn (Y x) nhp α n RCTE αn (Y x) 1 d N 0, 2(1 ) γ(x))γ2 (x) K γ(x) g(x) ( ) ( RCVaR λ,αn (Y x) nhp α n RCVaR λ,αn (Y x) 1 d N 0, γ2 (x)(λ λ 2γ(x)) 1 2γ(x) The RCTV αn (Y x) estimator involves the computation of a second order moment, it requires the stronger condition 0 < γ(x) < 1/4, ( ) ( ) RCTV αn (Y x) nhp α n RCTV αn (Y x) 1 d K 2 2 N 0, V γ(x), g(x) ) K 2 2 g(x) where V γ(x) = 8(1 γ(x))(1 2γ(x))(1 + 2γ(x) + 3γ2 (x)). (1 3γ(x))(1 4γ(x)) 25 / 41
26 Plan 1 Introduction to risk measures 2 Introduction to Extreme Value Analysis 3 Estimation of risk measures 4 Extrapolation 5 Application 26 / 41
27 A Weissman type estimator In Theorem 1, the condition nh p nα n provides a lower bound on the level of the risk measure to estimate. In consequence, α n must be an order of an extreme quantile within the sample. This restriction is a consequence of the use of a kernel estimator which cannot extrapolate beyond the maximum observation in the ball B(x, h n). Definition Let us consider (α n) n 1 and (β n) n 1 two positive sequences such that α n 0, β n 0 and 0 < β n < α n. A kernel adaptation of Weissman s estimator [1978] is given by ( ) aˆγn(x) αn RCTM W a,β n (Y x) = RCTM a,αn (Y x) β n }{{} extrapolation 27 / 41
28 Idea Recall that RCTM a,αn (Y x) = α aγ(x) n l a(α 1 x). n Thus, RCTM a,βn (Y x) RCTM a,αn (Y x) = ( αn β n ) aγ(x) l a(βn 1 x) l a(αn 1 x). Since l a(. x) is a slowly varying function, RCTM a,βn (Y x) RCTM a,αn (Y x) ( αn β n ) aγ(x). Replacing RCTM a,αn (Y x) and γ(x) by their estimators lead to the Weissmann estimator. 28 / 41
29 Asymptotic result Theorem 2 : Suppose the assumptions of Theorem 1 hold together with (F.3). Let ˆγ n(x) be an estimator of the conditional tail index such that ( ) nh p n α n(ˆγ n(x) γ(x)) d N 0, v 2 (x), with v(x) > 0. If, moreover (β n) n 1 is a positive sequence such that β n 0 and β n/α n 0 as n, then nh p n α n RCTM W a,β n (Y x) ( log(α n/β n) RCTM a,βn (Y x) 1 d N 0, (av(x)) 2). The condition β n/α n 0 allows us to extrapolate and choose a level β n arbitrarily small. 29 / 41
30 Estimation of the conditional tail index Without covariate : Hill [1975] Let Z 1,n Z n,n be the order statistics associated with i.i.d. heavy-tailed random variables Z 1,..., Z n and (α n) n 1 (0, 1). The Hill estimator is given by ˆγ (H) n, nα = 1 nα n n log Z n j+1,n log Z n nαn,n. nα n j=1 With a covariate : Z n j+1,n RVaR τj α n (Y x) A kernel version of the Hill estimator is given by / J J ˆγ n,αn (x) = (log RVaR τj α n (Y x) log RVaR τ1 α n (Y x)) log(τ 1/τ j ), j=1 j=1 where J 1 and (τ j ) j 1 is a decreasing sequence of weights. 30 / 41
31 Extrapolation The asymptotic normality of ˆγ n,αn (x) has been established by Daouia et al. [2011]. Suppose the assumptions of Theorem 1 hold together with (F.3), ( ) nh p n α n(ˆγ n,αn (x) γ(x)) d N 0, K 2 2γ 2 (x)v J /g(x), with V J = J j=1 ( ) /( 2(J j) + 1 J J 2 log(1/τ j ) τ j j=1 ) 2 In particular, if β n 0 and β n/α n 0, we obtain nh p n α n RCTE W β n (Y x) ( ) RCTE βn (Y x) 1 d N 0, K 2 2γ 2 (x)v J /g(x). 31 / 41
32 Plan 1 Introduction to risk measures 2 Introduction to Extreme Value Analysis 3 Estimation of risk measures 4 Extrapolation 5 Application 32 / 41
33 Daily rainfalls in the Cévennes-Vivarais region The Cévennes-Vivarais region raingauge stations Daily rainfalls (in mm) measured at N = 523 stations from 1958 to / 41
34 Objective The loss variable Y is the daily rainfall (in mm). The random value Y is observed with a covariate X = (latitude,longitude,altitude,etc...). Goal : for all x, estimate the following risk measures : RVaR β (Y x) with β = 1/( T ) T -year return level in hydrology. RCTE β (Y x) with β = 1/( T ) average rainfall over the T -year return level in hydrology. We consider the case T = 100 years (extrapolation needed). 34 / 41
35 A cross validation procedure to choose h n and α n : Step 1 Double loop on H = {h i ; i = 1,..., M} and on A = {α j ; j = 1,..., R}. Loop on all raingauge stations {x t; t = 1,..., N}. 35 / 41
36 A cross validation procedure to choose h n and α n : Step 1 Double loop on H = {h i ; i = 1,..., M} and on A = {α j ; j = 1,..., R}. Loop on all raingauge stations {x t; t = 1,..., N}. Consider one raingauge station x t 35 / 41
37 A cross validation procedure to choose h n and α n : Step 1 Double loop on H = {h i ; i = 1,..., M} and on A = {α j ; j = 1,..., R}. Loop on all raingauge stations {x t; t = 1,..., N}. Remove all other raingauge stations 35 / 41
38 A cross validation procedure to choose h n and α n : Step 1 Double loop on H = {h i ; i = 1,..., M} and on A = {α j ; j = 1,..., R}. Loop on all raingauge stations {x t; t = 1,..., N}. Remove all other raingauge stations Estimate γ > 0 using the classical Hill estimator. It only depends on α j. 35 / 41
39 A cross validation procedure to choose h n and α n : Step 1 Double loop on H = {h i ; i = 1,..., M} and on A = {α j ; j = 1,..., R}. Loop on all raingauge stations {x t; t = 1,..., N}. Remove all other raingauge stations Estimate γ > 0 using the classical Hill estimator. It only depends on α j. = We obtain ˆγ (H) n,t,α j 35 / 41
40 A cross validation procedure to choose h n and α n : Step 2 Remove the station x t 36 / 41
41 A cross validation procedure to choose h n and α n : Step 2 Work in B(x t, h i ) \ {x t} 36 / 41
42 A cross validation procedure to choose h n and α n : Step 2 Work in B(x t, h i ) \ {x t} Estimate γ(x) > 0 using the kernel version of the Hill estimator. It depends on α j and on h i. 36 / 41
43 A cross validation procedure to choose h n and α n : Step 2 Work in B(x t, h i ) \ {x t} Estimate γ(x) > 0 using the kernel version of the Hill estimator. It depends on α j and on h i. = We obtain ˆγ n,hi,α j (x t) 36 / 41
44 A cross validation procedure to choose h n and α n : Step 2 Work in B(x t, h i ) \ {x t} Estimate γ(x) > 0 using the kernel version of the Hill estimator. It depends on α j and on h i. = We obtain ˆγ n,hi,α j (x t) (h emp, α emp) = arg min median{(ˆγ (H) n,t,α j ˆγ n,hi,α j (x t)) 2, t {1,..., N}}. (h i,α j ) H A 36 / 41
45 Computation of RVaR W and RCTE W 523 Stations Regular grid : Two dimensional covariate X =(latitude, longitude). Bi-quadratic kernel : K(x) (1 x 2 ) 2 I { x 1}. Harmonic sequence of weights : (τ j ) j {1,...,9} = 1/j. Results of the procedure (h emp, α emp) = (24, 1/( )). 37 / 41
46 Estimated conditional tail index 38 / 41
47 RVaR W β (Y x) : 100-year return level 39 / 41
48 RCTE W β (Y x) above the 100-year return level 40 / 41
49 References Daouia, A., Gardes, L., Girard, S. and Lekina, A. (2011). Kernel estimators of extreme level curves, Test, 20, Gardes, L. and Girard, S. (2008). A moving window approach for nonparametric estimation of the conditional tail index, Journal of Multivariate Analysis, 99, Hill, B.M. (1975). A simple general approach to inference about the tail of a distribution, The Annals of Statistics, 3, Rockafellar, R.T. and Uryasev, S. (2000). Optimization of conditional value-at-risk. Journal of Risk, 2, Valdez, E.A. (2005). Tail conditional variance for elliptically contoured distributions, Belgian Actuarial Bulletin, 5, Weissman, I. (1978). Estimation of parameters and large quantiles based on the k largest observations, Journal of the American Statistical Association, 73, / 41
Nonparametric estimation of tail risk measures from heavy-tailed distributions
Nonparametric estimation of tail risk measures from heavy-tailed distributions Jonthan El Methni, Laurent Gardes & Stéphane Girard 1 Tail risk measures Let Y R be a real random loss variable. The Value-at-Risk
More informationEstimation of risk measures for extreme pluviometrical measurements
Estimation of risk measures for extreme pluviometrical measurements by Jonathan EL METHNI in collaboration with Laurent GARDES & Stéphane GIRARD 26th Annual Conference of The International Environmetrics
More informationFrontier estimation based on extreme risk measures
Frontier estimation based on extreme risk measures by Jonathan EL METHNI in collaboration with Ste phane GIRARD & Laurent GARDES CMStatistics 2016 University of Seville December 2016 1 Risk measures 2
More informationEstimation of the functional Weibull-tail coefficient
1/ 29 Estimation of the functional Weibull-tail coefficient Stéphane Girard Inria Grenoble Rhône-Alpes & LJK, France http://mistis.inrialpes.fr/people/girard/ June 2016 joint work with Laurent Gardes,
More informationNon-parametric estimation of extreme risk measures from conditional heavy-tailed distributions
Non-parametric estimation of extreme risk measures from conditional heavy-tailed distributions Jonathan El Methni, Laurent Gardes, Stephane Girard To cite this version: Jonathan El Methni, Laurent Gardes,
More informationNONPARAMETRIC ESTIMATION OF THE CONDITIONAL TAIL INDEX
NONPARAMETRIC ESTIMATION OF THE CONDITIONAL TAIL INDE Laurent Gardes and Stéphane Girard INRIA Rhône-Alpes, Team Mistis, 655 avenue de l Europe, Montbonnot, 38334 Saint-Ismier Cedex, France. Stephane.Girard@inrialpes.fr
More informationNonparametric estimation of extreme risks from conditional heavy-tailed distributions
Nonparametric estimation of extreme risks from conditional heavy-tailed distributions Jonathan El Methni, Laurent Gardes, Stephane Girard To cite this version: Jonathan El Methni, Laurent Gardes, Stephane
More informationSpatial analysis of extreme rainfalls In the Cévennes-Vivarais region
Spatial analysis of extreme rainfalls In the Cévennes-Vivarais region using a nearest neighbor approach Caroline Bernard-Michel (1), Laurent Gardes (1), Stéphane Girard (1), Gilles Molinié (2) (1) INRIA
More informationExtreme L p quantiles as risk measures
1/ 27 Extreme L p quantiles as risk measures Stéphane GIRARD (Inria Grenoble Rhône-Alpes) joint work Abdelaati DAOUIA (Toulouse School of Economics), & Gilles STUPFLER (University of Nottingham) December
More informationEstimation de mesures de risques à partir des L p -quantiles
1/ 42 Estimation de mesures de risques à partir des L p -quantiles extrêmes Stéphane GIRARD (Inria Grenoble Rhône-Alpes) collaboration avec Abdelaati DAOUIA (Toulouse School of Economics), & Gilles STUPFLER
More informationContributions to extreme-value analysis
Contributions to extreme-value analysis Stéphane Girard INRIA Rhône-Alpes & LJK (team MISTIS). 655, avenue de l Europe, Montbonnot. 38334 Saint-Ismier Cedex, France Stephane.Girard@inria.fr February 6,
More informationExtreme Value Analysis and Spatial Extremes
Extreme Value Analysis and Department of Statistics Purdue University 11/07/2013 Outline Motivation 1 Motivation 2 Extreme Value Theorem and 3 Bayesian Hierarchical Models Copula Models Max-stable Models
More informationQuantile prediction of a random eld extending the gaussian setting
Quantile prediction of a random eld extending the gaussian setting 1 Joint work with : Véronique Maume-Deschamps 1 and Didier Rullière 2 1 Institut Camille Jordan Université Lyon 1 2 Laboratoire des Sciences
More informationOn the estimation of the second order parameter in extreme-value theory
On the estimation of the second order parameter in extreme-value theory by El Hadji DEME (1,3), Laurent Gardes (2) and Stéphane Girard (3) (1) LERSTAD, Université Gaston Berger, Saint-Louis, Sénégal. (2)
More informationMFM Practitioner Module: Quantitiative Risk Management. John Dodson. October 14, 2015
MFM Practitioner Module: Quantitiative Risk Management October 14, 2015 The n-block maxima 1 is a random variable defined as M n max (X 1,..., X n ) for i.i.d. random variables X i with distribution function
More informationThe high order moments method in endpoint estimation: an overview
1/ 33 The high order moments method in endpoint estimation: an overview Gilles STUPFLER (Aix Marseille Université) Joint work with Stéphane GIRARD (INRIA Rhône-Alpes) and Armelle GUILLOU (Université de
More informationEstimation of the second order parameter for heavy-tailed distributions
Estimation of the second order parameter for heavy-tailed distributions Stéphane Girard Inria Grenoble Rhône-Alpes, France joint work with El Hadji Deme (Université Gaston-Berger, Sénégal) and Laurent
More informationExtreme geometric quantiles
1/ 30 Extreme geometric quantiles Stéphane GIRARD (Inria Grenoble Rhône-Alpes) Joint work with Gilles STUPFLER (Aix Marseille Université) ERCIM, Pisa, Italy, December 2014 2/ 30 Outline Extreme multivariate
More informationAnalysis methods of heavy-tailed data
Institute of Control Sciences Russian Academy of Sciences, Moscow, Russia February, 13-18, 2006, Bamberg, Germany June, 19-23, 2006, Brest, France May, 14-19, 2007, Trondheim, Norway PhD course Chapter
More informationIntroduction to Algorithmic Trading Strategies Lecture 10
Introduction to Algorithmic Trading Strategies Lecture 10 Risk Management Haksun Li haksun.li@numericalmethod.com www.numericalmethod.com Outline Value at Risk (VaR) Extreme Value Theory (EVT) References
More informationBayesian Modelling of Extreme Rainfall Data
Bayesian Modelling of Extreme Rainfall Data Elizabeth Smith A thesis submitted for the degree of Doctor of Philosophy at the University of Newcastle upon Tyne September 2005 UNIVERSITY OF NEWCASTLE Bayesian
More informationStatistical inference on Lévy processes
Alberto Coca Cabrero University of Cambridge - CCA Supervisors: Dr. Richard Nickl and Professor L.C.G.Rogers Funded by Fundación Mutua Madrileña and EPSRC MASDOC/CCA student workshop 2013 26th March Outline
More informationTime Series and Forecasting Lecture 4 NonLinear Time Series
Time Series and Forecasting Lecture 4 NonLinear Time Series Bruce E. Hansen Summer School in Economics and Econometrics University of Crete July 23-27, 2012 Bruce Hansen (University of Wisconsin) Foundations
More informationAnalysis methods of heavy-tailed data
Institute of Control Sciences Russian Academy of Sciences, Moscow, Russia February, 13-18, 2006, Bamberg, Germany June, 19-23, 2006, Brest, France May, 14-19, 2007, Trondheim, Norway PhD course Chapter
More informationEstimation of the Bivariate and Marginal Distributions with Censored Data
Estimation of the Bivariate and Marginal Distributions with Censored Data Michael Akritas and Ingrid Van Keilegom Penn State University and Eindhoven University of Technology May 22, 2 Abstract Two new
More informationLifetime Dependence Modelling using a Generalized Multivariate Pareto Distribution
Lifetime Dependence Modelling using a Generalized Multivariate Pareto Distribution Daniel Alai Zinoviy Landsman Centre of Excellence in Population Ageing Research (CEPAR) School of Mathematics, Statistics
More informationA proposal of a bivariate Conditional Tail Expectation
A proposal of a bivariate Conditional Tail Expectation Elena Di Bernardino a joint works with Areski Cousin b, Thomas Laloë c, Véronique Maume-Deschamps d and Clémentine Prieur e a, b, d Université Lyon
More informationPitfalls in Using Weibull Tailed Distributions
Pitfalls in Using Weibull Tailed Distributions Alexandru V. Asimit, Deyuan Li & Liang Peng First version: 18 December 2009 Research Report No. 27, 2009, Probability and Statistics Group School of Mathematics,
More informationEstimation of the extreme value index and high quantiles under random censoring
Estimation of the extreme value index and high quantiles under random censoring Jan Beirlant () & Emmanuel Delafosse (2) & Armelle Guillou (2) () Katholiee Universiteit Leuven, Department of Mathematics,
More informationExperience Rating in General Insurance by Credibility Estimation
Experience Rating in General Insurance by Credibility Estimation Xian Zhou Department of Applied Finance and Actuarial Studies Macquarie University, Sydney, Australia Abstract This work presents a new
More informationPractical conditions on Markov chains for weak convergence of tail empirical processes
Practical conditions on Markov chains for weak convergence of tail empirical processes Olivier Wintenberger University of Copenhagen and Paris VI Joint work with Rafa l Kulik and Philippe Soulier Toronto,
More informationESTIMATING BIVARIATE TAIL
Elena DI BERNARDINO b joint work with Clémentine PRIEUR a and Véronique MAUME-DESCHAMPS b a LJK, Université Joseph Fourier, Grenoble 1 b Laboratoire SAF, ISFA, Université Lyon 1 Framework Goal: estimating
More informationQuantile-quantile plots and the method of peaksover-threshold
Problems in SF2980 2009-11-09 12 6 4 2 0 2 4 6 0.15 0.10 0.05 0.00 0.05 0.10 0.15 Figure 2: qqplot of log-returns (x-axis) against quantiles of a standard t-distribution with 4 degrees of freedom (y-axis).
More information41903: Introduction to Nonparametrics
41903: Notes 5 Introduction Nonparametrics fundamentally about fitting flexible models: want model that is flexible enough to accommodate important patterns but not so flexible it overspecializes to specific
More informationCopula Regression RAHUL A. PARSA DRAKE UNIVERSITY & STUART A. KLUGMAN SOCIETY OF ACTUARIES CASUALTY ACTUARIAL SOCIETY MAY 18,2011
Copula Regression RAHUL A. PARSA DRAKE UNIVERSITY & STUART A. KLUGMAN SOCIETY OF ACTUARIES CASUALTY ACTUARIAL SOCIETY MAY 18,2011 Outline Ordinary Least Squares (OLS) Regression Generalized Linear Models
More informationSGN Advanced Signal Processing: Lecture 8 Parameter estimation for AR and MA models. Model order selection
SG 21006 Advanced Signal Processing: Lecture 8 Parameter estimation for AR and MA models. Model order selection Ioan Tabus Department of Signal Processing Tampere University of Technology Finland 1 / 28
More informationEstimating a frontier function using a high-order moments method
1/ 16 Estimating a frontier function using a high-order moments method Gilles STUPFLER (University of Nottingham) Joint work with Stéphane GIRARD (INRIA Rhône-Alpes) and Armelle GUILLOU (Université de
More informationIntroduction to Empirical Processes and Semiparametric Inference Lecture 02: Overview Continued
Introduction to Empirical Processes and Semiparametric Inference Lecture 02: Overview Continued Michael R. Kosorok, Ph.D. Professor and Chair of Biostatistics Professor of Statistics and Operations Research
More informationBickel Rosenblatt test
University of Latvia 28.05.2011. A classical Let X 1,..., X n be i.i.d. random variables with a continuous probability density function f. Consider a simple hypothesis H 0 : f = f 0 with a significance
More informationMultivariate Distributions
IEOR E4602: Quantitative Risk Management Spring 2016 c 2016 by Martin Haugh Multivariate Distributions We will study multivariate distributions in these notes, focusing 1 in particular on multivariate
More informationInference on distributions and quantiles using a finite-sample Dirichlet process
Dirichlet IDEAL Theory/methods Simulations Inference on distributions and quantiles using a finite-sample Dirichlet process David M. Kaplan University of Missouri Matt Goldman UC San Diego Midwest Econometrics
More informationFinancial Econometrics and Volatility Models Extreme Value Theory
Financial Econometrics and Volatility Models Extreme Value Theory Eric Zivot May 3, 2010 1 Lecture Outline Modeling Maxima and Worst Cases The Generalized Extreme Value Distribution Modeling Extremes Over
More informationA Framework for Daily Spatio-Temporal Stochastic Weather Simulation
A Framework for Daily Spatio-Temporal Stochastic Weather Simulation, Rick Katz, Balaji Rajagopalan Geophysical Statistics Project Institute for Mathematics Applied to Geosciences National Center for Atmospheric
More informationMax stable Processes & Random Fields: Representations, Models, and Prediction
Max stable Processes & Random Fields: Representations, Models, and Prediction Stilian Stoev University of Michigan, Ann Arbor March 2, 2011 Based on joint works with Yizao Wang and Murad S. Taqqu. 1 Preliminaries
More informationAsymptotic Statistics-III. Changliang Zou
Asymptotic Statistics-III Changliang Zou The multivariate central limit theorem Theorem (Multivariate CLT for iid case) Let X i be iid random p-vectors with mean µ and and covariance matrix Σ. Then n (
More informationHeavy Tailed Time Series with Extremal Independence
Heavy Tailed Time Series with Extremal Independence Rafa l Kulik and Philippe Soulier Conference in honour of Prof. Herold Dehling Bochum January 16, 2015 Rafa l Kulik and Philippe Soulier Regular variation
More informationThe Convergence Rate for the Normal Approximation of Extreme Sums
The Convergence Rate for the Normal Approximation of Extreme Sums Yongcheng Qi University of Minnesota Duluth WCNA 2008, Orlando, July 2-9, 2008 This talk is based on a joint work with Professor Shihong
More informationIntroduction to Rare Event Simulation
Introduction to Rare Event Simulation Brown University: Summer School on Rare Event Simulation Jose Blanchet Columbia University. Department of Statistics, Department of IEOR. Blanchet (Columbia) 1 / 31
More informationNotes, March 4, 2013, R. Dudley Maximum likelihood estimation: actual or supposed
18.466 Notes, March 4, 2013, R. Dudley Maximum likelihood estimation: actual or supposed 1. MLEs in exponential families Let f(x,θ) for x X and θ Θ be a likelihood function, that is, for present purposes,
More informationCross-Validation with Confidence
Cross-Validation with Confidence Jing Lei Department of Statistics, Carnegie Mellon University WHOA-PSI Workshop, St Louis, 2017 Quotes from Day 1 and Day 2 Good model or pure model? Occam s razor We really
More informationBregman superquantiles. Estimation methods and applications
Bregman superquantiles. Estimation methods and applications Institut de mathématiques de Toulouse 2 juin 2014 Joint work with F. Gamboa, A. Garivier (IMT) and B. Iooss (EDF R&D). 1 Coherent measure of
More informationAdditive functionals of infinite-variance moving averages. Wei Biao Wu The University of Chicago TECHNICAL REPORT NO. 535
Additive functionals of infinite-variance moving averages Wei Biao Wu The University of Chicago TECHNICAL REPORT NO. 535 Departments of Statistics The University of Chicago Chicago, Illinois 60637 June
More informationCross-Validation with Confidence
Cross-Validation with Confidence Jing Lei Department of Statistics, Carnegie Mellon University UMN Statistics Seminar, Mar 30, 2017 Overview Parameter est. Model selection Point est. MLE, M-est.,... Cross-validation
More information2. Variance and Covariance: We will now derive some classic properties of variance and covariance. Assume real-valued random variables X and Y.
CS450 Final Review Problems Fall 08 Solutions or worked answers provided Problems -6 are based on the midterm review Identical problems are marked recap] Please consult previous recitations and textbook
More informationFall 2017 STAT 532 Homework Peter Hoff. 1. Let P be a probability measure on a collection of sets A.
1. Let P be a probability measure on a collection of sets A. (a) For each n N, let H n be a set in A such that H n H n+1. Show that P (H n ) monotonically converges to P ( k=1 H k) as n. (b) For each n
More informationGoodness-of-fit testing and Pareto-tail estimation
Goodness-of-fit testing and Pareto-tail estimation Yuri Goegebeur Department of Statistics, University of Southern Denmar, e-mail: yuri.goegebeur@stat.sdu.d Jan Beirlant University Center for Statistics,
More informationA Conditional Approach to Modeling Multivariate Extremes
A Approach to ing Multivariate Extremes By Heffernan & Tawn Department of Statistics Purdue University s April 30, 2014 Outline s s Multivariate Extremes s A central aim of multivariate extremes is trying
More informationBayesian Semiparametric GARCH Models
Bayesian Semiparametric GARCH Models Xibin (Bill) Zhang and Maxwell L. King Department of Econometrics and Business Statistics Faculty of Business and Economics xibin.zhang@monash.edu Quantitative Methods
More informationQuantile Regression for Extraordinarily Large Data
Quantile Regression for Extraordinarily Large Data Shih-Kang Chao Department of Statistics Purdue University November, 2016 A joint work with Stanislav Volgushev and Guang Cheng Quantile regression Two-step
More informationBayesian Semiparametric GARCH Models
Bayesian Semiparametric GARCH Models Xibin (Bill) Zhang and Maxwell L. King Department of Econometrics and Business Statistics Faculty of Business and Economics xibin.zhang@monash.edu Quantitative Methods
More informationMath 576: Quantitative Risk Management
Math 576: Quantitative Risk Management Haijun Li lih@math.wsu.edu Department of Mathematics Washington State University Week 11 Haijun Li Math 576: Quantitative Risk Management Week 11 1 / 21 Outline 1
More informationModeling Real Estate Data using Quantile Regression
Modeling Real Estate Data using Semiparametric Quantile Regression Department of Statistics University of Innsbruck September 9th, 2011 Overview 1 Application: 2 3 4 Hedonic regression data for house prices
More informationOverview of Extreme Value Theory. Dr. Sawsan Hilal space
Overview of Extreme Value Theory Dr. Sawsan Hilal space Maths Department - University of Bahrain space November 2010 Outline Part-1: Univariate Extremes Motivation Threshold Exceedances Part-2: Bivariate
More informationTail negative dependence and its applications for aggregate loss modeling
Tail negative dependence and its applications for aggregate loss modeling Lei Hua Division of Statistics Oct 20, 2014, ISU L. Hua (NIU) 1/35 1 Motivation 2 Tail order Elliptical copula Extreme value copula
More informationNonparametric Regression. Badr Missaoui
Badr Missaoui Outline Kernel and local polynomial regression. Penalized regression. We are given n pairs of observations (X 1, Y 1 ),...,(X n, Y n ) where Y i = r(x i ) + ε i, i = 1,..., n and r(x) = E(Y
More informationMultivariate Normal-Laplace Distribution and Processes
CHAPTER 4 Multivariate Normal-Laplace Distribution and Processes The normal-laplace distribution, which results from the convolution of independent normal and Laplace random variables is introduced by
More informationExponential tail inequalities for eigenvalues of random matrices
Exponential tail inequalities for eigenvalues of random matrices M. Ledoux Institut de Mathématiques de Toulouse, France exponential tail inequalities classical theme in probability and statistics quantify
More informationNonparametric Econometrics
Applied Microeconometrics with Stata Nonparametric Econometrics Spring Term 2011 1 / 37 Contents Introduction The histogram estimator The kernel density estimator Nonparametric regression estimators Semi-
More informationNonparametric Inference via Bootstrapping the Debiased Estimator
Nonparametric Inference via Bootstrapping the Debiased Estimator Yen-Chi Chen Department of Statistics, University of Washington ICSA-Canada Chapter Symposium 2017 1 / 21 Problem Setup Let X 1,, X n be
More informationIntroduction to Regression
Introduction to Regression p. 1/97 Introduction to Regression Chad Schafer cschafer@stat.cmu.edu Carnegie Mellon University Introduction to Regression p. 1/97 Acknowledgement Larry Wasserman, All of Nonparametric
More informationNonparametric regression with martingale increment errors
S. Gaïffas (LSTA - Paris 6) joint work with S. Delattre (LPMA - Paris 7) work in progress Motivations Some facts: Theoretical study of statistical algorithms requires stationary and ergodicity. Concentration
More informationSTA 2201/442 Assignment 2
STA 2201/442 Assignment 2 1. This is about how to simulate from a continuous univariate distribution. Let the random variable X have a continuous distribution with density f X (x) and cumulative distribution
More informationRobust Backtesting Tests for Value-at-Risk Models
Robust Backtesting Tests for Value-at-Risk Models Jose Olmo City University London (joint work with Juan Carlos Escanciano, Indiana University) Far East and South Asia Meeting of the Econometric Society
More informationQuantile Processes for Semi and Nonparametric Regression
Quantile Processes for Semi and Nonparametric Regression Shih-Kang Chao Department of Statistics Purdue University IMS-APRM 2016 A joint work with Stanislav Volgushev and Guang Cheng Quantile Response
More informationStable Process. 2. Multivariate Stable Distributions. July, 2006
Stable Process 2. Multivariate Stable Distributions July, 2006 1. Stable random vectors. 2. Characteristic functions. 3. Strictly stable and symmetric stable random vectors. 4. Sub-Gaussian random vectors.
More informationLecture 3. Probability - Part 2. Luigi Freda. ALCOR Lab DIAG University of Rome La Sapienza. October 19, 2016
Lecture 3 Probability - Part 2 Luigi Freda ALCOR Lab DIAG University of Rome La Sapienza October 19, 2016 Luigi Freda ( La Sapienza University) Lecture 3 October 19, 2016 1 / 46 Outline 1 Common Continuous
More informationChange point in trending regression
Charles University, Prague ROBUST 2010 joint work with A. Aue, L.Horváth, J. Picek (Charles University, Prague) 1 (Charles University, Prague) 1 2 model-formulation (Charles University, Prague) 1 2 model-formulation
More informationApproximating the Integrated Tail Distribution
Approximating the Integrated Tail Distribution Ants Kaasik & Kalev Pärna University of Tartu VIII Tartu Conference on Multivariate Statistics 26th of June, 2007 Motivating example Let W be a steady-state
More informationInterval Estimation for AR(1) and GARCH(1,1) Models
for AR(1) and GARCH(1,1) Models School of Mathematics Georgia Institute of Technology 2010 This talk is based on the following papers: Ngai Hang Chan, Deyuan Li and (2010). Toward a Unified of Autoregressions.
More informationPattern Recognition and Machine Learning. Bishop Chapter 2: Probability Distributions
Pattern Recognition and Machine Learning Chapter 2: Probability Distributions Cécile Amblard Alex Kläser Jakob Verbeek October 11, 27 Probability Distributions: General Density Estimation: given a finite
More informationGeneralized Multivariate Rank Type Test Statistics via Spatial U-Quantiles
Generalized Multivariate Rank Type Test Statistics via Spatial U-Quantiles Weihua Zhou 1 University of North Carolina at Charlotte and Robert Serfling 2 University of Texas at Dallas Final revision for
More informationCan we do statistical inference in a non-asymptotic way? 1
Can we do statistical inference in a non-asymptotic way? 1 Guang Cheng 2 Statistics@Purdue www.science.purdue.edu/bigdata/ ONR Review Meeting@Duke Oct 11, 2017 1 Acknowledge NSF, ONR and Simons Foundation.
More informationModelling Risk on Losses due to Water Spillage for Hydro Power Generation. A Verster and DJ de Waal
Modelling Risk on Losses due to Water Spillage for Hydro Power Generation. A Verster and DJ de Waal Department Mathematical Statistics and Actuarial Science University of the Free State Bloemfontein ABSTRACT
More informationDoes k-th Moment Exist?
Does k-th Moment Exist? Hitomi, K. 1 and Y. Nishiyama 2 1 Kyoto Institute of Technology, Japan 2 Institute of Economic Research, Kyoto University, Japan Email: hitomi@kit.ac.jp Keywords: Existence of moments,
More informationNonparametric Estimation of the Dependence Function for a Multivariate Extreme Value Distribution
Nonparametric Estimation of the Dependence Function for a Multivariate Extreme Value Distribution p. /2 Nonparametric Estimation of the Dependence Function for a Multivariate Extreme Value Distribution
More informationResearch Article Strong Convergence Bound of the Pareto Index Estimator under Right Censoring
Hindawi Publishing Corporation Journal of Inequalities and Applications Volume 200, Article ID 20956, 8 pages doi:0.55/200/20956 Research Article Strong Convergence Bound of the Pareto Index Estimator
More informationThreshold estimation in marginal modelling of spatially-dependent non-stationary extremes
Threshold estimation in marginal modelling of spatially-dependent non-stationary extremes Philip Jonathan Shell Technology Centre Thornton, Chester philip.jonathan@shell.com Paul Northrop University College
More information5 Operations on Multiple Random Variables
EE360 Random Signal analysis Chapter 5: Operations on Multiple Random Variables 5 Operations on Multiple Random Variables Expected value of a function of r.v. s Two r.v. s: ḡ = E[g(X, Y )] = g(x, y)f X,Y
More information12 - Nonparametric Density Estimation
ST 697 Fall 2017 1/49 12 - Nonparametric Density Estimation ST 697 Fall 2017 University of Alabama Density Review ST 697 Fall 2017 2/49 Continuous Random Variables ST 697 Fall 2017 3/49 1.0 0.8 F(x) 0.6
More informationConcentration Inequalities for Random Matrices
Concentration Inequalities for Random Matrices M. Ledoux Institut de Mathématiques de Toulouse, France exponential tail inequalities classical theme in probability and statistics quantify the asymptotic
More informationBeyond tail median and conditional tail expectation: extreme risk estimation using tail Lp optimisation
Beyond tail median and conditional tail expectation: extreme ris estimation using tail Lp optimisation Laurent Gardes, Stephane Girard, Gilles Stupfler To cite this version: Laurent Gardes, Stephane Girard,
More informationAn Introduction to Bayesian Linear Regression
An Introduction to Bayesian Linear Regression APPM 5720: Bayesian Computation Fall 2018 A SIMPLE LINEAR MODEL Suppose that we observe explanatory variables x 1, x 2,..., x n and dependent variables y 1,
More informationError distribution function for parametrically truncated and censored data
Error distribution function for parametrically truncated and censored data Géraldine LAURENT Jointly with Cédric HEUCHENNE QuantOM, HEC-ULg Management School - University of Liège Friday, 14 September
More informationErnesto Mordecki 1. Lecture III. PASI - Guanajuato - June 2010
Optimal stopping for Hunt and Lévy processes Ernesto Mordecki 1 Lecture III. PASI - Guanajuato - June 2010 1Joint work with Paavo Salminen (Åbo, Finland) 1 Plan of the talk 1. Motivation: from Finance
More informationLecture 1: August 28
36-705: Intermediate Statistics Fall 2017 Lecturer: Siva Balakrishnan Lecture 1: August 28 Our broad goal for the first few lectures is to try to understand the behaviour of sums of independent random
More informationInvestigation of an Automated Approach to Threshold Selection for Generalized Pareto
Investigation of an Automated Approach to Threshold Selection for Generalized Pareto Kate R. Saunders Supervisors: Peter Taylor & David Karoly University of Melbourne April 8, 2015 Outline 1 Extreme Value
More informationLARGE SAMPLE BEHAVIOR OF SOME WELL-KNOWN ROBUST ESTIMATORS UNDER LONG-RANGE DEPENDENCE
LARGE SAMPLE BEHAVIOR OF SOME WELL-KNOWN ROBUST ESTIMATORS UNDER LONG-RANGE DEPENDENCE C. LÉVY-LEDUC, H. BOISTARD, E. MOULINES, M. S. TAQQU, AND V. A. REISEN Abstract. The paper concerns robust location
More informationFinal Examination Statistics 200C. T. Ferguson June 11, 2009
Final Examination Statistics 00C T. Ferguson June, 009. (a) Define: X n converges in probability to X. (b) Define: X m converges in quadratic mean to X. (c) Show that if X n converges in quadratic mean
More informationGaussian Processes. Le Song. Machine Learning II: Advanced Topics CSE 8803ML, Spring 2012
Gaussian Processes Le Song Machine Learning II: Advanced Topics CSE 8803ML, Spring 01 Pictorial view of embedding distribution Transform the entire distribution to expected features Feature space Feature
More informationIntroduction to Regression
Introduction to Regression Chad M. Schafer May 20, 2015 Outline General Concepts of Regression, Bias-Variance Tradeoff Linear Regression Nonparametric Procedures Cross Validation Local Polynomial Regression
More information