Bayesian Model Averaging for Multivariate Extreme Values
|
|
- Cecily Goodman
- 6 years ago
- Views:
Transcription
1 Bayesian Model Averaging for Multivariate Extreme Values Philippe Naveau Laboratoire des Sciences du Climat et l Environnement (LSCE) Gif-sur-Yvette, France joint work with A. Sabourin and A-L. Fougères FP7-ACQWA, GIS-PEPER, MIRACLE & ANR-McSim, MOPERA 14 novembre 2011
2 Environmental Data Analysis Extreme Value Theory Bayesian Model Averaging
3 Environmental Data Analysis Extreme Value Theory Bayesian Model Averaging
4 Air pollutants (Leeds, UK, winter 94-98, daily max) NO vs. PM10 (left), SO2 vs. PM10 (center), and SO2 vs. NO (right) Heffernan& Tawn 2004, Boldi & Davison, 2007, Cooley, Davis, Naveau, PM10 NO PM10 SO NO SO2
5 Typical question What is the probability of observing data in the blue box?
6 100 largest extremes NO SO PM PM NO S02
7 Environmental Data Analysis Extreme Value Theory Bayesian Model Averaging
8 Multivariate Extreme Value Theory (de Haan, Resnick and others) Maxima Max-stability G t (tz) =G(z) Regularly varying High quantiles Scaling property Λ(tA z )=t 1 Λ(A z ) Tail behavior Counting exceedances
9 Siméon Denis Poisson ( ) Counting excesses As a sum of random binary events, the variable N n that counts the number of events above the threshold u n has mean n Pr(X > u n) Poisson s theorem in 1837 If u n such that lim n Pr(X > un) =λ (0, ). n then N n follows approximately a Poisson variable N.
10 Still counting P(max(X 1,...,X n)/a n x, max(y 1,...,Y n)/a n y) =P(N n(a) =0) Y_i/a_n A y * * * * * * * * * * * * * * x X_i/a_n
11 Still counting P(max(X 1,...,X n)/a n x, max(y 1,...,Y n)/a n y) =P(N n(a) =0) Poisson again If then lim E(Nn(A)) = Λ(A), n lim P(Nn(A) =0) =P(N(A) =0) =exp( Λ(A)) n
12 Still counting P(max(X 1,...,X n)/a n x, max(y 1,...,Y n)/a n y) =P(N n(a) =0) Poisson again If then lim E(Nn(A)) = Λ(A), n lim P(Nn(A) =0) =P(N(A) =0) =exp( Λ(A)) n Two questions What is the sequence a n? What are its properties of Λ(A)?
13 Back to univariate case : Fréchet margins We know for the univariate GEV case with heavy-tailed and lim P(max(X 1,...,X n n)/a n x) =exp( x α ) with a n such that P(X > a n)=1/n Poisson condition with lim np(x/an Ax) =Λx(Ax) n Λ x(a x)=x α, for A x =[x, )
14 Scaling property Univariate case with Λ x(a x)=x α Λ x(ta x)=t α Λ x(a x) Multivariate case Λ(tA) =t α Λ(A)
15 Scaling property : an essential property of inference Λ(tA) =t α Λ(A) t A t α Λ t α A Area with data points
16 Interpreting the scaling property Λ(tA) =t α Λ(A) with α = 1 and y = y 1 + y 2 t y 1 + y 2 = t { {y : y/ y B and y >t} 1 y 1 + y 2 =1 B 1 t
17 Interpreting the scaling property Λ(tA) =t 1 Λ(A) A special case A = {x : x/r B and r > 1} where x, r = x and B any set belonging to the unit sphere A surprising property ta = {tx : x/r B and r > 1}, = {u : u/ u B and u > t}, with u = tx. This implies Λ({u : u/ u B and u > t}) =t 1 H (B) where H(.) spectral measure restricted to the unit sphere
18 Interpreting the scaling property Λ(tA) =t 1 Λ(A) A special case A = {x : x/r B and r > 1} where x, r = x and B any set belonging to the unit sphere A surprising property This implies ta = {tx : x/r B and r > 1}, = {u : u/ u B and u > t}, with u = tx. Λ({u : u/ u B and u > t}) =t 1 H (B) where H(.) spectral measure restricted to the unit sphere Independence between the radius r = x and the spectral measure The dependence among extremes is only captured by the spectral measure
19 Polar coordinates in 3D Radius r = x 1 + x 2 + x 3 and angle vector : w 1 = x 1 r, w 2 = x 2 r, w 3 = x 3 r
20 100 largest extremes NO SO PM PM NO S02
21 Dependence among the 100 angles W =(W 1, W 2, W 3 ) NO SO 2 PM 10
22 Our main problems How to find appropriate models to describe the dependence over the simplex? How to infer the parameters of our models? How to combine competing models?
23 An unique moment constraint the spectral measure H R Simplex w idh(w) = 1 d
24 An unique moment constraint the spectral measure H R Simplex w idh(w) = 1 d Non-parametric versus parametric In theory, there is no difference between theory and practice. But, in practice, there is. Jan L. A. van de Snepscheut or Yogi Berra
25 Environmental Data Analysis Extreme Value Theory Bayesian Model Averaging
26 Bayesian model averaging (BMA) Model 1 Model 2 S collection likelihood of distribution functions likelihood on S: j {h 1 (. θ 1 ), θ j Θ 1 }. {h 2 (. θ 2 ), θ 2 Θ 2 } verage consists in adding a prior layer to the m
27 Bayesian model averaging (BMA) Model 1 Model 2 S collection likelihood of distribution functions likelihood on S: j {h 1 (. θ 1 ), θ j Θ 1 }. {h 2 (. θ 2 ), θ 2 Θ 2 } verage consists in adding a prior layer to the m Objective Compute the posterior predictive density of the quantity of interest h(w data) =p(model 1 data) h 1 (w data)+p(model 2 data) h 2 (w data)
28 Bayesian model averaging (BMA) Model 1 Model 2 S collection likelihood of distribution functions likelihood on S: j {h 1 (. θ 1 ), θ j Θ 1 }. {h 2 (. θ 2 ), θ 2 Θ 2 } verage consists in adding a prior layer to the m Objective Compute the posterior predictive density of the quantity of interest h(w data) =p(model 1 data) h 1 (w data)+p(model 2 data) h 2 (w data) BMA : Sloughter et al. (2010), Raftery et al. (2005), Hoeting etal. (1999)
29 Bayesian model averaging (BMA) Model 1 Model 2 S collection likelihood of distribution functions likelihood on S: j {h 1 (. θ 1 ), θ j Θ 1 }. {h 2 (. θ 2 ), θ 2 Θ 2 } verage consists in adding a prior layer to the m Objective Compute the posterior predictive density of the quantity of interest h(w data) =p(model 1 data) h 1 (w data) + p(model 2 data) h 2 (w data) Z h 1 (w data) = h 1 (w data,θ 1 )[posterior of θ 1 data]
30 Bayesian model averaging (BMA) Model 1 Model 2 S collection likelihood of distribution functions likelihood on S: j {h 1 (. θ 1 ), θ j Θ 1 }. {h 2 (. θ 2 ), θ 2 Θ 2 } verage consists in adding a prior layer to the m Objective Compute the posterior predictive density of the quantity of interest h(w data) =p(model 1 data) h 1 (w data)+p(model 2 data) h 2 (w data) p(model 1 data) = p(data model 1) p(model 1) p(data) with Z p(data model 1) = marginal likelihood wrt θ 1 = h 1 (data θ 1 ) prior(θ 1 )
31 Model 1 Model 2 S collection likelihood of distribution functions likelihood on S: j {h 1 (. θ 1 ), θ j Θ 1 }. {h 2 (. θ 2 ), θ 2 Θ 2 } verage consists in adding a prior layer to the m Priors Priors Averaged Model } likelihood h(w (θ 1, θ 2 )) = p 1 h 1 (w θ 1 )+p 2 h 2 (w θ 2 ) Priors
32 Model 1 Marginal likelihood Z p(data model 1) = h 1(data θ 1) prior(θ 1) Computationally hard posterior weights p(model 1 data) = p(data model 1) p(model 1) p(data) Computationally Z easy Computationally easy Averaging Model 1 and 2 h(w data) =p(model 1 data) h 1 (w data)+p(model 2 data) h 2 (w data)
33 Environmental Data Analysis Extreme Value Theory Bayesian Model Averaging
34 Multivariate Extreme Value Theory and BMA A mixture of max-stable distributions is not max -stable A mixture of spectral measures is still a valid spectral measure, i.e. R Simplex w idh(w) = 1 d for dh = P p j dh j
35 Multivariate Extreme Value Theory and BMA Proposition 3.1. Let H 1,...H J be J angular spectral measures associated to J max-stable measures ν j (.). Let(p 1,...p J ) be a vector of positive weights summing to one. If H corresponds to their weighted average H = J j=1 p jh j, then (i) H is a valid spectral measure for a multivariate max-stable random vector M with unit-fréchet margins, exponent measure ν([0, x] c )= J p j ν j ([0, x] c ) j=1 (ii) M has max-combination representation J M = d p j M j j=1
36 Choosing two spectral parametric densities Model 1 PB Pairwise Beta Cooley, Davis, Naveau, 2010 Model 2 NL Nested Asym Logistic Gumbel, 1960, Tawn, 1990
37 Motivation Data EVT BMA BMA+EVT Wrapping up Simulation from two spectral parametric densities Model 1 PB w2 Model 2 NL w w alpha = 0.9 beta[1] = 15 beta[2] = 8 beta[3] = 0.5 w w3 w1
38 Our two spectral parametric densities where with Model 1 PB hpb(w α, β) = 1 i<j d hi,j(w α, βi,j) hi,j(w α, βi,j) = Kd(α) wij 2α 1 (1 wij) (d 2)α d+2 Γ(2βij) Γ 2 (βij) wβi,j 1 i/ij w βi,j 1 j/ij wi wij = wi + wj, w i/ij = wi + wj where hnl(w1w2) = 1 α 3α u v Model 2 2 α α NL (1 w12) 1 α 1 (w1w2) 1 α 12 α 1 u2(α12 1) v α α12 α12α uα12 2 v α 2 1 α = w 12 α 1 + w2 1 α 12 α = u α12 +(1 (w1 + w2)) 1 α. S
39 Simulation from two spectral parametric densities Algorithm 1. Model 1 PB (i) Choose uniformly a pair (i <j) (ii) Generate independently R ij Beta(2α +1, (d 2)α) Θ ij Beta(β i,j, β i,j) S Dirichlet d 2(1,...,1) (iii) Change variables back to define W via Model 2 NL Algorithm 2. (Stephenson, 2003) (i) Generate independently S PS(α) and S12 PS(α12). (ii) Simulate three independent standard exponentials E1,E2,E3 (iii) Set for i {1, 2}, Xi = S12S 1/α αα12 12 Ei and X3 = S α. E3 Then, X =(X1,X2,X3) has the desired distribution. Proof. If is generated according to the above algorithm, the cond W i = R ijθ ij W j = R ij(1 Θ ij) W [ (i,j)] =(1 R)S
40 Metropolis-Hasting at work Model 1 Marginal likelihood Z p(data model 1) = h 1(data θ 1) prior(θ 1) Computationally hard posterior weights p(model 1 data) = p(data model 1) p(model 1) p(data) Computationally Z easy Computationally easy Averaging Model 1 and 2 h(w data) =p(model 1 data) h 1 (w data)+p(model 2 data) h 2 (w data)
41 A few simulation results (with parameters tuned to our Leeds data) True model = PB Model Nsim p(data model) stdev(p(data model)) p(model data) PB NL
42 A few simulation results (with parameters tuned to our Leeds data) True model = PB Model Nsim p(data model) stdev(p(data model)) p(model data) PB NL True model = NL Model Nsim p(data model) stdev(p(data model)) p(model data) PB NL
43 A few simulation results (with parameters tuned to our Leeds data) True model = Mixture (PB + NL)/2 Model Nsim p(data model) stdev(p(data model)) p(model data) PB NL Kullback-Leibler divergence (small = good) KL(h true, h PB) KL(h true, h NL) KL(h true, h BMA) 0.014
44 Back to our example Environmental Data Analysis Extreme Value Theory Bayesian Model Averaging
45 Back to our example : priors (dotted) and posteriors (histo) for each parameter ModelPB Model NL log( ) log( [1]) alpha alpha log( [2]) log( [3])
46 Motivation Data EVT BMA BMA+EVT Wrapping up Back to our example : posteriors spectral densities ModelPB Model NL w2 PB w2 NL w3 w w3 w1
47 Back to our example ModelPB Model NL
48 Back to our example BMA verdict Model Nsim p(data model) stdev(p(data model)) p(model data) PB NL
49 Take home messages Feasibility of implementing BMA for multivariate extremes (in low dimensions) Computations can quickly become intensive The choice and number of parametric models are important Asymmetric nested logistic well tailored to represent bridges Pairwise beta model is flexible and be generalized (see Ballani and Schlater s extensions) More research needed to extend BMA to mixtures Going fully Bayesian non-parametric (Segers and colleagues, Boldi and Davison)
Bayesian 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 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 informationBayesian model mergings for multivariate extremes Application to regional predetermination of oods with incomplete data
Bayesian model mergings for multivariate extremes Application to regional predetermination of oods with incomplete data Anne Sabourin PhD defense University of Lyon 1 Committee: Anthony Davison, Johan
More informationNon parametric modeling of multivariate extremes with Dirichlet mixtures
Non parametric modeling of multivariate extremes with Dirichlet mixtures Anne Sabourin Institut Mines-Télécom, Télécom ParisTech, CNRS-LTCI Joint work with Philippe Naveau (LSCE, Saclay), Anne-Laure Fougères
More informationNew Classes of Multivariate Survival Functions
Xiao Qin 2 Richard L. Smith 2 Ruoen Ren School of Economics and Management Beihang University Beijing, China 2 Department of Statistics and Operations Research University of North Carolina Chapel Hill,
More informationOverview of Extreme Value Theory. Dr. Sawsan Hilal space
Overview of Extreme Value Theory Dr. Sawsan Hilal space Maths Department - University of Bahrain space November 2010 Outline Part-1: Univariate Extremes Motivation Threshold Exceedances Part-2: Bivariate
More 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 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 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 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 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 informationMultivariate Non-Normally Distributed Random Variables
Multivariate Non-Normally Distributed Random Variables An Introduction to the Copula Approach Workgroup seminar on climate dynamics Meteorological Institute at the University of Bonn 18 January 2008, Bonn
More informationExtreme value statistics: from one dimension to many. Lecture 1: one dimension Lecture 2: many dimensions
Extreme value statistics: from one dimension to many Lecture 1: one dimension Lecture 2: many dimensions The challenge for extreme value statistics right now: to go from 1 or 2 dimensions to 50 or more
More 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 informationBayesian Multivariate Extreme Value Thresholding for Environmental Hazards
Bayesian Multivariate Extreme Value Thresholding for Environmental Hazards D. Lupton K. Abayomi M. Lacer School of Industrial and Systems Engineering Georgia Institute of Technology Institute for Operations
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 informationApproximating the Conditional Density Given Large Observed Values via a Multivariate Extremes Framework, with Application to Environmental Data
Approximating the Conditional Density Given Large Observed Values via a Multivariate Extremes Framework, with Application to Environmental Data Daniel Cooley 1 Richard A. Davis 2 Philippe Naveau 3 1 Department
More informationSparse Representation of Multivariate Extremes with Applications to Anomaly Ranking
Sparse Representation of Multivariate Extremes with Applications to Anomaly Ranking Nicolas Goix Anne Sabourin Stéphan Clémençon Abstract Capturing the dependence structure of multivariate extreme events
More informationSemi-parametric estimation of non-stationary Pickands functions
Semi-parametric estimation of non-stationary Pickands functions Linda Mhalla 1 Joint work with: Valérie Chavez-Demoulin 2 and Philippe Naveau 3 1 Geneva School of Economics and Management, University of
More informationStatistics of Extremes
Statistics of Extremes Anthony Davison c 211 http://stat.epfl.ch Multivariate Extremes 19 Componentwise maxima.................................................. 194 Standardization........................................................
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 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 informationMadogram and asymptotic independence among maxima
Madogram and asymptotic independence among maxima Armelle Guillou, Philippe Naveau, Antoine Schorgen To cite this version: Armelle Guillou, Philippe Naveau, Antoine Schorgen. Madogram and asymptotic independence
More informationPrediction for Max-Stable Processes via an Approximated Conditional Density
Prediction for Max-Stable Processes via an Approximated Conditional Density Daniel Cooley 1 Richard A. Davis 2 Philippe Naveau 3 1 Department of Statistics Colorado State University, Fort Collins, CO USA
More 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 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 informationBernstein polynomial angular densities of multivariate extreme value distributions
Bernstein polynomial angular densities of multivariate extreme value distributions Timothy E. Hanson Department of Statistics, University of South Carolina, Columbia, SC 29208, USA email: hansont@stat.sc.edu
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 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 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 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 informationCPSC 540: Machine Learning
CPSC 540: Machine Learning MCMC and Non-Parametric Bayes Mark Schmidt University of British Columbia Winter 2016 Admin I went through project proposals: Some of you got a message on Piazza. No news is
More informationarxiv: v1 [stat.ap] 28 Nov 2014
COMBINING REGIONAL ESTIMATION AND HISTORICAL FLOODS: A MULTIVARIATE SEMI-PARAMETRIC PEAKS-OVER-THRESHOLD MODEL WITH CENSORED DATA Anne Sabourin 1, Benjamin Renard 2 arxiv:1411.7782v1 [stat.ap] 28 Nov 2014
More informationIntroduction to Probabilistic Graphical Models: Exercises
Introduction to Probabilistic Graphical Models: Exercises Cédric Archambeau Xerox Research Centre Europe cedric.archambeau@xrce.xerox.com Pascal Bootcamp Marseille, France, July 2010 Exercise 1: basics
More informationMultivariate Heavy Tails, Asymptotic Independence and Beyond
Multivariate Heavy Tails, endence and Beyond Sidney Resnick School of Operations Research and Industrial Engineering Rhodes Hall Cornell University Ithaca NY 14853 USA http://www.orie.cornell.edu/ sid
More informationA conditional approach for multivariate extreme values
J. R. Statist. Soc. B (2004) 66, Part 3, pp. 497 546 A conditional approach for multivariate extreme values Janet E. Heffernan and Jonathan A. Tawn Lancaster University, UK [Read before The Royal Statistical
More informationStatistics of Extremes
Statistics of Extremes Anthony Davison c 29 http://stat.epfl.ch Multivariate Extremes 184 Europe............................................................. 185 Componentwise maxima..................................................
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 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 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 informationFoundations of Nonparametric Bayesian Methods
1 / 27 Foundations of Nonparametric Bayesian Methods Part II: Models on the Simplex Peter Orbanz http://mlg.eng.cam.ac.uk/porbanz/npb-tutorial.html 2 / 27 Tutorial Overview Part I: Basics Part II: Models
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 informationPROBABILITY DISTRIBUTIONS. J. Elder CSE 6390/PSYC 6225 Computational Modeling of Visual Perception
PROBABILITY DISTRIBUTIONS Credits 2 These slides were sourced and/or modified from: Christopher Bishop, Microsoft UK Parametric Distributions 3 Basic building blocks: Need to determine given Representation:
More informationExtreme Value for Discrete Random Variables Applied to Avalanche Counts
Extreme Value for Discrete Random Variables Applied to Avalanche Counts Pascal Alain Sielenou IRSTEA / LSCE / CEN (METEO FRANCE) PostDoc (MOPERA and ECANA projects) Supervised by Nicolas Eckert and Philippe
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 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 informationBayes methods for categorical data. April 25, 2017
Bayes methods for categorical data April 25, 2017 Motivation for joint probability models Increasing interest in high-dimensional data in broad applications Focus may be on prediction, variable selection,
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 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 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 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 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 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 informationPREPRINT 2005:38. Multivariate Generalized Pareto Distributions HOLGER ROOTZÉN NADER TAJVIDI
PREPRINT 2005:38 Multivariate Generalized Pareto Distributions HOLGER ROOTZÉN NADER TAJVIDI Department of Mathematical Sciences Division of Mathematical Statistics CHALMERS UNIVERSITY OF TECHNOLOGY GÖTEBORG
More 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 informationWorkshop Copulas and Extremes
Workshop Copulas and Extremes Organization : Marianne Clausel & Stéphane Girard November, 19-20th, 2013 Tuesday, 9:00 9:30, Welcome coffee. Tuesday, 9:30 12:00, Session I, chair Anne-Catherine Favre. Some
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 informationStat260: Bayesian Modeling and Inference Lecture Date: February 10th, Jeffreys priors. exp 1 ) p 2
Stat260: Bayesian Modeling and Inference Lecture Date: February 10th, 2010 Jeffreys priors Lecturer: Michael I. Jordan Scribe: Timothy Hunter 1 Priors for the multivariate Gaussian Consider a multivariate
More informationBayesian Modeling of Air Pollution Extremes Using Nested Multivariate Max-Stable Processes
Bayesian Modeling of Air Pollution Extremes Using Nested Multivariate Max-Stable Processes arxiv:184.4588v1 [stat.ap] 18 Mar 218 Sabrina Vettori 1, Raphaël Huser 1, and Marc G. Genton 1 April 13, 218 Abstract
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 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 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 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 informationCS Lecture 19. Exponential Families & Expectation Propagation
CS 6347 Lecture 19 Exponential Families & Expectation Propagation Discrete State Spaces We have been focusing on the case of MRFs over discrete state spaces Probability distributions over discrete spaces
More informationCombining regional estimation and historical floods: a multivariate semi-parametric peaks-over-threshold model with censored data
Combining regional estimation and historical floods: a multivariate semi-parametric peaks-over-threshold model with censored data Anne Sabourin, Benjamin Renard To cite this version: Anne Sabourin, Benjamin
More informationExtreme Values on Spatial Fields p. 1/1
Extreme Values on Spatial Fields Daniel Cooley Department of Applied Mathematics, University of Colorado at Boulder Geophysical Statistics Project, National Center for Atmospheric Research Philippe Naveau
More informationMachine Learning Basics Lecture 7: Multiclass Classification. Princeton University COS 495 Instructor: Yingyu Liang
Machine Learning Basics Lecture 7: Multiclass Classification Princeton University COS 495 Instructor: Yingyu Liang Example: image classification indoor Indoor outdoor Example: image classification (multiclass)
More informationCMPE 58K Bayesian Statistics and Machine Learning Lecture 5
CMPE 58K Bayesian Statistics and Machine Learning Lecture 5 Multivariate distributions: Gaussian, Bernoulli, Probability tables Department of Computer Engineering, Boğaziçi University, Istanbul, Turkey
More informationCurve Fitting Re-visited, Bishop1.2.5
Curve Fitting Re-visited, Bishop1.2.5 Maximum Likelihood Bishop 1.2.5 Model Likelihood differentiation p(t x, w, β) = Maximum Likelihood N N ( t n y(x n, w), β 1). (1.61) n=1 As we did in the case of the
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 informationNonparametric Bayesian Methods - Lecture I
Nonparametric Bayesian Methods - Lecture I Harry van Zanten Korteweg-de Vries Institute for Mathematics CRiSM Masterclass, April 4-6, 2016 Overview of the lectures I Intro to nonparametric Bayesian statistics
More informationExponential Families
Exponential Families David M. Blei 1 Introduction We discuss the exponential family, a very flexible family of distributions. Most distributions that you have heard of are in the exponential family. Bernoulli,
More informationOn the occurrence times of componentwise maxima and bias in likelihood inference for multivariate max-stable distributions
On the occurrence times of componentwise maxima and bias in likelihood inference for multivariate max-stable distributions J. L. Wadsworth Department of Mathematics and Statistics, Fylde College, Lancaster
More informationNonparametric Bayesian modeling for dynamic ordinal regression relationships
Nonparametric Bayesian modeling for dynamic ordinal regression relationships Athanasios Kottas Department of Applied Mathematics and Statistics, University of California, Santa Cruz Joint work with Maria
More informationBayesian Learning. HT2015: SC4 Statistical Data Mining and Machine Learning. Maximum Likelihood Principle. The Bayesian Learning Framework
HT5: SC4 Statistical Data Mining and Machine Learning Dino Sejdinovic Department of Statistics Oxford http://www.stats.ox.ac.uk/~sejdinov/sdmml.html Maximum Likelihood Principle A generative model for
More informationA Bayesian Spatial Model for Exceedances Over a Threshold
A Bayesian Spatial odel for Exceedances Over a Threshold Fernando Ferraz do Nascimento and Bruno Sansó June 2, 2017 Abstract Extreme value theory focuses on the study of rare events and uses asymptotic
More informationA weather generator for simulating multivariate climatic series
A weather generator for simulating multivariate climatic series Denis Allard, with Nadine Brisson (AgroClim, INRA), Cédric Flecher (MetNext) and Philippe Naveau (LSCE, CNRS) Biostatistics and Spatial Processes
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 informationComputational statistics
Computational statistics Markov Chain Monte Carlo methods Thierry Denœux March 2017 Thierry Denœux Computational statistics March 2017 1 / 71 Contents of this chapter When a target density f can be evaluated
More informationLearning Bayesian network : Given structure and completely observed data
Learning Bayesian network : Given structure and completely observed data Probabilistic Graphical Models Sharif University of Technology Spring 2017 Soleymani Learning problem Target: true distribution
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 informationTail dependence in bivariate skew-normal and skew-t distributions
Tail dependence in bivariate skew-normal and skew-t distributions Paola Bortot Department of Statistical Sciences - University of Bologna paola.bortot@unibo.it Abstract: Quantifying dependence between
More informationModelling extreme-value dependence in high dimensions using threshold exceedances
Modelling extreme-value dependence in high dimensions using threshold exceedances Anna Kiriliouk A thesis submitted to the Université catholique de Louvain in partial fulfillment of the requirements for
More informationSupplementary Material for On the evaluation of climate model simulated precipitation extremes
Supplementary Material for On the evaluation of climate model simulated precipitation extremes Andrea Toreti 1 and Philippe Naveau 2 1 European Commission, Joint Research Centre, Ispra (VA), Italy 2 Laboratoire
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 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 information13 : Variational Inference: Loopy Belief Propagation and Mean Field
10-708: Probabilistic Graphical Models 10-708, Spring 2012 13 : Variational Inference: Loopy Belief Propagation and Mean Field Lecturer: Eric P. Xing Scribes: Peter Schulam and William Wang 1 Introduction
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 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 informationMath 576: Quantitative Risk Management
Math 576: Quantitative Risk Management Haijun Li lih@math.wsu.edu Department of Mathematics Washington State University Week 11 Haijun Li Math 576: Quantitative Risk Management Week 11 1 / 21 Outline 1
More informationBayesian model selection for computer model validation via mixture model estimation
Bayesian model selection for computer model validation via mixture model estimation Kaniav Kamary ATER, CNAM Joint work with É. Parent, P. Barbillon, M. Keller and N. Bousquet Outline Computer model validation
More informationFundamentals. CS 281A: Statistical Learning Theory. Yangqing Jia. August, Based on tutorial slides by Lester Mackey and Ariel Kleiner
Fundamentals CS 281A: Statistical Learning Theory Yangqing Jia Based on tutorial slides by Lester Mackey and Ariel Kleiner August, 2011 Outline 1 Probability 2 Statistics 3 Linear Algebra 4 Optimization
More informationConditioned limit laws for inverted max-stable processes
Conditioned limit laws for inverted max-stable processes Ioannis Papastathopoulos * and Jonathan A. Tawn * School of Mathematics and Maxwell Institute, University of Edinburgh, Edinburgh, EH9 3FD The Alan
More informationIntelligent Systems:
Intelligent Systems: Undirected Graphical models (Factor Graphs) (2 lectures) Carsten Rother 15/01/2015 Intelligent Systems: Probabilistic Inference in DGM and UGM Roadmap for next two lectures Definition
More informationBayesian non-parametric model to longitudinally predict churn
Bayesian non-parametric model to longitudinally predict churn Bruno Scarpa Università di Padova Conference of European Statistics Stakeholders Methodologists, Producers and Users of European Statistics
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 informationInvariant HPD credible sets and MAP estimators
Bayesian Analysis (007), Number 4, pp. 681 69 Invariant HPD credible sets and MAP estimators Pierre Druilhet and Jean-Michel Marin Abstract. MAP estimators and HPD credible sets are often criticized in
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 information