Conditional simulation of max-stable processes
|
|
- John Robinson
- 5 years ago
- Views:
Transcription
1 Biometrika (1),,, pp 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, Université de Poitiers, Téléport, BP 3179, F-8696 Futuroscope-Chasseneuil cede, France 8 clement.dombry@math.univ-poitiers.fr, frederic.eyi.minko@math.univ-poitiers.fr AND M. RIBATET Department of Mathematics, Université Montpellier, 4 place Eugène Bataillon, 3495 cede Montpellier, France mathieu.ribatet@math.univ-montp.fr SUMMARY Since many environmental processes are spatial in etent, a single etreme event may affect several locations, and their spatial dependence has to be appropriately taken into account. This paper proposes a framework for conditional simulation of ma-stable processes and gives closed forms for the regular conditional distributions of Brown Resnick and Schlather processes. We test the method on simulated data and give an application to etreme rainfall around Zurich and etreme temperatures. The proposed framework provides accurate conditional simulations and can handle real-sized problems. Some key words: Conditional simulation; Markov chain Monte Carlo; Ma-stable process; Precipitation; Regular conditional distribution; Temperature
2 C. DOMBRY, F. ÉYI-MINKO AND M. RIBATET 1. INTRODUCTION Ma-stable processes arise naturally when studying etremes of stochastic processes and therefore play a major role in the statistical modelling of spatial etremes (Buishand et al, 8; Padoan et al., 1; Davison et al., 1). Although a different spectral characterization of ma-stable processes eists (de Haan, 1984), for our purposes the most useful representation is (Schlather, ) Z() = ma i 1 ζ iy i (), R d, (1) where {ζ i } i 1 are the points of a Poisson process on (, ) with intensity dλ(ζ) = ζ dζ and the Y i are independent replicates of a non-negative stochastic process Y such that E{Y ()} = 1 for all R d. It is well known that Z is a ma-stable process on R d with unit Fréchet margins (Schlather, ; de Haan & Fereira, 6, page 37). Although (1) takes the pointwise maimum over an infinite number of points {ζ i } and processes Y i, it is possible to get approimate realizations from Z (Schlather, ; Oesting et al., 1). Based on (1) several parametric ma-stable models have been proposed (Schlather, ; Brown & Resnick, 1977; Kabluchko et al., 9; Davison et al., 1) that share the same finitedimensional distribution functions [ { pr{z( 1 ) z 1,..., Z( k ) z k } = ep E where k N, z 1,..., z k > and 1,..., k R d. ma j=1,...,k }] Y ( j ), z j Apart from the Smith model (Genton et al., 11), only the bivariate cumulative distribution functions are eplicitly known. To bypass this hurdle to inference based on ma-stable processes, de Haan & Pereira (6) propose a semi-parametric estimator and Padoan et al. (1) suggest the use of the maimum pairwise likelihood estimator
3 Conditional simulation of ma-stable processes 3 Paralleling the use of the variogram in classical geostatistics, the etremal coefficient function (Schlather & Tawn, 3; Cooley et al., 6) θ( 1 ) = z log pr{z( 1 ) z, Z( ) z} is widely used to summarize the spatial dependence of etremes for stationary processes. It takes values in the interval [1, ]; the lower bound indicates complete dependence, and the upper bound independence. The last decade has seen many advances to develop geostatistics for etremes, and software is available to practitioners (Wang, 1; Schlather, 1; Ribatet, 11). However an important tool currently missing is conditional simulation of ma-stable processes. In classical geostatistic based on Gaussian models, conditional simulation is well established (Chilès & Delfiner, 1999) and provides a framework to assess the distribution of a Gaussian random field given values observed at fied locations. For eample, conditional simulations of Gaussian processes have been used to model land topography (Mandelbrot, 198). Conditional simulation of ma-stable processes is a long-standing problem (Davis & Resnick, 1989, 1993). Wang & Stoev (11) provide a first solution, but their framework is limited to processes having a discrete spectral measure and thus may be too restrictive to appropriately model spatial dependence in comple situations. Based on the recent developments on the regular conditional distribution of ma-infinitely divisible processes, the aim of this paper is to provide a methodology for conditional simulation of ma-stable processes with continuous spectral measures. More formally for a study region X R d, our goal is to derive an algorithm to sample from the regular conditional distribution of Z {Z( 1 ) = z 1,..., Z( k ) = z k } for some z 1,..., z k > and k conditioning locations 1,..., k X
4 C. DOMBRY, F. ÉYI-MINKO AND M. RIBATET. CONDITIONAL SIMULATION OF MAX-STABLE PROCESSES 1. General framework This section reviews some key results of an unpublished paper available from the first author with a particular emphasis on ma-stable processes. Our goal is to give a more practical interpretation of their results from a simulation perspective. To this aim, we recall two key results and propose a procedure to generate conditional realizations of ma-stable processes. Let R X be the space of real-valued functions on X R d and let Φ = {ϕ i } i 1 be a Poisson point process on R X where ϕ i () = ζ i Y i () (i = 1,,...) with ζ i and Y i as in (1). We write f() = {f( 1 ),..., f( k )} for all random functions f : X R and = ( 1,..., k ) X k. It is not difficult to show that for all Borel sets A R k, the Poisson process {ϕ i ()} i 1 defined on R k has intensity measure Λ (A) = pr {ζy () A} ζ dζ. The point process Φ is called regular if the intensity measure Λ has an intensity function λ, i.e., Λ (dz) = λ (z)dz, for all X k. The first key result is that provided the Poisson process Φ is regular, the intensity function λ and the conditional intensity function λ s,z (u) = λ (s,)(u, z), (s, ) X m+k, u R m, z (, + ) k, () λ (z) drive how the conditioning terms {Z( j ) = z j } (j = 1,..., k) are met, see Steps 1 and in Theorem 1. The second key result is that, conditionally on Z() = z, the Poisson process Φ can be decomposed into two independent point processes, say Φ = Φ Φ +, where Φ = {ϕ Φ: ϕ( i ) < z i, i = 1,..., k}, Φ + = k {ϕ Φ: ϕ( i ) = z i }. i=
5 Conditional simulation of ma-stable processes 5 Before introducing a procedure to get conditional realizations of ma-stable processes, we introduce notation and make connections with the pioneering work of Wang & Stoev (11). A function ϕ Φ + such that ϕ( i ) = z i for some i {1,..., k} is called an etremal function associated to i and denoted by ϕ + i. It is easy to show that there eists almost surely a unique etremal function associated to i. Although Φ + = {ϕ + 1,..., ϕ + k } almost surely, it might happen that a single etremal function contributes to the random vector Z() at several locations i, e.g., ϕ + 1 = ϕ +. To take such repetitions into account, we define a random partition θ = (θ 1,..., θ l ) of the set { 1,..., k } into l = θ blocks, and etremal functions (ϕ + 1,..., ϕ+ l ) such that Φ+ = {ϕ + 1,..., ϕ+ l } and ϕ+ j ( i) = z i if i θ j and ϕ + j ( i) < z i if i / θ j (i = 1,..., k; j = 1,..., l). Wang & Stoev (11) call the partition θ the hitting scenario. The set of all possible partitions of { 1,..., k }, denoted P k, identifies all possible hitting scenarios. From a simulation perspective, the fact that Φ and Φ + are independent given Z() = z is especially convenient and suggests a three step procedure to sample from the conditional distribution of Z given Z() = z. THEOREM 1. Suppose that the point process Φ is regular and let (, s) X k+m. For τ = (τ 1,..., τ l ) P k and j = 1,..., l, define I j = {i: i τ j }, τj = ( i ) i Ij, z τj = (z i ) i Ij, τ c j = ( i ) i/ Ij and z τ c j = (z i ) i/ Ij. Consider the three step procedure: Step 1. Draw a random partition θ P k with distribution π (z, τ) = pr {θ = τ Z() = z} = where the normalization constant is C(, z) = τ τ P k j=1 1 C(, z) τ j=1 λ τj (z τj ) λ τ c τj,z τj (u j )du j, {u j <z τ c } j j λ τj (z θj ) λ τ c (u τj,z τj j )du j. {u j <z τ c } j j
6 C. DOMBRY, F. ÉYI-MINKO AND M. RIBATET Step. Given τ = (τ 1,..., τ l ), draw l independent random vectors ϕ + 1 (s),..., ϕ+ l (s) from the distribution } pr {ϕ + j (s) dv Z() = z, θ = τ = 1 { C j where 1 { } is the indicator function and C j = 1 {u<zτ c j }λ (s,τ c j ) τj,z τj (v, u)dudv, and define the random vector Z + (s) = ma j=1,...,l ϕ + j (s). } 1 {u<zτ c j }λ (s,τ c ) τj,z τj (v, u)du dv, j Step 3. Independently draw a Poisson point process {ζ i } i 1 on (, ) with intensity ζ dζ and {Y i } i 1 independent copies of Y, and define the random vector Z (s) = ma i 1 ζ iy i (s)1 {ζi Y i ()<z}. Then the random vector Z(s) = ma {Z + (s), Z (s)} follows the conditional distribution of Z(s) given Z() = z. The corresponding conditional cumulative distribution function is τ pr{z(s) a, Z() z} pr {Z(s) a Z() = z} = π (z, τ) F τ,j (a), (3) pr{z() z} τ P k j=1 where } F τ,j (a) = pr {ϕ + j (s) a Z() = z, θ = τ = {u<z τ c,v<a} λ (s, τ c ) τj,z τj (v, u)dudv j j {u<z τ c } λ. τ c τj,z τj (u)du j j The first term in the right hand side of (3) corresponds to Steps 1 and while the ratio is a consequence of Step 3. It is clear from (3) that the conditional random field Z {Z() = z} is not ma-stable.. Distribution of the etremal functions In this section we derive closed forms for the intensity function λ (z) and the conditional intensity function λ s,z (u) for two widely used ma-stable processes; the Brown Resnick (Brown
7 Conditional simulation of ma-stable processes 7 & Resnick, 1977; Kabluchko et al., 9) and the Schlather () processes. Details of the derivations are given in the Appendi. The Brown Resnick process corresponds to the case where Y () = ep{w () γ()} ( R d ) in (1) with W a centered Gaussian process with stationary increments, semivariogram γ and such that W (o) = almost surely where o denotes the origin of R d. For X k and provided the covariance matri Σ of the random vector W () is positive definite, the intensity function is ( λ (z) = C ep 1 ) k log zt Q log z + L log z zi 1, z (, ) k, with 1 k = (1) i=1,...,k, γ = {γ( i )} i=1,...,k, Q = Σ 1 Σ 1 1 k 1 T k Σ 1 1 T, L = k Σ 1 1 k C = (π) (1 k)/ Σ 1/ (1 T k Σ 1 1 k ) 1/ ep ( 1 T k Σ 1 γ 1 1 T k Σ 1 i=1 1 k γ 1 k { 1 ) T Σ 1, } (1 T k Σ 1 γ 1) 1 T 1 k Σ 1 1 k γt Σ 1 γ, and for all (s, ) X m+k, (u, z) (, ) m+k and provided the covariance matri Σ (s,) is positive definite, the conditional intensity function corresponds to a multivariate log-normal probability density function { λ s,z (u) = (π) m/ Σ s 1/ ep 1 } (log u µ s,z) T Σ 1 s (log u µ m s,z) u 1 i, where µ s,z R m and Σ s are the mean and covariance matri of the underlying normal distribution and are given by with } Σ 1 s = J m,k T Q (s,)j m,k, µ s,z = {L (s,) J m,k log z T T J m,k Q (s,) J m,k Σ s, J m,k = I m k,m, Jm,k = m,k I k, where I k denotes the k k identity matri and m,k the m n null matri. i=
8 C. DOMBRY, F. ÉYI-MINKO AND M. RIBATET The Schlather process considers the case Y () = (π) 1/ ma{, ε()} ( R d ) in (1) with ε a standard Gaussian process with correlation function ρ. The associated point process Φ is not regular and it is more convenient to consider the equivalent representation where Y () = (π) 1/ ε() ( R d ). For X k and provided the covariance matri Σ of the random vector ε() is positive definite, the intensity function is ( ) k + 1 λ (z) = π (k 1)/ Σ 1/ a (z) (k+1)/ Γ, z R k, where a (z) = z T Σ 1 z. For (s, ) X m+k, (u, z) R m+k and provided that the covariance matri Σ (s,) is positive definite, the conditional intensity function λ s,z (u) corresponds to the density of a multivariate Student distribution with k + 1 degrees of freedom, location parameter µ = Σ s: Σ 1 z, and scale matri Σ = a (z) ( Σs Σ s: Σ 1 ) Σ :s, Σ(s,) = k + 1 Σ s Σ :s 3. MARKOV CHAIN MONTE CARLO SAMPLER Section introduced a procedure to get realizations from the regular conditional distribution of ma-stable processes. This sampling scheme amounts to sampling from a discrete distribution whose state space corresponds to all possible partitions of the set of conditioning points, see Step 1 of Theorem 1. Hence, even for a moderate number k of conditioning locations, the state space P k becomes very large and the distribution π (z, ) cannot be computed eactly. A Gibbs sampler is especially convenient for Monte Carlo sampling from π (z, ). For τ P k, let τ j be the restriction of τ to the set { 1,..., k } \ { j }. Our goal is to simulate from the conditional distribution Σ s: Σ. pr(θ θ j = τ j ), (4)
9 Conditional simulation of ma-stable processes 9 where θ P k is a random partition which follows the target distribution π (z, ). Since the number of possible updates is always less than k, a combinatorial eplosion is avoided. Indeed for τ P k of size l, the number of partitions τ P k such that τ j = τ j for some j {1,..., k} is l if { j } is a partitioning set of τ, b + = l + 1 if { j } is not a partitioning set of τ, since the point j may be reallocated to any partitioning set of τ j or to a new one. To illustrate consider the set { 1,, 3 } and let τ = ({ 1, }, { 3 }). Then the possible partitions τ such that τ = τ are ({ 1, }, { 3 }), ({ 1 }, { }, { 3 }), ({ 1 }, {, 3 }), while there eist only two partitions such that τ 3 = τ 3, i.e., ({ 1, }, { 3 }), ({ 1,, 3 }). The distribution (4) has nice properties. Since for all τ P k such that τ j = τ j we have where pr (θ = τ θ j = τ j ) = π (z, τ ) τ P k π (z, τ)1 { τ j =τ j } w τ,j = λ τj (z τj ) λ τ c τj,z τj (u)du. {u<z τ c } j j τ j=1 w τ,j τ j=1 w, (5) τ,j Since many factors cancel out on the right hand side of (5), the Gibbs sampler is especially convenient. The most computationally demanding part of (5) is the evaluation of the integral {u<z τ c j } λ τ c j τj,z τj (u)du. For the Brown Resnick and Schlather processes, we follow the lines of Genz (199) and compute these probabilities using a separation of variables method which provides a transformation of the original integration problem to the unit hyper-cube. A quasi Monte Carlo scheme and antithetic variables are used to improve efficiency
10 C. DOMBRY, F. ÉYI-MINKO AND M. RIBATET Since it is not obvious how to implement a Gibbs sampler whose target distribution has support P k, the remainder of this section gives practical details. For any fied locations 1,..., k X, we first describe how each partition of { 1,..., k } is stored. To illustrate this consider the set { 1,, 3 } and the partition ({ 1, }; { 3 }). This partition is defined as (1, 1, ), indicating that 1 and belong to the partitioning set labeled 1 and 3 belongs to the partitioning set. There eist several equivalent notations for this partition: for eample one can use (,, 1) or (1, 1, 3). Since there is a one-one mapping between P k and the set } Pk {(a = 1,..., a k ): i {,..., k}, 1 = a 1 a i ma a j + 1, a i Z, 1 j<i we shall restrict our attention to the partitions that live in Pk, and going back to our eample we see that (1, 1, ) is valid but (,, 1) and (1, 1, 3) are not. For τ P k of size l, let r 1 = k i=1 1 {τ i =a j } and r = k i=1 1 {τ i =b}, i.e., the number of conditioning locations that belong to the partitioning sets a j and b where b {1,..., b + } with l (r 1 = 1), b + = l + 1 (r 1 1). Then the conditional probability distribution (5) satisfies pr(τ j = b τ i = a i, i j) 1 (b = a j ), (6a) w τ,b/(w τ,b w τ,aj ) w τ,bw τ,a j /(w τ,b w τ,aj ) w τ,bw τ,a j /w τ,aj (r 1 = 1, r, b a j ), (6b) (r 1 1, r, b a j ), (6c) (r 1 1, r =, b a j ), (6d) where τ = (a 1,..., a j 1, b, a j+1,..., a k ). Although τ may not belong to Pk, it corresponds to a unique partition of P k and we can use the bijection P k P k to recode τ into an element of P k. In (6a) (6d) the event {r 1 = 1, r =, b a j } is missing since {r 1 = 1, r = } implies
11 Conditional simulation of ma-stable processes 11 Table 1. Sample path properties of the ma-stable models. For the Brown Resnick model, the variogram parameters are set to ensure that the etremal coefficient function satisfies θ(115) = 1 7 while for the Schlather model the correlation function parameters are set to ensure that θ(1) = 1 5. Brown Resnick: γ(h) = (h/λ) κ Schlather: ρ(h) = ep{ (h/λ) κ } γ 1: Very wiggly γ : Wiggly γ 3: Smooth ρ 1: Very wiggly ρ : Wiggly ρ 3: Smooth λ κ that τ = τ, where the equality has to be understood in terms of elements of P k, and this case has been already taken into account with (6a). Once these conditional weights have been computed, the Gibbs sampler proceeds by updating each element of τ successively. We use a random scan implementation (Liu et al., 1995). One iteration selects an element of τ at random according to a given distribution, say p = (p 1,..., p k ), and then updates this element. Throughout this paper we will use the uniform random scan Gibbs sampler for which the selection distribution is a discrete uniform distribution, i.e., p = (k 1,..., k 1 ) SIMULATION STUDY In this section we check if our algorithm is able to produce realistic conditional simulations of Brown Resnick and Schlather processes. For each model, we consider three different sample path properties, as summarized in Table 1. These configurations were chosen such that the spatial dependence structures are similar to our applications in Section 5. In order to check if our sampling procedure is accurate and given a single conditional event {Z() = z} for each configuration, we generated 1 conditional realizations with standard
12 59 1 C. DOMBRY, F. ÉYI-MINKO AND M. RIBATET Z() 3 1 Z() 3 1 Z() Z() 1 Z() 1 Z() Fig. 1. Pointwise sample quantiles estimated from 1 conditional simulations of ma-stable processes with standard Gumbel margins with k = 5, 1, 15 conditioning locations. The top row shows results for the Brown Resnick models with semivariograms γ 3, γ, γ 1, from left to right. The bottom row shows results for the Schlather models with correlation functions ρ 3, ρ, ρ 1, from left to right. The solid black lines shows the pointwise 5, 5, 975 sample quantiles and the dashed grey lines that of a standard Gumbel distribution. The squares show the conditional points {( i, z i)} i=1,...,k. The solid grey lines show the simulated paths used to get the conditioning events. Gumbel margins. Figure 1 shows the pointwise sample quantiles obtained from these 1 simulated paths and compares them to unit Gumbel quantiles. As epected the conditional sample paths inherit the regularity driven by the shape parameter κ and there is less variability in regions close to conditioning locations. Since the Brown Resnick processes considered are ergodic
13 Conditional simulation of ma-stable processes 13 Table. Timings for conditional simulations of ma-stable processes on a 5 5 grid defined on the square [, 1 1/ ] for a varying number k of conditioning locations uniformly distributed over the region. The times, in seconds, are mean values over 1 conditional simulations; standard deviations are reported in brackets. Brown Resnick: γ(h) = (h/5) 5 Schlather: ρ(h) = ep { (h/8) 5} Step 1 Step Step 3 Overall Step 1 Step Step 3 Overall k = 5 1 ( 1) 49 (11) 1 4 ( 1) 5 (11) 1 4 ( ) 1 9 ( 7) 9 ( 3) 4 ( 8) k = 1 8 () 76 (18) 1 4 ( 1) 85 (19) 1 (4) 4 ( 8) 1 ( 3) 15 (4) k = 5 95 (38) 117 (3) 1 4 ( 1) 14 (61) 86 (4) 4 (1) 1 ( 3) 9 (43) k = (313) 348 (391) 1 5 ( 1) 931 (535) 367 () 6 (113) 1 ( 3) 43 (6) Conditional simulations with k = 5 do not use a Gibbs sampler (Kabluchko and Schlather, 1), for regions far away from any conditioning location the sample quantiles converges to that of a standard Gumbel distribution, indicating that the conditional event has no influence. This is not the case for the non-ergodic Schlather process. Most of the time the sample paths used to get the conditional events belong to the 95% pointwise confidence intervals, corroborating that our sampling procedure seems to be accurate. Table gives timings for conditional simulations of ma-stable processes on a 5 5 grid with a varying number of conditioning locations. Due to the combinatorial compleity of the partition set P k, the timings increase rapidly with respect to the number of conditioning points k. It is however reassuring that the algorithm is tractable when k 5; this covers many practical situations
14 C. DOMBRY, F. ÉYI-MINKO AND M. RIBATET (m) (km) Maimum rainfall in Zurich (mm) θ(h) Years Fig.. Left: Map of Switzerland showing the stations of the 4 rainfall gauges used for the analysis, with an insert showing the altitude. The station marked with a triangle corresponds to Zurich. Middle: Summer maimum rainfall values for at Zurich. Right: Comparison between the pairwise etremal coefficient estimates for the 51 original weather stations and the etremal coefficient function derived from a fitted Brown Resnick processes having semi variogram γ(h) = (h/λ) κ. The grey points are pairwise estimates; the black ones are binned estimates and the curve is the fitted etremal coefficient function h (km) APPLICATION 5 1. Etreme precipitation around Zurich The data considered here were previously analyzed by Davison et al. (1), who showed that Brown Resnick processes were among the best models for these data. The data are summer maimum rainfall for the years at 51 weather stations in the Plateau region of Switzerland, provided by the national meteorological service, MeteoSuisse. To ensure strong dependence between the conditioning locations, we consider as conditional locations the 4 weather stations that are at most 3km distant from Zurich and set as the conditional values the rainfall amounts recorded in the year, the year of the largest precipitation event recorded there between
15 Conditional simulation of ma-stable processes 15 Table 3. Distribution of the partition size for the rainfall and temperature data estimated from simulated Markov chains of length 15 and 1 respectively 675 Rainfall Temperature Partition size Frequencies (%) < , see Figure. The largest and smallest distances between the conditioning locations are around 55km and just over 4km respectively. A Brown Resnick process having semivariogram γ(h) = (h/λ) κ was fitted using the maimum pairwise likelihood estimator introduced by Padoan et al. (1) to simultaneously estimate the marginal parameters and the spatial dependence parameters λ and κ. As in Davison et al. (1), the marginal parameters were described by η() = β,η + β 1,η lon() + β,η lat(), σ() = β,σ + β 1,σ lon() + β,σ lat(), ξ() = β,ξ, where η(), σ(), ξ() are the location, scale and shape parameters of the generalized etreme value distribution and lon(), lat() denote the longitude and latitude of the stations at which the data are observed. The maimum pairwise likelihood estimates and standard errors for λ and κ are 38 (14) and 69 ( 7) and give a practical etremal range, i.e., the distance h + such that θ(h + ) = 1 7, of around 115km; see the right panel of Figure. Table 3 shows the distribution of the partition size estimated from a Markov chain of length 15. Around 65% of the time the summer maima observed at the 4 conditioning locations belonged to a single etremal function, i.e., only one storm event, and around 3% of the time to two different storms. Since the simulated Markov chain keeps a trace of all the simulated partitions, we looked at the partitions of size two and saw that around 65% of the time, at least one of the four up north conditioning locations was impacted by one etremal function while the remaining locations were always influenced by another one
16 C. DOMBRY, F. ÉYI-MINKO AND M. RIBATET Fig. 3. From left to right, maps on a 5 5 grid of the pointwise 5, 5 and 975 sample quantiles for rainfall (mm) obtained from 1 conditional simulations of Brown Resnick processes having semi variogram γ(h) = (h/38) 69. The rightmost panel plots the ratio of the pointwise confidence intervals with and without taking estimation uncertainties into account. The squares show the conditional locations Figure 3 plots the pointwise 5, 5 and 975 sample quantiles obtained from 1 conditional simulations of our fitted Brown Resnick process. The conditional median provides a point estimate for the rainfall at an ungauged location and the 5 and 975 conditional quantiles a 95% pointwise confidence interval. As indicated by Figure 1, the shape parameter κ has a major impact on the regularity of paths and on the width of the confidence interval. The value ˆκ 69 corresponds to very wiggly sample paths and wider confidence intervals. To assess the impact of parameter uncertainties on conditional simulations, the ratio of the width of the confidence intervals with and without parameter uncertainty is shown in the right panel of Figure 3. The uncertainties were taken into account by sampling from the asymptotic distribution of the maimum composite likelihood estimator, and drawing one conditional simulation for each realization. These ratios show no clear spatial pattern and the width of the confidence interval is increased by at most 1%
17 Basel (316) Zurich (556) Santis (49) Oeschberg (483) Neuchatel (485) Bad Ragaz (496) Engelberg (135) Bern (565) Chateau d Oe (985) Montreu (45) Montana (158) Gd St Bernard (47) (m) Conditional simulation of ma-stable processes 17 Locarno Monti (366) Lugano (73) Davos (159) Arosa (184) (km) Fig. 4. Left: Topographical map of Switzerland showing the sites and altitudes in metres above sea level of 16 weather stations for which annual maima temperature data are available. Right: Map of temperature anomalies ( C), i.e., the difference between the pointwise medians obtained from 1 conditional simulations and unconditional medians estimated from the fitted Schlather process. 5. Etreme temperatures in Switzerland The data considered here are annual maimum temperatures recorded at 16 sites in Switzerland during the period , see Figure 4. Following Davison and Gholamrezaee (1), we fit a Schlather process with an isotropic powered eponential correlation function and trend surfaces η() = β,η + β 1,η alt(), σ() = β,σ, ξ() = β,ξ + β 1,ξ alt(), where alt() denotes the altitude above mean sea level in kilometres and {η(), σ(), ξ()} are the location, scale and shape parameters of the generalized etreme value distribution at location. The spatial dependence parameter estimates and their standard errors are ˆλ = 6 (149) and ˆκ =.5 (.1) and yield a fitted etremal coefficient function similar to our test case ρ 3 in Section 4. In year 3, western Europe endured a heatwave believed to be most severe recorded since 154 (Beniston, 4). Switzerland was severely impacted by this event: the nationwide record temperature of 41 5 C was recorded that year in Grono, Graubünden, near Lugano. For our analysis we condition on the maimum temperatures observed in 3. Based on the fitted Schlather ( C)
18 C. DOMBRY, F. ÉYI-MINKO AND M. RIBATET model, we simulate a Markov chain of effective length 1 with a burn-in period of length 5 and a thinning lag of 1 iterations. The distribution of the partition size estimated from these Markov chains is shown in Table 3. Around 9% of the time the conditional realizations were a consequence of at most three etremal functions. Since our original observations were not summer maima but maimum daily values, a close inspection of the times series in year 3 reveals that the hottest temperatures occurred between the 11th and 13th of August and, to some etent, corroborates the distribution of Table 3. The right panel of Figure 4 shows the spatial distribution of temperature anomalies, i.e., the difference between the pointwise conditional medians obtained from 1 conditional simulations and the pointwise unconditional medians estimated from the fitted Schlather model. As epected, the largest deviations occur in the plateau region of Switzerland while appreciably smaller values are found in the Alps. The differences range between 5 C and 4 75 C and are consistent with the values reported by climatologists for mean values (Beniston, 4) ACKNOWLEDGEMENTS M. Ribatet was partly funded by the MIRACCLE-GICC and McSim ANR projects. The authors thank MeteoSuisse and Dr. S. A. Padoan for providing the precipitation data, Prof. A. C. Davison and Dr. M. M. Gholamrezaee for providing the temperature data and Dr. M. Oesting for pointing out an error in the computations for the Brown Resnick model APPENDIX Brown Resnick model For all X k and Borel set A R k Λ (A) = pr [ζ ep{w () γ()} A] ζ dζ = R k 1 {ζ ep{y γ()} A} f (y)dyζ dζ,
19 Conditional simulation of ma-stable processes 19 where f denotes the density of the random vector W (), i.e., a centered Gaussian random vector with covariance matri Σ and variance γ(). The change of variables z = ζ ep{y γ()} and r = log ζ yields with Since with Λ (A) = A λ (z) = f {log z r + γ()} k i=1 k i=1 z 1 i dze r dr = λ (z)dz A z 1 i f {log z r + γ()}e r dr. f {log z r + γ()}e r = (π) k/ Σ 1/ ep { 1 } P (r), P (r) = r 1 T k Σ 1 1 k r [ 1 T k Σ 1 {log z + γ()} 1 ] + {log z + γ()} T Σ 1 {log z + γ()}, standard computations for Gaussian integrals give The conditional intensity function is λ s,z (u) = C (s,) C ( λ (z) = C ep 1 ) log k zt Q log z + L log z z 1 i. { ep 1 log (u, z)t Q (s,) log(u, z) + L (s,) log(u, z) + 1 } log m zt Q log z L log z u 1 i, and since log(u, z) = J m,k log u + J m,k log z, it is not difficult to show that λ s,z (u) = C (s,) C i=1 { ep 1 } (log u µ s,z) T Σ 1 s (log u µ m s,z) u 1 i. Finally, the relation C (s,) /C = (π) m/ Σ s 1/ is a simple consequence of the normalization λs,z (u)du = 1. For all X k and Borel sets A R k Λ (A) = Schlather model { } pr (π) 1/ ζε() A ζ dζ = i=1 R k 1 {(π) 1/ ζy A} f (y)dyζ dζ, i=
20 C. DOMBRY, F. ÉYI-MINKO AND M. RIBATET where f denotes the density of the random vector ε(), i.e., a centered Gaussian random vector with covariance matri Σ. The change of variable z = (π) 1/ ζy gives { } Λ (A) = (π) k/ z f ζ (k+) dzdζ A (π) 1/ ζ ( = (π) k Σ 1/ ep 1 ) A 4πζ zt Σ 1 z ζ (k+) dzdζ { ( ) 1/ = (π) k Σ 1/ π 4π A z T Σ 1 z E[Xk 1 ]dz, X Weibull z T Σ 1, } z ( ) (k 1)/ ( ) = (π) k Σ 1/ π 4π k + 1 A z T Σ 1 z z T Σ 1 Γ dz z = λ (z)dz, A where λ (z) = π (k 1)/ Σ 1/ a (z) (k+1)/ Γ {(k + 1)/} and a (z) = z T Σ 1 z. For all u R m the conditional intensity function is λ s,z (u) = π m/ Σ (s,) 1/ { } (m+k+1)/ a(s,) (u, z) a Σ 1/ (z) m/ Γ ( ) m+k+1 a (z) Γ ( ). k+1 We start by focusing on the ratio a (s,) (u, z)/a (z). Since Σ s Σ :s Σ s: Σ 1 = Σ 1 straightforward algebra shows that a (s,) (u, z) a (z) (Σ s Σ s: Σ 1 Σ :s ) 1 (Σ s Σ s: Σ 1 Σ :s ) 1 Σ s: Σ 1, Σ :s (Σ s Σ s: Σ 1 Σ :s ) 1 Σ 1 Σ :s (Σ s Σ s: Σ 1 Σ :s ) 1 Σ s: Σ 1 + Σ 1 = 1 + (u µ)t Σ 1 (u µ), µ = Σ s: Σ 1 a (z) ( z, Σ = Σs Σ s: Σ 1 ) Σ :s. k + 1 k + 1 We now simplify the ratio Σ (s,) / Σ. Using the fact that Σ (s,) = Σ s Σ :s Σ s: Σ combined with some more algebra yields Σ (s,) Σ = I m Σ s: Σ s Σ s: Σ 1 Σ :s m,k, k,m Σ Σ 1 Σ :s I k = Σ s Σ s: Σ 1 Σ :s = { } m k + 1 Σ. a (z) Using the two previous results it is easily found that { } (m+k+1)/ λ s,z (u) = π m/ (k + 1) m/ Σ 1/ 1 + (u ( µ)t Σ 1 (u µ) Γ m+k+1 ) k + 1 Γ ( ), k
21 961 Conditional simulation of ma-stable processes 1 which corresponds to the density of a multivariate Student distribution with the stated parameters REFERENCES BENISTON, M. (4). The 3 heat wave in Europe: A shape of things to come? An analysis based on Swiss climatological data and model simulations. Geophys. Res. Lett. 31, L. BROWN, B. M. & RESNICK, S. (1977). Etreme values of independent stochastic processes. J. Appl. Prob. 14, BUISHAND, T. A., DE HAAN, L. & ZHOU, C. (8). On spatial etremes: With application to a rainfall problem. Annals of Applied Statistics, CHILÈS, J.-P. & DELFINER, P. (1999). Geostatistics: Modelling Spatial Uncertainty. New York: Wiley. COOLEY, D., NAVEAU, P. & PONCET, P. (6). Variograms for spatial ma-stable random fields. In Dependence in Probability and Statistics, P. Bertail, P. Soulier, P. Doukhan, P. Bickel, P. Diggle, S. Fienberg, U. Gather, I. Olkin & S. Zeger, eds., vol. 187 of Lecture Notes in Statistics. Springer New York, pp DAVIS, R. & RESNICK, S. (1989). Basic properties and prediction of ma-arma processes. Advances in Applied Probability 1, DAVIS, R. & RESNICK, S. (1993). Prediction of stationary ma-stable processes. Annals Of Applied Probability 3, DAVISON, A. C. & GHOLAMREZAEE, M. M. (1). Geostatistics of etremes. Proceedings of the Royal Society A: Mathematical, Physical and Engineering Science 468, DAVISON, A. C., PADOAN, S. A. & RIBATET, M. (1). Statistical modelling of spatial etremes. Statistical Science 7, DE HAAN, L. (1984). A spectral representation for ma-stable processes. The Annals of Probability 1, DE HAAN, L. & PEREIRA, T. T. (6). Spatial etremes: Models for the stationary case. The Annals of Statistics 34, DE HAAN, L. & FEREIRA, A. (6). Etreme Value Theory: An Introduction. Springer, New-York. GENTON, M. G., MA, Y. & SANG, H. (11). On the likelihood function of Gaussian ma-stable processes. Biometrika 98, GENZ, A. (199). Numerical computation of multivariate normal probabilities. J. Comp. Graph Stat 1, KABLUCHKO, Z., SCHLATHER, M. & DE HAAN, L. (9). Stationary ma-stable fields associated to negative definite functions. Ann. Prob. 37,
22 C. DOMBRY, F. ÉYI-MINKO AND M. RIBATET KABLUCHKO, Z. & SCHLATHER, M. (1). Ergodic properties of ma-infinitely divisible processes. Stochastic Processes and their Applications, 1, LIU, J., WONG, W. & KONG, A. (1995). Correlation structure and convergence rate of the Gibbs sampler with various scans. Journal of the Royal Statistical Society Series B 57, MANDELBROT, B. (198). The Fractal Geometry of Nature. New York: W. H. Freeman. OESTING, M., KABLUCHKO, Z. & SCHLATHER, M. (1). Simulation of Brown Resnick processes. Etremes 15, PADOAN, S. A., RIBATET, M. & SISSON, S. (1). Likelihood-based inference for ma-stable processes. Journal of the American Statistical Association 15, R DEVELOPMENT CORE TEAM (1). R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing, Vienna, Austria. ISBN RIBATET, M. (11). SpatialEtremes: Modelling Spatial Etremes. R package version SCHLATHER, M. (). Models for stationary ma-stable random fields. Etremes 5, SCHLATHER, M. (1). RandomFields: Simulation and Analysis of Random Fields. R package version..54. SCHLATHER, M. & TAWN, J. (3). A dependence measure for multivariate and spatial etremes: Properties and inference. Biometrika 9, WANG, Y. (1). malinear: Conditional sampling for ma-linear models. R package version 1.. WANG, Y. & STOEV, S. A. (11). Conditional sampling for spectrally discrete ma-stable random fields. Advances in Applied Probability 43,
Conditional simulations of max-stable processes
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é
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 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 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 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 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 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 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 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 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 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 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 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 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 informationJournal de la Société Française de Statistique. Spatial extremes: Max-stable processes at work
Journal de la Société Française de Statistique Submission Spatial etremes: Ma-stable processes at work Titre: Etrêmes spatiau : les processus ma-stables en pratique Ribatet Mathieu 1 Abstract: Since many
More informationStatistics of Extremes
Statistics of Extremes Anthony Davison c 29 http://stat.epfl.ch Multivariate Extremes 184 Europe............................................................. 185 Componentwise maxima..................................................
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 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 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 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 EXTREME VALUE ANALYSIS OF A SPATIAL PROCESS
REVSTAT Statistical Journal Volume 6, Number 1, March 008, 71 81 ON EXTREME VALUE ANALYSIS OF A SPATIAL PROCESS Authors: Laurens de Haan Erasmus University Rotterdam and University Lisbon, The Netherlands
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 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 informationST 740: Markov Chain Monte Carlo
ST 740: Markov Chain Monte Carlo Alyson Wilson Department of Statistics North Carolina State University October 14, 2012 A. Wilson (NCSU Stsatistics) MCMC October 14, 2012 1 / 20 Convergence Diagnostics:
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 informationSPATIAL ASSESSMENT OF EXTREME SIGNIFICANT WAVES HEIGHTS IN THE GULF OF LIONS
SPATIAL ASSESSMENT OF EXTREME SIGNIFICANT WAVES HEIGHTS IN THE GULF OF LIONS Romain Chailan 1, Gwladys Toulemonde 2, Frederic Bouchette 3, Anne Laurent 4, Florence Sevault 5 and Heloise Michaud 6 In the
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 informationReview. DS GA 1002 Statistical and Mathematical Models. Carlos Fernandez-Granda
Review DS GA 1002 Statistical and Mathematical Models http://www.cims.nyu.edu/~cfgranda/pages/dsga1002_fall16 Carlos Fernandez-Granda Probability and statistics Probability: Framework for dealing with
More informationThe Poisson storm process Extremal coefficients and simulation
The Poisson storm process Extremal coefficients and simulation C. Lantuéjoul 1, J.N. Bacro 2 and L. Bel 3 1 MinesParisTech 2 Université de Montpellier II 3 AgroParisTech 1 Storm process Introduced by Smith
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 informationMobile Robot Localization
Mobile Robot Localization 1 The Problem of Robot Localization Given a map of the environment, how can a robot determine its pose (planar coordinates + orientation)? Two sources of uncertainty: - observations
More informationSub-kilometer-scale space-time stochastic rainfall simulation
Picture: Huw Alexander Ogilvie Sub-kilometer-scale space-time stochastic rainfall simulation Lionel Benoit (University of Lausanne) Gregoire Mariethoz (University of Lausanne) Denis Allard (INRA Avignon)
More informationOn prediction and density estimation Peter McCullagh University of Chicago December 2004
On prediction and density estimation Peter McCullagh University of Chicago December 2004 Summary Having observed the initial segment of a random sequence, subsequent values may be predicted by calculating
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 HIERARCHICAL MAX-STABLE SPATIAL MODEL FOR EXTREME PRECIPITATION
Submitted to the Annals of Applied Statistics arxiv: arxiv:0000.0000 A HIERARCHICAL MAX-STABLE SPATIAL MODEL FOR EXTREME PRECIPITATION By Brian J. Reich and Benjamin A. Shaby North Carolina State University
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 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 informationBasics of Point-Referenced Data Models
Basics of Point-Referenced Data Models Basic tool is a spatial process, {Y (s), s D}, where D R r Chapter 2: Basics of Point-Referenced Data Models p. 1/45 Basics of Point-Referenced Data Models Basic
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 informationCRPS M-ESTIMATION FOR MAX-STABLE MODELS WORKING MANUSCRIPT DO NOT DISTRIBUTE. By Robert A. Yuen and Stilian Stoev University of Michigan
arxiv: math.pr/ CRPS M-ESTIMATION FOR MAX-STABLE MODELS WORKING MANUSCRIPT DO NOT DISTRIBUTE By Robert A. Yuen and Stilian Stoev University of Michigan Max-stable random fields provide canonical models
More informationHigh-Order Composite Likelihood Inference for Max-Stable Distributions and Processes
Supplementary materials for this article are available online. Please go to www.tandfonline.com/r/jcgs High-Order Composite Likelihood Inference for Max-Stable Distributions and Processes Stefano CASTRUCCIO,
More informationUsing spatial and spatial-extreme statistics to characterize snow avalanche cycles
Available online at www.sciencedirect.com Procedia Environmental Sciences 7 7 (011) 4 9 1 11 Spatial Statistics 011 - Mapping Global Change Using spatial and spatial-extreme statistics to characterize
More informationLecture 3: September 10
CS294 Markov Chain Monte Carlo: Foundations & Applications Fall 2009 Lecture 3: September 10 Lecturer: Prof. Alistair Sinclair Scribes: Andrew H. Chan, Piyush Srivastava Disclaimer: These notes have not
More informationMax-stable processes and annual maximum snow depth
Max-stable processes and annual maximum snow depth Juliette Blanchet blanchet@slf.ch 6th International Conference on Extreme Value Analysis June 23-26, 2009 Outline Motivation Max-stable process - Schlather
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 informationA space-time skew-t model for threshold exceedances
BIOMETRICS 000, 1 20 DOI: 000 000 0000 A space-time skew-t model for threshold exceedances Samuel A Morris 1,, Brian J Reich 1, Emeric Thibaud 2, and Daniel Cooley 2 1 Department of Statistics, North Carolina
More informationNonstationary spatial process modeling Part II Paul D. Sampson --- Catherine Calder Univ of Washington --- Ohio State University
Nonstationary spatial process modeling Part II Paul D. Sampson --- Catherine Calder Univ of Washington --- Ohio State University this presentation derived from that presented at the Pan-American Advanced
More informationA Practitioner s Guide to Generalized Linear Models
A Practitioners Guide to Generalized Linear Models Background The classical linear models and most of the minimum bias procedures are special cases of generalized linear models (GLMs). GLMs are more technically
More informationA Bayesian Hierarchical Model for Spatial Extremes with Multiple Durations
A Bayesian Hierarchical Model for Spatial Extremes with Multiple Durations Yixin Wang a, Mike K. P. So a a The Hong Kong University of Science and Technology, Clear Water Bay, Kowloon, Hong Kong Abstract
More informationApproximate inference, Sampling & Variational inference Fall Cours 9 November 25
Approimate inference, Sampling & Variational inference Fall 2015 Cours 9 November 25 Enseignant: Guillaume Obozinski Scribe: Basile Clément, Nathan de Lara 9.1 Approimate inference with MCMC 9.1.1 Gibbs
More informationSpatio-temporal precipitation modeling based on time-varying regressions
Spatio-temporal precipitation modeling based on time-varying regressions Oleg Makhnin Department of Mathematics New Mexico Tech Socorro, NM 87801 January 19, 2007 1 Abstract: A time-varying regression
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 informationPUBLICATIONS. Water Resources Research. Assessing the performance of the independence method in modeling spatial extreme rainfall
PUBLICATIONS Water Resources Research RESEARCH ARTICLE Key Points: The utility of the independence method is explored in modeling spatial extremes Four spatial models are compared in estimating marginal
More informationIntroduction to Geostatistics
Introduction to Geostatistics Abhi Datta 1, Sudipto Banerjee 2 and Andrew O. Finley 3 July 31, 2017 1 Department of Biostatistics, Bloomberg School of Public Health, Johns Hopkins University, Baltimore,
More informationModelling trends in the ocean wave climate for dimensioning of ships
Modelling trends in the ocean wave climate for dimensioning of ships STK1100 lecture, University of Oslo Erik Vanem Motivation and background 2 Ocean waves and maritime safety Ships and other marine structures
More informationRESEARCH REPORT. Estimation of sample spacing in stochastic processes. Anders Rønn-Nielsen, Jon Sporring and Eva B.
CENTRE FOR STOCHASTIC GEOMETRY AND ADVANCED BIOIMAGING www.csgb.dk RESEARCH REPORT 6 Anders Rønn-Nielsen, Jon Sporring and Eva B. Vedel Jensen Estimation of sample spacing in stochastic processes No. 7,
More informationAlternative implementations of Monte Carlo EM algorithms for likelihood inferences
Genet. Sel. Evol. 33 001) 443 45 443 INRA, EDP Sciences, 001 Alternative implementations of Monte Carlo EM algorithms for likelihood inferences Louis Alberto GARCÍA-CORTÉS a, Daniel SORENSEN b, Note a
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 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 informationSupplement to A Hierarchical Approach for Fitting Curves to Response Time Measurements
Supplement to A Hierarchical Approach for Fitting Curves to Response Time Measurements Jeffrey N. Rouder Francis Tuerlinckx Paul L. Speckman Jun Lu & Pablo Gomez May 4 008 1 The Weibull regression model
More 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 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 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 informationDefault Priors and Effcient Posterior Computation in Bayesian
Default Priors and Effcient Posterior Computation in Bayesian Factor Analysis January 16, 2010 Presented by Eric Wang, Duke University Background and Motivation A Brief Review of Parameter Expansion Literature
More informationCSC 2541: Bayesian Methods for Machine Learning
CSC 2541: Bayesian Methods for Machine Learning Radford M. Neal, University of Toronto, 2011 Lecture 3 More Markov Chain Monte Carlo Methods The Metropolis algorithm isn t the only way to do MCMC. We ll
More informationA Statistical Input Pruning Method for Artificial Neural Networks Used in Environmental Modelling
A Statistical Input Pruning Method for Artificial Neural Networks Used in Environmental Modelling G. B. Kingston, H. R. Maier and M. F. Lambert Centre for Applied Modelling in Water Engineering, School
More informationAvailable online at ScienceDirect. Procedia Environmental Sciences 27 (2015 ) Spatial Statistics 2015: Emerging Patterns
Available online at www.sciencedirect.com ScienceDirect Procedia Environmental Sciences 27 (2015 ) 75 82 Spatial Statistics 2015: Emerging Patterns Inferring spatio-temporal patterns in extreme snowfall
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 informationA User s Guide to the SpatialExtremes Package
A User s Guide to the SpatialExtremes Package Mathieu Ribatet Copyright 2009 Chair of Statistics École Polytechnique Fédérale de Lausanne Switzerland Contents Introduction 1 1 An Introduction to Max-Stable
More informationMonte Carlo Studies. The response in a Monte Carlo study is a random variable.
Monte Carlo Studies The response in a Monte Carlo study is a random variable. The response in a Monte Carlo study has a variance that comes from the variance of the stochastic elements in the data-generating
More informationModels for spatial data (cont d) Types of spatial data. Types of spatial data (cont d) Hierarchical models for spatial data
Hierarchical models for spatial data Based on the book by Banerjee, Carlin and Gelfand Hierarchical Modeling and Analysis for Spatial Data, 2004. We focus on Chapters 1, 2 and 5. Geo-referenced data arise
More informationSTATISTICS 174: APPLIED STATISTICS FINAL EXAM DECEMBER 10, 2002
Time allowed: 3 HOURS. STATISTICS 174: APPLIED STATISTICS FINAL EXAM DECEMBER 10, 2002 This is an open book exam: all course notes and the text are allowed, and you are expected to use your own calculator.
More informationApplications of Randomized Methods for Decomposing and Simulating from Large Covariance Matrices
Applications of Randomized Methods for Decomposing and Simulating from Large Covariance Matrices Vahid Dehdari and Clayton V. Deutsch Geostatistical modeling involves many variables and many locations.
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 informationSupplemental Information Likelihood-based inference in isolation-by-distance models using the spatial distribution of low-frequency alleles
Supplemental Information Likelihood-based inference in isolation-by-distance models using the spatial distribution of low-frequency alleles John Novembre and Montgomery Slatkin Supplementary Methods To
More informationA Gaussian state-space model for wind fields in the North-East Atlantic
A Gaussian state-space model for wind fields in the North-East Atlantic Julie BESSAC - Université de Rennes 1 with Pierre AILLIOT and Valï 1 rie MONBET 2 Juillet 2013 Plan Motivations 1 Motivations 2 Context
More information14.30 Introduction to Statistical Methods in Economics Spring 2009
MIT OpenCourseWare http://ocw.mit.edu 4.0 Introduction to Statistical Methods in Economics Spring 009 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms.
More informationAsymptotic inference for a nonstationary double ar(1) model
Asymptotic inference for a nonstationary double ar() model By SHIQING LING and DONG LI Department of Mathematics, Hong Kong University of Science and Technology, Hong Kong maling@ust.hk malidong@ust.hk
More informationState Space and Hidden Markov Models
State Space and Hidden Markov Models Kunsch H.R. State Space and Hidden Markov Models. ETH- Zurich Zurich; Aliaksandr Hubin Oslo 2014 Contents 1. Introduction 2. Markov Chains 3. Hidden Markov and State
More informationIntroduction to Spatial Data and Models
Introduction to Spatial Data and Models Sudipto Banerjee 1 and Andrew O. Finley 2 1 Department of Forestry & Department of Geography, Michigan State University, Lansing Michigan, U.S.A. 2 Biostatistics,
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 informationThreshold modeling of extreme spatial rainfall
WATER RESOURCES RESEARCH, VOL.???, XXXX, DOI:10.1029/, Threshold modeling of extreme spatial rainfall E. Thibaud, 1 R. Mutzner, 2 A.C. Davison, 1 Abstract. We propose an approach to spatial modeling of
More informationSTA 4273H: Statistical Machine Learning
STA 4273H: Statistical Machine Learning Russ Salakhutdinov Department of Computer Science! Department of Statistical Sciences! rsalakhu@cs.toronto.edu! h0p://www.cs.utoronto.ca/~rsalakhu/ Lecture 7 Approximate
More informationE(x i ) = µ i. 2 d. + sin 1 d θ 2. for d < θ 2 0 for d θ 2
1 Gaussian Processes Definition 1.1 A Gaussian process { i } over sites i is defined by its mean function and its covariance function E( i ) = µ i c ij = Cov( i, j ) plus joint normality of the finite
More informationA spatiotemporal auto-regressive moving average model for solar radiation
A spatiotemporal auto-regressive moving average model for solar radiation C.A. Glasbey and D.J. Allcroft Biomathematics and Statistics Scotland King s Buildings, Edinburgh, EH9 3JZ, Scotland January 11,
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 informationSpatial smoothing using Gaussian processes
Spatial smoothing using Gaussian processes Chris Paciorek paciorek@hsph.harvard.edu August 5, 2004 1 OUTLINE Spatial smoothing and Gaussian processes Covariance modelling Nonstationary covariance modelling
More informationUsing Estimating Equations for Spatially Correlated A
Using Estimating Equations for Spatially Correlated Areal Data December 8, 2009 Introduction GEEs Spatial Estimating Equations Implementation Simulation Conclusion Typical Problem Assess the relationship
More informationIntroduction to Spatial Data and Models
Introduction to Spatial Data and Models Sudipto Banerjee 1 and Andrew O. Finley 2 1 Biostatistics, School of Public Health, University of Minnesota, Minneapolis, Minnesota, U.S.A. 2 Department of Forestry
More informationModelling Dependence in Space and Time with Vine Copulas
Modelling Dependence in Space and Time with Vine Copulas Benedikt Gräler, Edzer Pebesma Abstract We utilize the concept of Vine Copulas to build multi-dimensional copulas out of bivariate ones, as bivariate
More informationNUMERICAL COMPUTATION OF THE CAPACITY OF CONTINUOUS MEMORYLESS CHANNELS
NUMERICAL COMPUTATION OF THE CAPACITY OF CONTINUOUS MEMORYLESS CHANNELS Justin Dauwels Dept. of Information Technology and Electrical Engineering ETH, CH-8092 Zürich, Switzerland dauwels@isi.ee.ethz.ch
More informationSolution of Partial Differential Equations
Solution of Partial Differential Equations Introduction and Problem Statement Combination of Variables R. Shankar Subramanian We encounter partial differential equations routinely in transport phenomena.
More informationMULTIDIMENSIONAL COVARIATE EFFECTS IN SPATIAL AND JOINT EXTREMES
MULTIDIMENSIONAL COVARIATE EFFECTS IN SPATIAL AND JOINT EXTREMES Philip Jonathan, Kevin Ewans, David Randell, Yanyun Wu philip.jonathan@shell.com www.lancs.ac.uk/ jonathan Wave Hindcasting & Forecasting
More 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 informationBayesian inference & process convolution models Dave Higdon, Statistical Sciences Group, LANL
1 Bayesian inference & process convolution models Dave Higdon, Statistical Sciences Group, LANL 2 MOVING AVERAGE SPATIAL MODELS Kernel basis representation for spatial processes z(s) Define m basis functions
More informationBayesian Linear Regression
Bayesian Linear Regression Sudipto Banerjee 1 Biostatistics, School of Public Health, University of Minnesota, Minneapolis, Minnesota, U.S.A. September 15, 2010 1 Linear regression models: a Bayesian perspective
More informationLecture 2 APPLICATION OF EXREME VALUE THEORY TO CLIMATE CHANGE. Rick Katz
1 Lecture 2 APPLICATION OF EXREME VALUE THEORY TO CLIMATE CHANGE Rick Katz Institute for Study of Society and Environment National Center for Atmospheric Research Boulder, CO USA email: rwk@ucar.edu Home
More 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 informationMultistate Modeling and Applications
Multistate Modeling and Applications Yang Yang Department of Statistics University of Michigan, Ann Arbor IBM Research Graduate Student Workshop: Statistics for a Smarter Planet Yang Yang (UM, Ann Arbor)
More information