Extreme Value Analysis and Spatial Extremes
|
|
- Meghan Neal
- 5 years ago
- Views:
Transcription
1 Extreme Value Analysis and Department of Statistics Purdue University 11/07/2013
2 Outline Motivation 1 Motivation 2 Extreme Value Theorem and 3 Bayesian Hierarchical Models Copula Models Max-stable Models
3 What are (Spatial) Extremes? Extreme events are defined to be rare and unexpected 100 year flood the level of flood water expected to be equaled or exceeded every 100 years on average Many extreme events of interest are spatio-temporal in natural
4 What are (Spatial) Extremes? Extreme events are defined to be rare and unexpected 100 year flood the level of flood water expected to be equaled or exceeded every 100 years on average Many extreme events of interest are spatio-temporal in natural
5 What are (Spatial) Extremes? Extreme events are defined to be rare and unexpected 100 year flood the level of flood water expected to be equaled or exceeded every 100 years on average Many extreme events of interest are spatio-temporal in natural
6 What are (Spatial) Extremes? Extreme events are defined to be rare and unexpected 100 year flood the level of flood water expected to be equaled or exceeded every 100 years on average Many extreme events of interest are spatio-temporal in natural
7 Why study extremes? Although infrequent, extremes have large human impact European heat wave example. Around 70,000 were killed!
8 Why study extreme? "There is always going to be an element of doubt, as one is extrapolating into into areas one does not know about. But what EVT is doing is making the best use of whatever data you have about extreme phenomena." Richard Smith
9 Usual vs Extremes Relies on asymptotic theory to provide models for the tail Uses only the "extreme" observations to fit the model
10 History Motivation 1920 s: Foundations of asymptotic argument developed by Fisher and Tippett 1940 s: Asymptotic theory unified and extended by Gnedenko and von Mises 1950 s: Use of asymptotic distributions for statistical modelling by Gumbel and Jenkinson 1970 s: Classic limit laws generalized by Pickands
11 History Motivation 1980 s: Leadbetter (and others) extend theory to stationary processes 1990 s: Multivariate and other techniques explored as a means to improve inference 2000 s: Interest in spatial and spatio-temporal applications, and in finance
12 Probability Framework Let X 1,, X n iid F and define Mn = max{x 1,, X n } Then the distribution function of M n is P(M n x) = P(X 1 x,, X n x) = P(X 1 x) P(X n x) = F n (x) Remark F n (x) n = { 0 if F(x) < 1 1 if F(x) = 1
13 Classical Limit Laws Motivation Recall the Central Limit Theorem: X n µ σ n d N(0, 1) rescaling is the key to obtain a non-degenerate distribution Question: Can we get the limiting distribution of M n b n a n for suitable sequence {a n } > 0 and {b n }?
14 Classical Limit Laws Motivation Recall the Central Limit Theorem: X n µ σ n d N(0, 1) rescaling is the key to obtain a non-degenerate distribution Question: Can we get the limiting distribution of M n b n a n for suitable sequence {a n } > 0 and {b n }?
15 Theorem (Fisher Tippett Gnedenko theorem) If there exist sequences of constants a n > 0 and b n such that, as n P( M n b n x) d G(x) a n for some non-degenerate distribution G, then G belongs to either the Gumbel, the Fréchet or the Weibull family Gumbel : G(x) = exp(exp( x)) < x < ; { 0 x 0, Fréchet : G(x) = exp( x α ) x > 0, α > 0; { exp( ( x) Weibull : G(x) = α ) x < 0, α > 0, 1 x 0;
16 Example: Exponential Maxima X Exp(1) F n (x + log n) = (1 exp( x log n)) n = (1 1 n exp( x))n n exp{ exp( x)}
17 Generalized Extreme Value Distribution (GEV) This family encompasses all three extreme value limit families: where x + = max(x, 0) { [ G(x) = exp 1 + ξ( x µ σ ] 1 } ) ξ + µ and σ are location and scale parameters ξ is a shape parameter determining the rate of tail decay, with ξ > 0 giving the heavy-tailed (Fréchet) case ξ = 0 giving the light-tailed (Gumbel) case ξ < 0 giving the bounded-tailed (Weibull) case
18 Quantiles and Return Levels In terms of quantiles, take 0 < p < 1 and define x p such that: { [ G(x p ) = exp 1 + ξ( x p µ ] 1 } ξ ) = 1 p σ + x p = µ σ [ 1 { log(1 p) ξ }] ξ In the extreme value terminology, x p is the return level associated with the return period 1 p
19 Max-Stability Motivation Definition A distribution G is said to be max-stable if G k (a k x + b k ) = G(x), k N for some constants a k > 0 and b k Taking powers of G results only in a change of location and scale A distribution is max-stable it is a GEV distribution
20 Some Remarks Motivation There has been some work on the convergence rate of M n to the limiting regime, which depends on the underling distribution For statisticians, we use the GEV as an approximate distribution for sample maximal for "finite" n assess the fit empirically Direct use the GEV rather than three types separately
21 Block Maximum Approach Determine the block size and compute maximal for blocks usually annual maximal Fit the GEV to the maximal and assess fit usually via likelihood based techniques Perform inference for return levels, probabilities, etc
22 Diagnostics Motivation
23 Point Process Approach Motivation The block maximum method ignores much of the data which may also relevant to extreme we would like to use the data more efficient. Alternatives: peaks over thresholds r-largest order statistics Both are special cases of a point process representation, under which we approximate the exceedances over a threshold by a two-dimensional Poisson process
24
25 Point Process Limit: Basic Idea iid Suppose X 1,, X n F and { M n b n a n } converges to GEV distribution Construct a sequence of point processes on R 2 by { i P n = ( n + 1, X } i b n ) : i = 1,, n a n P n n P, where P is a Poisson process
26 Formally Linking Extremes and Point Processes Given: P( Mn bn a n { } x) exp (1 + ξx) 1 ξ log P( M n b n a n log P n ( X b n a n x) (1 + ξx) 1 ξ x) (1 + ξx) 1 ξ n log(1 P( X b n a n > x)) (1 + ξx) 1 ξ np( X b n a n > x) (1 + ξx) 1 ξ
27 Extremes and Point Processes Cond t np( X b n a n > x) (1 + ξx) 1 ξ Therefore, creating a series of point processes { } Xi b n : i = 1,, n n a n these will converge to an inhomogeneous Poisson process with measure ν((x, )) = (1 + ξx) 1 ξ for sets bounded away from zero (exceedance example).
28 Generalized Pareto Distribution (GPD) for Exceedances P(X i > x + u X i > u) = np(x i > x + u) np(x i > u) ( 1 + ξ x+u b n ) 1 ξ a n 1 + ξ x bn a n = ( 1 + ) 1 ξx ξ a n + ξ(u b n ) Survival function of generalized Pareto distribution
29 Theorem (Pickands Balkema de Haan theorem) Let X 1, iid F, and let F u be their conditional excess distribution function. Pickands (1975), Balkema and de Haan (1974) posed that for a large class of underlying distribution functions F, and large u, F u is well approximated by the generalized Pareto distribution GPD. That is: F u (y) GPD ξ,σ(u),u (y) u where GPD ξ,σ(u),u (y) = { 1 (1 + ξy σ(u) ) 1 ξ ξ 0, 1 exp( y σ(u) ) ξ = 0;
30 Fort Collins Data: Block Maximum vs. Threshold-exceedance GPD has lower standard error for ξ, lower estimate as well. GPD has narrower confidence interval. Have not yet discussed threshold selection procedure.
31 Threshold Selection Motivation Bias variance trade off: threshold too low bias because of the model asymptotics being invalid; threshold too high variance is large due to few data points Figure: Mean residual life plot(mrl): MRL is linear when GPD holds
32 Temporal Dependence Question: Is the GEV still the limiting distribution for block maxima of a stationary (but not independent) sequence {X i }? Answer: Yes, so long as mixing conditions hold. (Leadbetter et al., 1983) What does this mean for inference? Block maximum approach: GEV still correct for marginal. Since block maximum data likely have negligible dependence, proceed as usual Threshold exceedance approach: GPD is correct for the marginal. If extremes occur in clusters, estimation affected as likelihood assumes independence of threshold exceedances
33 Temporal Dependence Question: Is the GEV still the limiting distribution for block maxima of a stationary (but not independent) sequence {X i }? Answer: Yes, so long as mixing conditions hold. (Leadbetter et al., 1983) What does this mean for inference? Block maximum approach: GEV still correct for marginal. Since block maximum data likely have negligible dependence, proceed as usual Threshold exceedance approach: GPD is correct for the marginal. If extremes occur in clusters, estimation affected as likelihood assumes independence of threshold exceedances
34 Temporal Dependence Question: Is the GEV still the limiting distribution for block maxima of a stationary (but not independent) sequence {X i }? Answer: Yes, so long as mixing conditions hold. (Leadbetter et al., 1983) What does this mean for inference? Block maximum approach: GEV still correct for marginal. Since block maximum data likely have negligible dependence, proceed as usual Threshold exceedance approach: GPD is correct for the marginal. If extremes occur in clusters, estimation affected as likelihood assumes independence of threshold exceedances
35 Temporal Dependence Question: Is the GEV still the limiting distribution for block maxima of a stationary (but not independent) sequence {X i }? Answer: Yes, so long as mixing conditions hold. (Leadbetter et al., 1983) What does this mean for inference? Block maximum approach: GEV still correct for marginal. Since block maximum data likely have negligible dependence, proceed as usual Threshold exceedance approach: GPD is correct for the marginal. If extremes occur in clusters, estimation affected as likelihood assumes independence of threshold exceedances
36 Remarks on To estimate the tail, EVT uses only extreme observations Tail parameter ξ is extremely important but hard to estimate GPD is the limiting distribution for threshold exceedances Threshold exceedance approaches allow the user to retain more data than block-maximum approaches, thereby reducing the uncertainty with parameter estimates Temporal dependence in the data is more of an issue in threshold exceedance models. One can either decluster, or alternatively, adjust inference
37 Remarks on To estimate the tail, EVT uses only extreme observations Tail parameter ξ is extremely important but hard to estimate GPD is the limiting distribution for threshold exceedances Threshold exceedance approaches allow the user to retain more data than block-maximum approaches, thereby reducing the uncertainty with parameter estimates Temporal dependence in the data is more of an issue in threshold exceedance models. One can either decluster, or alternatively, adjust inference
38 Remarks on To estimate the tail, EVT uses only extreme observations Tail parameter ξ is extremely important but hard to estimate GPD is the limiting distribution for threshold exceedances Threshold exceedance approaches allow the user to retain more data than block-maximum approaches, thereby reducing the uncertainty with parameter estimates Temporal dependence in the data is more of an issue in threshold exceedance models. One can either decluster, or alternatively, adjust inference
39 Remarks on To estimate the tail, EVT uses only extreme observations Tail parameter ξ is extremely important but hard to estimate GPD is the limiting distribution for threshold exceedances Threshold exceedance approaches allow the user to retain more data than block-maximum approaches, thereby reducing the uncertainty with parameter estimates Temporal dependence in the data is more of an issue in threshold exceedance models. One can either decluster, or alternatively, adjust inference
40 Remarks on To estimate the tail, EVT uses only extreme observations Tail parameter ξ is extremely important but hard to estimate GPD is the limiting distribution for threshold exceedances Threshold exceedance approaches allow the user to retain more data than block-maximum approaches, thereby reducing the uncertainty with parameter estimates Temporal dependence in the data is more of an issue in threshold exceedance models. One can either decluster, or alternatively, adjust inference
41 Examples A central aim of multivariate extremes is trying to find an appropriate structure to describe tail dependence
42 What is a Multivariate Extreme? Let Z m = (Z m,1,, Z m,d ) T, m N be an iid sequence of random vectors. Want to extract a subset of data considered "extreme" Block-maximum: Mn = ( n m=1 Z m,1,, n Leads to modeling with multivariate max-stable distributions m=1 Z m,d) T Marginal-exceedance: For each marginal i = 1,, d, find an appropriate threshold u i and retain data where Z m,i > u i. Leads to multivariate generalized Pareto distribution Norm-exceedance: For a given norm retain data where Z m > z. Leads to description by multivariate regular variation
43 What is a Multivariate Extreme?
44 Multivariate Models Motivation Let X 1, X 2, iid d-dimensional random vectors with distribution F Our interest is in the (non degenerate) limiting distribution { } Xi b i n P max x G(x) i=1,,n a n for some sequences a n > 0 and b n R d. G is called a multivariate extreme value distribution Transform to common marginals Z with unit Fréchet i.e., Z (x,,, ) = = Z (,,, x) = exp( 1 x ) to model the dependence structure
45 Multivariate Models cond t P(Z 1 { z 1,, Z d z} d ) = exp V (z 1,, z d ) z 1,, z d > 0 The exponent measure V (z 1,, z d ) satisfies V (tz 1,, tz d ) = t 1 V (z 1,, z d ) = V is homogeneous of order -1 V (z 1,, z d ) = { 1 z z D 1 min(z 1,,z D ) if Z 1,, Z d are independent if Z 1,, Z d are entirely dependent
46 Spectral Representations Pickands, 1981 V (z 1,, z d ) = S D max( w 1 z 1,, w D z D ) dm(w 1,, w D ) where M is a measure on the D-dimensional simplex S D wd dm(w 1,, w D ) = 1 for each d No simple parametric forms for V
47 measure of extremal dependence { P(Z 1 z,, Z D z) = exp V (1,,1) z } exp( θ D z ) z > 0 θ D is the extrmal coefficient θ D = 1 fully dependent θ D = D independent In bivariate case lim z P(Z 2 > z Z 1 > z) = 2 θ D When D = 2 madogram [Cooley, Naveau and Poncet, 2006] µ F = 1 2 E{ F(Z 1 F(Z 2 )) } θ 2 = 1 + 2ν F 1 2ν F provide a good estimator of θ 2
48 Motivation Developed originally to predict probability distributions of ore grades for mining operations is based largely on the theory of Gaussian random processes Y (x) = µ(x) + e(x) + ε(x), x D R 2 E(Y (x)) = µ(x), Cov(Y (x), Y (y)) = C( x y ) The main propose of geostatistics is interpolation
49 Figure: A realization of Gaussian random processes
50 Bayesian Hierarchical Models Copula Models Max-stable Models Three-Level Spatial Hierarchical Model Data level: Likelihood which characterizes the distribution of the observed data given the parameters at the process level Y i (x d ) {µ(x d ), σ(x d ), ξ(x d )} GEV {µ(x d ), σ(x d ), ξ(x d )} i = 1,, n, d = 1,, D Process level: Latent process captured by spatial model for the data level parameters Prior level: Prior distributions put on the parameters
51 Bayesian Hierarchical Models Copula Models Max-stable Models Three-Level Spatial Hierarchical Model Suppose the response variables {Y (x)} are independent conditionally on an unobserved latent process {S(x)} Assume the {S(x)} follows a Gaussian process and induce spatial dependence in {Y (x)} by integration over the latent process Common approach in geostatistics with non-normal response [Diggle et al. 1998,2007]. Usually performed in a Bayesian setting using Markov chain Monte Carlo (MCMC)
52 Pros and Cons Motivation Bayesian Hierarchical Models Copula Models Max-stable Models The quantile surfaces are realistic After averaging over S(x) the marginal distribution of {Y (x)} is NOT GEV The spatial dependence is ignored because of conditional independence
53 Pros and Cons Motivation Bayesian Hierarchical Models Copula Models Max-stable Models Figure: One realisation of the latent variable model, showing the lack of local spatial structure
54 Copula Motivation Bayesian Hierarchical Models Copula Models Max-stable Models The D dimensional joint distribution F of any random variable (Y 1,, Y D ) may be written as F (y 1,, y n ) = C{F 1 (y 1 ),, F D (y D )} Where F 1,, F D are univariate marginal distributions of Y 1,, Y D and C is a copula The copula is uniquely determined for distributions F with absolutely continuous marginal With continuous and strictly increasing marginals, the copula corresponds to F(y 1,, y n ) = C{F 1 (y 1 ),, F D (y D )} is as follows: { } C(u 1,, u D ) = F F 1 1 (u 1 ),, F 1 D (u D )
55 Copula Motivation Bayesian Hierarchical Models Copula Models Max-stable Models The D dimensional joint distribution F of any random variable (Y 1,, Y D ) may be written as F (y 1,, y n ) = C{F 1 (y 1 ),, F D (y D )} Where F 1,, F D are univariate marginal distributions of Y 1,, Y D and C is a copula The copula is uniquely determined for distributions F with absolutely continuous marginal With continuous and strictly increasing marginals, the copula corresponds to F(y 1,, y n ) = C{F 1 (y 1 ),, F D (y D )} is as follows: { } C(u 1,, u D ) = F F 1 1 (u 1 ),, F 1 D (u D )
56 Copula Motivation Bayesian Hierarchical Models Copula Models Max-stable Models The D dimensional joint distribution F of any random variable (Y 1,, Y D ) may be written as F (y 1,, y n ) = C{F 1 (y 1 ),, F D (y D )} Where F 1,, F D are univariate marginal distributions of Y 1,, Y D and C is a copula The copula is uniquely determined for distributions F with absolutely continuous marginal With continuous and strictly increasing marginals, the copula corresponds to F(y 1,, y n ) = C{F 1 (y 1 ),, F D (y D )} is as follows: { } C(u 1,, u D ) = F F 1 1 (u 1 ),, F 1 D (u D )
57 Example Motivation Bayesian Hierarchical Models Copula Models Max-stable Models Gaussian copula: Y 1,, Y D N d (0, Ω) { C(u 1,, u d ) = Φ Φ 1 (u 1 ),, Φ 1 (u d ) : Ω} Student t copula: { C(u 1,, u d ) = T ν Extremal Copula: Tν 1 (u 1 ),, Tν 1 (u D ) : Ω} Y 1,, Y D multivariate GEV By max-stability C(u m 1,, um D ) = Cm (u 1,, u D ), 0 < u 1,, u d < 1, m N
58 Pros and Cons Motivation Bayesian Hierarchical Models Copula Models Max-stable Models On the "copula scale : the dependence seems more or less OK But this is no longer true at the original, i.e., extremal scale. Figure: One simulation from the fitted Gaussian copula model
59 Max-stable Processes Bayesian Hierarchical Models Copula Models Max-stable Models Definition Let Y m (x), x D R 2, m = 1,, n be independent copies of Y (x), and let M n (x) = max Y m (x). Y (x) is termed max-stable if there exist {a n (x)} and {b n (x)} such that P( M n(x) b n (x) a n (x) y(x)) n P(G(x) y(x))
60 Spectral Representations Bayesian Hierarchical Models Copula Models Max-stable Models Theorem (Schlather, 2002) Z (x), x D R 2 is max-stable with unit Fréchet marginals There exist iid positive stochastic processes V 1 (x), V 2 (x), with E[V i (x)] = 1 x D and E[sup x D V (x)] < and an independent point process {s i } i=1 on R + with intensity measure r 2 dr such that Z (x) d = max i N s i V i
61 Some Models Motivation Bayesian Hierarchical Models Copula Models Max-stable Models Smith: Y i (x) = φ(x i U i )(), {U i } i 1 points of a homogeneous Poisson process on R 2 Schlather: Y i (x) = 2πε i (x), ε i ( ) standard Gaussian process Geometric: Y i (x) = exp{σε i (x) σ2 2 } Brown Res: Y i (x) = exp{ε i (x) γ(x)}, ε i ( ) intrinsically stationary Gaussian process with (semi) variogram γ Figure: One realization of max-stable processes
62 Inference Motivation Bayesian Hierarchical Models Copula Models Max-stable Models For the max-stable process models, only the bivariate distributions are known Composite likelihoods [Lindsay, 1988] are used to obtain estimation Since we have the bivariate distributions we will use the pair wise likelihood l p (θ, y) = n K 1 m=1 i=1 j=i+1 k log f (ym, i fm; j θ) Not a true likelihood Over uses the data each observation appears K-1 times
63 Pros and Cons Motivation Bayesian Hierarchical Models Copula Models Max-stable Models Justified by extreme value theory Able to describe asymptotic dependence Describe everything that is asymptotically independent as exactly independent Only suitable for observations that are annual maximal at this point
64 Motivation Although infrequent, extremes have large human impact Uses only the "extreme" observations to fit the model No role for normal distribution! Spatial extremes modeling is challenge both in theoretical and computational aspects
65
Overview 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 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 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 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 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 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 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 informationA Conditional Approach to Modeling Multivariate Extremes
A Approach to ing Multivariate Extremes By Heffernan & Tawn Department of Statistics Purdue University s April 30, 2014 Outline s s Multivariate Extremes s A central aim of multivariate extremes is trying
More informationBayesian 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 informationStatistics of Extremes
Statistics of Extremes Anthony Davison c 211 http://stat.epfl.ch Multivariate Extremes 19 Componentwise maxima.................................................. 194 Standardization........................................................
More informationLecture Quantitative Finance Spring Term 2015
on bivariate Lecture Quantitative Finance Spring Term 2015 Prof. Dr. Erich Walter Farkas Lecture 07: April 2, 2015 1 / 54 Outline on bivariate 1 2 bivariate 3 Distribution 4 5 6 7 8 Comments and conclusions
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 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 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 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 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 informationSTATISTICAL MODELS FOR QUANTIFYING THE SPATIAL DISTRIBUTION OF SEASONALLY DERIVED OZONE STANDARDS
STATISTICAL MODELS FOR QUANTIFYING THE SPATIAL DISTRIBUTION OF SEASONALLY DERIVED OZONE STANDARDS Eric Gilleland Douglas Nychka Geophysical Statistics Project National Center for Atmospheric Research Supported
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 informationOverview of Extreme Value Analysis (EVA)
Overview of Extreme Value Analysis (EVA) Brian Reich North Carolina State University July 26, 2016 Rossbypalooza Chicago, IL Brian Reich Overview of Extreme Value Analysis (EVA) 1 / 24 Importance of extremes
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 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 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 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 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 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 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 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 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 informationLikelihood inference for threshold exceedances
31 Spatial Extremes Anthony Davison École Polytechnique Fédérale de Lausanne, Switzerland Raphaël Huser King Abdullah University, Thuwal, Saudi Arabia Emeric Thibaud École Polytechnique Fédérale de Lausanne,
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 informationMax-stable processes: Theory and Inference
1/39 Max-stable processes: Theory and Inference Mathieu Ribatet Institute of Mathematics, EPFL Laboratory of Environmental Fluid Mechanics, EPFL joint work with Anthony Davison and Simone Padoan 2/39 Outline
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 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 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 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 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 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 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 informationExtremes and Atmospheric Data
Extremes and Atmospheric Data Eric Gilleland Research Applications Laboratory National Center for Atmospheric Research 2007-08 Program on Risk Analysis, Extreme Events and Decision Theory, opening workshop
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 informationHigh-frequency data modelling using Hawkes processes
Valérie Chavez-Demoulin joint work with High-frequency A.C. Davison data modelling and using A.J. Hawkes McNeil processes(2005), J.A EVT2013 McGill 1 /(201 High-frequency data modelling using Hawkes processes
More informationHigh-frequency data modelling using Hawkes processes
High-frequency data modelling using Hawkes processes Valérie Chavez-Demoulin 1 joint work J.A McGill 1 Faculty of Business and Economics, University of Lausanne, Switzerland Boulder, April 2016 Boulder,
More informationBayesian Spatial Modeling of Extreme Precipitation Return Levels
Bayesian Spatial Modeling of Extreme Precipitation Return Levels Daniel Cooley 1,2 Doug Nychka 2 Philippe Naveau 3,4 1 Department of Statistics Colorado State University, Fort Collins, CO 2 Geophysical
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 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 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 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 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 informationVaR vs. Expected Shortfall
VaR vs. Expected Shortfall Risk Measures under Solvency II Dietmar Pfeifer (2004) Risk measures and premium principles a comparison VaR vs. Expected Shortfall Dependence and its implications for risk measures
More informationPeaks-Over-Threshold Modelling of Environmental Data
U.U.D.M. Project Report 2014:33 Peaks-Over-Threshold Modelling of Environmental Data Esther Bommier Examensarbete i matematik, 30 hp Handledare och examinator: Jesper Rydén September 2014 Department of
More informationEfficient Estimation of Distributional Tail Shape and the Extremal Index with Applications to Risk Management
Journal of Mathematical Finance, 2016, 6, 626-659 http://www.scirp.org/journal/jmf ISSN Online: 2162-2442 ISSN Print: 2162-2434 Efficient Estimation of Distributional Tail Shape and the Extremal Index
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 informationWei-han Liu Department of Banking and Finance Tamkang University. R/Finance 2009 Conference 1
Detecting Structural Breaks in Tail Behavior -From the Perspective of Fitting the Generalized Pareto Distribution Wei-han Liu Department of Banking and Finance Tamkang University R/Finance 2009 Conference
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 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 informationModelling Dependence with Copulas and Applications to Risk Management. Filip Lindskog, RiskLab, ETH Zürich
Modelling Dependence with Copulas and Applications to Risk Management Filip Lindskog, RiskLab, ETH Zürich 02-07-2000 Home page: http://www.math.ethz.ch/ lindskog E-mail: lindskog@math.ethz.ch RiskLab:
More informationSTA 2201/442 Assignment 2
STA 2201/442 Assignment 2 1. This is about how to simulate from a continuous univariate distribution. Let the random variable X have a continuous distribution with density f X (x) and cumulative distribution
More 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 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 informationStat 5101 Lecture Notes
Stat 5101 Lecture Notes Charles J. Geyer Copyright 1998, 1999, 2000, 2001 by Charles J. Geyer May 7, 2001 ii Stat 5101 (Geyer) Course Notes Contents 1 Random Variables and Change of Variables 1 1.1 Random
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 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 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 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 informationLosses Given Default in the Presence of Extreme Risks
Losses Given Default in the Presence of Extreme Risks Qihe Tang [a] and Zhongyi Yuan [b] [a] Department of Statistics and Actuarial Science University of Iowa [b] Smeal College of Business Pennsylvania
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 informationExtreme Event Modelling
Extreme Event Modelling Liwei Wu, SID: 52208712 Department of Mathematics City University of Hong Kong Supervisor: Dr. Xiang Zhou March 31, 2014 Contents 1 Introduction 4 2 Theory and Methods 5 2.1 Asymptotic
More informationModelling spatial extremes using max-stable processes
Chapter 1 Modelling spatial extremes using max-stable processes Abstract This chapter introduces the theory related to the statistical modelling of spatial extremes. The presented modelling strategy heavily
More informationA spatio-temporal model for extreme precipitation simulated by a climate model
A spatio-temporal model for extreme precipitation simulated by a climate model Jonathan Jalbert Postdoctoral fellow at McGill University, Montréal Anne-Catherine Favre, Claude Bélisle and Jean-François
More informationLifetime Dependence Modelling using a Generalized Multivariate Pareto Distribution
Lifetime Dependence Modelling using a Generalized Multivariate Pareto Distribution Daniel Alai Zinoviy Landsman Centre of Excellence in Population Ageing Research (CEPAR) School of Mathematics, Statistics
More informationMarkov Chain Monte Carlo methods
Markov Chain Monte Carlo methods Tomas McKelvey and Lennart Svensson Signal Processing Group Department of Signals and Systems Chalmers University of Technology, Sweden November 26, 2012 Today s learning
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 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 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 informationAccommodating measurement scale uncertainty in extreme value analysis of. northern North Sea storm severity
Introduction Model Analysis Conclusions Accommodating measurement scale uncertainty in extreme value analysis of northern North Sea storm severity Yanyun Wu, David Randell, Daniel Reeve Philip Jonathan,
More informationBayesian Methods for Machine Learning
Bayesian Methods for Machine Learning CS 584: Big Data Analytics Material adapted from Radford Neal s tutorial (http://ftp.cs.utoronto.ca/pub/radford/bayes-tut.pdf), Zoubin Ghahramni (http://hunch.net/~coms-4771/zoubin_ghahramani_bayesian_learning.pdf),
More informationQuantitative Modeling of Operational Risk: Between g-and-h and EVT
: Between g-and-h and EVT Paul Embrechts Matthias Degen Dominik Lambrigger ETH Zurich (www.math.ethz.ch/ embrechts) Outline Basel II LDA g-and-h Aggregation Conclusion and References What is Basel II?
More informationAnalysis of Bivariate Extreme Values
Analysis of Bivariate Extreme Values Chris Egeland Busuttil Master of Science in Physics and Mathematics Submission date: February 2015 Supervisor: Arvid Næss, MATH Norwegian University of Science and
More informationSTA205 Probability: Week 8 R. Wolpert
INFINITE COIN-TOSS AND THE LAWS OF LARGE NUMBERS The traditional interpretation of the probability of an event E is its asymptotic frequency: the limit as n of the fraction of n repeated, similar, and
More informationTail negative dependence and its applications for aggregate loss modeling
Tail negative dependence and its applications for aggregate loss modeling Lei Hua Division of Statistics Oct 20, 2014, ISU L. Hua (NIU) 1/35 1 Motivation 2 Tail order Elliptical copula Extreme value copula
More 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 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 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 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 informationSYSM 6303: Quantitative Introduction to Risk and Uncertainty in Business Lecture 4: Fitting Data to Distributions
SYSM 6303: Quantitative Introduction to Risk and Uncertainty in Business Lecture 4: Fitting Data to Distributions M. Vidyasagar Cecil & Ida Green Chair The University of Texas at Dallas Email: M.Vidyasagar@utdallas.edu
More information