Extreme value statistics: from one dimension to many. Lecture 1: one dimension Lecture 2: many dimensions
|
|
- Heather Perkins
- 5 years ago
- Views:
Transcription
1 Extreme value statistics: from one dimension to many Lecture 1: one dimension Lecture 2: many dimensions
2 The challenge for extreme value statistics right now: to go from 1 or 2 dimensions to 50 or more The challenge for computation: To maximize fifty-dimensional likelihoods with hundreds of parameters The challenge for probability: To construct parametric models which makes this possible and to understand the models (and then model validation, asymptotics, ) Extreme rainfall at many catchments, flooding of many dykes, storm insurance, landslide risk assessment, consistent risk estimation for financial portfolios, extreme influenza epidemics,
3 Notation Bold symbols are d-variate vectors. For instance, γγ = (γγ 1,..., γγ dd ) and 00 = (0,..., 0) RR dd. Operations and relations are componentwise, with shorter vectors recycled. For instance aaaa + bb = (aa 1 xx 1 + bb 1,..., aa dd xx dd + bb dd ), xx yy iiii xx jj yy jj for jj = 1,..., dd, and tt γγ = (tt γγ 1,, tt γγ dd)
4 dd-dimensional generalized extreme value (GEV) distributions No finite-dimensional parametrization Marginal distributions GEV: GG ii xx = ee 1+xx μμ ii σσ ii + Can be expressed through spectral representations. 1/γγ ii Can be expressed in terms of point processes, see later slides
5 Let ηη, dd be the vector of lower endpoints of the marginal distributions of GG The GEV distributions are the limit distributions of maxima: Let XX 1, XX 2, be an i.i.d. sequence of dd-dimensional vectors with cdf FF. If for some scaling and location sequences aa nn > 00 and bb nn PPPP aa 1 nn nn ii=1 XX ii bb n ηη xx dd GG xx, as nn where GG has non-degenerate margins, then GG is a GEV distribution. Conversely all GEV distributions can be obtained in this way.
6 The GEV distributions are the max-stable distributions: Let XX 1, XX 2, be an i.i.d. sequence with nondegenerate marginal cdf-s GG ii. If for some scaling and location sequences αα nn > 0 and ββ nn 1 nn Pr αα nn XX ii ββ nn xx = GG xx, for nn 1, 2, ii=1 then GG is a GEV distribution. Conversely all GEV distributions can be obtained in this way.
7 Think of XX 1, XX 2,, XX n as points in RR dd. Subtract aa nn and divide by bb nn to get a point process with points XX ii aa nn bb nn ηη; ii = 1, nn. It converges as nn iff XX DD(GG), to a limit point process xx 2 xx 1 {WW ii RR ii γγ } where (in sequel think of γγ > 00) WW ii are i.i.d. copies of some ( any ) positive stochastic dd-dimensional vector and the RR ii are points of a unit rate Poisson process on RR + WW 1 /RR 1 γγ WW 2 /RR 2 γγ WW 3 /RR 3 γγ 1 2 3
8 All d-dimensional GEV distributions may be obtained as γγ G xx = PPPP(max{WW ii RR ii ii } xx) Often the vector WW = WW (dd) is thought of as obtained by sampling from a process WW(ss) with ss RR dd (as in figure on previous slide) Often G is defined on a standardized scale, say γγ = 11, and model then is transformed to the real scale Smith process: WW is a non-random Gaussian density function Brown-Resnick process: WW(ss) = exp{xx ss σσ 2 /2}, with XX a Gaussian process with variance σσ 2 Schlater model
9 The block maxima method As in one dimension: Get observations XX 1,, XX nn of block maxima; assume the observations are i.i.d and follow a GEV distribution; use XX 1,, XX nn to estimate the parameters of the GEV distribution Likelihood estimation: Distribution functions often are possible to calculate, but differentiation with respect to the dd component variables often lead to combinatorial explosion of number of terms (because often GG xx = exp{ l(xx)}) Composite likelihood: Use sum of likelihoods of pairs, or triplets, or Stephenson-Tawn likelihoods: in between block maxima and PoT A. C. Davison, S. A. Padoan and M. Ribatet Statistical Modeling of Spatial Extremes, Statistical Science 2012 R. Huser, A. C. Davison, M. G. Genton Likelihood estimators for multivariate extremes, Extremes 2016
10 Copula modelling (i.e. modeling after transforming marginal distributions to uniformity) is often used Popular copulas: logistic = Gumbel; asymmetric logistic; inverted logistic; Gaussian; Husler-Reiss; bilogistic; Dirichlet; distributions with same copula, but different margins: Uniform; Fréchet; Gumbel; and Laplace
11 pointwise 25-year return levels for rainfall (mm) obtained from latent variable and max-stable models. Top and bottom rows: lower and upper bounds of 95% pointwise credible/confidence intervals. Middle row: predictive pointwise posterior mean and pointwise estimates. Left column: latent variable model. Middle column: Another latent variable model. Right column: Extremal t copula model. Copied from Davison, Padoan & Ribatet paper
12 Multivariate peaks over thresholds modelling and likelihood inference Or: News from the exciting exploration of GP-territory Joint with Anna Kiriliouk, Johan Segers, Maud Thomas, Jennifer Wadsworth
13 The philosophy is simple Extreme episodes are often quite different from ordinary everyday behavior, and ordinary behavior then has little to say about extremes, so that only other extreme events give useful information about future extreme events Operationally: choose thresholds uu 1,, uu dd, say that an extreme episode occurs if at least one of the components of the observation YY = (YY 1,, Y dd ) exceeds its threshold, only model times of occurence and undershoots and overshoots, XX = YY uu, of extreme episodes Theory: times of occurrence follows a Poisson process, XX is a generalized Pareto (GP) distribution
14 Generalized Pareto(GP) distributions HH xx = 1 log GG(00) log GG(xx) GG(xx 00), where GG is a GEV cdf with GG 00 > 0
15 The GP distributions are the limit distributions threshold excesses: Let XX~ FF. If there exist continuous threshold and scaling functions uu tt and ss tt > 0, with FF uu tt < 1 and FF uu tt 1 as tt, such that PPPP ss tt 1 XX uu tt xx XX uu tt dd HH xx, as nn, where HH has non-degenerate margins, then HH is a GP distribution. Conversely all GP distributions can be obtained in this way. maxima of FF are in the domain of attraction of a GEV cdf if and only if threshold excesses of FF are in the threshold domain of attraction of a GP cdf HHH
16 The GP distributions are the threshold-stable distribution: Let XX~HH, with HH nondegenerate. If there exist continuous threshold and scaling functions uu tt, ss tt 00 and with FF uu tt < 1 and FF uu tt 1 as tt, such that PPPP ss tt 1 XX uu tt xx XX > uu tt = HH xx, for all tt 1, then HH is a GP distribution. Conversely all GP distributions has this property.
17 Closed under Properties of the class of GP distributions Scaling Increase of level Taking limits Mixing Going to conditional margins Not closed under Location changes Going to (unconditional) margins Conditional 1-d marginals GP: HH jj + xx = γγ jj σσ jj xx (= 1 ee xx/σσ jj if γγ jj = 0) 1/γγ jj
18 Multivariate PoT modelling Should ideally contain three components: 1) A model for the behavior of extreme episodes in the physical world 2) Understanding of how conditioning on threshold exceedance changes the distribution obtained in 1) 3) The possibility to incorporate trends in the models ( likelihood inference) Inference should be for the original observations, and not after (approximate) transformation of margins of observations to some standard form, say standard Pareto or standard GP
19 Think of XX 1, XX 2,, XX n as points in RR dd. Subtract aa nn and divide by bb nn to get a point process with points XX ii aa nn bb nn ηη; ii = 1, nn. It converges as nn iff XX DD(GG), to a limit point process xx 2 xx 1 {WW ii RR ii γγ } where (in sequel think of γγ > 00) WW ii are i.i.d. copies of some ( any ) positive stochastic dd-dimensional vector and the RR ii are points of a unit rate Poisson process on RR + WW 1 /RR 1 γγ WW 2 /RR 2 γγ WW 3 /RR 3 γγ 1 2 3
20 A typical point (positive shape) RR TT γγ RR an arbitrary dd-dimensional random vector, subject to EE RR jj <, and RR jj 0 if γγ jj > 0 TT improper random variable, has density = 1 with respect to Lebesgue measure on RR + RR and TT are mutually independent If γγ jj > 0 then TT γγ jj has density xx 1 γγ jj 1 with respect to Lebesgue measure on RR + If γγ jj = 0 then RR jj TT γγ jj is defined to mean RR jj log TT and log TT has density ee xx with respect to Lebesgue measure on RR
21 Conditioning on threshold exceedance FF d.f. of RR in typical point RR TT γγ. XX uu = XX uu, PPPP XX uu xx XX uu 00 = PPrr XX uu xx PPrr XX uu xx 00 PPrr XX uu 00 Set XX = WW RR γγ and uu = σσ γγ and condition (formally) on TT uu 2 xx 2 1 not here 1 uu 1 xx 1 (R) HH xx = 0 FF tt γγ xx+ σσ ddtt FF tt γγ xx 00+ σσ γγ γγ 0 FF(tt γγσσ γγ )ddtt is a GP and all GP-s can be obtained this way dddd
22 (R) HH rr xx = Three representations 0 FFRR tt γγ xx+ σσ ddtt FF γγ RR tt γγ xx 00+ σσ dddd γγ 0 FF RR (tt γγσσ γγ )ddtt is a GP and all GP-s can be obtained this way (R ss ) HH ss xx = 0 FFSS 1 γγ log xx+σσ γγ +log tt ddtt FF SS( 1 γγ log xx 00+σσ γγ 0 FF SS (log tt)ddtt is a GP and all GP-s can be obtained this way +log tt) dddd (R ss ) HH SS xx = 0 1 FF SS = FFSS 0 1 γγ log xx + σσ γγ + log tt dddd 1 γγ log xx + σσ γγ tt ee tt dddd for SS = SS max SS jj and all GP-s can be obtained this way Ferreira de Haan (2014) representation
23 Densities Explicit formulas for densities for all three representations Formula for (R ss ) density simplest!! Censored likelihood contributions computable
24 Estimation Likelihood for (R): NN dd 0 tt jj=1 γγ jj ff tt γγ xx ii + γγ σσ ii=1 0 FF(tt γγγγ σσ )ddtt dddd One-dimensional integrals, so relatively easy to compute if ff, FF are easy to compute. In some cases integrals can be computed analytically Censoring makes computation harder Parenthesis: Coefficient of asymptotic dependence: Pr(FF 1 (XX 1 ) > qq,, FF dd (XX dd ) > qq) χχ dd : = lim. qq 1 1 qq GEV and GP modelling aimed at asymptotic dependence, i.e. χχ > 0 for some subset of components. Substantial literature (Heffernan, Tawn, Resnick, Wadsworth, ) on modelling under asymptotic independence, when χχ = 0
25 Parametrization and identifiability Understanding is on the (R) scale Choice between (R), (R ss ), and (R ss ) a choice of parametrization σσ and γγ parameters of conditional margins Replacing RR by ZZ γγ RR doesn t change HH rr. Similar properties of HH ss and HH ss (R ss ) and (R ss ) continuous at γγ ii = 0. Not so for (R) Much left to explore
26 Probabilities and conditional probabilities Probabilities of general events given by same formula as HH rr Conditional densities same type of formulas Conditional probabilities same type of formulas (prediction) if γγ 1 = = γγ dd =: γγ then weighted sums of components, conditioned to be positive, are also GP
27 If γγ 1 = = γγ dd γγ then weighted sums of components, conditioned to be positive, are also GP RR Pf. If σσ = TT γγ γγ (RR 1 σσ 1, RR 2 σσ 2 ) then the sum of the TT γγ γγ TT γγ γγ components equals RR 1+RR 2 σσ 1+σσ 2 =: RR σσ, and by the TT γγ γγ TT γγ γγ onedimensional (R) representation, the conditional distribution of RR σσ given that it is positive is a one-dimensional GP TT γγ γγ distribution with parameters σσ and γγ
28 Simulation 1. Direct for (R ss ) 2. Change of measure + rejection sampling for (R) and (R ss ) 3. Change of measure + MCMC for (R) and (R ss ) 4. Approximate method for (R) and (R ss )
29 Before Gudrun After Gudrun SEK 2,8 billion loss (LF) Windstorm Gudrun, January % of total loss Remeber: analysis based on data : 1% chance biggest damage next 15-year period will be more than SEK 2,5 billion??? Bivariate logistic PoT analysis basered on data : 10% chance biggest damage next 15-year period will be more than SEK 7 billion (this time we didn t miss damage to forrest but did we miss something else?) Brodin & Rootzén, Insurance 2009
30 Landslides Landslide caused by one day with very extreme rainfall, or by two consequtive days with extreme rainfall. Simplest model: XX 1, XX 2, independent positive variables, γγ 1 = γγ 2 = γγ RR TT γγ = XX 1 XX 2 TT γγ, XX 1 + XX 2 TT γγ = XX 1 XX 2 TT γγ, XX 1 XX 2 + XX 1 XX 2 TT γγ = rainiest day, rainest day + rainest adjacent day
31 Slide: Anna Kiriliouk
32 Slides: Jennifer Wadsworth
33 Modelling strategy
34
35
36 Prediction of influenza epidemics in France yy 1 =#cases week 3,, yy 8 =#cases week 11 γγ = 00 (from data), RR ii NN(mm ii, ss ii 2 ), mm 1 = 0 Model: conditional distribution of RR σσ log TT given that RR 1 > σσ log TT Total #cases GP with γγ = 0, σσ = σσ 1 + σσ 8 LR-test of parabolic form for the mm ii, equality of the ss ii 2 #cases per week normal epidemics extreme epidemics
37 Challenges Model construction Parametrization Time series Prediction Asymptotics H. Rootzén, J. Segers, and J.L. Wadsworth Multivariate peaks over thresholds models, ArXive 2016
Modelling 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 informationVariations. ECE 6540, Lecture 02 Multivariate Random Variables & Linear Algebra
Variations ECE 6540, Lecture 02 Multivariate Random Variables & Linear Algebra Last Time Probability Density Functions Normal Distribution Expectation / Expectation of a function Independence Uncorrelated
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 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 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 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 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 informationStatistics of Extremes
Statistics of Extremes Anthony Davison c 211 http://stat.epfl.ch Multivariate Extremes 19 Componentwise maxima.................................................. 194 Standardization........................................................
More informationMultivariate generalized Pareto distributions
Multivariate generalized Pareto distributions Holger Rootzén and Nader Tajvidi Abstract Statistical inference for extremes has been a subject of intensive research during the past couple of decades. One
More informationLecture 3. STAT161/261 Introduction to Pattern Recognition and Machine Learning Spring 2018 Prof. Allie Fletcher
Lecture 3 STAT161/261 Introduction to Pattern Recognition and Machine Learning Spring 2018 Prof. Allie Fletcher Previous lectures What is machine learning? Objectives of machine learning Supervised and
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 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 informationMultivariate peaks over thresholds models
Extremes DOI 1.17/s1687-17-294-4 Multivariate peaks over thresholds models Holger Rootzén 1 Johan Segers 2 Jennifer L. Wadsworth 3 Received: 18 October 216 / Revised: 2 April 217 / Accepted: 21 April 217
More informationGeneral Strong Polarization
General Strong Polarization Madhu Sudan Harvard University Joint work with Jaroslaw Blasiok (Harvard), Venkatesan Gurswami (CMU), Preetum Nakkiran (Harvard) and Atri Rudra (Buffalo) December 4, 2017 IAS:
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 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 informationClassical Extreme Value Theory - An Introduction
Chapter 1 Classical Extreme Value Theory - An Introduction 1.1 Introduction Asymptotic theory of functions of random variables plays a very important role in modern statistics. The objective of the asymptotic
More 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 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 informationWorksheets for GCSE Mathematics. Algebraic Expressions. Mr Black 's Maths Resources for Teachers GCSE 1-9. Algebra
Worksheets for GCSE Mathematics Algebraic Expressions Mr Black 's Maths Resources for Teachers GCSE 1-9 Algebra Algebraic Expressions Worksheets Contents Differentiated Independent Learning Worksheets
More informationMultivariate Pareto distributions: properties and examples
Multivariate Pareto distributions: properties and examples Ana Ferreira 1, Laurens de Haan 2 1 ISA UTL and CEAUL, Portugal 2 Erasmus Univ Rotterdam and CEAUL EVT2013 Vimeiro, September 8 11 Univariate
More 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 informationGeneral Linear Models. with General Linear Hypothesis Tests and Likelihood Ratio Tests
General Linear Models with General Linear Hypothesis Tests and Likelihood Ratio Tests 1 Background Linear combinations of Normals are Normal XX nn ~ NN μμ, ΣΣ AAAA ~ NN AAμμ, AAAAAA A sum of squared, standardized
More informationBayesian Model Averaging for Multivariate Extreme Values
Bayesian Model Averaging for Multivariate Extreme Values Philippe Naveau naveau@lsce.ipsl.fr Laboratoire des Sciences du Climat et l Environnement (LSCE) Gif-sur-Yvette, France joint work with A. Sabourin
More informationBayesian 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 informationFinancial Econometrics and Volatility Models Extreme Value Theory
Financial Econometrics and Volatility Models Extreme Value Theory Eric Zivot May 3, 2010 1 Lecture Outline Modeling Maxima and Worst Cases The Generalized Extreme Value Distribution Modeling Extremes Over
More informationModels 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 informationNon-Parametric Weighted Tests for Change in Distribution Function
American Journal of Mathematics Statistics 03, 3(3: 57-65 DOI: 0.593/j.ajms.030303.09 Non-Parametric Weighted Tests for Change in Distribution Function Abd-Elnaser S. Abd-Rabou *, Ahmed M. Gad Statistics
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 informationSimulation of Max Stable Processes
Simulation of Max Stable Processes Whitney Huang Department of Statistics Purdue University November 6, 2013 1 / 30 Outline 1 Max-Stable Processes 2 Unconditional Simulations 3 Conditional Simulations
More informationEstimate by the L 2 Norm of a Parameter Poisson Intensity Discontinuous
Research Journal of Mathematics and Statistics 6: -5, 24 ISSN: 242-224, e-issn: 24-755 Maxwell Scientific Organization, 24 Submied: September 8, 23 Accepted: November 23, 23 Published: February 25, 24
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 informationECE 6540, Lecture 06 Sufficient Statistics & Complete Statistics Variations
ECE 6540, Lecture 06 Sufficient Statistics & Complete Statistics Variations Last Time Minimum Variance Unbiased Estimators Sufficient Statistics Proving t = T(x) is sufficient Neyman-Fischer Factorization
More informationGeneral Strong Polarization
General Strong Polarization Madhu Sudan Harvard University Joint work with Jaroslaw Blasiok (Harvard), Venkatesan Gurswami (CMU), Preetum Nakkiran (Harvard) and Atri Rudra (Buffalo) May 1, 018 G.Tech:
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 informationSupport Vector Machines. CSE 4309 Machine Learning Vassilis Athitsos Computer Science and Engineering Department University of Texas at Arlington
Support Vector Machines CSE 4309 Machine Learning Vassilis Athitsos Computer Science and Engineering Department University of Texas at Arlington 1 A Linearly Separable Problem Consider the binary classification
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 informationPeaks over thresholds modelling with multivariate generalized Pareto distributions
Peaks over thresholds modelling with multivariate generalized Pareto distributions Anna Kiriliouk Johan Segers Université catholique de Louvain Institut de Statistique, Biostatistique et Sciences Actuarielles
More informationExpectation Propagation performs smooth gradient descent GUILLAUME DEHAENE
Expectation Propagation performs smooth gradient descent 1 GUILLAUME DEHAENE In a nutshell Problem: posteriors are uncomputable Solution: parametric approximations 2 But which one should we choose? Laplace?
More informationGradient expansion formalism for generic spin torques
Gradient expansion formalism for generic spin torques Atsuo Shitade RIKEN Center for Emergent Matter Science Atsuo Shitade, arxiv:1708.03424. Outline 1. Spintronics a. Magnetoresistance and spin torques
More informationMultiple Regression Analysis: Inference ECONOMETRICS (ECON 360) BEN VAN KAMMEN, PHD
Multiple Regression Analysis: Inference ECONOMETRICS (ECON 360) BEN VAN KAMMEN, PHD Introduction When you perform statistical inference, you are primarily doing one of two things: Estimating the boundaries
More informationModule 7 (Lecture 27) RETAINING WALLS
Module 7 (Lecture 27) RETAINING WALLS Topics 1.1 RETAINING WALLS WITH METALLIC STRIP REINFORCEMENT Calculation of Active Horizontal and vertical Pressure Tie Force Factor of Safety Against Tie Failure
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 informationWork, Energy, and Power. Chapter 6 of Essential University Physics, Richard Wolfson, 3 rd Edition
Work, Energy, and Power Chapter 6 of Essential University Physics, Richard Wolfson, 3 rd Edition 1 With the knowledge we got so far, we can handle the situation on the left but not the one on the right.
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 informationLecture 11. Kernel Methods
Lecture 11. Kernel Methods COMP90051 Statistical Machine Learning Semester 2, 2017 Lecturer: Andrey Kan Copyright: University of Melbourne This lecture The kernel trick Efficient computation of a dot product
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 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 informationLecture 6. Notes on Linear Algebra. Perceptron
Lecture 6. Notes on Linear Algebra. Perceptron COMP90051 Statistical Machine Learning Semester 2, 2017 Lecturer: Andrey Kan Copyright: University of Melbourne This lecture Notes on linear algebra Vectors
More informationLecture 3. Linear Regression
Lecture 3. Linear Regression COMP90051 Statistical Machine Learning Semester 2, 2017 Lecturer: Andrey Kan Copyright: University of Melbourne Weeks 2 to 8 inclusive Lecturer: Andrey Kan MS: Moscow, PhD:
More informationOn The Cauchy Problem For Some Parabolic Fractional Partial Differential Equations With Time Delays
Journal of Mathematics and System Science 6 (216) 194-199 doi: 1.17265/2159-5291/216.5.3 D DAVID PUBLISHING On The Cauchy Problem For Some Parabolic Fractional Partial Differential Equations With Time
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 informationTerms of Use. Copyright Embark on the Journey
Terms of Use All rights reserved. No part of this packet may be reproduced, stored in a retrieval system, or transmitted in any form by any means - electronic, mechanical, photo-copies, recording, or otherwise
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 informationMath 171 Spring 2017 Final Exam. Problem Worth
Math 171 Spring 2017 Final Exam Problem 1 2 3 4 5 6 7 8 9 10 11 Worth 9 6 6 5 9 8 5 8 8 8 10 12 13 14 15 16 17 18 19 20 21 22 Total 8 5 5 6 6 8 6 6 6 6 6 150 Last Name: First Name: Student ID: Section:
More informationCS Lecture 8 & 9. Lagrange Multipliers & Varitional Bounds
CS 6347 Lecture 8 & 9 Lagrange Multipliers & Varitional Bounds General Optimization subject to: min ff 0() R nn ff ii 0, h ii = 0, ii = 1,, mm ii = 1,, pp 2 General Optimization subject to: min ff 0()
More informationWorksheets for GCSE Mathematics. Quadratics. mr-mathematics.com Maths Resources for Teachers. Algebra
Worksheets for GCSE Mathematics Quadratics mr-mathematics.com Maths Resources for Teachers Algebra Quadratics Worksheets Contents Differentiated Independent Learning Worksheets Solving x + bx + c by factorisation
More informationA new procedure for sensitivity testing with two stress factors
A new procedure for sensitivity testing with two stress factors C.F. Jeff Wu Georgia Institute of Technology Sensitivity testing : problem formulation. Review of the 3pod (3-phase optimal design) procedure
More informationModule 7 (Lecture 25) RETAINING WALLS
Module 7 (Lecture 25) RETAINING WALLS Topics Check for Bearing Capacity Failure Example Factor of Safety Against Overturning Factor of Safety Against Sliding Factor of Safety Against Bearing Capacity Failure
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 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 informationRevision : Thermodynamics
Revision : Thermodynamics Formula sheet Formula sheet Formula sheet Thermodynamics key facts (1/9) Heat is an energy [measured in JJ] which flows from high to low temperature When two bodies are in thermal
More informationChapter 22 : Electric potential
Chapter 22 : Electric potential What is electric potential? How does it relate to potential energy? How does it relate to electric field? Some simple applications What does it mean when it says 1.5 Volts
More informationTranslator Design Lecture 16 Constructing SLR Parsing Tables
SLR Simple LR An LR(0) item (item for short) of a grammar G is a production of G with a dot at some position of the right se. Thus, production A XYZ yields the four items AA XXXXXX AA XX YYYY AA XXXX ZZ
More informationCDS 101/110: Lecture 6.2 Transfer Functions
CDS 11/11: Lecture 6.2 Transfer Functions November 2, 216 Goals: Continued study of Transfer functions Review Laplace Transform Block Diagram Algebra Bode Plot Intro Reading: Åström and Murray, Feedback
More informationUncertain Compression & Graph Coloring. Madhu Sudan Harvard
Uncertain Compression & Graph Coloring Madhu Sudan Harvard Based on joint works with: (1) Adam Kalai (MSR), Sanjeev Khanna (U.Penn), Brendan Juba (WUStL) (2) Elad Haramaty (Harvard) (3) Badih Ghazi (MIT),
More informationPHY103A: Lecture # 9
Semester II, 2017-18 Department of Physics, IIT Kanpur PHY103A: Lecture # 9 (Text Book: Intro to Electrodynamics by Griffiths, 3 rd Ed.) Anand Kumar Jha 20-Jan-2018 Summary of Lecture # 8: Force per unit
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 informationJasmin Smajic1, Christian Hafner2, Jürg Leuthold2, March 23, 2015
Jasmin Smajic, Christian Hafner 2, Jürg Leuthold 2, March 23, 205 Time Domain Finite Element Method (TD FEM): Continuous and Discontinuous Galerkin (DG-FEM) HSR - University of Applied Sciences of Eastern
More informationA brief summary of the chapters and sections that we cover in Math 141. Chapter 2 Matrices
A brief summary of the chapters and sections that we cover in Math 141 Chapter 2 Matrices Section 1 Systems of Linear Equations with Unique Solutions 1. A system of linear equations is simply a collection
More informationReview for Exam Hyunse Yoon, Ph.D. Adjunct Assistant Professor Department of Mechanical Engineering, University of Iowa
Review for Exam2 11. 13. 2015 Hyunse Yoon, Ph.D. Adjunct Assistant Professor Department of Mechanical Engineering, University of Iowa Assistant Research Scientist IIHR-Hydroscience & Engineering, University
More informationMultivariate generalized Pareto distributions
Bernoulli 12(5), 2006, 917 930 Multivariate generalized Pareto distributions HOLGER ROOTZÉN 1 and NADER TAJVIDI 2 1 Chalmers University of Technology, S-412 96 Göteborg, Sweden. E-mail rootzen@math.chalmers.se
More 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 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 informationRandom Vectors Part A
Random Vectors Part A Page 1 Outline Random Vectors Measurement of Dependence: Covariance Other Methods of Representing Dependence Set of Joint Distributions Copulas Common Random Vectors Functions of
More informationHaar Basis Wavelets and Morlet Wavelets
Haar Basis Wavelets and Morlet Wavelets September 9 th, 05 Professor Davi Geiger. The Haar transform, which is one of the earliest transform functions proposed, was proposed in 90 by a Hungarian mathematician
More informationChapter 5: Spectral Domain From: The Handbook of Spatial Statistics. Dr. Montserrat Fuentes and Dr. Brian Reich Prepared by: Amanda Bell
Chapter 5: Spectral Domain From: The Handbook of Spatial Statistics Dr. Montserrat Fuentes and Dr. Brian Reich Prepared by: Amanda Bell Background Benefits of Spectral Analysis Type of data Basic Idea
More informationGeneral Strong Polarization
General Strong Polarization Madhu Sudan Harvard University Joint work with Jaroslaw Blasiok (Harvard), Venkatesan Guruswami (CMU), Preetum Nakkiran (Harvard) and Atri Rudra (Buffalo) Oct. 5, 018 Stanford:
More informationVariable Selection in Predictive Regressions
Variable Selection in Predictive Regressions Alessandro Stringhi Advanced Financial Econometrics III Winter/Spring 2018 Overview This chapter considers linear models for explaining a scalar variable when
More informationQuestion 1 carries a weight of 25%; question 2 carries 25%; question 3 carries 20%; and question 4 carries 30%.
UNIVERSITY OF EAST ANGLIA School of Economics Main Series PGT Examination 2017-18 FINANCIAL ECONOMETRIC THEORY ECO-7024A Time allowed: 2 HOURS Answer ALL FOUR questions. Question 1 carries a weight of
More informationFermi Surfaces and their Geometries
Fermi Surfaces and their Geometries Didier Ndengeyintwali Physics Department, Drexel University, Philadelphia, Pennsylvania 19104, USA (Dated: May 17, 2010) 1. Introduction The Pauli exclusion principle
More informationStatistics of Extremes
Statistics of Extremes Anthony Davison c 29 http://stat.epfl.ch Multivariate Extremes 184 Europe............................................................. 185 Componentwise maxima..................................................
More informationOn the spatial dependence of extreme ocean storm seas
Accepted by Ocean Engineering, August 2017 On the spatial dependence of extreme ocean storm seas Emma Ross a, Monika Kereszturi b, Mirrelijn van Nee c, David Randell d, Philip Jonathan a, a Shell Projects
More information1. The graph of a function f is given above. Answer the question: a. Find the value(s) of x where f is not differentiable. Ans: x = 4, x = 3, x = 2,
1. The graph of a function f is given above. Answer the question: a. Find the value(s) of x where f is not differentiable. x = 4, x = 3, x = 2, x = 1, x = 1, x = 2, x = 3, x = 4, x = 5 b. Find the value(s)
More informationGeneral Strong Polarization
General Strong Polarization Madhu Sudan Harvard University Joint work with Jaroslaw Blasiok (Harvard), Venkatesan Guruswami (CMU), Preetum Nakkiran (Harvard) and Atri Rudra (Buffalo) Oct. 8, 018 Berkeley:
More informationExam Programme VWO Mathematics A
Exam Programme VWO Mathematics A The exam. The exam programme recognizes the following domains: Domain A Domain B Domain C Domain D Domain E Mathematical skills Algebra and systematic counting Relationships
More informationSpatial data analysis
Spatial data analysis Global Health Sciences Global Health Group Data Science Africa 2016 Ricardo Andrade Outline The geographic context Geostatistics Non-linear models Discrete processes Time interactions
More informationSolving Fuzzy Nonlinear Equations by a General Iterative Method
2062062062062060 Journal of Uncertain Systems Vol.4, No.3, pp.206-25, 200 Online at: www.jus.org.uk Solving Fuzzy Nonlinear Equations by a General Iterative Method Anjeli Garg, S.R. Singh * Department
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 information" = Y(#,$) % R(r) = 1 4& % " = Y(#,$) % R(r) = Recitation Problems: Week 4. a. 5 B, b. 6. , Ne Mg + 15 P 2+ c. 23 V,
Recitation Problems: Week 4 1. Which of the following combinations of quantum numbers are allowed for an electron in a one-electron atom: n l m l m s 2 2 1! 3 1 0 -! 5 1 2! 4-1 0! 3 2 1 0 2 0 0 -! 7 2-2!
More informationComparing Different Estimators of Reliability Function for Proposed Probability Distribution
American Journal of Mathematics and Statistics 21, (2): 84-94 DOI: 192/j.ajms.212. Comparing Different Estimators of Reliability Function for Proposed Probability Distribution Dhwyia S. Hassun 1, Nathie
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 informationOn the Estimation and Application of Max-Stable Processes
On the Estimation and Application of Max-Stable Processes Zhengjun Zhang Department of Statistics University of Wisconsin Madison, WI 53706 USA Richard L Smith Department of Statistics University of North
More 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 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 informationMax-stable Processes for Threshold Exceedances in Spatial Extremes
Max-stable Processes for Threshold Exceedances in Spatial Extremes Soyoung Jeon A dissertation submitted to the faculty of the University of North Carolina at Chapel Hill in partial fulfillment of the
More informationData Preprocessing. Jilles Vreeken IRDM 15/ Oct 2015
Data Preprocessing Jilles Vreeken 22 Oct 2015 So, how do you pronounce Jilles Yill-less Vreeken Fray-can Okay, now we can talk. Questions of the day How do we preprocess data before we can extract anything
More informationSECTION 7: FAULT ANALYSIS. ESE 470 Energy Distribution Systems
SECTION 7: FAULT ANALYSIS ESE 470 Energy Distribution Systems 2 Introduction Power System Faults 3 Faults in three-phase power systems are short circuits Line-to-ground Line-to-line Result in the flow
More informationSummary. Forward modeling of the magnetic fields using rectangular cells. Introduction
Yue Zhu, University of Utah, Michael S. Zhdanov, University of Utah, TechnoImaging and MIPT, and Martin Čuma, University of Utah and TechnoImaging Summary Conventional 3D magnetic inversion methods are
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