Modelling Multivariate Peaks-over-Thresholds using Generalized Pareto Distributions
|
|
- Norman Price
- 5 years ago
- Views:
Transcription
1 Modelling Multivariate Peaks-over-Thresholds using Generalized Pareto Distributions Anna Kiriliouk 1 Holger Rootzén 2 Johan Segers 1 Jennifer L. Wadsworth 3 1 Université catholique de Louvain (BE) 2 Chalmers University of Technology (S) 3 Lancaster University (UK) Risk, Extremes and Contagion Université de Nanterre, May 25 26, 2016 Kiriliouk/Rootzén/Segers/Wadsworth Multivariate Generalized Pareto Distributions Nanterre, 26 May / 37
2 Multivariate Peaks-over-Thresholds ξ exceeds a multivariate threshold u: ξ u j = 1,..., d : ξ j > u j ξ 2 u 2 Distribution of excess ξ u given ξ u? ξ 1 u 1 Other exceedances (Thibaud and Opitz, 2013; Dombry and Ribatet, 2015): j = 1,..., d : ξ j > u j ξ ξ d > u u d... Kiriliouk/Rootzén/Segers/Wadsworth Multivariate Generalized Pareto Distributions Nanterre, 26 May / 37
3 Multivariate Generalized Pareto Distributions Maxima and exceedances Parametrization Representations and densities Stability Inference Kiriliouk/Rootzén/Segers/Wadsworth Multivariate Generalized Pareto Distributions Nanterre, 26 May / 37
4 Generalized Extreme-Value distributions Random vector ξ in R d, cdf F ξ (x) = P(ξ x). Assumption: Max-domain of attraction There exist scaling a n > 0 and location b n R d sequences such that F n ξ (a nx + b n ) G(x), n F n is cdf of vector of componentwise maxima of iid sample from F. Limit G is necessarily Generalized Extreme-Value = Max-Stable Kiriliouk/Rootzén/Segers/Wadsworth Multivariate Generalized Pareto Distributions Nanterre, 26 May / 37
5 Multivariate Generalized Pareto distribution If F n ξ (a nx + b n ) G(x) as n, then, provided 0 < G(0) < 1 (w.l.o.g.), L ( an 1 ) d (ξ b n ) η ξ b n Multivariate Generalized Pareto H where η j [, ) is the lower endpoint of G j. Cumulative distribution function H = GP(G): for x such that G(x) > 0, H(x) = log G(x 0) ( log G(x)) log G(0) Beirlant et al. (2004) Rootzén and Tajvidi (2006) Falk et al. (2010) Ferreira and de Haan (2014) Kiriliouk/Rootzén/Segers/Wadsworth Multivariate Generalized Pareto Distributions Nanterre, 26 May / 37
6 Independent maxima: single-component exceedances Generalized Extreme-Value G: ( G(x 1, x 2 ) = exp 1 x ) x x 1 > 1, x 2 > 1 Generalized Pareto H: law of { ( 1, Y) wp 1/2, (X 1, X 2 ) = (Y, 1) wp 1/2 X with P(Y t) = y for y 0 (univariate Generalized Pareto). X 1 Kiriliouk/Rootzén/Segers/Wadsworth Multivariate Generalized Pareto Distributions Nanterre, 26 May / 37
7 Completely dependent maxima and exceedances Generalized Extreme-Value G: ( ) 1 G(x 1, x 2 ) = exp min(x 1, x 2 ) + 1 x 1 > 1, x 2 > 1 Generalized Pareto H: law of (Y, Y) X with P(Y y) = y for y 0 (univariate Generalized Pareto). X 1 Kiriliouk/Rootzén/Segers/Wadsworth Multivariate Generalized Pareto Distributions Nanterre, 26 May / 37
8 1 Example: logistic Generalized Extreme-Value G(x 1, x 2 ): exp[ { 2 j=1 (1 + x j) 1/θ } θ ] x 1 > 1, x 2 > 1 Generalized Pareto H: pdf h(x 1, x 2 ): 2 θ 2 x 1 x 2 { 2 j=1 (1+x j) 1/θ } θ Support: x 1 > 1 x 2 > 1 x 1 x 2 > θ = Kiriliouk/Rootzén/Segers/Wadsworth Multivariate Generalized Pareto Distributions Nanterre, 26 May / 37
9 Multivariate Generalized Pareto Distributions Maxima and exceedances Parametrization Representations and densities Stability Inference Kiriliouk/Rootzén/Segers/Wadsworth Multivariate Generalized Pareto Distributions Nanterre, 26 May / 37
10 Generalized Extreme-Value distributions Margins: shape γ j R, location µ j R, scale α j (0, ) exp[ {1 + γ j (x j µ j )/α j } 1/γ j ] if γ j 0, G j (x j ) = exp[ exp{ (x j µ j )/α j }] if γ j = 0 Stable tail dependence function (stdf) l : [0, ) d [0, ) G(x) = exp[ l{ log G 1 (x 1 ),..., log G d (x d )}] Notation: G = GEV(γ, µ, α, l) Kiriliouk/Rootzén/Segers/Wadsworth Multivariate Generalized Pareto Distributions Nanterre, 26 May / 37
11 Identifiability issue Max-stability: G = GEV(γ, µ, α, l) = G t = GEV(γ, µ(t), α(t), l), 0 < t < with µ(t) =... and α(t) =... (exercise). However: GP(G t ) = GP(G). Therefore: not all GEV parameters are identifiable from H = GP(G). Solution? Kiriliouk/Rootzén/Segers/Wadsworth Multivariate Generalized Pareto Distributions Nanterre, 26 May / 37
12 Cumulative distribution function Let G = GEV(γ, µ, α, l). Suppose σ = α γµ (0, ) Then 0 < G(0) < 1. Let H = GP(G). For x R such that σ + γx > 0, ( H(x) = l π (1 + γ(x 0)/σ) 1/γ) l (π (1 + γx/σ) 1/γ) where, in case also x 0, H j (0) = π j H j (x j ) H j (0) = (1 + γ jx j /σ j ) 1/γ j H(x) = l ( H 1 (x 1 ),..., H d (x d ) ) Parametrization: H = GP(γ, σ, π, l). Constraint: l(π) = 1. Kiriliouk/Rootzén/Segers/Wadsworth Multivariate Generalized Pareto Distributions Nanterre, 26 May / 37
13 Multivariate Generalized Pareto Distributions Maxima and exceedances Parametrization Representations and densities Stability Inference Kiriliouk/Rootzén/Segers/Wadsworth Multivariate Generalized Pareto Distributions Nanterre, 26 May / 37
14 Standardization We have X GP(γ, σ, π, l) if and only if we have and Z GP(0, 1, π, l). X = σ eγz 1 γ The support of Z is contained in [, ) \ [, 0] and P[Z z] = l ( πe (z 0)) l ( πe z) How to construct Z? Kiriliouk/Rootzén/Segers/Wadsworth Multivariate Generalized Pareto Distributions Nanterre, 26 May / 37
15 Spectral representation A random vector S in [, 0] d is a spectral random vector if: P[max(S 1,..., S d ) = 0] = 1 P[S j > ] > 0, j = 1,..., d Let E be a unit exponential random variable independent of S. Then where Z := S + E GP(0, 1, π, l) π j = E[e S j ] [ l(y) = E max j=1,...,d { e S }] j y j π j H(z) = 1 E[1 e max(s z) ] Conversely: given π and l with l(π) = 1, there exists a unique (in law) spectral random vector S such that the above holds. Kiriliouk/Rootzén/Segers/Wadsworth Multivariate Generalized Pareto Distributions Nanterre, 26 May / 37
16 Constructing spectral random vectors (1) Let T be a random vector in [, ) d such that Then 0 < E[e T j ] <, j = 1,..., d P[max(T 1,..., T d ) > ] = 1 S = T max(t) is a spectral random vector. The associated GP(0, 1, π, l) = GP S (0, 1, L(S)) distribution is determined by π j = E[e Tj max(t) ] [ { l(y) = E max j=1,...,d y j }] e T j max(t) E[e Tj max(t) ] H(z) = 1 E[1 e max(t z) max(t) ] Kiriliouk/Rootzén/Segers/Wadsworth Multivariate Generalized Pareto Distributions Nanterre, 26 May / 37
17 Density (1) Suppose T has support included in R d and is absolutely continuous with Lebesgue density f T. Then the GP T (0, 1, L(T)) distribution has Lebesgue density h given by h(z) = 1(z 0) 1 e max(z) 0 f T (z + log t) t 1 dt. Kiriliouk/Rootzén/Segers/Wadsworth Multivariate Generalized Pareto Distributions Nanterre, 26 May / 37
18 Example: max-normalized log-gaussian generators (1) Pdf of H = GP T (0, 1, N d (0, Σ)): h(z) = (2π)(1 d)/2 Σ 1/2 (1 T Σ 1 1) 1/2 e max(z) exp { 12 (Σ zt 1 Σ 1 11 T Σ 1 ) } 1 T z Σ 1 1 z R d \ (, 0] d Kiriliouk/Rootzén/Segers/Wadsworth Multivariate Generalized Pareto Distributions Nanterre, 26 May / 37
19 Example: max-normalized log-gaussian generators (2) pdf of (Z 1, Z 2 ) GP T with γ 1 = γ 2 = 0 σ 1 = σ 2 = 0 (T 1, T 2 ) N 2 (0, ( )) 1 ρ ρ ρ =.5 Kiriliouk/Rootzén/Segers/Wadsworth Multivariate Generalized Pareto Distributions Nanterre, 26 May / 37
20 Example: max-normalized log-gaussian generators (2) pdf of (Z 1, Z 2 ) GP T with γ 1 = γ 2 = 0 σ 1 = σ 2 = 0 (T 1, T 2 ) N 2 (0, ( )) 1 ρ ρ ρ = 0 Kiriliouk/Rootzén/Segers/Wadsworth Multivariate Generalized Pareto Distributions Nanterre, 26 May / 37
21 Example: max-normalized log-gaussian generators (2) pdf of (Z 1, Z 2 ) GP T with γ 1 = γ 2 = 0 σ 1 = σ 2 = 0 (T 1, T 2 ) N 2 (0, ( )) 1 ρ ρ ρ =.5 Kiriliouk/Rootzén/Segers/Wadsworth Multivariate Generalized Pareto Distributions Nanterre, 26 May / 37
22 Constructing spectral random vectors (2) Let U be a random vector in [, ) d such that 0 < E[e U j ] < for all j. Let S be the spectral random vector defined in distribution by P[S ] = E [ 1{U max(u) } e max(u)] E[e max(u) ] The associated GP(0, 1, π, l) = GP S (0, 1, L(S)) distribution is given by π j = E[eU j ] E[e max(u) ] [ l(y) = E max j=1,...,d { y j e U }] j E[e U j] H(z) = 1 E[emax(U) e max(u z) ] E[e max(u) ] Kiriliouk/Rootzén/Segers/Wadsworth Multivariate Generalized Pareto Distributions Nanterre, 26 May / 37
23 Density (2) Suppose U has support included in R d and is absolutely continuous with Lebesgue density f T. Then the GP U (0, 1, L(U)) distribution has Lebesgue density h given by h(z) = 1(z 0) 1 E[e max(u) ] 0 f U (z + log t) dt. Kiriliouk/Rootzén/Segers/Wadsworth Multivariate Generalized Pareto Distributions Nanterre, 26 May / 37
24 Example: independent Fréchet generators e U j independent Fréchet shape α > 0 scales λ j > 0 GP U (0, 1, L(U)) density: h(z) e α d j=1 z j ( d j=1 ( ez j λ j ) α ) d 1/α z R d \ (, 0] d explicit proportionality constant margins dependence γ 1 = γ 2 = 0 α = 2 σ 1 = σ 2 = 1 λ 1 = 2, λ 2 = 1 Kiriliouk/Rootzén/Segers/Wadsworth Multivariate Generalized Pareto Distributions Nanterre, 26 May / 37
25 Example: independent Beta generators e U j independent Beta shapes (α j, 1) scales λ j GP U (0, 1, L(U)) density: h(z) z R d \ (, 0] d d j=1 ez j/α j max(λe z ) d j=1 1/α j explicit proportionality constant margins dependence γ 1 = γ 2 = 0 α 1 = 2, α 2 = 3 σ 1 = σ 2 = 1 λ 1 = 2, λ 2 = 1 Kiriliouk/Rootzén/Segers/Wadsworth Multivariate Generalized Pareto Distributions Nanterre, 26 May / 37
26 Point process representation iid U, U 1, U 2,... in [, ) d such that 0 < E[e U j ] < for all j unit-rate Poisson process 0 < R 1 < R 2 <... (U i ) i and (R i ) i are independent Point process consisting of points ξ i (i = 1, 2,...) where ξ i = U i log R i The law, G, of max i ξ i is GEV (de Haan, 1984): log G(z) = E[e max(u z) ] Associated GP: H = GP(G) = GP U (0, 1, L(U)) Kiriliouk/Rootzén/Segers/Wadsworth Multivariate Generalized Pareto Distributions Nanterre, 26 May / 37
27 Multivariate Generalized Pareto Distributions Maxima and exceedances Parametrization Representations and densities Stability Inference Kiriliouk/Rootzén/Segers/Wadsworth Multivariate Generalized Pareto Distributions Nanterre, 26 May / 37
28 Lower-dimensional margins Let X be GP and let J {1,..., d}. The law of X J = (X j ) j J is not GP. The law of X J given that max j J X j > 0 is GP: X GP(γ, σ, π, l) = L(X J X J 0) = GP (γ J, σ J, (P[X j > u j X J 0]) j J, l J ) X GP U (γ, σ, U) = L(X J X J 0) = GP U (γ J, σ J, L(U J )) If J = {j}, we find that L(X j X j > 0) = GP(γ j, σ j ). Kiriliouk/Rootzén/Segers/Wadsworth Multivariate Generalized Pareto Distributions Nanterre, 26 May / 37
29 Threshold stability Univariate GP distributions are threshold stable. What if multivariate? Let X GP(γ, σ, π, l). Let u [0, ) d be such that P(X j > u j ) > 0 for all j. Then L(X u X u) = GP (γ, σ + γu, (P[X j > u j X u]) j, l) Change of marginal parameters as in univariate case Change from π j = P[X j > 0] to P[X j > u j X u] Same stdf l Kiriliouk/Rootzén/Segers/Wadsworth Multivariate Generalized Pareto Distributions Nanterre, 26 May / 37
30 Linear combinations Let X GP S (γ, σ, π, l) be such that γ 1 =... = γ d =: γ. Let a [0, ) d and write a X = d j=1 a jx j. If P[a X > 0] > 0, then L(a X a X > 0) = GP(γ, a σ). Kiriliouk/Rootzén/Segers/Wadsworth Multivariate Generalized Pareto Distributions Nanterre, 26 May / 37
31 Linear transformations Let X GP S (γ, σ, L(S)) be such that γ 1 =... = γ d =: γ. Let A = (a i,j ) i,j [0, ) m d be such that P[A ix > 0] > 0 for all i. For x R m such that A iσ + γx i > 0 for all i = 1,..., m, [ ] P[AX x] = E 1 max {(1 + γx i/a iσ) 1/γ e U i } i=1,...,m where U = (U 1,..., U m ) is given by γ 1 log ( d j=1 U i = p i,j e γs ) j if γ 0, d j=1 p i,j S j if γ = 0, where p i,j = a i,j σ j /A iσ j. As a consequence, L(AX AX 0) = GP U (γ, Aσ, L(U)). Kiriliouk/Rootzén/Segers/Wadsworth Multivariate Generalized Pareto Distributions Nanterre, 26 May / 37
32 Multivariate Generalized Pareto Distributions Maxima and exceedances Parametrization Representations and densities Stability Inference Kiriliouk/Rootzén/Segers/Wadsworth Multivariate Generalized Pareto Distributions Nanterre, 26 May / 37
33 Inference problem Data ξ 1,..., ξ n. Threshold u. Threshold exceedances: {ξ i ξ i u}. Model: L(ξ ξ u) = u + σ eγz 1 γ GP(0, 1, π, l) Z GP T (0, 1, L(T)) GP U (0, 1, L(U)) Parametric model for l or L(T) or L(U), parameter α = Parametric model for f Z = Likelihood inference on θ = (γ, σ, α)? Nonparametric inference: empirical stdf ˆl n? Kiriliouk/Rootzén/Segers/Wadsworth Multivariate Generalized Pareto Distributions Nanterre, 26 May / 37
34 Data points with some components far beneath the threshold Exceedance x u (on original scale): excess x u 0. Model by Generalized Pareto? Ill-justified if x j u j for some j. Censor components that are too low. X 1 X 2 u 1 u 2 Kiriliouk/Rootzén/Segers/Wadsworth Multivariate Generalized Pareto Distributions Nanterre, 26 May / 37
35 Censored likelihood 1. Choose second, lower threshold v < u. 2. If x j v j, replace x j by v j : censoring from below. 3. Likelihood contribution of a point x u but x v: h(w C u C, x D u D ) dw C (,v C ] C = {j = 1,..., d : x j v j } censored variables; D = {j = 1,..., d : x j > v j } uncensored variables. Comparative study of likelihood-based estimators: Huser et al. (2014). Kiriliouk/Rootzén/Segers/Wadsworth Multivariate Generalized Pareto Distributions Nanterre, 26 May / 37
36 Censored likelihood: Proof of concept X GP U (γ, σ, log Beta) α 1 α 2 α 3 λ 1 λ 2 σ 1 σ 2 γ 1 γ 2 Boxplots of normalized parameter estimates, 30 repetitions Censoring at v = u = 0 Kiriliouk/Rootzén/Segers/Wadsworth Multivariate Generalized Pareto Distributions Nanterre, 26 May / 37
37 Overall summary Modelling excess X = ξ u conditionally on ξ u: = Multivariate Generalized Pareto distribution H Q How does H look like? A X = σ eγz 1 γ with Z GP(0, 1, π, l) Q How to construct parametric models for it? A via spectral representations Q How to fit such models? A Censored likelihood Kiriliouk/Rootzén/Segers/Wadsworth Multivariate Generalized Pareto Distributions Nanterre, 26 May / 37
38 Outlook Computational challenges Likelihood: proportionality constant E[e max(u) ]? Simulation: change-of-measure f U (u) e max(u) f U (u)? Likelihood optimisation in case many parameters?... Statistical challenges Graceful blending with models for the bulk of the distribution? Covariates? Threshold choice? Model construction: parsimony vs flexibility?... Thank you! Kiriliouk/Rootzén/Segers/Wadsworth Multivariate Generalized Pareto Distributions Nanterre, 26 May / 37
39 Bibliography Beirlant, J., Y. Goegebeur, J. Segers, and J. Teugels (2004). Statistics of Extremes: Theory and Applications. John Wiley & Sons. de Haan, L. (1984). A spectral representation for max-stable processes. The Annals of Probability 12(4), Dombry, C. and M. Ribatet (2015). Functional regular variations, pareto processes and peaks over threshold. Statistics and Its Interface 8, Falk, M., J. Hüsler, and R.-D. Reiss (2010). Laws of small numbers: extremes and rare events. Springer Science & Business Media. Ferreira, A. and L. de Haan (2014). The generalized Pareto process; with a view towards application and simulation. Bernoulli 20(4), Huser, R., A. C. Davison, and M. G. Genton (2014). A comparative study of likelihood estimators for multivariate extremes. Extremes (arxiv: ). Rootzén, H. and N. Tajvidi (2006). Multivariate generalized Pareto distributions. Bernoulli 12(5), Thibaud, E. and T. Opitz (2013). Efficient inference and simulation for elliptical Pareto processes. arxiv preprint arxiv: Kiriliouk/Rootzén/Segers/Wadsworth Multivariate Generalized Pareto Distributions Nanterre, 26 May / 37
Extreme value statistics: from one dimension to many. Lecture 1: one dimension Lecture 2: many dimensions
Extreme value statistics: from one dimension to many Lecture 1: one dimension Lecture 2: many dimensions The challenge for extreme value statistics right now: to go from 1 or 2 dimensions to 50 or more
More informationThe extremal elliptical model: Theoretical properties and statistical inference
1/25 The extremal elliptical model: Theoretical properties and statistical inference Thomas OPITZ Supervisors: Jean-Noel Bacro, Pierre Ribereau Institute of Mathematics and Modeling in Montpellier (I3M)
More informationPeaks over thresholds modelling with multivariate generalized Pareto distributions
Peaks over thresholds modelling with multivariate generalized Pareto distributions Anna Kiriliouk Johan Segers Université catholique de Louvain Institut de Statistique, Biostatistique et Sciences Actuarielles
More informationMultivariate Pareto distributions: properties and examples
Multivariate Pareto distributions: properties and examples Ana Ferreira 1, Laurens de Haan 2 1 ISA UTL and CEAUL, Portugal 2 Erasmus Univ Rotterdam and CEAUL EVT2013 Vimeiro, September 8 11 Univariate
More informationMultivariate peaks over thresholds models
Extremes DOI 1.17/s1687-17-294-4 Multivariate peaks over thresholds models Holger Rootzén 1 Johan Segers 2 Jennifer L. Wadsworth 3 Received: 18 October 216 / Revised: 2 April 217 / Accepted: 21 April 217
More informationEstimation of spatial max-stable models using threshold exceedances
Estimation of spatial max-stable models using threshold exceedances arxiv:1205.1107v1 [stat.ap] 5 May 2012 Jean-Noel Bacro I3M, Université Montpellier II and Carlo Gaetan DAIS, Università Ca Foscari -
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 informationModelling extreme-value dependence in high dimensions using threshold exceedances
Modelling extreme-value dependence in high dimensions using threshold exceedances Anna Kiriliouk A thesis submitted to the Université catholique de Louvain in partial fulfillment of the requirements for
More informationPREPRINT 2005:38. Multivariate Generalized Pareto Distributions HOLGER ROOTZÉN NADER TAJVIDI
PREPRINT 2005:38 Multivariate Generalized Pareto Distributions HOLGER ROOTZÉN NADER TAJVIDI Department of Mathematical Sciences Division of Mathematical Statistics CHALMERS UNIVERSITY OF TECHNOLOGY GÖTEBORG
More informationStatistics of Extremes
Statistics of Extremes Anthony Davison c 211 http://stat.epfl.ch Multivariate Extremes 19 Componentwise maxima.................................................. 194 Standardization........................................................
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 informationMultivariate generalized Pareto distributions
Multivariate generalized Pareto distributions Holger Rootzén and Nader Tajvidi Abstract Statistical inference for extremes has been a subject of intensive research during the past couple of decades. One
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 informationOn the occurrence times of componentwise maxima and bias in likelihood inference for multivariate max-stable distributions
On the occurrence times of componentwise maxima and bias in likelihood inference for multivariate max-stable distributions J. L. Wadsworth Department of Mathematics and Statistics, Fylde College, Lancaster
More informationModels and estimation.
Bivariate generalized Pareto distribution practice: Models and estimation. Eötvös Loránd University, Budapest, Hungary 7 June 2011, ASMDA Conference, Rome, Italy Problem How can we properly estimate the
More informationAssessing Dependence in Extreme Values
02/09/2016 1 Motivation Motivation 2 Comparison 3 Asymptotic Independence Component-wise Maxima Measures Estimation Limitations 4 Idea Results Motivation Given historical flood levels, how high should
More informationAbstract: In this short note, I comment on the research of Pisarenko et al. (2014) regarding the
Comment on Pisarenko et al. Characterization of the Tail of the Distribution of Earthquake Magnitudes by Combining the GEV and GPD Descriptions of Extreme Value Theory Mathias Raschke Institution: freelancer
More informationBivariate generalized Pareto distribution
Bivariate generalized Pareto distribution in practice Eötvös Loránd University, Budapest, Hungary Minisymposium on Uncertainty Modelling 27 September 2011, CSASC 2011, Krems, Austria Outline Short summary
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 informationExtreme Value Theory and Applications
Extreme Value Theory and Deauville - 04/10/2013 Extreme Value Theory and Introduction Asymptotic behavior of the Sum Extreme (from Latin exter, exterus, being on the outside) : Exceeding the ordinary,
More informationMultivariate generalized Pareto distributions
Bernoulli 12(5), 2006, 917 930 Multivariate generalized Pareto distributions HOLGER ROOTZÉN 1 and NADER TAJVIDI 2 1 Chalmers University of Technology, S-412 96 Göteborg, Sweden. E-mail rootzen@math.chalmers.se
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 informationBayesian inference for multivariate extreme value distributions
Bayesian inference for multivariate extreme value distributions Sebastian Engelke Clément Dombry, Marco Oesting Toronto, Fields Institute, May 4th, 2016 Main motivation For a parametric model Z F θ of
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 informationarxiv: v2 [stat.me] 25 Sep 2012
Estimation of Hüsler-Reiss distributions and Brown-Resnick processes arxiv:107.6886v [stat.me] 5 Sep 01 Sebastian Engelke Institut für Mathematische Stochastik, Georg-August-Universität Göttingen, Goldschmidtstr.
More informationVariable inspection plans for continuous populations with unknown short tail distributions
Variable inspection plans for continuous populations with unknown short tail distributions Wolfgang Kössler Abstract The ordinary variable inspection plans are sensitive to deviations from the normality
More informationMAS223 Statistical Inference and Modelling Exercises
MAS223 Statistical Inference and Modelling Exercises The exercises are grouped into sections, corresponding to chapters of the lecture notes Within each section exercises are divided into warm-up questions,
More informationNonlinear Time Series Modeling
Nonlinear Time Series Modeling Part II: Time Series Models in Finance Richard A. Davis Colorado State University (http://www.stat.colostate.edu/~rdavis/lectures) MaPhySto Workshop Copenhagen September
More informationarxiv: v4 [math.pr] 7 Feb 2012
THE MULTIVARIATE PIECING-TOGETHER APPROACH REVISITED STEFAN AULBACH, MICHAEL FALK, AND MARTIN HOFMANN arxiv:1108.0920v4 math.pr] 7 Feb 2012 Abstract. The univariate Piecing-Together approach (PT) fits
More informationThe Generalized Pareto process; with application
The Generalized Pareto process; with application Ana Ferreira ISA, Universidade Técnica de Lisboa and CEAUL Tapada da Ajuda 1349-017 Lisboa, Portugal Laurens de Haan Erasmus University Rotterdam and CEAUL
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 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 informationBayesian nonparametrics for multivariate extremes including censored data. EVT 2013, Vimeiro. Anne Sabourin. September 10, 2013
Bayesian nonparametrics for multivariate extremes including censored data Anne Sabourin PhD advisors: Anne-Laure Fougères (Lyon 1), Philippe Naveau (LSCE, Saclay). Joint work with Benjamin Renard, IRSTEA,
More informationSharp statistical tools Statistics for extremes
Sharp statistical tools Statistics for extremes Georg Lindgren Lund University October 18, 2012 SARMA Background Motivation We want to predict outside the range of observations Sums, averages and proportions
More informationBridging Asymptotic Independence and Dependence in Spatial Extremes Using Gaussian Scale Mixtures
Bridging Asymptotic Independence and Dependence in Spatial Extremes Using Gaussian Scale Mixtures Raphaël Huser 1, Thomas Opitz 2 and Emeric Thibaud 3 May 24, 217 Abstract Gaussian scale mixtures are constructed
More informationA Bayesian Spatial Model for Exceedances Over a Threshold
A Bayesian Spatial odel for Exceedances Over a Threshold Fernando Ferraz do Nascimento and Bruno Sansó June 2, 2017 Abstract Extreme value theory focuses on the study of rare events and uses asymptotic
More informationBayesian Model Averaging for Multivariate Extreme Values
Bayesian Model Averaging for Multivariate Extreme Values Philippe Naveau naveau@lsce.ipsl.fr Laboratoire des Sciences du Climat et l Environnement (LSCE) Gif-sur-Yvette, France joint work with A. Sabourin
More informationMULTIDIMENSIONAL COVARIATE EFFECTS IN SPATIAL AND JOINT EXTREMES
MULTIDIMENSIONAL COVARIATE EFFECTS IN SPATIAL AND JOINT EXTREMES Philip Jonathan, Kevin Ewans, David Randell, Yanyun Wu philip.jonathan@shell.com www.lancs.ac.uk/ jonathan Wave Hindcasting & Forecasting
More informationChapter 5 continued. Chapter 5 sections
Chapter 5 sections Discrete univariate distributions: 5.2 Bernoulli and Binomial distributions Just skim 5.3 Hypergeometric distributions 5.4 Poisson distributions Just skim 5.5 Negative Binomial distributions
More informationBivariate generalized Pareto distribution in practice: models and estimation
Bivariate generalized Pareto distribution in practice: models and estimation November 9, 2011 Pál Rakonczai Department of Probability Theory and Statistics, Eötvös University, Hungary e-mail: paulo@cs.elte.hu
More informationChange Point Analysis of Extreme Values
Change Point Analysis of Extreme Values TIES 2008 p. 1/? Change Point Analysis of Extreme Values Goedele Dierckx Economische Hogeschool Sint Aloysius, Brussels, Belgium e-mail: goedele.dierckx@hubrussel.be
More informationSimulation of Max Stable Processes
Simulation of Max Stable Processes Whitney Huang Department of Statistics Purdue University November 6, 2013 1 / 30 Outline 1 Max-Stable Processes 2 Unconditional Simulations 3 Conditional Simulations
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 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 informationBrief Review of Probability
Maura Department of Economics and Finance Università Tor Vergata Outline 1 Distribution Functions Quantiles and Modes of a Distribution 2 Example 3 Example 4 Distributions Outline Distribution Functions
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 informationModeling Extremal Events Is Not Easy: Why the Extreme Value Theorem Cannot Be As General As the Central Limit Theorem
University of Texas at El Paso DigitalCommons@UTEP Departmental Technical Reports (CS) Department of Computer Science 5-2015 Modeling Extremal Events Is Not Easy: Why the Extreme Value Theorem Cannot Be
More informationarxiv: v2 [math.pr] 28 Jul 2011
ON EXTREME VALUE PROCESSES AND THE FUNCTIONAL D-NORM arxiv:1107.5136v2 [math.pr] 28 Jul 2011 STEFAN AULBACH, MICHAEL FALK AND MARTIN HOFMANN Abstract. We introduce some mathematical framework for functional
More informationHIERARCHICAL MODELS IN EXTREME VALUE THEORY
HIERARCHICAL MODELS IN EXTREME VALUE THEORY Richard L. Smith Department of Statistics and Operations Research, University of North Carolina, Chapel Hill and Statistical and Applied Mathematical Sciences
More informationMULTIVARIATE EXTREMES AND RISK
MULTIVARIATE EXTREMES AND RISK Richard L. Smith Department of Statistics and Operations Research University of North Carolina Chapel Hill, NC 27599-3260 rls@email.unc.edu Interface 2008 RISK: Reality Durham,
More informationMax-stable Processes for Threshold Exceedances in Spatial Extremes
Max-stable Processes for Threshold Exceedances in Spatial Extremes Soyoung Jeon A dissertation submitted to the faculty of the University of North Carolina at Chapel Hill in partial fulfillment of the
More informationSTATISTICAL METHODS FOR RELATING TEMPERATURE EXTREMES TO LARGE-SCALE METEOROLOGICAL PATTERNS. Rick Katz
1 STATISTICAL METHODS FOR RELATING TEMPERATURE EXTREMES TO LARGE-SCALE METEOROLOGICAL PATTERNS Rick Katz Institute for Mathematics Applied to Geosciences National Center for Atmospheric Research Boulder,
More informationBias-corrected goodness-of-fit tests for Pareto-type behavior
Bias-corrected goodness-of-fit tests for Pareto-type behavior Yuri Goegebeur University of Southern Denmark, Department of Statistics JB Winsløws Ve 9B DK5000 Odense C, Denmark E-mail: yurigoegebeur@statsdudk
More informationSparse Representation of Multivariate Extremes with Applications to Anomaly Ranking
Sparse Representation of Multivariate Extremes with Applications to Anomaly Ranking Nicolas Goix Anne Sabourin Stéphan Clémençon Abstract Capturing the dependence structure of multivariate extreme events
More informationApproximate Bayesian computation for spatial extremes via open-faced sandwich adjustment
Approximate Bayesian computation for spatial extremes via open-faced sandwich adjustment Ben Shaby SAMSI August 3, 2010 Ben Shaby (SAMSI) OFS adjustment August 3, 2010 1 / 29 Outline 1 Introduction 2 Spatial
More informationBayesian Point Process Modeling for Extreme Value Analysis, with an Application to Systemic Risk Assessment in Correlated Financial Markets
Bayesian Point Process Modeling for Extreme Value Analysis, with an Application to Systemic Risk Assessment in Correlated Financial Markets Athanasios Kottas Department of Applied Mathematics and Statistics,
More informationSome conditional extremes of a Markov chain
Some conditional extremes of a Markov chain Seminar at Edinburgh University, November 2005 Adam Butler, Biomathematics & Statistics Scotland Jonathan Tawn, Lancaster University Acknowledgements: Janet
More information01 Probability Theory and Statistics Review
NAVARCH/EECS 568, ROB 530 - Winter 2018 01 Probability Theory and Statistics Review Maani Ghaffari January 08, 2018 Last Time: Bayes Filters Given: Stream of observations z 1:t and action data u 1:t Sensor/measurement
More informationAPPROXIMATING THE GENERALIZED BURR-GAMMA WITH A GENERALIZED PARETO-TYPE OF DISTRIBUTION A. VERSTER AND D.J. DE WAAL ABSTRACT
APPROXIMATING THE GENERALIZED BURR-GAMMA WITH A GENERALIZED PARETO-TYPE OF DISTRIBUTION A. VERSTER AND D.J. DE WAAL ABSTRACT In this paper the Generalized Burr-Gamma (GBG) distribution is considered to
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 informationarxiv: v1 [math.pr] 24 Aug 2018
On a class of norms generated by nonnegative integrable distributions arxiv:180808200v1 [mathpr] 24 Aug 2018 Michael Falk a and Gilles Stupfler b a Institute of Mathematics, University of Würzburg, Würzburg,
More informationOn the Estimation and Application of Max-Stable Processes
On the Estimation and Application of Max-Stable Processes Zhengjun Zhang Department of Statistics University of Wisconsin Madison, WI 53706, USA Co-author: Richard Smith EVA 2009, Fort Collins, CO Z. Zhang
More informationSampling Distributions
Sampling Distributions In statistics, a random sample is a collection of independent and identically distributed (iid) random variables, and a sampling distribution is the distribution of a function of
More informationAn M-estimator of Spatial Tail Dependence Einmahl, John; Kiriliouk, A.; Krajina, A.; Segers, J.J.J.
Tilburg University An M-estimator of Spatial Tail Dependence Einmahl, John; Kiriliou, A.; Krajina, A.; Segers, J.J.J. Document version: Early version, also nown as pre-print Publication date: 2014 Lin
More informationPhysically-Based Statistical Models of Extremes arising from Extratropical Cyclones
Lancaster University STOR603: PhD Proposal Physically-Based Statistical Models of Extremes arising from Extratropical Cyclones Author: Paul Sharkey Supervisors: Jonathan Tawn Jenny Wadsworth Simon Brown
More informationGeneralized Logistic Distribution in Extreme Value Modeling
Chapter 3 Generalized Logistic Distribution in Extreme Value Modeling 3. Introduction In Chapters and 2, we saw that asymptotic or approximated model for large values is the GEV distribution or GP distribution.
More informationA PRACTICAL WAY FOR ESTIMATING TAIL DEPENDENCE FUNCTIONS
Statistica Sinica 20 2010, 365-378 A PRACTICAL WAY FOR ESTIMATING TAIL DEPENDENCE FUNCTIONS Liang Peng Georgia Institute of Technology Abstract: Estimating tail dependence functions is important for applications
More informationGeneralized additive modelling of hydrological sample extremes
Generalized additive modelling of hydrological sample extremes Valérie Chavez-Demoulin 1 Joint work with A.C. Davison (EPFL) and Marius Hofert (ETHZ) 1 Faculty of Business and Economics, University of
More informationAnalysing geoadditive regression data: a mixed model approach
Analysing geoadditive regression data: a mixed model approach Institut für Statistik, Ludwig-Maximilians-Universität München Joint work with Ludwig Fahrmeir & Stefan Lang 25.11.2005 Spatio-temporal regression
More informationOn heavy tailed time series and functional limit theorems Bojan Basrak, University of Zagreb
On heavy tailed time series and functional limit theorems Bojan Basrak, University of Zagreb Recent Advances and Trends in Time Series Analysis Nonlinear Time Series, High Dimensional Inference and Beyond
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 informationLecture 2: Repetition of probability theory and statistics
Algorithms for Uncertainty Quantification SS8, IN2345 Tobias Neckel Scientific Computing in Computer Science TUM Lecture 2: Repetition of probability theory and statistics Concept of Building Block: Prerequisites:
More informationA Closer Look at the Hill Estimator: Edgeworth Expansions and Confidence Intervals
A Closer Look at the Hill Estimator: Edgeworth Expansions and Confidence Intervals Erich HAEUSLER University of Giessen http://www.uni-giessen.de Johan SEGERS Tilburg University http://www.center.nl EVA
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 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 informationBayesian Inference for Clustered Extremes
Newcastle University, Newcastle-upon-Tyne, U.K. lee.fawcett@ncl.ac.uk 20th TIES Conference: Bologna, Italy, July 2009 Structure of this talk 1. Motivation and background 2. Review of existing methods Limitations/difficulties
More informationStat 542: Item Response Theory Modeling Using The Extended Rank Likelihood
Stat 542: Item Response Theory Modeling Using The Extended Rank Likelihood Jonathan Gruhl March 18, 2010 1 Introduction Researchers commonly apply item response theory (IRT) models to binary and ordinal
More informationTwo practical tools for rainfall weather generators
Two practical tools for rainfall weather generators Philippe Naveau naveau@lsce.ipsl.fr Laboratoire des Sciences du Climat et l Environnement (LSCE) Gif-sur-Yvette, France FP7-ACQWA, GIS-PEPER, MIRACLE
More informationDiscussion on Human life is unlimited but short by Holger Rootzén and Dmitrii Zholud
Extremes (2018) 21:405 410 https://doi.org/10.1007/s10687-018-0322-z Discussion on Human life is unlimited but short by Holger Rootzén and Dmitrii Zholud Chen Zhou 1 Received: 17 April 2018 / Accepted:
More informationMultivariate Non-Normally Distributed Random Variables
Multivariate Non-Normally Distributed Random Variables An Introduction to the Copula Approach Workgroup seminar on climate dynamics Meteorological Institute at the University of Bonn 18 January 2008, Bonn
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 informationRandom Variables and Their Distributions
Chapter 3 Random Variables and Their Distributions A random variable (r.v.) is a function that assigns one and only one numerical value to each simple event in an experiment. We will denote r.vs by capital
More informationApproximating the Conditional Density Given Large Observed Values via a Multivariate Extremes Framework, with Application to Environmental Data
Approximating the Conditional Density Given Large Observed Values via a Multivariate Extremes Framework, with Application to Environmental Data Daniel Cooley 1 Richard A. Davis 2 Philippe Naveau 3 1 Department
More informationModelação de valores extremos e sua importância na
Modelação de valores extremos e sua importância na segurança e saúde Margarida Brito Departamento de Matemática FCUP (FCUP) Valores Extremos - DemSSO 1 / 12 Motivation Consider the following events Occurance
More informationA Note on Tail Behaviour of Distributions. the max domain of attraction of the Frechét / Weibull law under power normalization
ProbStat Forum, Volume 03, January 2010, Pages 01-10 ISSN 0974-3235 A Note on Tail Behaviour of Distributions in the Max Domain of Attraction of the Frechét/ Weibull Law under Power Normalization S.Ravi
More informationStatistics of Extremes
Statistics of Extremes Anthony Davison c 29 http://stat.epfl.ch Multivariate Extremes 184 Europe............................................................. 185 Componentwise maxima..................................................
More informationSkew Generalized Extreme Value Distribution: Probability Weighted Moments Estimation and Application to Block Maxima Procedure
Skew Generalized Extreme Value Distribution: Probability Weighted Moments Estimation and Application to Block Maxima Procedure Pierre Ribereau ( ), Esterina Masiello ( ), Philippe Naveau ( ) ( ) Université
More informationChapter 3 : Likelihood function and inference
Chapter 3 : Likelihood function and inference 4 Likelihood function and inference The likelihood Information and curvature Sufficiency and ancilarity Maximum likelihood estimation Non-regular models EM
More informationProbability and Distributions
Probability and Distributions What is a statistical model? A statistical model is a set of assumptions by which the hypothetical population distribution of data is inferred. It is typically postulated
More informationChapter 5. Chapter 5 sections
1 / 43 sections Discrete univariate distributions: 5.2 Bernoulli and Binomial distributions Just skim 5.3 Hypergeometric distributions 5.4 Poisson distributions Just skim 5.5 Negative Binomial distributions
More informationEVA Tutorial #2 PEAKS OVER THRESHOLD APPROACH. Rick Katz
1 EVA Tutorial #2 PEAKS OVER THRESHOLD APPROACH Rick Katz Institute for Mathematics Applied to Geosciences National Center for Atmospheric Research Boulder, CO USA email: rwk@ucar.edu Home page: www.isse.ucar.edu/staff/katz/
More informationA Few Special Distributions and Their Properties
A Few Special Distributions and Their Properties Econ 690 Purdue University Justin L. Tobias (Purdue) Distributional Catalog 1 / 20 Special Distributions and Their Associated Properties 1 Uniform Distribution
More informationSTAT 992 Paper Review: Sure Independence Screening in Generalized Linear Models with NP-Dimensionality J.Fan and R.Song
STAT 992 Paper Review: Sure Independence Screening in Generalized Linear Models with NP-Dimensionality J.Fan and R.Song Presenter: Jiwei Zhao Department of Statistics University of Wisconsin Madison April
More informationIntroduction to Machine Learning
Introduction to Machine Learning Introduction to Probabilistic Methods Varun Chandola Computer Science & Engineering State University of New York at Buffalo Buffalo, NY, USA chandola@buffalo.edu Chandola@UB
More informationA space-time skew-t model for threshold exceedances
BIOMETRICS 000, 1 20 DOI: 000 000 0000 A space-time skew-t model for threshold exceedances Samuel A Morris 1,, Brian J Reich 1, Emeric Thibaud 2, and Daniel Cooley 2 1 Department of Statistics, North Carolina
More informationIntroduction to Systems Analysis and Decision Making Prepared by: Jakub Tomczak
Introduction to Systems Analysis and Decision Making Prepared by: Jakub Tomczak 1 Introduction. Random variables During the course we are interested in reasoning about considered phenomenon. In other words,
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 informationSpatial extreme value theory and properties of max-stable processes Poitiers, November 8-10, 2012
Spatial extreme value theory and properties of max-stable processes Poitiers, November 8-10, 2012 November 8, 2012 15:00 Clement Dombry Habilitation thesis defense (in french) 17:00 Snack buet November
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 informationBIAS CORRECTION IN MULTIVARIATE EXTREMES. By Anne-Laure Fougères Laurens de Haan and Cécile Mercadier Université Lyon 1 and Erasmus University
Submitted to the Annals of Statistics BIAS CORRECTION IN MULTIVARIATE EXTREMES By Anne-Laure Fougères Laurens de Haan and Cécile Mercadier Université Lyon 1 and Erasmus University The estimation of the
More information