Conditional simulations of max-stable processes
|
|
- Janice Quinn
- 5 years ago
- Views:
Transcription
1 Conditional simulations of max-stable processes C. Dombry, F. Éyi-Minko, M. Ribatet Laboratoire de Mathématiques et Application, Université de Poitiers Institut de Mathématiques et de Modélisation, Université Montpellier 2 Conditional simulations of max-stable processes Mathieu Ribatet 1 / 24
2 Goal Max-stable processes are widely used for modelling spatial extremes since they arise as the only possible (non degenerate) limit of pointwise maxima over independent replicates, i.e., max i=1,...,n X i (x) b n (x) a n (x) Z (x), n, x X R d, for some normalizing functions a n > 0 and b n and where X i are independent copies of a stochastic process X. Conditional simulations of max-stable processes Mathieu Ribatet 2 / 24
3 Goal Max-stable processes are widely used for modelling spatial extremes since they arise as the only possible (non degenerate) limit of pointwise maxima over independent replicates, i.e., max i=1,...,n X i (x) b n (x) a n (x) Z (x), n, x X R d, for some normalizing functions a n > 0 and b n and where X i are independent copies of a stochastic process X. Can we get a procedure for conditional simulations of max-stable processes (with continuous spectral measure)? Conditional simulations of max-stable processes Mathieu Ribatet 2 / 24
4 Setup Recall that any max-stable process has the spectral characterization where Z ( )=max i 1 ζ i Y i ( ), Y i ( ) are independent copies of a non negative stochastic process such that E[Y (x)]=1 for all x X ; {ζ i } i 1 are the points of a Poisson process on (0, ) with intensity dλ(ζ)=ζ 2 dζ. Given a study region X R d, we want to sample from Z ( ) {Z (x 1 )=z 1,..., Z (x k )=z k }, for some z 1,..., z k > 0 and k conditioning locations x 1,..., x k X. Conditional simulations of max-stable processes Mathieu Ribatet 3 / 24
5 Random partitions Sampling scheme of max-stable processes Conditional simulations of max-stable processes Mathieu Ribatet 4 / 24
6 Decomposition of Φ Random partitions Sampling scheme Let Φ a point process on (0, ) k whose atoms are ϕ i (x)=ζ i Y i (x), x=(x 1,..., x k ). Consider the two following point processes Φ = { ϕ Φ: ϕ(x i )<z i, for all i {1,...,k} }, (sub-extremal functions) Φ + = { ϕ Φ: ϕ(x i )=z i, for some i {1,...,k} }.(extremal functions) Clearly Φ=Φ Φ + and Φ + = {ϕ + 1,...,ϕ+ k }={ϕ+ 1,...,ϕ+ l }, a.s. (1 l k). Conditional simulations of max-stable processes Mathieu Ribatet 5 / 24
7 Decomposition of Φ Random partitions Sampling scheme Let Φ a point process on (0, ) k whose atoms are ϕ i (x)=ζ i Y i (x), x=(x 1,..., x k ). Consider the two following point processes Φ = { ϕ Φ: ϕ(x i )<z i, for all i {1,...,k} }, (sub-extremal functions) Φ + = { ϕ Φ: ϕ(x i )=z i, for some i {1,...,k} }.(extremal functions) Clearly Φ=Φ Φ + and Φ + = {ϕ + 1,...,ϕ+ k }={ϕ+ 1,...,ϕ+ l }, a.s. (1 l k). Key point #1: Conditionally on Z (x)=z, Φ and Φ + are independent. Conditional simulations of max-stable processes Mathieu Ribatet 5 / 24
8 Conditional intensity function Random partitions Sampling scheme Z (x)=max i 1 ζ i Y i (x)=max i 1 ϕ i (x) The Poisson point process {ϕ i (x)} i 1 has intensity measure Λ x (A)= 0 Pr{ζY (x) A}ζ 2 dζ, Borel set A R k. We assume that Φ is regular, i.e., Λ x (dz)=λ x (z)dz, for all x X k. Conditional simulations of max-stable processes Mathieu Ribatet 6 / 24
9 Conditional intensity function Random partitions Sampling scheme Z (x)=max i 1 ζ i Y i (x)=max i 1 ϕ i (x) The Poisson point process {ϕ i (x)} i 1 has intensity measure Λ x (A)= 0 Pr{ζY (x) A}ζ 2 dζ, Borel set A R k. We assume that Φ is regular, i.e., Λ x (dz)=λ x (z)dz, for all x X k. λ s x,z (u)= λ (s,x)(u, z) λ x (z) Key point #2: The conditional intensity function is the (regular) conditional distribution of Z (x) if we integrate w.r.t. all possible partitions of x. But not that of Z ( )!!! Conditional simulations of max-stable processes Mathieu Ribatet 6 / 24
10 Random partitions? Random partitions Sampling scheme Z(x) x 2 x 5 x 1 x 3 x 4 x Conditional simulations of max-stable processes Mathieu Ribatet 7 / 24
11 Random partitions? Random partitions Sampling scheme Z(x) ϕ 1 + x 2 x 5 x 1 x 3 x 4 x Conditional simulations of max-stable processes Mathieu Ribatet 7 / 24
12 Random partitions? Random partitions Sampling scheme Z(x) ϕ 1 + ϕ 2 + x 2 x 5 x 1 x 3 x 4 x Conditional simulations of max-stable processes Mathieu Ribatet 7 / 24
13 Random partitions? Random partitions Sampling scheme Z(x) ϕ 1 + ϕ 2 + ϕ 3 + x 2 x 5 x 1 x 3 x 4 x Conditional simulations of max-stable processes Mathieu Ribatet 7 / 24
14 Random partitions? Random partitions Sampling scheme Z(x) ϕ 1 + ϕ 2 + ϕ 3 + ϕ 4 + x 2 x 5 x 1 x 3 x 4 x Conditional simulations of max-stable processes Mathieu Ribatet 7 / 24
15 Random partitions? Random partitions Sampling scheme Z(x) ϕ 1 + ϕ 2 + ϕ 3 + ϕ 4 + x 2 x 5 x 1 x 3 x 4 Here the set {x 1,..., x 5 } is partitioned into ({x 1, x 3 },{x 2 },{x 4 },{x 5 }) x Conditional simulations of max-stable processes Mathieu Ribatet 7 / 24
16 Random partitions? Random partitions Sampling scheme Z(x) ϕ 1 + ϕ 2 + ϕ 3 + ϕ 4 + x 2 x 5 x 1 x 3 x 4 Here the set {x 1,..., x 5 } is partitioned into ({x 1, x 3 },{x 2 },{x 4 },{x 5 }) The hitting bounds {z i } i=1,...,k might be reached by several extremal functions, i.e., Φ + = {ϕ + 1,...,ϕ+ k }={ϕ+ 1,...,ϕ+ l } a.s., 1 l k. So we need to take into account all possible ways these hitting bounds are reached: the hitting scenarios x Conditional simulations of max-stable processes Mathieu Ribatet 7 / 24
17 Why should we bother about Φ? Random partitions Sampling scheme Z(x) x 2 x 5 x 1 x 3 x 4 x Conditional simulations of max-stable processes Mathieu Ribatet 8 / 24
18 Why should we bother about Φ? Random partitions Sampling scheme Z(x) max Φ + x 2 x 5 x 1 x 3 x 4 x Conditional simulations of max-stable processes Mathieu Ribatet 8 / 24
19 Why should we bother about Φ? Random partitions Sampling scheme Z(x) max Φ + max Φ x 2 x 5 x 1 x 3 x 4 x Conditional simulations of max-stable processes Mathieu Ribatet 8 / 24
20 Why should we bother about Φ? Random partitions Sampling scheme Z(x) max Φ + max Φ x 2 x 5 x 1 x 3 x 4 x Conditional simulations of max-stable processes Mathieu Ribatet 8 / 24
21 Why should we bother about Φ? Random partitions Sampling scheme Z(x) max Φ + max Φ x 2 x 5 x 1 x 3 x 4 The atoms of Φ + are only of interest if we restrict our attention to the conditioning points x; But most often one would like to get realizations at s x. The atoms of Φ are needed since it is likely that maxφ (s)>maxφ + (s)! x Conditional simulations of max-stable processes Mathieu Ribatet 8 / 24
22 A three step procedure Random partitions Sampling scheme The above key points suggest a three step sampling scheme: Step 1 Draw a random partition τ, i.e., a hitting scenario; Step 2 Given τ of size l, draw the extremal functions ϕ + 1,...,ϕ+ l independently; Step 3 Independently from Steps 1 & 2, draw the sub-extremal functions ϕ i, i 1. Conditional simulations of max-stable processes Mathieu Ribatet 9 / 24
23 Step 1: The random partitions Random partitions Sampling scheme Let P k the set of all possible partitions of the set {x 1,..., x k }. Draw a random partition τ P k with distribution π x (z,τ)= 1 C (x, z) τ j=1 λ xτ j (z τj ) }{{} density that some bounds are reached, i.e., the zτ j {u<z τ c j } λ xτ c x τ j j,z τ j (u)du, }{{} probability to lie below the remaining bounds, i.e., below the z τ c j where the normalization constant C (x, z) is given by C (x, z)= θ P k θ j=1 λ xθj (z θj ) λ xθ c x θj,z θj (u)du, {u<z θ c } j j and τ is the size of the partition τ. Conditional simulations of max-stable processes Mathieu Ribatet 10 / 24
24 Step 2: The extremal functions Random partitions Sampling scheme Given τ=(τ 1,...,τ l ), draw l independent random vectors ϕ + 1 (s),...,ϕ+ (s) from the distribution l { } [ ] Pr ϕ + j (s) dv j = 1 C j 1 {u<zτ c }λ (s,xτ c ) x τ j j j,z τ j (v j, u) du }{{} density of ϕ Φ such that ϕ(x τj )= z τj dv j, where 1 { } is the indicator function and C j = 1 {u<zτ c j }λ (s,xτ c j ) x τ j,z τ j (v j, u)dudv j. Define the random vector Z + (s)= max j=1,...,l ϕ+ j (s), s Xm. Conditional simulations of max-stable processes Mathieu Ribatet 11 / 24
25 Step 3: The sub-extremal functions Random partitions Sampling scheme Independently draw {ζ i } i 1 a Poisson point process on (0, ) with intensity ζ 2 dζ and {Y i ( )} i 1 independent copies of Y ( ) Define the random vector Z (s)=max i 1 ζ i Y i (s)1 {ζi Y i (x)<z}, s X m. Conditional simulations of max-stable processes Mathieu Ribatet 12 / 24
26 Step 3: The sub-extremal functions Random partitions Sampling scheme Independently draw {ζ i } i 1 a Poisson point process on (0, ) with intensity ζ 2 dζ and {Y i ( )} i 1 independent copies of Y ( ) Define the random vector Z (s)=max i 1 ζ i Y i (s)1 {ζi Y i (x)<z}, s X m. Then provided Φ is regular, the random vector Z (s)=max { Z + (s), Z (s) } follows the conditional distribution of Z (s) given Z (x) = z. Conditional simulations of max-stable processes Mathieu Ribatet 12 / 24
27 As an aside Random partitions Sampling scheme The conditional cumulative distribution function is { } τ Pr[Z (s) a, Z (x) z] Pr {Z (s) a Z (x)=z}= π x (z,τ) F τ, j (a), τ P k j=1 Pr[Z (x) z] }{{}}{{} Steps 1 & 2 Step 3 where F τ, j (a)= {y<z τ c,u<a} λ (s,x τ c ) x τ j j j,z τ j (u, y)dydu {y<z τ c } λ. t τ c x τ j j j,z τ j (y)dy Conditional simulations of max-stable processes Mathieu Ribatet 13 / 24
28 As an aside Random partitions Sampling scheme The conditional cumulative distribution function is { } τ Pr[Z (s) a, Z (x) z] Pr {Z (s) a Z (x)=z}= π x (z,τ) F τ, j (a), τ P k j=1 Pr[Z (x) z] }{{}}{{} Steps 1 & 2 Step 3 where F τ, j (a)= {y<z τ c,u<a} λ (s,x τ c ) x τ j j j,z τ j (u, y)dydu {y<z τ c } λ. t τ c x τ j j j,z τ j (y)dy Remark. It is clear that Z ( ) {Z (x) = z} is not max-stable. Conditional simulations of max-stable processes Mathieu Ribatet 13 / 24
29 Brown Resnick processes Random partitions Sampling scheme Example 1. If Z is a Brown Resnick process, i.e., Z (x)=max i 1 ζ i exp{ε i (x) γ(x)}, x X, then the intensity function is ( λ x (z)= C x exp 1 ) k 2 logzt Q x logz+l x logz z 1 i, z (0, ) k, i=1 and the conditional intensity function is λ s x,z (u)=(2π) m/2 Σ s x 1/2 exp { 1 } m 2 (logu µ s x,z) T Σ 1 s x (logu µ s x,z) i.e., the extremal functions are log-normal processes. i=1 u 1 i, Conditional simulations of max-stable processes Mathieu Ribatet 14 / 24
30 Schlather processes Random partitions Sampling scheme Example 2. If Z is a Schlather process, i.e., Z (x)= 2πmax i 1 ζ i max{0,ε i (x)}, x X, then the intensity function is ( ) k+ 1 λ x (z)=π (k 1)/2 Σ x 1/2 a x (z) (k+1)/2 Γ, z R k, 2 where a x (z)=z T Σ 1 x z, and the conditional intensity function is λ s x,z (u)=π m/2 (k+ 1) m/2 Σ 1/2 { 1+ (u µ)t Σ 1 (u µ) k+ 1 i.e., the extremal functions are Student processes. } ) (m+k+1)/2 Γ( m+k+1 Γ 2 ( k+1 2 ), Conditional simulations of max-stable processes Mathieu Ribatet 15 / 24
31 Computational burden Full conditional The full conditional are nice! 2. Markov chain Monte Carlo sampler (for Step 1) Conditional simulations of max-stable processes Mathieu Ribatet 16 / 24
32 Do you recognize these numbers? Computational burden Full conditional The full conditional are nice! Conditional simulations of max-stable processes Mathieu Ribatet 17 / 24
33 Do you recognize these numbers? Computational burden Full conditional The full conditional are nice! These are the first 20 Bell numbers. Remark. Recall that Bell(k) is the number of partitions of a set with k elements. Hence with our terminology we have # hitting scenarios=card(p k )=Bell(k). Conditional simulations of max-stable processes Mathieu Ribatet 17 / 24
34 Computational burden Computational burden Full conditional The full conditional are nice! In Step 1, we need to sample from a discrete distribution whose state space is P k. Combinatorial explosion Hence we cannot compute C (x, z) in π x (z,τ)= 1 C (x, z) τ j=1 λ xτ j (z τj ) λ xτ c x τ {u<z τ c } j j,z τ j (u)du. j Conditional simulations of max-stable processes Mathieu Ribatet 18 / 24
35 Computational burden Computational burden Full conditional The full conditional are nice! In Step 1, we need to sample from a discrete distribution whose state space is P k. Combinatorial explosion Hence we cannot compute C (x, z) in π x (z,τ)= 1 C (x, z) τ j=1 λ xτ j (z τj ) λ xτ c x τ {u<z τ c } j j,z τ j (u)du. j Use of MCMC samplers to sample from the target π x (z, ). Remark. We will use a Gibbs sampler since the full conditional are especially convenient. Conditional simulations of max-stable processes Mathieu Ribatet 18 / 24
36 Full conditional Computational burden Full conditional The full conditional are nice! Want to sample from Pr[θ θ j = τ j ], θ π x (z, ) where τ j is the restriction of τ P k to the set {x 1,..., x k } \ {x j }. The number of possible states for θ is b + = { l if {x j } is a partitioning set of τ, l+1 otherwise. Example 3. For τ=({x 1, x 2 },{x 3 }) we have l=2 and Restriction Possible states τ 2 = τ 2 τ 3 = τ 3 ({x 1, x 2 },{x 3 }) ({x 1, x 2 },{x 3 }) ({x 1 },{x 2, x 3 }) ({x 1, x 2, x 3 }) ({x 1 },{x 2 },{x 3 }) Conditional simulations of max-stable processes Mathieu Ribatet 19 / 24
37 The full conditional are nice! Computational burden Full conditional The full conditional are nice! For all τ P k such that τ j = τ j, Pr[θ= τ π x (z,τ ) θ j = τ j ]= π x (z, τ)1 { τ j =τ j } τ P k where w τ, j = λ xτ j (z τj ) {u<z τ c } λ x τ c x τ j j j,z τ j (u)du. τ j=1 w τ, j τ j=1 w τ, j, In particular at most 4 weights w, need to be evaluated and the Gibbs sampler is especially convenient! Conditional simulations of max-stable processes Mathieu Ribatet 20 / 24
38 Conditional simulations of max-stable processes Mathieu Ribatet 21 / 24
39 Precipitation around Zurich in Figure 1: From left to right, maps on a grid of the pointwise 0 025, 0 5 and sample quantiles for rainfall (mm) obtained from conditional simulations of Brown Resnick processes having semi variogram γ(h)=(h/38) The rightmost panel plots the ratio of the width of the pointwise confidence intervals with and without taking estimation uncertainties into account. The squares show the conditional locations. Conditional simulations of max-stable processes Mathieu Ribatet 22 / 24
40 Temperature anomalies for the 2003 European heatwave Basel (316) Zurich (556) Santis (2490) Oeschberg (483) Neuchatel (485) Bad Ragaz (496) Engelberg (1035) Bern (565) Chateau d Oex (985) Montreux (405) Montana (1508) Gd St Bernard (2472) (m) Locarno Monti (366) Lugano (273) Davos (1590) Arosa (1840) (km) ( C) Figure 2: Left: Topographical map of Switzerland showing the sites and altitudes in metres above sea level of 16 weather stations for which annual maxima temperature data are available. Right: Map of temperature anomalies ( C), i.e., the difference between the pointwise medians obtained from conditional simulations and unconditional medians estimated from the fitted Schlather process. As expected the largest deviations occur in the plateau region of Switzerland The differences range between 2 5 C and 4 75 C Conditional simulations of max-stable processes Mathieu Ribatet 23 / 24
41 THANK YOU! ANY QUESTIONS? Dombry, C. Éyi-Minko, F. and Ribatet, M. (2012) Conditional simulation of max-stable processes To appear in Biometrika.
Conditional simulation of max-stable processes
Biometrika (1),,, pp. 1 1 3 4 5 6 7 C 7 Biometrika Trust Printed in Great Britain Conditional simulation of ma-stable processes BY C. DOMBRY, F. ÉYI-MINKO, Laboratoire de Mathématiques et Application,
More informationarxiv: v1 [stat.me] 16 Dec 2011
Biometrika (2011),,, pp. 1 15 C 2007 Biometrika Trust Printed in Great Britain Conditional simulations of Brown Resnick processes arxiv:1112.3891v1 [stat.me] 16 Dec 2011 BY C. DOMBRY, F. ÉYI-MINKO, Laboratoire
More informationGeostatistics of Extremes
Geostatistics of Extremes Anthony Davison Joint with Mohammed Mehdi Gholamrezaee, Juliette Blanchet, Simone Padoan, Mathieu Ribatet Funding: Swiss National Science Foundation, ETH Competence Center for
More informationStatistics of Extremes
Statistics of Extremes Anthony Davison c 29 http://stat.epfl.ch Multivariate Extremes 184 Europe............................................................. 185 Componentwise maxima..................................................
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 informationModelling spatial extremes using max-stable processes
Chapter 1 Modelling spatial extremes using max-stable processes Abstract This chapter introduces the theory related to the statistical modelling of spatial extremes. The presented modelling strategy heavily
More informationStatistical modelling of spatial extremes using the SpatialExtremes package
Statistical modelling of spatial extremes using the SpatialExtremes package M. Ribatet University of Montpellier TheSpatialExtremes package Mathieu Ribatet 1 / 39 Rationale for thespatialextremes package
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 informationSoutenance d Habilitation à Diriger des Recherches. Théorie spatiale des extrêmes et propriétés des processus max-stables
Soutenance d Habilitation à Diriger des Recherches Théorie spatiale des extrêmes et propriétés des processus max-stables CLÉMENT DOMBRY Laboratoire de Mathématiques et Applications, Université de Poitiers.
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 informationConditional Simulation of Max-Stable Processes
1 Conditional Simulation of Ma-Stable Processes Clément Dombry UniversitédeFranche-Comté Marco Oesting INRA, AgroParisTech Mathieu Ribatet UniversitéMontpellier2andISFA,UniversitéLyon1 CONTENTS Abstract...
More informationGeostatistics of extremes
468, 581 608 doi:10.1098/rspa.2011.0412 Published online 12 October 2011 Geostatistics of extremes BY A. C. DAVISON 1, * AND M. M. GHOLAMREZAEE 2 1 Ecole Polytechnique Fédérale de Lausanne, EPFL-FSB-MATHAA-STAT,
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 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 informationNeatest and Promptest Manner. E d i t u r ami rul)lihher. FOIt THE CIIILDIIES'. Trifles.
» ~ $ ) 7 x X ) / ( 8 2 X 39 ««x» ««! «! / x? \» «({? «» q «(? (?? x! «? 8? ( z x x q? ) «q q q ) x z x 69 7( X X ( 3»«! ( ~«x ««x ) (» «8 4 X «4 «4 «8 X «x «(» X) ()»» «X «97 X X X 4 ( 86) x) ( ) z z
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 informationSpatial extreme value theory and properties of max-stable processes Poitiers, November 8-10, 2012
Spatial extreme value theory and properties of max-stable processes Poitiers, November 8-10, 2012 November 8, 2012 15:00 Clement Dombry Habilitation thesis defense (in french) 17:00 Snack buet November
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 informationEXACT SIMULATION OF BROWN-RESNICK RANDOM FIELDS
EXACT SIMULATION OF BROWN-RESNICK RANDOM FIELDS A.B. DIEKER AND T. MIKOSCH Abstract. We propose an exact simulation method for Brown-Resnick random fields, building on new representations for these stationary
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 information' Liberty and Umou Ono and Inseparablo "
3 5? #< q 8 2 / / ) 9 ) 2 ) > < _ / ] > ) 2 ) ) 5 > x > [ < > < ) > _ ] ]? <
More informationMANY BILLS OF CONCERN TO PUBLIC
- 6 8 9-6 8 9 6 9 XXX 4 > -? - 8 9 x 4 z ) - -! x - x - - X - - - - - x 00 - - - - - x z - - - x x - x - - - - - ) x - - - - - - 0 > - 000-90 - - 4 0 x 00 - -? z 8 & x - - 8? > 9 - - - - 64 49 9 x - -
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 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 informationHidden Markov Models for precipitation
Hidden Markov Models for precipitation Pierre Ailliot Université de Brest Joint work with Peter Thomson Statistics Research Associates (NZ) Page 1 Context Part of the project Climate-related risks for
More informationA. H. Hall, 33, 35 &37, Lendoi
7 X x > - z Z - ----»»x - % x x» [> Q - ) < % - - 7»- -Q 9 Q # 5 - z -> Q x > z»- ~» - x " < z Q q»» > X»? Q ~ - - % % < - < - - 7 - x -X - -- 6 97 9
More informationLOWELL WEEKLY JOURNAL
Y -» $ 5 Y 7 Y Y -Y- Q x Q» 75»»/ q } # ]»\ - - $ { Q» / X x»»- 3 q $ 9 ) Y q - 5 5 3 3 3 7 Q q - - Q _»»/Q Y - 9 - - - )- [ X 7» -» - )»? / /? Q Y»» # X Q» - -?» Q ) Q \ Q - - - 3? 7» -? #»»» 7 - / Q
More informationMCMC Sampling for Bayesian Inference using L1-type Priors
MÜNSTER MCMC Sampling for Bayesian Inference using L1-type Priors (what I do whenever the ill-posedness of EEG/MEG is just not frustrating enough!) AG Imaging Seminar Felix Lucka 26.06.2012 , MÜNSTER Sampling
More informationELEC633: Graphical Models
ELEC633: Graphical Models Tahira isa Saleem Scribe from 7 October 2008 References: Casella and George Exploring the Gibbs sampler (1992) Chib and Greenberg Understanding the Metropolis-Hastings algorithm
More informationMonte Carlo conditioning on a sufficient statistic
Seminar, UC Davis, 24 April 2008 p. 1/22 Monte Carlo conditioning on a sufficient statistic Bo Henry Lindqvist Norwegian University of Science and Technology, Trondheim Joint work with Gunnar Taraldsen,
More informationClimate Change: the Uncertainty of Certainty
Climate Change: the Uncertainty of Certainty Reinhard Furrer, UZH JSS, Geneva Oct. 30, 2009 Collaboration with: Stephan Sain - NCAR Reto Knutti - ETHZ Claudia Tebaldi - Climate Central Ryan Ford, Doug
More informationComputer Vision Group Prof. Daniel Cremers. 10a. Markov Chain Monte Carlo
Group Prof. Daniel Cremers 10a. Markov Chain Monte Carlo Markov Chain Monte Carlo In high-dimensional spaces, rejection sampling and importance sampling are very inefficient An alternative is Markov Chain
More informationProbabilistic Machine Learning
Probabilistic Machine Learning Bayesian Nets, MCMC, and more Marek Petrik 4/18/2017 Based on: P. Murphy, K. (2012). Machine Learning: A Probabilistic Perspective. Chapter 10. Conditional Independence Independent
More informationBayesian spatial hierarchical modeling for temperature extremes
Bayesian spatial hierarchical modeling for temperature extremes Indriati Bisono Dr. Andrew Robinson Dr. Aloke Phatak Mathematics and Statistics Department The University of Melbourne Maths, Informatics
More informationIntroduction to Probability and Statistics (Continued)
Introduction to Probability and Statistics (Continued) Prof. icholas Zabaras Center for Informatics and Computational Science https://cics.nd.edu/ University of otre Dame otre Dame, Indiana, USA Email:
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 informationMarkov Chain Monte Carlo, Numerical Integration
Markov Chain Monte Carlo, Numerical Integration (See Statistics) Trevor Gallen Fall 2015 1 / 1 Agenda Numerical Integration: MCMC methods Estimating Markov Chains Estimating latent variables 2 / 1 Numerical
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 informationPart III. A Decision-Theoretic Approach and Bayesian testing
Part III A Decision-Theoretic Approach and Bayesian testing 1 Chapter 10 Bayesian Inference as a Decision Problem The decision-theoretic framework starts with the following situation. We would like to
More informationName of research institute or organization: Federal Office of Meteorology and Climatology MeteoSwiss
Name of research institute or organization: Federal Office of Meteorology and Climatology MeteoSwiss Title of project: The weather in 2016 Report by: Stephan Bader, Climate Division MeteoSwiss English
More informationFin285a:Computer Simulations and Risk Assessment Section 2.3.2:Hypothesis testing, and Confidence Intervals
Fin285a:Computer Simulations and Risk Assessment Section 2.3.2:Hypothesis testing, and Confidence Intervals Overview Hypothesis testing terms Testing a die Testing issues Estimating means Confidence intervals
More informationPattern Recognition and Machine Learning. Bishop Chapter 11: Sampling Methods
Pattern Recognition and Machine Learning Chapter 11: Sampling Methods Elise Arnaud Jakob Verbeek May 22, 2008 Outline of the chapter 11.1 Basic Sampling Algorithms 11.2 Markov Chain Monte Carlo 11.3 Gibbs
More informationBayesian GLMs and Metropolis-Hastings Algorithm
Bayesian GLMs and Metropolis-Hastings Algorithm We have seen that with conjugate or semi-conjugate prior distributions the Gibbs sampler can be used to sample from the posterior distribution. In situations,
More informationQ SON,' (ESTABLISHED 1879L
( < 5(? Q 5 9 7 00 9 0 < 6 z 97 ( # ) $ x 6 < ( ) ( ( 6( ( ) ( $ z 0 z z 0 ) { ( % 69% ( ) x 7 97 z ) 7 ) ( ) 6 0 0 97 )( 0 x 7 97 5 6 ( ) 0 6 ) 5 ) 0 ) 9%5 z» 0 97 «6 6» 96? 0 96 5 0 ( ) ( ) 0 x 6 0
More informationMarkov random fields. The Markov property
Markov random fields The Markov property Discrete time: (X k X k!1,x k!2,... = (X k X k!1 A time symmetric version: (X k! X!k = (X k X k!1,x k+1 A more general version: Let A be a set of indices >k, B
More informationControl Variates for Markov Chain Monte Carlo
Control Variates for Markov Chain Monte Carlo Dellaportas, P., Kontoyiannis, I., and Tsourti, Z. Dept of Statistics, AUEB Dept of Informatics, AUEB 1st Greek Stochastics Meeting Monte Carlo: Probability
More informationA Fully Nonparametric Modeling Approach to. BNP Binary Regression
A Fully Nonparametric Modeling Approach to Binary Regression Maria Department of Applied Mathematics and Statistics University of California, Santa Cruz SBIES, April 27-28, 2012 Outline 1 2 3 Simulation
More informationThe Mixture Approach for Simulating New Families of Bivariate Distributions with Specified Correlations
The Mixture Approach for Simulating New Families of Bivariate Distributions with Specified Correlations John R. Michael, Significance, Inc. and William R. Schucany, Southern Methodist University The mixture
More informationF denotes cumulative density. denotes probability density function; (.)
BAYESIAN ANALYSIS: FOREWORDS Notation. System means the real thing and a model is an assumed mathematical form for the system.. he probability model class M contains the set of the all admissible models
More informationClimatic change in the Alps
Climatic change in the Alps Prof. Martin Beniston Martin.Beniston@unige.ch Wengen-2006 Workshop Overview Introduction Current and future climate in the Alps Potential impacts Conclusions 1 Introduction
More informationSTAT 425: Introduction to Bayesian Analysis
STAT 425: Introduction to Bayesian Analysis Marina Vannucci Rice University, USA Fall 2017 Marina Vannucci (Rice University, USA) Bayesian Analysis (Part 2) Fall 2017 1 / 19 Part 2: Markov chain Monte
More informationLecture 7 and 8: Markov Chain Monte Carlo
Lecture 7 and 8: Markov Chain Monte Carlo 4F13: Machine Learning Zoubin Ghahramani and Carl Edward Rasmussen Department of Engineering University of Cambridge http://mlg.eng.cam.ac.uk/teaching/4f13/ Ghahramani
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 informationMonte Carlo integration
Monte Carlo integration Eample of a Monte Carlo sampler in D: imagine a circle radius L/ within a square of LL. If points are randoml generated over the square, what s the probabilit to hit within circle?
More informationComputer Vision Group Prof. Daniel Cremers. 11. Sampling Methods: Markov Chain Monte Carlo
Group Prof. Daniel Cremers 11. Sampling Methods: Markov Chain Monte Carlo Markov Chain Monte Carlo In high-dimensional spaces, rejection sampling and importance sampling are very inefficient An alternative
More informationMarkov and Gibbs Random Fields
Markov and Gibbs Random Fields Bruno Galerne bruno.galerne@parisdescartes.fr MAP5, Université Paris Descartes Master MVA Cours Méthodes stochastiques pour l analyse d images Lundi 6 mars 2017 Outline The
More informationBayesian Regression Linear and Logistic Regression
When we want more than point estimates Bayesian Regression Linear and Logistic Regression Nicole Beckage Ordinary Least Squares Regression and Lasso Regression return only point estimates But what if we
More informationTerminology Suppose we have N observations {x(n)} N 1. Estimators as Random Variables. {x(n)} N 1
Estimation Theory Overview Properties Bias, Variance, and Mean Square Error Cramér-Rao lower bound Maximum likelihood Consistency Confidence intervals Properties of the mean estimator Properties of the
More informationBisection Ideas in End-Point Conditioned Markov Process Simulation
Bisection Ideas in End-Point Conditioned Markov Process Simulation Søren Asmussen and Asger Hobolth Department of Mathematical Sciences, Aarhus University Ny Munkegade, 8000 Aarhus C, Denmark {asmus,asger}@imf.au.dk
More informationII&Ij <Md Tmlaiiiiiit, aad once in Ihe y a w Teataa m i, the vmb thatalmta oot Uiaapirit world. into as abode or wotld by them- CooTBOtioa
382 4 7 q X
More informationPOSTERIOR PROPRIETY IN SOME HIERARCHICAL EXPONENTIAL FAMILY MODELS
POSTERIOR PROPRIETY IN SOME HIERARCHICAL EXPONENTIAL FAMILY MODELS EDWARD I. GEORGE and ZUOSHUN ZHANG The University of Texas at Austin and Quintiles Inc. June 2 SUMMARY For Bayesian analysis of hierarchical
More informationLikelihood inference for threshold exceedances
31 Spatial Extremes Anthony Davison École Polytechnique Fédérale de Lausanne, Switzerland Raphaël Huser King Abdullah University, Thuwal, Saudi Arabia Emeric Thibaud École Polytechnique Fédérale de Lausanne,
More informationMixture models. Mixture models MCMC approaches Label switching MCMC for variable dimension models. 5 Mixture models
5 MCMC approaches Label switching MCMC for variable dimension models 291/459 Missing variable models Complexity of a model may originate from the fact that some piece of information is missing Example
More informationIntroduction to Machine Learning CMU-10701
Introduction to Machine Learning CMU-10701 Markov Chain Monte Carlo Methods Barnabás Póczos & Aarti Singh Contents Markov Chain Monte Carlo Methods Goal & Motivation Sampling Rejection Importance Markov
More informationStat 535 C - Statistical Computing & Monte Carlo Methods. Arnaud Doucet.
Stat 535 C - Statistical Computing & Monte Carlo Methods Arnaud Doucet Email: arnaud@cs.ubc.ca 1 1.1 Outline Introduction to Markov chain Monte Carlo The Gibbs Sampler Examples Overview of the Lecture
More informationLOWELL JOURNAL. MUST APOLOGIZE. such communication with the shore as Is m i Boimhle, noewwary and proper for the comfort
- 7 7 Z 8 q ) V x - X > q - < Y Y X V - z - - - - V - V - q \ - q q < -- V - - - x - - V q > x - x q - x q - x - - - 7 -» - - - - 6 q x - > - - x - - - x- - - q q - V - x - - ( Y q Y7 - >»> - x Y - ] [
More informationStatistics & Data Sciences: First Year Prelim Exam May 2018
Statistics & Data Sciences: First Year Prelim Exam May 2018 Instructions: 1. Do not turn this page until instructed to do so. 2. Start each new question on a new sheet of paper. 3. This is a closed book
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 informationUniversal examples. Chapter The Bernoulli process
Chapter 1 Universal examples 1.1 The Bernoulli process First description: Bernoulli random variables Y i for i = 1, 2, 3,... independent with P [Y i = 1] = p and P [Y i = ] = 1 p. Second description: Binomial
More informationACCOUNTING FOR INPUT-MODEL AND INPUT-PARAMETER UNCERTAINTIES IN SIMULATION. <www.ie.ncsu.edu/jwilson> May 22, 2006
ACCOUNTING FOR INPUT-MODEL AND INPUT-PARAMETER UNCERTAINTIES IN SIMULATION Slide 1 Faker Zouaoui Sabre Holdings James R. Wilson NC State University May, 006 Slide From American
More informationBayesian Gaussian Process Regression
Bayesian Gaussian Process Regression STAT8810, Fall 2017 M.T. Pratola October 7, 2017 Today Bayesian Gaussian Process Regression Bayesian GP Regression Recall we had observations from our expensive simulator,
More informationMarkov Chain Monte Carlo
Markov Chain Monte Carlo Recall: To compute the expectation E ( h(y ) ) we use the approximation E(h(Y )) 1 n n h(y ) t=1 with Y (1),..., Y (n) h(y). Thus our aim is to sample Y (1),..., Y (n) from f(y).
More informationComputer Vision Group Prof. Daniel Cremers. 11. Sampling Methods
Prof. Daniel Cremers 11. Sampling Methods Sampling Methods Sampling Methods are widely used in Computer Science as an approximation of a deterministic algorithm to represent uncertainty without a parametric
More informationOn Markov Chain Monte Carlo
MCMC 0 On Markov Chain Monte Carlo Yevgeniy Kovchegov Oregon State University MCMC 1 Metropolis-Hastings algorithm. Goal: simulating an Ω-valued random variable distributed according to a given probability
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 informationST5212: Survival Analysis
ST51: Survival Analysis 8/9: Semester II Tutorial 1. A model for lifetimes, with a bathtub-shaped hazard rate, is the exponential power distribution with survival fumction S(x) =exp{1 exp[(λx) α ]}. (a)
More informationConditional Sampling for Max Stable Random Fields
Conditional Sampling for Max Stable Random Fields Yizao Wang Department of Statistics, the University of Michigan April 30th, 0 th GSPC at Duke University, Durham, North Carolina Joint work with Stilian
More informationHastings-within-Gibbs Algorithm: Introduction and Application on Hierarchical Model
UNIVERSITY OF TEXAS AT SAN ANTONIO Hastings-within-Gibbs Algorithm: Introduction and Application on Hierarchical Model Liang Jing April 2010 1 1 ABSTRACT In this paper, common MCMC algorithms are introduced
More informationMultimodal Nested Sampling
Multimodal Nested Sampling Farhan Feroz Astrophysics Group, Cavendish Lab, Cambridge Inverse Problems & Cosmology Most obvious example: standard CMB data analysis pipeline But many others: object detection,
More informationPrinciples of Bayesian Inference
Principles of Bayesian Inference Sudipto Banerjee University of Minnesota July 20th, 2008 1 Bayesian Principles Classical statistics: model parameters are fixed and unknown. A Bayesian thinks of parameters
More informationA spatio-temporal model for extreme precipitation simulated by a climate model
A spatio-temporal model for extreme precipitation simulated by a climate model Jonathan Jalbert Postdoctoral fellow at McGill University, Montréal Anne-Catherine Favre, Claude Bélisle and Jean-François
More informationComputer Vision Group Prof. Daniel Cremers. 14. Sampling Methods
Prof. Daniel Cremers 14. Sampling Methods Sampling Methods Sampling Methods are widely used in Computer Science as an approximation of a deterministic algorithm to represent uncertainty without a parametric
More informationIce Mass & Sea Level Change Unit 2: Temperature - A Global Trendsetter
Ice Mass & Sea Level Change Unit 2: Temperature - A Global Trendsetter Becca Walker and Leigh Stearns Global air temperatures are calculated using satellite observations, ground-based weather stations,
More informationBayesian Phylogenetics:
Bayesian Phylogenetics: an introduction Marc A. Suchard msuchard@ucla.edu UCLA Who is this man? How sure are you? The one true tree? Methods we ve learned so far try to find a single tree that best describes
More informationDivision of Informatics, University of Edinburgh
T E H U N I V E R S I T Y O H F R G Division of Informatics, University of Edinburgh E D I N B U Institute for Adaptive and Neural Computation An Analysis of Contrastive Divergence Learning in Gaussian
More informationBayesian Modeling of Conditional Distributions
Bayesian Modeling of Conditional Distributions John Geweke University of Iowa Indiana University Department of Economics February 27, 2007 Outline Motivation Model description Methods of inference Earnings
More informationMetropolis Hastings. Rebecca C. Steorts Bayesian Methods and Modern Statistics: STA 360/601. Module 9
Metropolis Hastings Rebecca C. Steorts Bayesian Methods and Modern Statistics: STA 360/601 Module 9 1 The Metropolis-Hastings algorithm is a general term for a family of Markov chain simulation methods
More informationP A L A C E P IE R, S T. L E O N A R D S. R a n n o w, q u a r r y. W WALTER CR O TC H, Esq., Local Chairman. E. CO O PER EVANS, Esq.,.
? ( # [ ( 8? [ > 3 Q [ ««> » 9 Q { «33 Q> 8 \ \ 3 3 3> Q»«9 Q ««« 3 8 3 8 X \ [ 3 ( ( Z ( Z 3( 9 9 > < < > >? 8 98 ««3 ( 98 < # # Q 3 98? 98 > > 3 8 9 9 ««««> 3 «>
More informationMonte Carlo-based statistical methods (MASM11/FMS091)
Monte Carlo-based statistical methods (MASM11/FMS091) Jimmy Olsson Centre for Mathematical Sciences Lund University, Sweden Lecture 12 MCMC for Bayesian computation II March 1, 2013 J. Olsson Monte Carlo-based
More informationBayesian Nonparametric Regression for Diabetes Deaths
Bayesian Nonparametric Regression for Diabetes Deaths Brian M. Hartman PhD Student, 2010 Texas A&M University College Station, TX, USA David B. Dahl Assistant Professor Texas A&M University College Station,
More informationVariational inference
Simon Leglaive Télécom ParisTech, CNRS LTCI, Université Paris Saclay November 18, 2016, Télécom ParisTech, Paris, France. Outline Introduction Probabilistic model Problem Log-likelihood decomposition EM
More informationBayesian Inference. Chapter 9. Linear models and regression
Bayesian Inference Chapter 9. Linear models and regression M. Concepcion Ausin Universidad Carlos III de Madrid Master in Business Administration and Quantitative Methods Master in Mathematical Engineering
More informationModeling with Itô Stochastic Differential Equations
Modeling with Itô Stochastic Differential Equations 2.4-2.6 E. Allen presentation by T. Perälä 27.0.2009 Postgraduate seminar on applied mathematics 2009 Outline Hilbert Space of Stochastic Processes (
More informationHybrid Dirichlet processes for functional data
Hybrid Dirichlet processes for functional data Sonia Petrone Università Bocconi, Milano Joint work with Michele Guindani - U.T. MD Anderson Cancer Center, Houston and Alan Gelfand - Duke University, USA
More informationName of research institute or organization: Federal Office of Meteorology and Climatology MeteoSwiss
Name of research institute or organization: Federal Office of Meteorology and Climatology MeteoSwiss Title of project: The weather in 2017 Report by: Stephan Bader, Climate Division MeteoSwiss English
More informationIn this chapter we study several functions that are useful in calculus and other areas of mathematics.
Calculus 5 7 Special functions In this chapter we study several functions that are useful in calculus and other areas of mathematics. 7. Hyperbolic trigonometric functions The functions we study in this
More informationBranching Brownian motion seen from the tip
Branching Brownian motion seen from the tip J. Berestycki 1 1 Laboratoire de Probabilité et Modèles Aléatoires, UPMC, Paris 09/02/2011 Joint work with Elie Aidekon, Eric Brunet and Zhan Shi J Berestycki
More informationMAS2317/3317: Introduction to Bayesian Statistics
MAS2317/3317: Introduction to Bayesian Statistics Case Study 2: Bayesian Modelling of Extreme Rainfall Data Semester 2, 2014 2015 Motivation Over the last 30 years or so, interest in the use of statistical
More informationMarkov Chain Monte Carlo Lecture 4
The local-trap problem refers to that in simulations of a complex system whose energy landscape is rugged, the sampler gets trapped in a local energy minimum indefinitely, rendering the simulation ineffective.
More information