MULTIVARIATE EXTREMES AND RISK
|
|
- Fay Clarke
- 6 years ago
- Views:
Transcription
1 MULTIVARIATE EXTREMES AND RISK Richard L. Smith Department of Statistics and Operations Research University of North Carolina Chapel Hill, NC Interface 2008 RISK: Reality Durham, N.C. May 22,
2 I. OVERVIEW OF UNIVARIATE EXTREME VALUE THEORY II. MULTIVARIATE EXTREME VALUE THEORY III. ALTERNATIVE FORMULATIONS IV. MAX-STABLE PROCESSES 2
3 I. OVERVIEW OF UNIVARIATE EXTREME VALUE THEORY 3
4 INSURANCE RISK EXAMPLE (Smith and Goodman, 2000) The data consist of all insurance claims experienced by a large international oil company over a threshold during a 15-year period a total of 393 claims. Total of all 393 claims: largest claims: 776.2, 268.0, 142.0, 131.0, 95.8, 56.8, 46.2, 45.2, 40.4,
5 Some plots of the insurance data. 5
6 Some problems: 1. What is the distribution of very large claims? 2. Is there any evidence of a change of the distribution over time? 3. What is the influence of the different types of claim? 4. How should one characterize the risk to the company? More precisely, what probability distribution can one put on the amount of money that the company will have to pay out in settlement of large insurance claims over a future time period of, say, three years? 6
7 EXTREME VALUE DISTRIBUTIONS X 1, X 2,..., i.i.d., F (x) = Pr{X i x}, M n = max(x 1,..., X n ), Pr{M n x} = F (x) n. For non-trivial results must renormalize: find a n > 0, b n such that { } Mn b n Pr x = F (a n x + b n ) n H(x). a n The Three Types Theorem (Fisher-Tippett, Gnedenko) asserts that if nondegenerate H exists, it must be one of three types: H(x) = exp( e x ), all x (Gumbel) { 0 x < 0 H(x) = exp( x α ) x > 0 (Fréchet) { exp( x H(x) = α ) x < 0 1 x > 0 (Weibull) In Fréchet and Weibull, α > 0. 7
8 The three types may be combined into a single generalized extreme value (GEV) distribution: ( H(x) = exp 1 + ξ x µ ) 1/ξ ψ, (y + = max(y, 0)) where µ is a location parameter, ψ > 0 is a scale parameter and ξ is a shape parameter. ξ 0 corresponds to the Gumbel distribution, ξ > 0 to the Fréchet distribution with α = 1/ξ, ξ < 0 to the Weibull distribution with α = 1/ξ. ξ > 0: long-tailed case, 1 F (x) x 1/ξ, ξ = 0: exponential tail ξ < 0: short-tailed case, finite endpoint at µ ξ/ψ + 8
9 EXCEEDANCES OVER THRESHOLDS Consider the distribution of X conditionally on exceeding some high threshold u: F u (y) = F (u + y) F (u). 1 F (u) As u ω F = sup{x : F (x) < 1}, often find a limit F u (y) G(y; σ u, ξ) where G is generalized Pareto distribution (GPD) G(y; σ, ξ) = 1 ( 1 + ξ y ) 1/ξ σ + Equivalence to three types theorem established by Pickands (1975).. 9
10 The Generalized Pareto Distribution G(y; σ, ξ) = 1 ( 1 + ξ y ) 1/ξ σ +. ξ > 0: long-tailed (equivalent to usual Pareto distribution), tail like x 1/ξ, ξ = 0: take limit as ξ 0 to get G(y; σ, 0) = 1 exp i.e. exponential distribution with mean σ, ( y ), σ ξ < 0: finite upper endpoint at σ/ξ. 10
11 POINT PROCESS APPROACH Two-dimensional plot of exceedance times and exceedance levels forms a nonhomogeneous Poisson process with Λ(A) = (t 2 t 1 )Ψ(y; µ, ψ, ξ) Ψ(y; µ, ψ, ξ) = ( 1 + ξ y µ ψ ) 1/ξ (1 + ξ(y µ)/ψ > 0). 11
12 Illustration of point process model. 12
13 APPLICATION TO INSURANCE DATA We apply the GPD and point process approaches to the 393 insurance claims described at the beginning of the talk. GPD fits to various thresholds: u N u Mean σ ξ Excess
14 Point process approach: u N u µ log ψ ξ (4.4) (0.24) (0.09) (5.2) (0.31) (0.16) (5.5) (0.31) (0.21) (5.7) (0.32) (0.25) (3.9) (0.46) (0.45) (5.7) (0.56) (0.53) (8.6) (0.66) (0.56) Standard errors are in parentheses 14
15 Predictive Distributions of Future Losses What is the probability distribution of future losses over a specific time period, say 1 year? Let Y be future total loss. Distribution function G(y; µ, ψ, ξ) in practice this must itself be simulated. 15
16 Traditional frequentist approach: where ˆµ, ˆψ, ˆξ are MLEs. Ĝ(y) = G(y; ˆµ, ˆψ, ˆξ) Bayesian: G(y) = G(y; µ, ψ, ξ)dπ(µ, ψ, ξ X) where π( X) denotes posterior density given data X. 16
17 Estimated posterior densities for the three parameters, and for the predictive distribution function. Four independent Monte Carlo runs are shown for each plot. 17
18 II. MULTIVARIATE EXTREME VALUE THEORY 18
19 Multivariate extreme value theory applies when we are interested in the joint distribution of extremes from several random variables. Examples: Winds and waves on an offshore structure Meteorological variables, e.g. temperature and precipitation Air pollution variables, e.g. ozone and sulfur dioxide Finance, e.g. price changes in several stocks or indices Spatial extremes, e.g. joint distributions of extreme precipitation at several locations 19
20 LIMIT THEOREMS FOR MULTIVARIATE SAMPLE MAXIMA Let Y i = (Y i1...y id ) T be i.i.d. D-dimensional vectors, i = 1, 2,... M nd = max{y 1d,..., Y nd }(1 d D) d th-component maximum Look for constants a nd, b nd such that { Mnd b Pr nd x d, d = 1,..., D a nd } G(x 1,..., x D ). Vector notation: Pr { Mn b n a n x } G(x). 20
21 Two general points: 1. If G is a multivariate extreme value distribution, then each of its marginal distributions must be one of the univariate extreme value distributions, and therefore can be represented in GEV form 2. The form of the limiting distribution is invariant under monotonic transformation of each component. Therefore, without loss of generality we can transform each marginal distribution into a specified form. For most of our applications, it is convenient to assume the Fréchet form: Pr{X d x} = exp( x α ), x > 0, d = 1,..., D. Here α > 0. The case α = 1 is called unit Fréchet. 21
22 Basics of Multivariate Regular Variation (following Resnick (2007), Chapter 6) After transformation of margins, lim t t Pr { Yi b(t) A } = ν(a) b regularly varying function of index α > 0 (i.e. t for fixed x > 0), ν a measure on the cone satisfying for any scalar t > 0. E = [0, ] D {0} ν(ta) = t α ν(a) b(tx) b(t) xα as 22
23 X2 ta 0 A X1 23
24 The last statement implies that ν can be decomposed into a product of radial and angular components. Define S D = {(x 1,..., x D ) : x 1 0,..., x D 0, x x D = 1}. Consider sets A of form A = x E : x > r, for some S S D, x = D j=1 x j. { } x x S Then for some measure H on S D. ν(a) = r α H(S) 24
25 First Interpretation: Consider i.i.d. MRV. vectors {X i, i = 1, 2,...} whose distribution is Let P n be a measure on [0, ] D consisting of the points { X 1 b(n),..., X } n b(n). Let A be a measurable set on E, then the expected number of points of P n in A is { } Xi n Pr b(n) A ν(a) as n. With some measure-theoretic formalities, this shows that P n converges vaguely to a nonhomogeneous Poisson process on E with intensity measure ν. 25
26 Second Interpretation: Fix x 1 0,..., x D 0, D d=1 x d > 0. Let A be the complement of Then [0, x 1 ] [0, x 2 ]... [0, x D ]. { Pr n Xi1 b(n) x 1,..., X } id b(n) x D is the probability that P n places no points in the set A. By Poisson limit theorem, this probability tends to e ν(a) as n. Therefore, the limit of (1) is where V (x) = ν(a). (1) G(x) = exp { V (x)} (2) Moreover, using the radial-spectral decomposition of ν, V (x) = S D max d=1,...,d ( ) wj x j dh(w). (3) 26
27 The function V (x) is called the exponent measure and formula (3) is the Pickands representation. If we fix d {1,..., D} with 0 < x d <, and define x d = + for d d, then so we must have V (x) = S D = 1 x d max d=1,...,d ( ) wd x d S D w d dh(w) dh(w) S D w d dh(w) = 1, d = 1,..., D, (4) to ensure that the marginal distributions are correct. 27
28 Note that kv (x) = V ( x k) (which is in fact another characterization of V ) so G k (x) = exp( kv (x)) = exp = G ( ( x k) V. ( x k)) Hence G is max-stable. In particular, if X 1,..., X k are i.i.d. from G, then max{x 1,..., X k } (vector of componentwise maxima) has the same distribution as kx 1. 28
29 EXAMPLES 29
30 Logistic (Gumbel and Goldstein, 1964) Check: 1. V (x/k) = kv (x) V (x) = D d=1 x r d 1/r 2. V ((+, +,..., x d,..., +, + ) = x 1 d 3. e V (x) is a valid c.d.f., r 1. Limiting cases: r = 1: independent components r : the limiting case when X i1 = X i2 =... = X id with probability 1. 30
31 Asymmetric logistic (Tawn 1990) V (x) = c C ( ) rc 1/r θi,c c, i c x i where C is the class of non-empty subsets of {1,..., D}, r c 1, θ i,c = 0 if i / c, θ i,c 0, c C θ i,c = 1 for each i. Negative logistic (Joe 1990) V (x) = 1 x j + c C: c 2 ( 1) c ( ) rc 1/r θi,c c, i c x i r c 0, θ i,c = 0 if i / c, θ i,c 0, c C( 1) c θ i,c 1 for each i. 31
32 Tilted Dirichlet (Coles and Tawn 1991) A general construction: Suppose h is an arbitrary positive function on S d with m d = S D u d h (u)du <, then define h(w) = ( m k w k ) (D+1) D d=1 ( m d h m1 w 1,..., m ) Dw D. mk w k mk w k h is density of positive measure H, S D u d dh(u) = 1 each d. As a special case of this, they considered Dirichlet density Leads to h(w) = D d=1 h (u) = Γ( α j ) d Γ(α d ) D d=1 α d Γ(α d ) Γ( αd + 1) ( α d w d ) D+1 u α d 1 d. D d=1 Disadvantage: need for numerical integration ( ) αd 1 αd w d. αk w k 32
33 ESTIMATION Parametric Non/semi-parametric Both approaches have problems, e.g. nonregular behavior of MLE even in finite-parameter problems; curse of dimensionality if D large. Typically proceed by transforming margins to unit Fréchet first, though there are advantages in joint estimation of marginal distributions and dependence structure (Shi 1995) 33
34 III. ALTERNATIVE FORMULATIONS OF MULTIVARIATE EXTREMES 34
35 Ledford-Tawn-Ramos approach Heffernan-Tawn approach However the Heffernan-Tawn approach will be the subject of a later talk so we concentrate here on the Ledford-Tawn-Ramos approach 35
36 The first paper to suggest that multivariate extreme value theory (as defined so far) might not be general enough was Ledford and Tawn (1996). Suppose (Z 1, Z 2 ) are a bivariate random vector with unit Fréchet margins. Traditional cases lead to { const. r 1 dependent cases Pr{Z 1 > r, Z 2 > r} const. r 2 exact independent The first case covers all bivariate extreme value distributions except the independent case. However, Ledford and Tawn showed by example that for a number of cases of practical interest, Pr{Z 1 > r, Z 2 > r} L(r)r 1/η, where L is a slowly varying function ( L(rx) L(r) η (0, 1]. 1 as r ) and Estimation: used fact that 1/η is Pareto index for min(z 1, Z 2 ). 36
37 More general case (Ledford and Tawn 1997): Pr{Z 1 > z 1, Z 2 > z 2, } = L(z 1, z 2 )z c 1 1 z c 2 2, 0 < η 1; c 1 + c 2 = 1 η ; L slowly varying in sense that L(tz g(z 1, z 2 ) = lim 1, tz 2 ). t L(t, t) ( ) They showed g(z 1, z 2 ) = g z1 z 1 +z but were unable to estimate 2 g directly needed to make parametric assumptions about this. More recently, Resnick and co-authors were able to make a more rigorous mathematical theory using concept of hidden regular variation (see e.g. Resnick 2002, Maulik and Resnick 2005, Heffernan and Resnick 2005; see also Section 9.4 of Resnick (2007)). 37
38 Ramos-Ledford (2007a, 2007b) approach: Pr{Z 1 > z 1, Z 2 > z 2, } = L(z 1, z 2 )(z 1 z 2 ) 1/(2η), L(ux lim 1, ux 2 ) = g(x u 1, x 2 ) L(u, u) ( ) x1 = g x 1 + x 2 Limiting joint survivor function Pr{X 1 > x 1, X 2 > x 2 } = 1 = η { 1 η 0 1 min ( w, 1 w )} 1/η dh x 1 x η(w), 2 0 {min (w, 1 w)}1/η dh η(w). 38
39 Multivariate generalization to D > 2: Pr{X 1 > x 1,..., X D > x D } = 1 = η S D η S D { min 1 d D { ( wd min 1 d D w j x d )} 1/η dh η (w), } 1/η dh η(w). Open problems: How to find sufficiently rich classes of H η(w) to derive parametric families suitable for real data How to do the subsequent estimation 39
40 IV. MAX-STABLE PROCESSES 40
41 Suppose Y t, t T is a stochastic process (e.g. spatial, temporal) Normalization condition: for each n 1, there exist constants a nt, b nt, t T such that, for any t 1,..., t m T, Pr n Y tj b ntj a ntj x tj, j = 1,..., m G t 1,...,t m (x t1,..., x tm ).(5) Then G t1,...,t m (x t1,..., x tm ) is a MEVD. When (5) holds in a stochastic process sense the limiting process is called max-stable. 41
42 For this section of the talk, I describe a particular class of max-stable processes called M4 processes (Smith and Weissman (1996),Zhang (2002),Smith (2003) Chamú Morales (2005),Zhang and Smith (2007)). Alternative approaches will be described in the talk by D. Cooley. 42
43 Suppose {X id, i = 0, ±1, ±2,..., d = 1,..., D} form a stationary time series in D dimensions. Suppose we transform each X id to a variable Y id which is unit Fréchet (P {Y y} = e 1/y ). For example, a common method would be to fit a GPD to all exceedances of X id above a high threshold, use the empirical distribution function for all values below the threshold, and apply the probability integral transformation. The process {Y id } is called a multivariate maxima of moving maxima process (M4 for short), if it is defined by Y id = max l max k a l,k,d Z l,i k. k a l,k,d = 1 for where a l,k,d 0 are coefficients such that l each d, and {Z l,i, l = 1,..., L, i = 0, ±1, ±2,..., } is a double array of independent unit Fréchet random variables. 43
44 Elementary Properties Pr{Y id y id for all i, d} = Pr{a l,k,d Z l,i k y id for all l, k, i, d} So = Pr{a l,k,d Z l,i y i+k,d for all l, k, i, d} = Pr{Z l,i min = exp l k,d i y i+k,d a l,k,d for all l, i} max k,d a l,k,d. y i+k,d Pr{Y id ny id for all i, d} n = Pr{Y id y id for all i, d} (6) and for a particular (i, d), Pr{Y id y id } = exp ( ) l k a l,k,d y i,d = exp ( 1 y i,d ). (7) (7) shows that the margins are unit Fréchet while (6) shows that the joint distributions are max-stable. 44
45 If we ignore the time series component and just treat this as a model for a single random vector (Y 1,..., Y D ), the equation simplifies to Pr{Y d y d, d = 1,..., D} = exp l Recall the formula for the spectral measure: Pr{Y d y d, d = 1,..., D} = exp { S D max j=1,...,d max d a l,d y d ( ) wj y j. (8) dh(w) where S D is simplex. Thus, (8) is equivalent to defining a spectral measure H that is discrete with mass at points (a l,1,..., a l,d ) d a l,d S D. The condition S D w d H(w) = 1, for d = 1,..., D, is equivalent to the statement that l a l,d = 1 for all d. } 45
46 Difficulties of Estimation Y id = max l max k a l,k,d Z l,i k. Problem of signature patterns. If a specific Z l,k is much larger than its neighbors, then for i close to k, Y id a l,i k,d. Under the model, this pattern will be replicated throughout the process. Such patterns cannot be expected in real data, so we must devise methods of estimation that take account of the fact that the M4 process is not an exact model. 46
47 Davis and Resnick (1989, 1993) developed many properties of max-arma processes but did not devise a good statistical approach. Hall, Peng and Yao (2002) proposed an alternative estimation scheme for max-arma processes, based on the empirical distribution distribution. Chamú Morales (PhD thesis, 2005) developed an approximate MCMC technique Zhang and Smith (2007) (but originally Zhang (2002)) generalized the Hall-Peng-Yao technique to M4 processes 47
48 APPLICATION OF M4 PROCESSES TO FINANCIAL DATA (Smith 2003) 48
49 We use neg daily returns from closing prices of stock prices in three companies, Pfizer, GE and Citibank. Univariate EVT models were fitted to threshold u =.02 with results: Series N u µ log ψ ξ (SE) (SE) (SE) Pfizer (.0029) (.132) (.051) GE (.0029) (.143) (.062) Citibank (.0036) (.119) (.012) However, diagnostic plots show a problem sign of volatility. (See in particular plot (b) QQ plot based on intervals between threshold exceedances these are clearly inconsistent with a simple Poisson processes.) 49
50 50
51 Estimate volatility by fitting a GARCH(1,1) model to each series. Standardize returns data, recompute threshold analysis (u = 1.2). Diagnostic plots now look OK (Pfizer shown, others similar) Series N u µ log ψ ξ (SE) (SE) (SE) Pfizer (.155) (.148) (.061) GE (.130) (.128) (.053) Citibank (.157) (.126) (.050) 51
52 52
53 53
54 After fitting a univariate extreme value model to each series, the exceedances over the threshold for each series are transformed to have marginal Fréchet distributions. On the transformed scale, the data in each series consist of all values in excess of threshold 1. On this transformed scale, pairwise scatterplots are shown of the three series against each other on the same day (top 3 plots), and against series on neighboring days. The two numbers on each plot show the expected number of joint exceedances based on an independence assumption, and the observed number of joint exceedances. 54
55 Pfizer day 0 GE day (35,107) Pfizer day 0 Citibank day (31,91) GE day 0 Citibank day (31,122) Pfizer day 0 GE day (35,36) Pfizer day 0 Citibank day (31,33) GE day 0 Citibank day (31,47) Pfizer day 0 GE day (35,43) Pfizer day 0 Citibank day (31,41) GE day 0 Citibank day (31,42) Pfizer day 0 Pfizer day (35,39) GE day 0 GE day (36,43) Citibank day 0 Citibank day (27,35) 55
56 Also shown is a plot of Fréchet exceedances for the three series on the same day, normalized to have total 1, plotted in barycentric coordinates. The three circles near the corner points P, G and C correspond to days for which that series along had an exceedance. 56
57 57
58 After fitting an M4 model, a new data set was simulated using the fitted model. The following figures correspond to the earlier ones, but using simulated data. This is a check on the realism of the fitted model. 58
59 59
60 Pfizer day 0 GE day (39,93) Pfizer day 0 Citibank day (32,62) GE day 0 Citibank day (32,91) Pfizer day 0 GE day (39,65) Pfizer day 0 Citibank day (32,47) GE day 0 Citibank day (32,46) Pfizer day 0 GE day (39,71) Pfizer day 0 Citibank day (32,46) GE day 0 Citibank day (32,49) Pfizer day 0 Pfizer day (38,61) GE day 0 GE day (40,68) Citibank day 0 Citibank day (26,42) 60
61 61
62 Finally, we attempt to validate the model by calibrating observed vs. expected probabilities of extreme events under the model. The extreme event considered is that there is at least one exceedance of a specific threshold u by one of the three series in one of the next 10 days after a given day. To make the comparison honest, the period of study is divided into four periods each of length just uder 5 years. The univariate and multivariate EV model is fitted to each of the first three 5- year period, and used to predict extreme events in the following period. The final plot shows observed (dashed lines) and expected (solid lines) counts for a sequence of thresholds u. There is excellent agreement between the two curves. 62
63 63
64 REFERENCES 64
65 Chamú Morales, F. (2005), Estimation of max-stable processes using Monte Carlo methods with applications to financial risk assessment. Ph.D. dissertation, Department of Statistics and Operations Research, University of North Carolina, Chapel Hill. Coles, S.G. (2001), An Introduction to Statistical Modeling of Extreme Values. Springer Verlag, New York. Coles, S.G. and Tawn, J.A. (1991), Modelling extreme multivariate events. J.R. Statist. Soc. B 53, Davis, R.A. and Resnick, S.I. (1989), Basic properties and prediction of max-arma processes. Advances in Applied Probability 21, Davis, R.A. and Resnick, S.I. (1993), Prediction of stationary max-stable processes. Annals of Applied Probability 3, Finkenstadt, B. and Rootzén, H. (editors) (2003), Extreme Values in Finance, Telecommunications and the Environment. Chapman and Hall/CRC Press, London. 65
66 Gumbel, E.J. and Goldstein, N. (1964), Analysis of empirical bivariate extremal distributions. J. Amer. Statist. Soc. 59, Hall, P., Peng, L. and Yao, Q. (2002), Moving-maximum models for extrema of time series. Journal of Statistical Planning and Inference 103, Heffernan, J.E. and Resnick, S. (2005), Hidden regular variation and the rank transform. Advances in Applied Probabilitty 37-2, Joe, H. (1990), Families of min-stable multivariate exponential and multivariate extreme value distributions. Statist. and Probab. Letters 9, Ledford, A.W. and Tawn, J.A. (1996), Modelling dependence within joint tail regions, J.R. Statist. Soc. B 59, Ledford, A.W. and Tawn, J.A. (1997), Statistics for near independence in multivariate extreme values. Biometrika 83,
67 Maulik, K. and Resnick, S. (2005), Characterizations and examples of hidden regular variation. Extremes 7, Ramos, A. and Ledford, A. (2007a), A new class of models for bivariate joint tails. Preprint. Ramos, A. and Ledford, A. (2007b), An alternative point process framework for modeling multivariate extreme values. Preprint. Resnick, S. (2002), Second order regular variation and asymptotic independence. Extremes 5, Resnick, S. (2007), Heavy-Tail Phenomena: Probabilistic and Statistical Modeling. Springer, New York. Shi, D. (1995), Fisher information for a multivariate extreme value distribution. Biometrika 82, Smith, R.L. (2003), Statistics of extremes, with applications in environment, insurance and finance. Chapter 1 of Extreme Values in Finance, Telecommunications and the Environment, edited by B. Finkenstädt and H. Rootzén, Chapman and Hall/CRC Press, London, pp
68 Smith, R.L. and Goodman, D. (2000), Bayesian risk analysis. Chapter 17 of Extremes and Integrated Risk Management, edited by P. Embrechts. Risk Books, London, Smith, R.L. and Weissman, I. (1996), Characterization and estimation of the multivariate extremal index. Unpublished paper, available from Tawn, J.A. (1990), Modelling multivariate extreme value distributions. Biometrika 77, Zhang, Z. (2002), Multivariate Extremes, Max-Stable Process Estimation and Dynamic Financial Modeling. Ph.D. dissertation, Department of Statistics, University of North Carolina, Chapel Hill. Zhang, Z. and Smith, R.L. (2007), On the estimation and application of max-stable processes. (Re-)Submitted for publication. 68
HIERARCHICAL 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 informationRISK AND EXTREMES: ASSESSING THE PROBABILITIES OF VERY RARE EVENTS
RISK AND EXTREMES: ASSESSING THE PROBABILITIES OF VERY RARE EVENTS Richard L. Smith Department of Statistics and Operations Research University of North Carolina Chapel Hill, NC 27599-3260 rls@email.unc.edu
More informationRISK ANALYSIS AND EXTREMES
RISK ANALYSIS AND EXTREMES Richard L. Smith Department of Statistics and Operations Research University of North Carolina Chapel Hill, NC 27599-3260 rls@email.unc.edu Opening Workshop SAMSI program on
More informationRichard L. Smith Department of Statistics and Operations Research University of North Carolina Chapel Hill, NC
EXTREME VALUE THEORY Richard L. Smith Department of Statistics and Operations Research University of North Carolina Chapel Hill, NC 27599-3260 rls@email.unc.edu AMS Committee on Probability and Statistics
More informationNew Classes of Multivariate Survival Functions
Xiao Qin 2 Richard L. Smith 2 Ruoen Ren School of Economics and Management Beihang University Beijing, China 2 Department of Statistics and Operations Research University of North Carolina Chapel Hill,
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 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 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 informationThe Behavior of Multivariate Maxima of Moving Maxima Processes
The Behavior of Multivariate Maxima of Moving Maxima Processes Zhengjun Zhang Department of Mathematics Washington University Saint Louis, MO 6313-4899 USA Richard L. Smith Department of Statistics University
More informationExtreme Precipitation: An Application Modeling N-Year Return Levels at the Station Level
Extreme Precipitation: An Application Modeling N-Year Return Levels at the Station Level Presented by: Elizabeth Shamseldin Joint work with: Richard Smith, Doug Nychka, Steve Sain, Dan Cooley Statistics
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 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 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 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 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 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 informationTail dependence in bivariate skew-normal and skew-t distributions
Tail dependence in bivariate skew-normal and skew-t distributions Paola Bortot Department of Statistical Sciences - University of Bologna paola.bortot@unibo.it Abstract: Quantifying dependence between
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 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 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 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 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 Richard L Smith Department of Statistics University of North
More informationClassical Extreme Value Theory - An Introduction
Chapter 1 Classical Extreme Value Theory - An Introduction 1.1 Introduction Asymptotic theory of functions of random variables plays a very important role in modern statistics. The objective of the asymptotic
More informationCONTAGION VERSUS FLIGHT TO QUALITY IN FINANCIAL MARKETS
EVA IV, CONTAGION VERSUS FLIGHT TO QUALITY IN FINANCIAL MARKETS Jose Olmo Department of Economics City University, London (joint work with Jesús Gonzalo, Universidad Carlos III de Madrid) 4th Conference
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 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 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 informationEmma Simpson. 6 September 2013
6 September 2013 Test What is? Beijing during periods of low and high air pollution Air pollution is composed of sulphur oxides, nitrogen oxides, carbon monoxide and particulates. Particulates are small
More informationStatistics of Extremes
Statistics of Extremes Anthony Davison c 211 http://stat.epfl.ch Multivariate Extremes 19 Componentwise maxima.................................................. 194 Standardization........................................................
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 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 informationTail Approximation of Value-at-Risk under Multivariate Regular Variation
Tail Approximation of Value-at-Risk under Multivariate Regular Variation Yannan Sun Haijun Li July 00 Abstract This paper presents a general tail approximation method for evaluating the Valueat-Risk of
More informationOn the Estimation and Application of Max-Stable Processes
On the Estimation and Application of Ma-Stable Processes Zhengjun Zhang Department of Mathematics Washington University Saint Louis, MO 63130-4899 USA Richard L Smith Department of Statistics University
More informationA New Class of Tail-dependent Time Series Models and Its Applications in Financial Time Series
A New Class of Tail-dependent Time Series Models and Its Applications in Financial Time Series Zhengjun Zhang Department of Mathematics, Washington University, Saint Louis, MO 63130-4899, USA Abstract
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 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 informationAPPLICATION OF EXTREMAL THEORY TO THE PRECIPITATION SERIES IN NORTHERN MORAVIA
APPLICATION OF EXTREMAL THEORY TO THE PRECIPITATION SERIES IN NORTHERN MORAVIA DANIELA JARUŠKOVÁ Department of Mathematics, Czech Technical University, Prague; jarus@mat.fsv.cvut.cz 1. Introduction The
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 informationEstimating Bivariate Tail: a copula based approach
Estimating Bivariate Tail: a copula based approach Elena Di Bernardino, Université Lyon 1 - ISFA, Institut de Science Financiere et d'assurances - AST&Risk (ANR Project) Joint work with Véronique Maume-Deschamps
More informationSpatial Extremes in Atmospheric Problems
Spatial Extremes in Atmospheric Problems Eric Gilleland Research Applications Laboratory (RAL) National Center for Atmospheric Research (NCAR), Boulder, Colorado, U.S.A. http://www.ral.ucar.edu/staff/ericg
More informationMultivariate Heavy Tails, Asymptotic Independence and Beyond
Multivariate Heavy Tails, endence and Beyond Sidney Resnick School of Operations Research and Industrial Engineering Rhodes Hall Cornell University Ithaca NY 14853 USA http://www.orie.cornell.edu/ sid
More informationZwiers FW and Kharin VV Changes in the extremes of the climate simulated by CCC GCM2 under CO 2 doubling. J. Climate 11:
Statistical Analysis of EXTREMES in GEOPHYSICS Zwiers FW and Kharin VV. 1998. Changes in the extremes of the climate simulated by CCC GCM2 under CO 2 doubling. J. Climate 11:2200 2222. http://www.ral.ucar.edu/staff/ericg/readinggroup.html
More informationModelling Multivariate Peaks-over-Thresholds using Generalized Pareto Distributions
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
More informationA conditional approach for multivariate extreme values
J. R. Statist. Soc. B (2004) 66, Part 3, pp. 497 546 A conditional approach for multivariate extreme values Janet E. Heffernan and Jonathan A. Tawn Lancaster University, UK [Read before The Royal Statistical
More informationThe estimation of M4 processes with geometric moving patterns
AISM (2008) 60:121 150 DOI 10.1007/s10463-006-0078-0 The estimation of M4 processes with geometric moving patterns Zhengjun Zhang Received: 16 September 2003 / Revised: 28 June 2006 / Published online:
More informationExtreme 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 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 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 informationEXTREMAL MODELS AND ENVIRONMENTAL APPLICATIONS. Rick Katz
1 EXTREMAL MODELS AND ENVIRONMENTAL APPLICATIONS Rick Katz Institute for Study of Society and Environment National Center for Atmospheric Research Boulder, CO USA email: rwk@ucar.edu Home page: www.isse.ucar.edu/hp_rick/
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 informationWhat Can We Infer From Beyond The Data? The Statistics Behind The Analysis Of Risk Events In The Context Of Environmental Studies
What Can We Infer From Beyond The Data? The Statistics Behind The Analysis Of Risk Events In The Context Of Environmental Studies Sibusisiwe Khuluse, Sonali Das, Pravesh Debba, Chris Elphinstone Logistics
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 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 informationShape of the return probability density function and extreme value statistics
Shape of the return probability density function and extreme value statistics 13/09/03 Int. Workshop on Risk and Regulation, Budapest Overview I aim to elucidate a relation between one field of research
More informationLecture 2 APPLICATION OF EXREME VALUE THEORY TO CLIMATE CHANGE. Rick Katz
1 Lecture 2 APPLICATION OF EXREME VALUE THEORY TO CLIMATE CHANGE Rick Katz Institute for Study of Society and Environment National Center for Atmospheric Research Boulder, CO USA email: rwk@ucar.edu Home
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 informationPrediction for Max-Stable Processes via an Approximated Conditional Density
Prediction for Max-Stable Processes via an Approximated Conditional Density Daniel Cooley 1 Richard A. Davis 2 Philippe Naveau 3 1 Department of Statistics Colorado State University, Fort Collins, CO USA
More informationUNIVERSITY OF CALGARY. Inference for Dependent Generalized Extreme Values. Jialin He A THESIS SUBMITTED TO THE FACULTY OF GRADUATE STUDIES
UNIVERSITY OF CALGARY Inference for Dependent Generalized Extreme Values by Jialin He A THESIS SUBMITTED TO THE FACULTY OF GRADUATE STUDIES IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE OF
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 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 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 informationSemi-parametric estimation of non-stationary Pickands functions
Semi-parametric estimation of non-stationary Pickands functions Linda Mhalla 1 Joint work with: Valérie Chavez-Demoulin 2 and Philippe Naveau 3 1 Geneva School of Economics and Management, University of
More informationSTATISTICS FOR CLIMATE SCIENCE
STATISTICS FOR CLIMATE SCIENCE Richard L Smith University of North Carolina and SAMSI VI-MSS Workshop on Environmental Statistics Kolkata, March 2-4, 2015 www.unc.edu/~rls/kolkata.html 1 1 In Memoriam
More informationOn the Application of the Generalized Pareto Distribution for Statistical Extrapolation in the Assessment of Dynamic Stability in Irregular Waves
On the Application of the Generalized Pareto Distribution for Statistical Extrapolation in the Assessment of Dynamic Stability in Irregular Waves Bradley Campbell 1, Vadim Belenky 1, Vladas Pipiras 2 1.
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 informationExtremogram and ex-periodogram for heavy-tailed time series
Extremogram and ex-periodogram for heavy-tailed time series 1 Thomas Mikosch University of Copenhagen Joint work with Richard A. Davis (Columbia) and Yuwei Zhao (Ulm) 1 Zagreb, June 6, 2014 1 2 Extremal
More informationStatistics for extreme & sparse data
Statistics for extreme & sparse data University of Bath December 6, 2018 Plan 1 2 3 4 5 6 The Problem Climate Change = Bad! 4 key problems Volcanic eruptions/catastrophic event prediction. Windstorms
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 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 informationStatistics of Extremes
Statistics of Extremes Anthony Davison c 29 http://stat.epfl.ch Multivariate Extremes 184 Europe............................................................. 185 Componentwise maxima..................................................
More informationNon parametric modeling of multivariate extremes with Dirichlet mixtures
Non parametric modeling of multivariate extremes with Dirichlet mixtures Anne Sabourin Institut Mines-Télécom, Télécom ParisTech, CNRS-LTCI Joint work with Philippe Naveau (LSCE, Saclay), Anne-Laure Fougères
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 informationA class of probability distributions for application to non-negative annual maxima
Hydrol. Earth Syst. Sci. Discuss., doi:.94/hess-7-98, 7 A class of probability distributions for application to non-negative annual maxima Earl Bardsley School of Science, University of Waikato, Hamilton
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 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 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 informationIntroduction A basic result from classical univariate extreme value theory is expressed by the Fisher-Tippett theorem. It states that the limit distri
The dependence function for bivariate extreme value distributions { a systematic approach Claudia Kluppelberg Angelika May October 2, 998 Abstract In this paper e classify the existing bivariate models
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 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 informationA simple graphical method to explore tail-dependence in stock-return pairs
A simple graphical method to explore tail-dependence in stock-return pairs Klaus Abberger, University of Konstanz, Germany Abstract: For a bivariate data set the dependence structure can not only be measured
More informationEstimation of the Angular Density in Multivariate Generalized Pareto Models
in Multivariate Generalized Pareto Models René Michel michel@mathematik.uni-wuerzburg.de Institute of Applied Mathematics and Statistics University of Würzburg, Germany 18.08.2005 / EVA 2005 The Multivariate
More informationExtremogram and Ex-Periodogram for heavy-tailed time series
Extremogram and Ex-Periodogram for heavy-tailed time series 1 Thomas Mikosch University of Copenhagen Joint work with Richard A. Davis (Columbia) and Yuwei Zhao (Ulm) 1 Jussieu, April 9, 2014 1 2 Extremal
More informationData. Climate model data from CMIP3
Data Observational data from CRU (Climate Research Unit, University of East Anglia, UK) monthly averages on 5 o x5 o grid boxes, aggregated to JJA average anomalies over Europe: spatial averages over 10
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 informationStatistical Methods for Clusters of Extreme Values
Statistical Methods for Clusters of Extreme Values Christopher A. T. Ferro, B.Sc. M.Sc. Submitted for the degree of Doctor of Philosophy at Lancaster University, September 2003. I declare that the work
More informationAccounting for Choice of Measurement Scale in Extreme Value Modeling
Accounting for Choice of Measurement Scale in Extreme Value Modeling J. L. Wadsworth Department of Mathematics and Statistics Lancaster University, LA1 4YF, UK J. A. Tawn Department of Mathematics and
More informationPENULTIMATE APPROXIMATIONS FOR WEATHER AND CLIMATE EXTREMES. Rick Katz
PENULTIMATE APPROXIMATIONS FOR WEATHER AND CLIMATE EXTREMES Rick Katz Institute for Mathematics Applied to Geosciences National Center for Atmospheric Research Boulder, CO USA Email: rwk@ucar.edu Web site:
More informationMultivariate extremes. Anne-Laure Fougeres. Laboratoire de Statistique et Probabilites. INSA de Toulouse - Universite Paul Sabatier 1
Multivariate extremes Anne-Laure Fougeres Laboratoire de Statistique et Probabilites INSA de Toulouse - Universite Paul Sabatier 1 1. Introduction. A wide variety of situations concerned with extreme events
More informationEstimation of Operational Risk Capital Charge under Parameter Uncertainty
Estimation of Operational Risk Capital Charge under Parameter Uncertainty Pavel V. Shevchenko Principal Research Scientist, CSIRO Mathematical and Information Sciences, Sydney, Locked Bag 17, North Ryde,
More informationTail Dependence of Multivariate Pareto Distributions
!#"%$ & ' ") * +!-,#. /10 243537698:6 ;=@?A BCDBFEHGIBJEHKLB MONQP RS?UTV=XW>YZ=eda gihjlknmcoqprj stmfovuxw yy z {} ~ ƒ }ˆŠ ~Œ~Ž f ˆ ` š œžÿ~ ~Ÿ œ } ƒ œ ˆŠ~ œ
More informationChallenges in implementing worst-case analysis
Challenges in implementing worst-case analysis Jon Danielsson Systemic Risk Centre, lse,houghton Street, London WC2A 2AE, UK Lerby M. Ergun Systemic Risk Centre, lse,houghton Street, London WC2A 2AE, UK
More information100-year and 10,000-year Extreme Significant Wave Heights How Sure Can We Be of These Figures?
100-year and 10,000-year Extreme Significant Wave Heights How Sure Can We Be of These Figures? Rod Rainey, Atkins Oil & Gas Jeremy Colman, Independent Consultant Wave crest elevations: BP s EI 322 A Statoil
More informationExtreme Value Theory as a Theoretical Background for Power Law Behavior
Extreme Value Theory as a Theoretical Background for Power Law Behavior Simone Alfarano 1 and Thomas Lux 2 1 Department of Economics, University of Kiel, alfarano@bwl.uni-kiel.de 2 Department of Economics,
More informationEXTREMAL QUANTILES OF MAXIMUMS FOR STATIONARY SEQUENCES WITH PSEUDO-STATIONARY TREND WITH APPLICATIONS IN ELECTRICITY CONSUMPTION ALEXANDR V.
MONTENEGRIN STATIONARY JOURNAL TREND WITH OF ECONOMICS, APPLICATIONS Vol. IN 9, ELECTRICITY No. 4 (December CONSUMPTION 2013), 53-63 53 EXTREMAL QUANTILES OF MAXIMUMS FOR STATIONARY SEQUENCES WITH PSEUDO-STATIONARY
More informationFORECAST VERIFICATION OF EXTREMES: USE OF EXTREME VALUE THEORY
1 FORECAST VERIFICATION OF EXTREMES: USE OF EXTREME VALUE THEORY Rick Katz Institute for Study of Society and Environment National Center for Atmospheric Research Boulder, CO USA Email: rwk@ucar.edu Web
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 informationModels for Spatial Extremes. Dan Cooley Department of Statistics Colorado State University. Work supported in part by NSF-DMS
Models for Spatial Extremes Dan Cooley Department of Statistics Colorado State University Work supported in part by NSF-DMS 0905315 Outline 1. Introduction Statistics for Extremes Climate and Weather 2.
More informationModels for the extremes of Markov chains
Biometrika (1998), 85, 4, pp. 851-867 Printed in Great Britain Models for the extremes of Markov chains BY PAOLA BORTOT Department of Statistical Sciences, University of Padova, 35121 Padova, Italy bortot@hal.stat.unipd.it
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 informationExtreme Value Analysis of Multivariate High Frequency Wind Speed Data
Extreme Value Analysis of Multivariate High Frequency Wind Speed Data Christina Steinkohl Center of Mathematical Sciences, Technische Universität München, D-85748 Garching, Germany Richard A. Davis Department
More information