Sharp statistical tools Statistics for extremes
|
|
- Alicia Horn
- 5 years ago
- Views:
Transcription
1 Sharp statistical tools Statistics for extremes Georg Lindgren Lund University October 18, 2012 SARMA
2 Background Motivation We want to predict outside the range of observations Sums, averages and proportions Normal distribution Extremes Central limit theorem Successes in large experiments Normal distribution inappropriate Bulk of data may be misleading Extremes are rare Sparsity of data may be compensated by proper (extreme value) theory
3 Background Overview of contents Topics Extreme value distributions Block maxima POT - Peaks Over Threshold Estimation, return period, uncertainties Extremes with cyclic or linear trend Extremes with other covariates CO2, NAO, AO,... Useful statistics packages The R environment - need packages extremes and evd for the computer experiments The WAFO package in Matlab and Python (Downloadable from code.google.com)
4 Background A typical example: annual maximum of daily precipitation 100 years of yearly maximum of daily precipitation in Fort Collins precipitation (in)
5 Background A second typical example Monthly maximum precipitation values has 12 times as many data but season Precipitation (in)
6 Background Some basic facts about statistical extremes Notation: M n = max(x 1, X 2,..., X n ) is the maximum of n observations of a random quantity. Example: X k = maximum precipitation year k Problem: Find the probability Prob(M n > x) for large n and reasonable x. In particular when x > largest observed data!
7 Background Return period return level X 1, X 2,... a sequence of random quantities with common distribution, e.g. one observation per year (a mean, a yearly maximum,..., anything) 100 year return value x 100 is the value that is exceeded on average once in 100 years. The same as Prob(X 1 > x 100 ) = 1/100 ( Prob(M 100 > x 100 ) = ) /e = The probability that the 100-year return level is exceeded at least one time during a 100 year period is 63%. The 1000-year value exceeded with the same probability at least once during 1000 years, etc. IF YEARS ARE INDEPENDENT OF EACH OTHER. USUALLY,
8 Background Wave data from the Pacific Ocean 582 measurements of H s = Significant wave-height at buoy in the Pacific. About 14 data per month, Dec-Feb, during 7 years. Hs t
9 Background Histogram of Pacific wave height Histogram of Hs Hs F(x) CDF of Hs exceedances over 7m CDF of exponential distribution +7m x Conclusion: the tail of the distribution has a more standardized distribution than the central part.
10 GEV and GPD Goals for statistical extreme value analysis One-dimensional: From a series of observations (often in time) estimate the probability of high values, of the order of the largest observation; OR HIGHER, even MUCH HIGHER Make a statement about how uncertain the prediction is Identify possible non-stationarity in extremes can be different from non-stationarity in averages or standard deviation Find covariates that explain the occurrence of extreme values Utilize data as much as possible balance between number of data available and bias Multi-dimensional: find probabilities of combinations of extreme (or not so extreme) values in two or more series (research is going on) Extremes in random fields (even more research is going on)
11 GEV and GPD The GEV and the GPD Block maxima: Take one extreme value per time unit (e.g. day, month, year,...). Ideal situation: stationary conditions, no trend. Exceedance analysis: Take all values over a certain threshold gives more data that are more representative for the global maximum that the bulk data. Statistical extreme value analysis relies on two mathematical facts Block maxima: The maximum of many observations has a GEV distribution. Exceedance analysis: The exceedances have a GPD distribution.
12 GEV and GPD The Extremal types theorem An almost true theorem: When n is large, the distribution of the maximum M n is a GEV = Generalized Extreme Value distribution: { ( Prob(M n x) exp 1 + ξ x µ ) } 1/ξ ψ where ξ = shape (type), ψ = scale, µ = location ξ < 0 = an upper limit exists, ξ > 0 = a lower limit exists, ξ = 0 = Gumbel distribution, no limits exist +
13 GEV and GPD Light, Gumbel, or Heavy tails Light tail ξ < 0 log( log(f)) Gumbel Probability Plot X Gumbel Probability Plot 5 10 Double tail ξ = Heavy 200 tail ξ > log( log(f)) log( log(f)) X Gumbel Probability Plot X
14 GEV and GPD The Tail theorem Another almost true theorem: Almost all statistical distribution have a tail that is GPD they have a Generalized Pareto Distribution. Take only observations that above a fixed threshold u and let Y = X u be the exceedance. ( Prob(Y y) ξ y ) 1/ξ σ + where ξ = shape (type), σ = scale, (µ = location = 0). ξ = 0 = exponential tail, ξ > 0 = heavy tail, ξ < 0 = limited tail
15 GEV and GPD GPD-tail in normal distribution Normal tail is GPD with ξ = Normal distribution The tail > 2 of a normal distribution 0.8 F(x) Red = empirical cdf of exceedances over 2 Blue = estimated GPD x 1.8
16 Block maxima Block maxima with block = one year Many observations over a year (hourly, semi-daily, daily) take the maximum. Assume a GEV distribution for yearly maximum. Estimate parameters in GEV by Maximum likelihood method in R-package extremes (or WAFO) fit <- gev.fit(data)
17 Block maxima 100 years of precipitation in Fort Collins, CO Fort Collins, CO, daily precipitation: odie/rain.html Time series of daily precipitation, Strong annual cycle wet in late spring, dry in winter No long-term trend Recent flood, July 28, 1997
18 Block maxima Annual daily maximum precipitation, Fort Collins precipitation (in) year
19 Block maxima Seasons can be handled by making separate analyses for each month Precipitation (in)
20 Block maxima Some R-code data <- read.table(file="ftcprec.prn", header=true) plot(data$mn,data$prec/100,ylab="precipitation (in)", xlab="month") YearlyMax = aggregate(prec Year,data=data,max) YearlyMax[,2]= YearlyMax[,2]/100 plot(yearlymax,type= l,col= blue ) library(extremes) ftcanmax <- read.table(file="ftcanmax.prn", header=true) fit <- gev.fit(ftcanmax$prec/100) gev.diag(fit)
21 Block maxima Estimated parameter values + standard errror Parameter Estimate Standard error Location µ Scale ψ Shape ξ
22 Block maxima Is the GEV a Gumbel distribution? Is the shape parameter ξ = > 0 really significantly different from 0? Can be tested by a likelihood ratio test: Dev = 2 log max likelihood under restricted model (ξ = 0) max likelihood under full model If restricted model is correct, this has a chi-squared distribution; d.f. = # parameters in full model - # parameters in restricted model Test by extremes (ξ = 0 = Gumbel): fit0 <- gum.fit(ftcanmax$prec/100) Dev <- 2*(fit0$nllh - fit$nllh) pval = pchisq( Dev, 1, lower.tail=f) (= 0.038)
23 Block maxima A 95% confidence interval for shape estimate is (0.009, 0.369) gev.profxi(fit, fit$mle[3] -0.25, fit$mle[3] +0.25) Profile Log likelihood
24 Block maxima Diagnostics and conclusions We want return values for N = 10, 100, 100, 1000 years in GEV: x N = µ ψ {1 ( ln(1 1/N)) ξ} ξ Estimated standard errors gives confidence limits for return levels. gev.diag(fit) An example for the exercises; fit.rl <- return.level(fit,rperiods=c(10,100), conf = 0.05)
25 Block maxima GEV summary Probability Plot Quantile Plot Model Empirical Empirical Model Return Level Plot Density Plot Return Level f(z) e 01 1e+00 1e+01 1e+02 1e+03 Return Period z
26 Block maxima Trends in extremes Water level in the Japan Sea 20 Water level 15 (m) Maximum 5 min water level (m) Time (h)
27 Block maxima Alt I: Normalize (simplistic) Take average and standard deviation for each 5-minute period, subtract and divide. Take maximum over 5-minute period, and fit a GPD. Diagnostic plot: 1 Probability plot Density plot F(x) x Residual Quantile Plot x Residual Probability Plot Model (gev) Model (gev) Empirical Fit method: PWM, Fit p value: Empirical
28 Block maxima Alt II: fit with trend yura5fit<-gev.fit(y5max[,1]) Est.: S.e.: yura5fit.mul<-gev.fit(y5max[,1], ydat=as.matrix(1:length(y5max[,1])),mul=1) Est S.e.: yura5fit.sigl<-gev.fit(y5max[,1], ydat=as.matrix(1:length(y5max[,1])),sigl=1) Est.: S.E.:
29 Exceedance analysis The dilemma of statistical extremes We want to predict events that have never been observed! From 20 years data can one say something about the 100 year return value?
30 Exceedance analysis 20 years of monthly data 5 20 years of monthly data
31 Exceedance analysis Use more high values Waste of data to use only yearly maximum 20 years of data = 240 monthly data The smallest yearly maximum (year 7) is X 7 = There are 42 monthly values greater than 1.67 Can one use all 42? Or all those 48 greater than 1.5? Or the 84 values greater than 1?
32 Exceedance analysis High values are rare and occur randomly Take a reasonably high level u try many! Estimate the rate by which this level is exceeded per time unit (e.g. per year), λ = λ u : λ = Observed number of exceedances over u Total observation time u = 1.5 gives λ 1.5 = 48/20 = 2.4 Reasonable to assume that N = the number of exceedances over 1.5 any year has a Poisson distribution with mean λ = 2.4: Prob(N = k) = e λ λ k /k!
33 Exceedance analysis Generaliserad Pareto fördelning - GPD Exceedances over a high threshold are more representative for the global extremes than the buld data Almost all distributions have a Generalized Pareto-tail, GPD With Y = X u = exceedance over level u: ( Prob(Y y) ξ y ) 1/ξ σ + Exponential tail: ξ = 0; Heavy tail: ξ > 0; Limited tail: ξ < 0
34 Exceedance analysis GPD-tail in normal distribution Normal tail is GPD with ξ = Normal distribution The tail > 2 of a normal distribution 0.8 F(x) Red = empirical cdf of exceedances over 2 Blue = estimated GPD x 1.8
35 Exceedance analysis Poisson + GPD = GEV N = number of exceedances Y j = X j u over u is (approximately) Poisson distributed with expectation λ The size of exceedances Y 1,..., Y N, have (approximately) a GPD If M = yearly maximum = u + max(y 1,..., Y N ), then for x > u: Prob(M x) = Prob(N = 0) + =... = exp { ( λ 1 + ξ x u σ Prob(N = n, Y 1,..., Y n x u) n=1 ) 1/ξ + } (1)
36 Exceedance analysis Poisson + GPD = GEV, contd (1) is a GEV distribution Prob(M x) = exp Translation from Poisson+GPD to GEV: { ψ = σ λ ξ ( 1 + ξ x µ ψ µ = u + ψ σ ξ ) 1/ξ To get maximum over n years, replace λ with nλ in (1) + }
37 Exceedance analysis Choice of threshold How to choose u? Note: assume GPD above threshold u Diagnostic: A GPD has linear mean excess E(X u X > u) = σ + ξu 1 ξ Plot the mean of all excesses over u as function of u. Take u as the smallest threshold for which the curve to the right is linear The slope is ξ/(1 ξ) om ξ < 1.
38 Exceedance analysis Mean excess plot Plot of E(X u X > u) for 20 years of monthly data: 1.2 Mean exceedance over threshold
39 Exceedance analysis Diagnostics in GPD analysis A plot of mean excess can be had to interpret Alternative: Estimate a full GPD for different thresholds If the tail above u 0 is GPD then all exceedances over u > u 0 are also GPD with the same shape parameter ξ but with different scal parameter σ u = σ u0 + ξ (u u 0 ) Modified scale = σ u ξ u should be constant, independent of u, when GPD is appropriate
40 Exceedance analysis Estimated CDF for yearly maximum From 20 yearly maxima: c = 0.14, µ = 2.81, a = Empirical and GEV estimated cdf (PWM method) True CDF for yearly maximum F(x) CDF for estimated GEV x
41 Exceedance analysis Estimated CDF by POT method 84 exceedances over u = 1 and GPD estimate gives ξ = 0.04, b = 2.38, a = Tail probability 10 1 Tail by POT method True CDF 10 2 Tail by direct GEV estimation
42 Exceedance analysis Fatalities in English coal mines Time and number of fatalities
43 Exceedance analysis GEV? We try GEV on all data (not really motivated!): Empirical and GEV estimated cdf (PWM method) F(x) x
44 Exceedance analysis Consider only big accidents POT There are 25 accidents with > 100 dead. Fit GPD to data > % of these exceed 350. CDF for deaths > 100 and GPD F(x) x
45 Exceedance analysis Vargas, Venezuela, catastrophe 1999 In December 15, 1999, the Vargas province was hit by 410 mm rain.
46 Exceedance analysis Rain in Venezuela - GEV or Gumbel? Gumbel distribution is GEV with shape parameter ξ = 0. Based on rain statistics the maximum daily rain distribution was estimated with GEV: scale = 19.9, location = 49.2, shape = The shape parameter ξ is not significantly different from 0, so Gumbel could perhaps be used instead of GEV. Estimated year return value for daily rain is x = 249 mm. In December 1999 the Vargas province was hit by 410 mm rain during one day. From the full GEV the estimate of the year value had been x = 468 mm, much closer to the real value. A 95% one-sided confidence interval is x < 1030 mm.
47 Multivariate extremes General recipe for bivariate extremes Assume, for a sequence of days, X 1, X 2,... is the wind speed and Y 1, Y 2,... is the wave height measured at a North Sea platform. The platform is designed for extreme 2 waves. It is also designed for extreme 1.8 Extremes come togehter winds. What about the combined 1.6 effect of wind and waves (and current)? Make marginal EVA and transform 1 Extremes just add x and y to standard scale. Plot 0.8 pairs of transformed extremes to examine joint extreme behavior. Three Extremes come alone main types of dependence/nondependence:
48 Some literature Some literature Beirlant, J., Goegebeur, Y., Segers, J., Teugels, J. (2004), Statistics of Extremes: Theory and Applications. Wiley. Coles S. (2001), An Introduction to Statistical Modeling of Extreme Values. Springer-Verlag. Gilleland, E., Katz, R., Young, G. (Feb. 2012), Package extremes. Nelsen, R.B. (2006), An Introduction to Copulas. Wiley. The WAFO group (2011), WAFO a MATLAB toolbox for analysis of random waves and loads. Lund univ. and
Overview 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 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 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 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 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 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 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 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 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 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 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 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 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 informationOn the modelling of extreme droughts
Modelling and Management of Sustainable Basin-scale Water Resource Systems (Proceedings of a Boulder Symposium, July 1995). IAHS Publ. no. 231, 1995. 377 _ On the modelling of extreme droughts HENRIK MADSEN
More informationChange Point Analysis of Extreme Values
Change Point Analysis of Extreme Values TIES 2008 p. 1/? Change Point Analysis of Extreme Values Goedele Dierckx Economische Hogeschool Sint Aloysius, Brussels, Belgium e-mail: goedele.dierckx@hubrussel.be
More informationEstimation of Quantiles
9 Estimation of Quantiles The notion of quantiles was introduced in Section 3.2: recall that a quantile x α for an r.v. X is a constant such that P(X x α )=1 α. (9.1) In this chapter we examine quantiles
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 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 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 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 informationR&D Research Project: Scaling analysis of hydrometeorological time series data
R&D Research Project: Scaling analysis of hydrometeorological time series data Extreme Value Analysis considering Trends: Methodology and Application to Runoff Data of the River Danube Catchment M. Kallache,
More informationIT S TIME FOR AN UPDATE EXTREME WAVES AND DIRECTIONAL DISTRIBUTIONS ALONG THE NEW SOUTH WALES COASTLINE
IT S TIME FOR AN UPDATE EXTREME WAVES AND DIRECTIONAL DISTRIBUTIONS ALONG THE NEW SOUTH WALES COASTLINE M Glatz 1, M Fitzhenry 2, M Kulmar 1 1 Manly Hydraulics Laboratory, Department of Finance, Services
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 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 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 informationJournal of Environmental Statistics
jes Journal of Environmental Statistics February 2010, Volume 1, Issue 3. http://www.jenvstat.org Exponentiated Gumbel Distribution for Estimation of Return Levels of Significant Wave Height Klara Persson
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 informationANALYZING SEASONAL TO INTERANNUAL EXTREME WEATHER AND CLIMATE VARIABILITY WITH THE EXTREMES TOOLKIT. Eric Gilleland and Richard W.
P2.15 ANALYZING SEASONAL TO INTERANNUAL EXTREME WEATHER AND CLIMATE VARIABILITY WITH THE EXTREMES TOOLKIT Eric Gilleland and Richard W. Katz Research Applications Laboratory, National Center for Atmospheric
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 informationUsing statistical methods to analyse environmental extremes.
Using statistical methods to analyse environmental extremes. Emma Eastoe Department of Mathematics and Statistics Lancaster University December 16, 2008 Focus of talk Discuss statistical models used to
More informationModeling daily precipitation in Space and Time
Space and Time SWGen - Hydro Berlin 20 September 2017 temporal - dependence Outline temporal - dependence temporal - dependence Stochastic Weather Generator Stochastic Weather Generator (SWG) is a stochastic
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 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 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 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 informationSTT 843 Key to Homework 1 Spring 2018
STT 843 Key to Homework Spring 208 Due date: Feb 4, 208 42 (a Because σ = 2, σ 22 = and ρ 2 = 05, we have σ 2 = ρ 2 σ σ22 = 2/2 Then, the mean and covariance of the bivariate normal is µ = ( 0 2 and Σ
More informationAbstract: In this short note, I comment on the research of Pisarenko et al. (2014) regarding the
Comment on Pisarenko et al. Characterization of the Tail of the Distribution of Earthquake Magnitudes by Combining the GEV and GPD Descriptions of Extreme Value Theory Mathias Raschke Institution: freelancer
More 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 informationHydrologic Design under Nonstationarity
Hydrologic Design under Nonstationarity Jayantha Obeysekera ( Obey ), SFWMD Jose D. Salas, Colorado State University Hydroclimatology & Engineering Adaptaton (HYDEA) Subcommittee Meeting May 30, 2017,
More informationClimate predictability beyond traditional climate models
Climate predictability beyond traditional climate models Rasmus E. Benestad & Abdelkader Mezghani Rasmus.benestad@met.no More heavy rain events? More heavy rain events? Heavy precipitation events with
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 informationGENERALIZED LINEAR MODELING APPROACH TO STOCHASTIC WEATHER GENERATORS
GENERALIZED LINEAR MODELING APPROACH TO STOCHASTIC WEATHER GENERATORS Rick Katz Institute for Study of Society and Environment National Center for Atmospheric Research Boulder, CO USA Joint work with Eva
More informationStatistical Tools for Making Predictions. LCPC Training Week 5 October 2010 Dr Colin Caprani
www.colincaprani.com 1 Introduction About me: Did PhD in Bridge Traffic Loading (BTL) in UCD, 2001-2005 Concentrated on statistical analysis of BTL Continuing research in the BTL (and new areas) About
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 informationChange Point Analysis of Extreme Values
Change Point Analysis of Extreme Values DGVFM Stuttgart 27 APRIL 2012 p. 1/3 Change Point Analysis of Extreme Values Goedele Dierckx Economische Hogeschool Sint Aloysius, Brussels, Belgium Jef L. Teugels
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 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 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 informationPrecipitation Extremes in the Hawaiian Islands and Taiwan under a changing climate
Precipitation Extremes in the Hawaiian Islands and Taiwan under a changing climate Pao-Shin Chu Department of Atmospheric Sciences University of Hawaii-Manoa Y. Ruan, X. Zhao, D.J. Chen, and P.L. Lin December
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 informationModelling Risk on Losses due to Water Spillage for Hydro Power Generation. A Verster and DJ de Waal
Modelling Risk on Losses due to Water Spillage for Hydro Power Generation. A Verster and DJ de Waal Department Mathematical Statistics and Actuarial Science University of the Free State Bloemfontein ABSTRACT
More informationSpatial and temporal extremes of wildfire sizes in Portugal ( )
International Journal of Wildland Fire 2009, 18, 983 991. doi:10.1071/wf07044_ac Accessory publication Spatial and temporal extremes of wildfire sizes in Portugal (1984 2004) P. de Zea Bermudez A, J. Mendes
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 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 informationAPPROXIMATING THE GENERALIZED BURR-GAMMA WITH A GENERALIZED PARETO-TYPE OF DISTRIBUTION A. VERSTER AND D.J. DE WAAL ABSTRACT
APPROXIMATING THE GENERALIZED BURR-GAMMA WITH A GENERALIZED PARETO-TYPE OF DISTRIBUTION A. VERSTER AND D.J. DE WAAL ABSTRACT In this paper the Generalized Burr-Gamma (GBG) distribution is considered to
More 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 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 informationColorado State University, Fort Collins, CO Weather Station Monthly Summary Report
Colorado State University, Fort Collins, CO Weather Station Monthly Summary Report Month: December Year: 2016 Temperature: Mean T max was 39.7 F which is 3.1 below the 1981-2010 normal for the month. This
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 informationProbability Methods in Civil Engineering Prof. Dr. Rajib Maity Department of Civil Engineering Indian Institute of Technology, Kharagpur
Probability Methods in Civil Engineering Prof. Dr. Rajib Maity Department of Civil Engineering Indian Institute of Technology, Kharagpur Lecture No. # 33 Probability Models using Gamma and Extreme Value
More informationMAS2317/3317: Introduction to Bayesian Statistics
MAS2317/3317: Introduction to Bayesian Statistics Case Study 2: Bayesian Modelling of Extreme Rainfall Data Semester 2, 2014 2015 Motivation Over the last 30 years or so, interest in the use of statistical
More informationStatistical Assessment of Extreme Weather Phenomena Under Climate Change
Statistical Assessment of Extreme Weather Phenomena Under Climate Change NCAR Advanced Study Program Summer Colloquium 2011, 624 June Practice Sets for the R tutorial on EVA in R 1 Fitting the GEV to data
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 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 informationGENERALIZED LINEAR MODELING APPROACH TO STOCHASTIC WEATHER GENERATORS
GENERALIZED LINEAR MODELING APPROACH TO STOCHASTIC WEATHER GENERATORS Rick Katz Institute for Study of Society and Environment National Center for Atmospheric Research Boulder, CO USA Joint work with Eva
More informationUncertainties in extreme surge level estimates from observational records
Uncertainties in extreme surge level estimates from observational records By H.W. van den Brink, G.P. Können & J.D. Opsteegh Royal Netherlands Meteorological Institute, P.O. Box 21, 373 AE De Bilt, The
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 informationFreeway rear-end collision risk for Italian freeways. An extreme value theory approach
XXII SIDT National Scientific Seminar Politecnico di Bari 14 15 SETTEMBRE 2017 Freeway rear-end collision risk for Italian freeways. An extreme value theory approach Gregorio Gecchele Federico Orsini University
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 informationChange Point Analysis of Extreme Values
Change Point Analysis of Extreme Values Lisboa 2013 p. 1/3 Change Point Analysis of Extreme Values Goedele Dierckx Economische Hogeschool Sint Aloysius, Brussels, Belgium Jef L. Teugels Katholieke Universiteit
More informationMULTIDIMENSIONAL COVARIATE EFFECTS IN SPATIAL AND JOINT EXTREMES
MULTIDIMENSIONAL COVARIATE EFFECTS IN SPATIAL AND JOINT EXTREMES Philip Jonathan, Kevin Ewans, David Randell, Yanyun Wu philip.jonathan@shell.com www.lancs.ac.uk/ jonathan Wave Hindcasting & Forecasting
More 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 informationR.Garçon, F.Garavaglia, J.Gailhard, E.Paquet, F.Gottardi EDF-DTG
Homogeneous samples and reliability of probabilistic models : using an atmospheric circulation patterns sampling for a better estimation of extreme rainfall probability R.Garçon, F.Garavaglia, J.Gailhard,
More informationINDIAN INSTITUTE OF SCIENCE STOCHASTIC HYDROLOGY. Lecture -27 Course Instructor : Prof. P. P. MUJUMDAR Department of Civil Engg., IISc.
INDIAN INSTITUTE OF SCIENCE STOCHASTIC HYDROLOGY Lecture -27 Course Instructor : Prof. P. P. MUJUMDAR Department of Civil Engg., IISc. Summary of the previous lecture Frequency factors Normal distribution
More informationTime Series Analysis -- An Introduction -- AMS 586
Time Series Analysis -- An Introduction -- AMS 586 1 Objectives of time series analysis Data description Data interpretation Modeling Control Prediction & Forecasting 2 Time-Series Data Numerical data
More informationTowards a more physically based approach to Extreme Value Analysis in the climate system
N O A A E S R L P H Y S IC A L S C IE N C E S D IV IS IO N C IR E S Towards a more physically based approach to Extreme Value Analysis in the climate system Prashant Sardeshmukh Gil Compo Cecile Penland
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 informationFuture extreme precipitation events in the Southwestern US: climate change and natural modes of variability
Future extreme precipitation events in the Southwestern US: climate change and natural modes of variability Francina Dominguez Erick Rivera Fernandez Hsin-I Chang Christopher Castro AGU 2010 Fall Meeting
More informationAutomated, Efficient, and Practical Extreme Value Analysis with Environmental Applications
Automated, Efficient, and Practical Extreme Value Analysis with Environmental Applications arxiv:1611.08261v1 [stat.me] 24 Nov 2016 Brian M. Bader, Ph.D. University of Connecticut, 2016 ABSTRACT Although
More informationMODELING EXTREME VALUES OF PROCESSES OBSERVED AT IRREGULAR TIME STEPS: APPLICATION TO SIGNIFICANT WAVE HEIGHT
The Annals of Applied Statistics 2014, Vol. 8, No. 1, 622 647 DOI: 10.1214/13-AOAS711 Institute of Mathematical Statistics, 2014 MODELING EXTREME VALUES OF PROCESSES OBSERVED AT IRREGULAR TIME STEPS: APPLICATION
More informationCopulas. MOU Lili. December, 2014
Copulas MOU Lili December, 2014 Outline Preliminary Introduction Formal Definition Copula Functions Estimating the Parameters Example Conclusion and Discussion Preliminary MOU Lili SEKE Team 3/30 Probability
More information5.2 Annual maximum sea levels in Venice
Chapter 5 Non stationary extremes 5.1 Introduction In the context of environmental processes, it is common to observe non stationarity for example, due to different seasons having different climate patterns,
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 informationRegional Estimation from Spatially Dependent Data
Regional Estimation from Spatially Dependent Data R.L. Smith Department of Statistics University of North Carolina Chapel Hill, NC 27599-3260, USA December 4 1990 Summary Regional estimation methods are
More informationModeling Extremal Events Is Not Easy: Why the Extreme Value Theorem Cannot Be As General As the Central Limit Theorem
University of Texas at El Paso DigitalCommons@UTEP Departmental Technical Reports (CS) Department of Computer Science 5-2015 Modeling Extremal Events Is Not Easy: Why the Extreme Value Theorem Cannot Be
More informationFin285a:Computer Simulations and Risk Assessment Section 6.2 Extreme Value Theory Daníelson, 9 (skim), skip 9.5
Fin285a:Computer Simulations and Risk Assessment Section 6.2 Extreme Value Theory Daníelson, 9 (skim), skip 9.5 Overview Extreme value distributions Generalized Pareto distributions Tail shapes Using power
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 informationSkew Generalized Extreme Value Distribution: Probability Weighted Moments Estimation and Application to Block Maxima Procedure
Skew Generalized Extreme Value Distribution: Probability Weighted Moments Estimation and Application to Block Maxima Procedure Pierre Ribereau ( ), Esterina Masiello ( ), Philippe Naveau ( ) ( ) Université
More informationPRELIMINARY DRAFT FOR DISCUSSION PURPOSES
Memorandum To: David Thompson From: John Haapala CC: Dan McDonald Bob Montgomery Date: February 24, 2003 File #: 1003551 Re: Lake Wenatchee Historic Water Levels, Operation Model, and Flood Operation This
More informationHomework 2. For the homework, be sure to give full explanations where required and to turn in any relevant plots.
Homework 2 1 Data analysis problems For the homework, be sure to give full explanations where required and to turn in any relevant plots. 1. The file berkeley.dat contains average yearly temperatures for
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 informationColorado State University, Fort Collins, CO Weather Station Monthly Summary Report
Colorado State University, Fort Collins, CO Weather Station Monthly Summary Report Month: January Year: 2019 Temperature: Mean T max was 47.2 F which is 3.0 above the 1981-2010 normal for the month. This
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 informationFalse Alarm Probability based on bootstrap and extreme-value methods for periodogram peaks
False Alarm Probability based on bootstrap and extreme-value methods for periodogram peaks M. Süveges ISDC Data Centre for Astrophysics, University of Geneva, Switzerland Abstract One of the most pertinent
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 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 informationFrequency Estimation of Rare Events by Adaptive Thresholding
Frequency Estimation of Rare Events by Adaptive Thresholding J. R. M. Hosking IBM Research Division 2009 IBM Corporation Motivation IBM Research When managing IT systems, there is a need to identify transactions
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 informationApplications of Tail Dependence II: Investigating the Pineapple Express. Dan Cooley Grant Weller Department of Statistics Colorado State University
Applications of Tail Dependence II: Investigating the Pineapple Express Dan Cooley Grant Weller Department of Statistics Colorado State University Joint work with: Steve Sain, Melissa Bukovsky, Linda Mearns,
More information2015: A YEAR IN REVIEW F.S. ANSLOW
2015: A YEAR IN REVIEW F.S. ANSLOW 1 INTRODUCTION Recently, three of the major centres for global climate monitoring determined with high confidence that 2015 was the warmest year on record, globally.
More information