Generalized additive modelling of hydrological sample extremes
|
|
- Leslie Benson
- 6 years ago
- Views:
Transcription
1 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 Lausanne, Switzerland October 31, 2013 MFEW01, Isaac Newton Institute Valérie Chavez-Demoulin Joint work with A.C. Generalized Davison (EPFL) additive and modelling Marius Hofert hydrological (ETHZ) sample (Lausanne) extremes 1 / 28
2 Hydrological time series Hydrological time series are the result of complex dynamical processes (precipitation, snow accumulation and melt, evatranspiration,...) Daily maximum flow (1923:2008) Day of year Valérie Chavez-Demoulin Joint work with A.C. Generalized Davison (EPFL) additive and modelling Marius Hofert hydrological (ETHZ) sample (Lausanne) extremes 2 / 28
3 Point changes like moves of a station, changes in measuring instruments and hydro-electric installations may lead to discontinuities Daily maximum flow Year Valérie Chavez-Demoulin Joint work with A.C. Generalized Davison (EPFL) additive and modelling Marius Hofert hydrological (ETHZ) sample (Lausanne) extremes 3 / 28
4 Presence of seasonality, trends,... Daily streamflow records and their extremes are often dependent and not identically distributed Non-stationarity within one year or over longer periods = In this case EVT is not directly applicable! Variation due to the specifications of the station may be summarized parametrically Changes in time do not need to have a specific parametric form Valérie Chavez-Demoulin Joint work with A.C. Generalized Davison (EPFL) additive and modelling Marius Hofert hydrological (ETHZ) sample (Lausanne) extremes 4 / 28
5 Aim We combine the point process for exceedances with smoothing methods to give a flexible exploratory approach to model changes in the high flow exceedances The data are not declustered as we aim to model both long term and short term dependence Uncertainty assessment is made through appropriate confidence intervals Valérie Chavez-Demoulin Joint work with A.C. Generalized Davison (EPFL) additive and modelling Marius Hofert hydrological (ETHZ) sample (Lausanne) extremes 5 / 28
6 The block maxima approach I Consider X 1,..., X q iid from F (x) Suppose there exists a point x 0 (perhaps + ) such that lim x x0 F (x) = 1. For any fixed x < x 0 we have P(max{X 1,..., X q } x) = P(X i x, i = 1,..., q) = F q (x) which tends to 0 as q Given suitable sequences {a q } and {b q } of normalizing constants leading to W q = aq 1 {max(x 1,..., X q ) b q }, the non-degenerate limiting distribution must be a generalized extreme value (GEV) distribution { ( H µ,ψ,ξ (w) = exp 1 + ξ w µ ) } 1/ξ ψ < µ <, ψ > 0, < ξ < Valérie Chavez-Demoulin Joint work with A.C. Generalized Davison (EPFL) additive and modelling Marius Hofert hydrological (ETHZ) sample (Lausanne) extremes 6 / 28 +
7 The block maxima approach II We fit the GEV distribution to the series of (typically) annual maximum data (the blocks) From m blocks of size q W = (M q (1),..., M q (m) ) Construct a log likelihood by assuming we have independent observations from a GEV with density h µ,ψ,ξ (w) ( m l(µ, ψ, ξ; W ) = log h µ,ψ,ξ (M q (i) i=1 = ξ, µ, ψ )1 (i) {1+ξ(M q µ)/ψ>0} ) Valérie Chavez-Demoulin Joint work with A.C. Generalized Davison (EPFL) additive and modelling Marius Hofert hydrological (ETHZ) sample (Lausanne) extremes 7 / 28
8 Peaks-over-Threshold Method (I) Let X t1,..., X tn denote the exceedances over a high threshold u with corresponding excesses Y ti = X ti u, i {1,..., n} 1) the number of exceedances N t approximately follows a Poisson process with intensity λ, that is, N t Poi(λ(t)) with integrated rate function Λ(t) = λt 2) the excesses Y t1,..., Y tn over u approximately follow (independently t of N t ) a generalized Pareto distribution (GPD), denoted by GPD(ξ, β) for ξ R, β > 0, with distribution function G ξ,β (x) = { 1 ( 1 + ξx/β ) 1/ξ, if ξ 0, 1 exp( x/β), if ξ = 0, for x 0, if ξ 0, and x [0, β/ξ], if ξ < 0 Valérie Chavez-Demoulin Joint work with A.C. Generalized Davison (EPFL) additive and modelling Marius Hofert hydrological (ETHZ) sample (Lausanne) extremes 8 / 28
9 Peaks-over-Threshold Method (II) Asymptotic independence between Poisson exceedance times and GPD excesses (λt )n n L(λ, ξ, β; Y ) = exp( λt ) g ξ,β (Y ti ), n! i=1 where Y = (Y t1,..., Y tn ) and g ξ,β is the density of G ξ,β = l(λ, ξ, β; Y ) = l(λ; Y ) + l(ξ, β; Y ), l(λ; Y ) = λt + n log(λ) + log(t n /n!) and l(ξ, β; Y ) = with l(ξ, β; y) = n l(ξ, β; Y ti ) i=1 { log(β) (1 + 1/ξ) log(1 + ξy/β), if ξ 0, log(β) y/β, if ξ = 0, = Maximization can be carried out separately Valérie Chavez-Demoulin Joint work with A.C. Generalized Davison (EPFL) additive and modelling Marius Hofert hydrological (ETHZ) sample (Lausanne) extremes 9 / 28
10 A dynamic EVT approach Let θ R r be the vector of r EVT model parameters (r = 3 in both the GEV model and POT representation) Let x be the vector of covariates; t be the time, then a general model } θ i = g i {x T η i + h i (t), i = 1,..., r g i is a link function; η i R p are the parameter vectors; h i is smoothed nonparametric function of t A where A R is the subset on which t is observed θ R r can be estimated by maximizing the penalized log-likelihood r ] l(θ; y) [γ i h i (t) 2 dt A i=1 where l(θ; y) is the log-likelihood based on the EVT model (either block-maxima or POT) Valérie Chavez-Demoulin Joint work with A.C. Generalized Davison (EPFL) additive and modelling Marius Hofert hydrological (ETHZ) sample (Lausanne) extremes 10 / 28
11 Dynamic POT Approach: Non-homogeneous Poisson model for number of exceedances For the number of exceedances, a non-homogeneous Poisson process rate λ = λ(x, t) = exp(x η λ + h λ (t)) This model is a standard generalized additive (GAM) model Embedded models are compared using LRS (parametric part) Degrees of freedom for the non-parametric part (smoothing spline) are chosen using AIC (cf below) Valérie Chavez-Demoulin Joint work with A.C. Generalized Davison (EPFL) additive and modelling Marius Hofert hydrological (ETHZ) sample (Lausanne) extremes 11 / 28
12 Dynamic POT Approach: GAM GPD for the exceedance sizes For the GPD parameters, we replace β by which is orthogonal to ξ ν = log((1 + ξ)β) The corresponding reparameterized log-likelihood l r for the excesses is thus l r (ξ, ν; Y ) = l(ξ, exp(ν)/(1 + ξ); Y ) We assume that ξ and ν are of the form ξ = ξ(x, t) = x η ξ + h ξ (t) ν = ν(x, t) = x η ν + h ν (t) Valérie Chavez-Demoulin Joint work with A.C. Generalized Davison (EPFL) additive and modelling Marius Hofert hydrological (ETHZ) sample (Lausanne) extremes 12 / 28
13 Penalized maximum likelihood estimator (I) In order to fit reasonably smooth functions h ξ, h ν we use the penalized likelihood l p (η ξ, h ξ, η ν, h ν ; z 1,..., z n ) = l r (ξ, ν; y) γ ξ T 0 h ξ (t)2 dt γ ν T 0 h ν(t) 2 dt where γ ξ, γ ν 0 denote smoothing parameters, y = (y t1,..., y tn ), and l r (ξ, ν; y) = n l r (ξ i, ν i ; y ti ) i=1 for l r (ξ i, ν i ; y ti ) = l(ξ i, exp(ν i )/(1 + ξ i ); y ti ) Larger values of the smoothing parameters lead to smoother fitted curves. A related quantity is the equivalent degrees of freedom Valérie Chavez-Demoulin Joint work with A.C. Generalized Davison (EPFL) additive and modelling Marius Hofert hydrological (ETHZ) sample (Lausanne) extremes 13 / 28
14 Penalized maximum likelihood estimator (II) Let 0 = s 0 < s 1 < < s m < s m+1 = T denote the ordered and distinct values among {t 1,..., t n } for a natural cubic spline h with knots s 1,..., s m T 0 h (t) 2 dt = h Kh, where h = (h s1,..., h sm ) = (h(s 1 ),..., h(s m )) and K is a symmetric (m, m)-matrix of rank m 2 only depending on the knots s 1,..., s m Valérie Chavez-Demoulin Joint work with A.C. Generalized Davison (EPFL) additive and modelling Marius Hofert hydrological (ETHZ) sample (Lausanne) extremes 14 / 28
15 Penalized maximum likelihood estimator (III) The penalized log-likelihood can thus be written as l p (η ξ, h ξ, η ν, h ν ; z 1,..., z n ) = l r (ξ, ν; y) γ ξ h ξ Kh ξ γ ν h ν Kh ν with h ξ = (h ξ (s 1 ),..., h ξ (s m )) and h ν = (h ν (s 1 ),..., h ν (s m )). Backfitting algorithm for estimating simultaneously ξ and ν (thus β) Confidence intervals calculated using post-blackend Bootstrap Valérie Chavez-Demoulin Joint work with A.C. Generalized Davison (EPFL) additive and modelling Marius Hofert hydrological (ETHZ) sample (Lausanne) extremes 15 / 28
16 Return level estimation Based on the fitted values ˆλ, ˆξ, and ˆβ for a fixed covariate vector x and time point t, one can compute (depending on x and t) estimates of the 1/p-year return level R 1 p = u + ˆβ (( ) ˆξ ) 1 ˆξ pˆλ Valérie Chavez-Demoulin Joint work with A.C. Generalized Davison (EPFL) additive and modelling Marius Hofert hydrological (ETHZ) sample (Lausanne) extremes 16 / 28
17 Demonstration for simulated data (I) Details are provided as a demo in the package QRM of R (with Marius Hofert, ETHZ) We generate a data set of exceedances over a time period of 10 years for two groups (Group A and Group B) The simulated losses are drawn from a (non-stationary) generalized Pareto distribution depending on the covariates year and group We then fit the (Poisson process) intensity λ and the two parameters ξ and β of the generalized Pareto distribution depending on year and Group We then calculate a 99.9% VaR (equivalent to the return revel) Valérie Chavez-Demoulin Joint work with A.C. Generalized Davison (EPFL) additive and modelling Marius Hofert hydrological (ETHZ) sample (Lausanne) extremes 17 / 28
18 Demonstration for simulated data (II) λ^ with pointwise asymptotic two sided 0.95% confidence intervals A B Year λ^ 0.95 CI true number of losses Valérie Chavez-Demoulin Joint work with A.C. Generalized Davison (EPFL) additive and modelling Marius Hofert hydrological (ETHZ) sample (Lausanne) extremes 18 / 28
19 Demonstration for simulated data (III) ξ^ with bootstrapped pointwise two sided 0.95% confidence intervals A B Year ξ^ 0.95 CI true ξ Valérie Chavez-Demoulin Joint work with A.C. Generalized Davison (EPFL) additive and modelling Marius Hofert hydrological (ETHZ) sample (Lausanne) extremes 19 / 28
20 Demonstration for simulated data (IV) A B Year β^ 0.95 CI true β β^ with bootstrapped pointwise two sided 0.95% confidence intervals Valérie Chavez-Demoulin Joint work with A.C. Davison (EPFL) and Marius Hofert (ETHZ) (Lausanne) Generalized additive modelling of hydrological sample extremes 20 / 28
21 Demonstration for simulated data (V) e+02 2e+03 5e+03 2e+04 5e+04 2e+05 A B Year VaR CI true VaR0.999 VaR0.999 with bootstrapped ptw. two sided 0.95% confidence intervals Valérie Chavez-Demoulin Joint work with A.C. Davison (EPFL) and Marius Hofert (ETHZ) (Lausanne) Generalized additive modelling of hydrological sample extremes 21 / 28
22 Application to Muota-Ingenbohl data: AIC Exploratory purpose: use high degrees of freedom Nonparametric estimate of the Poisson intensity for the dependence on the day of year (degrees of freedom = 20) and for the dependence on year (degrees of freedom = 10) lambda Jun (161) 19Jul (200) 27Oct (300) lambda Year Day of year Valérie Chavez-Demoulin Joint work with A.C. Generalized Davison (EPFL) additive and modelling Marius Hofert hydrological (ETHZ) sample (Lausanne) extremes 22 / 28
23 Application to Muota-Ingenbohl data: AIC Automatic procedure for selecting the smoothing parameter via AIC AIC for xi~s(dayofyear,df) AIC for nu~s(dayofyear,df) Df Df AIC for xi~s(year,df) AIC for nu~s(year,df) Df Df Valérie Chavez-Demoulin Joint work with A.C. Generalized Davison (EPFL) additive and modelling Marius Hofert hydrological (ETHZ) sample (Lausanne) extremes 23 / 28
24 Application to Muota-Ingenbohl data: Model and parameters estimates The selected model: log λ(d, t) = h (3) λ (3) (t) + g λ (d) ˆξ(d, t) = h (2) ξ (t) + g (2) ξ (d) ˆν(d, t) = h ν (2) (t) + g ν (2) (d) xi nu lambda Year Year Year Valérie Chavez-Demoulin Joint work with A.C. Generalized Davison (EPFL) additive and modelling Marius Hofert hydrological (ETHZ) sample (Lausanne) extremes 24 / 28
25 Application to Muota-Ingenbohl data: Goodness-of-fit If the excesses Y t1,..., Y tn (approximately) follow a GPD(ξ, β), then R i Exp(1), i {1,..., n}, where R i = log(1 + Y ti ξ i /β i )/ξ i, i {1,..., n} We can thus graphically check whether (approximately) r ti = log(1 + y ti ˆξ i / ˆβ i )/ˆξ i, i {1,..., n}, are distributed as independent standard exponential variables Residuals Quantiles of Exponential Valérie Chavez-Demoulin Joint work with A.C. Generalized Davison (EPFL) additive and modelling Marius Hofert hydrological (ETHZ) sample (Lausanne) extremes 25 / 28
26 Application to Muota-Ingenbohl data: 20-years return level years return level Valérie Chavez-Demoulin Joint work with A.C. Davison (EPFL) and Marius Hofert (ETHZ) (Lausanne) Generalized additive modelling of hydrological sample extremes 26 / 28
27 Thank you! Valérie Chavez-Demoulin Joint work with A.C. Generalized Davison (EPFL) additive and modelling Marius Hofert hydrological (ETHZ) sample (Lausanne) extremes 27 / 28
28 References Chavez-Demoulin, V. and Davison, A.C., (2005), Generalized additive models for sample extremes. Applied Statistics 54(1), Chavez-Demoulin V. and Embrechts P. (2004), Smooth extremal models in finance and insurance. Journal of Risk and Insurance, 71(2), Chavez-Demoulin V. and Embrechts P. and Hofert, M.(2013) An extreme value approach for modeling Operational Risk losses depending on covariates. Submitted. Yee, T.W. and Stephenson, A.G. (2007). Vector generalized linear and additive extreme value models. Extremes,9, Laurini, F. and Pauli, F. (2009). Smoothing sample extremes: The mixed model approach. Comp. Statist. Data Anal., 53, Pauli, F. and Coles, S.G. (2001). Penalized likelihood inference in extreme value analyses. J. Appl. Statist., 28, Valérie Chavez-Demoulin Joint work with A.C. Generalized Davison (EPFL) additive and modelling Marius Hofert hydrological (ETHZ) sample (Lausanne) extremes 28 / 28
High-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 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 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 informationModeling Operational Risk Depending on Covariates. An Empirical Investigation.
Noname manuscript No. (will be inserted by the editor) Modeling Operational Risk Depending on Covariates. An Empirical Investigation. Paul Embrechts Kamil J. Mizgier Xian Chen Received: date / Accepted:
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 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 informationExtreme Value Analysis and Spatial Extremes
Extreme Value Analysis and Department of Statistics Purdue University 11/07/2013 Outline Motivation 1 Motivation 2 Extreme Value Theorem and 3 Bayesian Hierarchical Models Copula Models Max-stable Models
More informationExtreme Value Theory and Applications
Extreme Value Theory and Deauville - 04/10/2013 Extreme Value Theory and Introduction Asymptotic behavior of the Sum Extreme (from Latin exter, exterus, being on the outside) : Exceeding the ordinary,
More 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 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 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 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 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 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 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 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 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 informationBayesian Point Process Modeling for Extreme Value Analysis, with an Application to Systemic Risk Assessment in Correlated Financial Markets
Bayesian Point Process Modeling for Extreme Value Analysis, with an Application to Systemic Risk Assessment in Correlated Financial Markets Athanasios Kottas Department of Applied Mathematics and Statistics,
More informationA Conditional Approach to Modeling Multivariate Extremes
A Approach to ing Multivariate Extremes By Heffernan & Tawn Department of Statistics Purdue University s April 30, 2014 Outline s s Multivariate Extremes s A central aim of multivariate extremes is trying
More 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 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 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 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 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 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 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 informationModelling geoadditive survival data
Modelling geoadditive survival data Thomas Kneib & Ludwig Fahrmeir Department of Statistics, Ludwig-Maximilians-University Munich 1. Leukemia survival data 2. Structured hazard regression 3. Mixed model
More informationTesting Hypothesis. Maura Mezzetti. Department of Economics and Finance Università Tor Vergata
Maura Department of Economics and Finance Università Tor Vergata Hypothesis Testing Outline It is a mistake to confound strangeness with mystery Sherlock Holmes A Study in Scarlet Outline 1 The Power Function
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 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 informationParameter Estimation
Parameter Estimation Chapters 13-15 Stat 477 - Loss Models Chapters 13-15 (Stat 477) Parameter Estimation Brian Hartman - BYU 1 / 23 Methods for parameter estimation Methods for parameter estimation Methods
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 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 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 informationConstruction of confidence intervals for extreme rainfall quantiles
Risk Analysis VIII 93 Construction of confidence intervals for extreme rainfall quantiles A. T. Silva 1, M. M. Portela 1, J. Baez & M. Naghettini 3 1 Instituto Superior Técnico, Portugal Universidad Católica
More informationLecture 26: Likelihood ratio tests
Lecture 26: Likelihood ratio tests Likelihood ratio When both H 0 and H 1 are simple (i.e., Θ 0 = {θ 0 } and Θ 1 = {θ 1 }), Theorem 6.1 applies and a UMP test rejects H 0 when f θ1 (X) f θ0 (X) > c 0 for
More informationBayesian covariate models in extreme value analysis
Bayesian covariate models in extreme value analysis David Randell, Philip Jonathan, Kathryn Turnbull, Mathew Jones EVA 2015 Ann Arbor Copyright 2015 Shell Global Solutions (UK) EVA 2015 Ann Arbor June
More informationNonparametric inference in hidden Markov and related models
Nonparametric inference in hidden Markov and related models Roland Langrock, Bielefeld University Roland Langrock Bielefeld University 1 / 47 Introduction and motivation Roland Langrock Bielefeld University
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 informationStatistics - Lecture One. Outline. Charlotte Wickham 1. Basic ideas about estimation
Statistics - Lecture One Charlotte Wickham wickham@stat.berkeley.edu http://www.stat.berkeley.edu/~wickham/ Outline 1. Basic ideas about estimation 2. Method of Moments 3. Maximum Likelihood 4. Confidence
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 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 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 informationGeneralized Additive Models
Generalized Additive Models The Model The GLM is: g( µ) = ß 0 + ß 1 x 1 + ß 2 x 2 +... + ß k x k The generalization to the GAM is: g(µ) = ß 0 + f 1 (x 1 ) + f 2 (x 2 ) +... + f k (x k ) where the functions
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 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 informationSTAT 135 Lab 6 Duality of Hypothesis Testing and Confidence Intervals, GLRT, Pearson χ 2 Tests and Q-Q plots. March 8, 2015
STAT 135 Lab 6 Duality of Hypothesis Testing and Confidence Intervals, GLRT, Pearson χ 2 Tests and Q-Q plots March 8, 2015 The duality between CI and hypothesis testing The duality between CI and hypothesis
More informationExtreme value modelling of rainfalls and
Universite de Paris Sud Master s Thesis Extreme value modelling of rainfalls and flows Author: Fan JIA Supervisor: Pr Elisabeth Gassiat Dr Elena Di Bernardino Pr Michel Bera A thesis submitted in fulfilment
More informationof the 7 stations. In case the number of daily ozone maxima in a month is less than 15, the corresponding monthly mean was not computed, being treated
Spatial Trends and Spatial Extremes in South Korean Ozone Seokhoon Yun University of Suwon, Department of Applied Statistics Suwon, Kyonggi-do 445-74 South Korea syun@mail.suwon.ac.kr Richard L. Smith
More informationVariable Selection and Model Choice in Survival Models with Time-Varying Effects
Variable Selection and Model Choice in Survival Models with Time-Varying Effects Boosting Survival Models Benjamin Hofner 1 Department of Medical Informatics, Biometry and Epidemiology (IMBE) Friedrich-Alexander-Universität
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 informationPh.D. Qualifying Exam Friday Saturday, January 6 7, 2017
Ph.D. Qualifying Exam Friday Saturday, January 6 7, 2017 Put your solution to each problem on a separate sheet of paper. Problem 1. (5106) Let X 1, X 2,, X n be a sequence of i.i.d. observations from a
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 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 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 informationarxiv: v2 [stat.me] 25 Sep 2012
Estimation of Hüsler-Reiss distributions and Brown-Resnick processes arxiv:107.6886v [stat.me] 5 Sep 01 Sebastian Engelke Institut für Mathematische Stochastik, Georg-August-Universität Göttingen, Goldschmidtstr.
More informationLikelihood Ratio Tests. that Certain Variance Components Are Zero. Ciprian M. Crainiceanu. Department of Statistical Science
1 Likelihood Ratio Tests that Certain Variance Components Are Zero Ciprian M. Crainiceanu Department of Statistical Science www.people.cornell.edu/pages/cmc59 Work done jointly with David Ruppert, School
More informationTREND AND VARIABILITY ANALYSIS OF RAINFALL SERIES AND THEIR EXTREME
TREND AND VARIABILITY ANALYSIS OF RAINFALL SERIES AND THEIR EXTREME EVENTS J. Abaurrea, A. C. Cebrián. Dpto. Métodos Estadísticos. Universidad de Zaragoza. Abstract: Rainfall series and their corresponding
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 informationNon stationary extremes
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, or
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 informationExtreme value theory and high quantile convergence
Journal of Operational Risk 51 57) Volume 1/Number 2, Summer 2006 Extreme value theory and high quantile convergence Mikhail Makarov EVMTech AG, Baarerstrasse 2, 6300 Zug, Switzerland In this paper we
More informationn =10,220 observations. Smaller samples analyzed here to illustrate sample size effect.
Chapter 7 Parametric Likelihood Fitting Concepts: Chapter 7 Parametric Likelihood Fitting Concepts: Objectives Show how to compute a likelihood for a parametric model using discrete data. Show how to compute
More informationQuantile-quantile plots and the method of peaksover-threshold
Problems in SF2980 2009-11-09 12 6 4 2 0 2 4 6 0.15 0.10 0.05 0.00 0.05 0.10 0.15 Figure 2: qqplot of log-returns (x-axis) against quantiles of a standard t-distribution with 4 degrees of freedom (y-axis).
More informationBayesian trend analysis for daily rainfall series of Barcelona
Adv. Geosci., 26, 71 76, 2010 doi:10.5194/adgeo-26-71-2010 Author(s) 2010. CC Attribution 3.0 License. Advances in Geosciences Bayesian trend analysis for daily rainfall series of Barcelona M. I. Ortego
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 informationModelling spatially-dependent non-stationary extremes with application to hurricane-induced wave heights
Modelling spatially-dependent non-stationary extremes with application to hurricane-induced wave heights Paul Northrop and Philip Jonathan March 5, 2010 Abstract In environmental applications it is found
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 informationModeling Real Estate Data using Quantile Regression
Modeling Real Estate Data using Semiparametric Quantile Regression Department of Statistics University of Innsbruck September 9th, 2011 Overview 1 Application: 2 3 4 Hedonic regression data for house prices
More informationWhat Can We Infer From Beyond The Data? The Statistics Behind The Analysis Of Risk Events In The Context Of Environmental Studies
What Can We Infer From Beyond The Data? The Statistics Behind The Analysis Of Risk Events In The Context Of Environmental Studies Sibusisiwe Khuluse, Sonali Das, Pravesh Debba, Chris Elphinstone Logistics
More informationOn Backtesting Risk Measurement Models
On Backtesting Risk Measurement Models Hideatsu Tsukahara Department of Economics, Seijo University e-mail address: tsukahar@seijo.ac.jp 1 Introduction In general, the purpose of backtesting is twofold:
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 informationIntroduction to Regression
Introduction to Regression p. 1/97 Introduction to Regression Chad Schafer cschafer@stat.cmu.edu Carnegie Mellon University Introduction to Regression p. 1/97 Acknowledgement Larry Wasserman, All of Nonparametric
More informationRobust and Efficient Estimation for the Generalized Pareto Distribution
Robust and Efficient Estimation for the Generalized Pareto Distribution Sergio F. Juárez Faculty of Statistics and Informatics Veracruzana University, Xalapa, Ver, México email: sejuarez@uv.mx and William
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 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 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 informationStat 710: Mathematical Statistics Lecture 12
Stat 710: Mathematical Statistics Lecture 12 Jun Shao Department of Statistics University of Wisconsin Madison, WI 53706, USA Jun Shao (UW-Madison) Stat 710, Lecture 12 Feb 18, 2009 1 / 11 Lecture 12:
More informationBIO5312 Biostatistics Lecture 13: Maximum Likelihood Estimation
BIO5312 Biostatistics Lecture 13: Maximum Likelihood Estimation Yujin Chung November 29th, 2016 Fall 2016 Yujin Chung Lec13: MLE Fall 2016 1/24 Previous Parametric tests Mean comparisons (normality assumption)
More informationPh.D. Qualifying Exam Friday Saturday, January 3 4, 2014
Ph.D. Qualifying Exam Friday Saturday, January 3 4, 2014 Put your solution to each problem on a separate sheet of paper. Problem 1. (5166) Assume that two random samples {x i } and {y i } are independently
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 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 informationMathematical statistics
October 1 st, 2018 Lecture 11: Sufficient statistic Where are we? Week 1 Week 2 Week 4 Week 7 Week 10 Week 14 Probability reviews Chapter 6: Statistics and Sampling Distributions Chapter 7: Point Estimation
More informationPayer, Küchenhoff: Modelling extreme wind speeds in the context of risk analysis for high speed trains
Payer, Küchenhoff: Modelling extreme wind speeds in the context of risk analysis for high speed trains Sonderforschungsbereich 386, Paper 295 (2002) Online unter: http://epub.ub.uni-muenchen.de/ Projektpartner
More informationL-momenty s rušivou regresí
L-momenty s rušivou regresí Jan Picek, Martin Schindler e-mail: jan.picek@tul.cz TECHNICKÁ UNIVERZITA V LIBERCI ROBUST 2016 J. Picek, M. Schindler, TUL L-momenty s rušivou regresí 1/26 Motivation 1 Development
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 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 informationNonparametric Inference In Functional Data
Nonparametric Inference In Functional Data Zuofeng Shang Purdue University Joint work with Guang Cheng from Purdue Univ. An Example Consider the functional linear model: Y = α + where 1 0 X(t)β(t)dt +
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 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 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 informationQualifying Exam CS 661: System Simulation Summer 2013 Prof. Marvin K. Nakayama
Qualifying Exam CS 661: System Simulation Summer 2013 Prof. Marvin K. Nakayama Instructions This exam has 7 pages in total, numbered 1 to 7. Make sure your exam has all the pages. This exam will be 2 hours
More informationSCHOOL OF MATHEMATICS AND STATISTICS. Linear and Generalised Linear Models
SCHOOL OF MATHEMATICS AND STATISTICS Linear and Generalised Linear Models Autumn Semester 2017 18 2 hours Attempt all the questions. The allocation of marks is shown in brackets. RESTRICTED OPEN BOOK EXAMINATION
More informationSupplement to A Hierarchical Approach for Fitting Curves to Response Time Measurements
Supplement to A Hierarchical Approach for Fitting Curves to Response Time Measurements Jeffrey N. Rouder Francis Tuerlinckx Paul L. Speckman Jun Lu & Pablo Gomez May 4 008 1 The Weibull regression model
More informationReliable Inference in Conditions of Extreme Events. Adriana Cornea
Reliable Inference in Conditions of Extreme Events by Adriana Cornea University of Exeter Business School Department of Economics ExISta Early Career Event October 17, 2012 Outline of the talk Extreme
More informationNormalized kernel-weighted random measures
Normalized kernel-weighted random measures Jim Griffin University of Kent 1 August 27 Outline 1 Introduction 2 Ornstein-Uhlenbeck DP 3 Generalisations Bayesian Density Regression We observe data (x 1,
More informationAn application of the GAM-PCA-VAR model to respiratory disease and air pollution data
An application of the GAM-PCA-VAR model to respiratory disease and air pollution data Márton Ispány 1 Faculty of Informatics, University of Debrecen Hungary Joint work with Juliana Bottoni de Souza, Valdério
More informationLecture 17: Likelihood ratio and asymptotic tests
Lecture 17: Likelihood ratio and asymptotic tests Likelihood ratio When both H 0 and H 1 are simple (i.e., Θ 0 = {θ 0 } and Θ 1 = {θ 1 }), Theorem 6.1 applies and a UMP test rejects H 0 when f θ1 (X) f
More informationSimultaneous Confidence Bands for the Coefficient Function in Functional Regression
University of Haifa From the SelectedWorks of Philip T. Reiss August 7, 2008 Simultaneous Confidence Bands for the Coefficient Function in Functional Regression Philip T. Reiss, New York University Available
More information