Max-stable processes: Theory and Inference
|
|
- Sharyl Harrell
- 5 years ago
- Views:
Transcription
1 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 2/39 Outline Introduction Max-stable processes Spatial dependence of max-stable processes Inference for max-stable processes Model Selection Application References
3 3/39 Outline Introduction Max-stable processes Spatial dependence of max-stable processes Inference for max-stable processes Model Selection Application References
4 4/39 EVT has been applied to many environmental processes - rainfall, temperature, river flow, wind, fire sizes,... However many of these processes are spatial Univariate/Multivariate analysis of extremes are not enough Problem Let Y 1,...,Y n be i.i.d. replicates of a stochastic process Y defined on K R d. We might be interested in estimating Pr[Y (x) > f (x),x K], f some deterministic function Pr[ K 1 Y (x)dx > y], K 1 K Pr[Y (x 2 ) > y 2 Y (x 1 ) > y 1 ], (x 1,x 2 ) K K
5 5/39 Outline Introduction Max-stable processes Spatial dependence of max-stable processes Inference for max-stable processes Model Selection Application References
6 Let X 1,X 2,... be independent replications of a stochastic processes on R d with continuous sample path We are interested in the limiting process { max i=1,...,n X i (x) b n (x) a n (x) } x K d {Y (x)} x R d (1) for some non-degenerate process Y, a n (x) > 0 and b n (x) continuous functions. de Haan shows that the class of limit process in (1) coincides with the class of max-stable processes. Definition A stochastic process Y with non-degenerate marginals is called max-stable if there are continuous function a n (x) > 0 and b n (x) such that if Y 1,Y 2,... are i.i.d. copies of Y then n = 1,2,... (2) Y i b n max i=1,...,n a n d = Y, 6/39
7 7/39 It is often more convenient to work with specified margins Most authors work with unit Fréchet margins i.e. Pr[Y (x) z] = exp( 1/z) Definition (Special case) A stochastic process Z is said max-stable with unit Fréchet margins if max i=1,...,n n 1 d Z i = Z, n = 1,2,... (3) Theorem (Spectral representation) Let {ξ i } i 1 be the points of a Poisson process on R + with intensity dλ(ξ) = ξ 2 dξ, and {Y i ( )} i 1 be i.i.d. replicates of a process on R d such that E[max{0,Y (x)}] = 1. Then Z(x) = max i 1 ξ i max{0,y i (x)} is a max-stable process with unit Fréchet margins.
8 8/39 Using the spectral representation, it is possible get parametric max-stable models through the definition of the process Y In the next slides, I will introduce some models
9 9/39 Smith s model [Smith, 1990] Let Y i (x) = ϕ(x X i ) where {X i } i 1 are the points of a homogeneous Poisson process and ϕ is the zero mean multivariate normal density with covariance matrix Σ, both on R d. Remark This model is often referred to as the storm profile model Lets focus on rainfall events. You can think about ξ i as the intensity of the i-th storm, {X i } as the storm centres ϕ as a deterministic function that drives how the amount of rain decreases from the storm centre. ξ iϕ(x Xi) X 3 X 1 X 2
10 10/39 Schlather s model [Schlather, 2002] Later Schlather proposed to take Y i (x) = 2πY,i (x) where Y,i are independent replications of a standard Gaussian process with correlation ρ Remark This model has no deterministic shape for the storms......but random ones Storm structure is driven by the correlation function ρ(h) But it has some drawbacks - see later ξ imax{0, Yi(x)}
11 11/39 Other models One can define various models in the spirit of the Schlather one The most popular are Geometric gaussian model Y i (x) = exp where ε is a standard Gaussian process. } {σε(x) σ2 2 Brown-Resnick processes { } Y i (x) = exp ε(x) σ2 (x) 2 where ε is a Gaussian process with stationary increments, σ 2 (x) = Var[ε(x)]. We ll see later what are the assets and drawbacks of these models
12 12/39 y y x x y y x Not yet implemented! x Figure: Top left: Smith. Top right: Schlather. Bottom left: Geometric. Bottom right: Brown Resnick. Gumbel margins. Simulations were performed using the SpatialExtremes package.
13 13/39 Outline Introduction Max-stable processes Spatial dependence of max-stable processes Inference for max-stable processes Model Selection Application References
14 14/39 Motivation The modelling of univariate extreme events is limited if the process is spatial by nature The knowledge of dependence between extreme events occuring at several locations is essential In this section, we will introduce some notions describing the spatial structure of extreme events
15 15/39 Any multivariate extreme value distribution has the form Pr[Z(x 1 ) z 1,...,Z(x k ) z k ] = exp { V (z 1,...,z k )} (4) where V has some homogeneity properties depending on the margins. If we suppose unit Fréchet margins, V is homogeneous of order 1 so that log (Pr[Z(x 1 ) z,...,z(x k ) z]) = V (z,...,z) hom. = θ k z where θ k = V (1,...,1) is called the extremal coefficient. θ k is a measure of dependence between extremes and satisfies 1 θ k k θ k = 1 complete dependence θ k = k independence For max-stable processes, an important measure of dependence is the extremal coefficient function θ(x 1 x 2 ) = z log Pr[Z(x 1 ) z,z(x 2 ) z] (5)
16 16/39 It is easy to see that the extremal coefficient functions for the models introduced are ( ) (x1 x Smith θ(x 1 x 2 ) = 2Φ 2 ) T Σ 1 (x 1 x 2 ) 2 1 ρ(x 1 x 2 ) Schlather θ(x 1 x 2 ) = 1 + ( 2 Geometric Gaussian θ(x 1 x 2 ) = 2Φ 2 ( ) γ(x1 x Brown Resnick θ(x 1 x 2 ) = 2Φ 2 ) σ 2 (1 ρ(x 1 x 2 )) The Schlather model has an upper bound of in R 2 and in R 3. In addition, as most of the correlation functions take only positive values, the upper bound is 1 + 1/2. Independence never reached! The geometric gaussian model has an upper bound of 2Φ(σ/ 2). Hence, if σ is large enough, independence is reached - at least numerically. 2 )
17 17/39 Variogram-based approaches In conventional geostatistics, the variogram determines the degree of spatial dependence γ(x 1 x 2 ) = 1 2 Var [Y (x 1) Y (x 2 )] = σ 2 {1 ρ(x 1 x 2 )} But for max-stable processes having unit Fréchet margins, E[Z(x)] = Var[Z(x)] = + Hence we need variogram-based approaches for max-stable random fields (6)
18 18/39 Let Z be a stationary max-stable random fields with unit Fréchet margins Let F(x) = exp( 1/z) The F-madogram [Cooley et al., 2006] is defined as Then we have ν F (x 1 x 2 ) = 1 2 E [ F {Z(x 1)} F {Z(x 2 )} ] (7) 2ν F (x 1 x 2 ) = θ(x 1 x 2 ) 1 θ(x 1 x 2 ) + 1 (8) or conversely θ(x 1 x 2 ) = 1 + 2ν F(x 1 x 2 ) 1 2ν F (x 1 x 2 ) (9) For the proof just note that x y = 2max{x,y} x y and you re set.
19 19/39 Outline Introduction Max-stable processes Spatial dependence of max-stable processes Inference for max-stable processes Model Selection Application References
20 20/39 A first statement It is almost impossible to get analytical CDF for Z( ) Indeed as Z(x) = max i ξ i Y + i (x), Y + i (x) = max{0,y i (x)} Hence to get the k-variate CDF we need to compute F(z 1,...,z k ) = Pr[max ξ i Y + i i (x j ) z j,j = 1,...,k] [ = Pr ξ i z ] j Y +, i,j = 1,...,k i (x j ) [ ( ) z j = exp I ξ > min R d R j y j + ξ 2 dξdp Y (y 1,...,y k ) [ = exp max R d j y + j z j dp Y (y 1,...,y k ) The number of possible cases becomes quickly intractable when k gets large Only the bivariate densities can easily be derived analytically ]
21 21/39 Bivariate CDF Smith [1990] j Pr[Z(x 1 ) z 1, Z(x 2 ) z 2 ] = exp 1 a Φ z a log z «2 1 a Φ z 1 z a log z «ff 1 z 2 where Φ is the standard normal CDF and a 2 = x T Σ 1 x. Schlather [2002] ( Pr[Z(x 1 ) z 1, Z(x 2 ) z 2 ] = exp «s!) z 1 z {ρ(h) + 1} 2 z 1 z 2 (z 1 + z 2 ) 2 where h = x 1 x 2. Geometric Gaussian and Brown Resnick processses The bivariate CDFs are the same as for the Smith model where Geom a 2 = 2σ 2 {1 ρ(h)} B. R. a 2 = γ(x 2 x 1 ), γ( ) is the variogram of ǫ( ) If γ(h) h 2, h 0 then B. R. Smith (isotropic).
22 22/39 Definition Composite likelihoods Let {f (y;θ),y Y,θ Θ} a parametric statistical model, where Y R K, Θ R p, K 1 and p 1. Consider a set of (marginal or conditional) events {A i : A i F,i I }, where I N and F is a σ-algebra on Y. A log-composite likelihood is defined as l c (θ;y) = i I w i log f (y A i ;θ) where f (y A i ;θ) = f ({y j Y : y j A i } ;θ), y = (y 1,...,y n ) and {w i,i I } is a set of suitable weights. In a nutshell, log-composite likelihood are just a linear combination of (smaller) log-likelihood entities
23 23/39 Why does it works? First, note that the full likelihood is a special case of composite likelihood For i being fixed, log f (y A i ;θ) is a valid log-likelihood Thus leading to an unbiased estimating equation log f (y A i ;θ) = 0 Finally l c (θ;y) = i I w i log f (y A i ;θ) = 0 is unbiased - as a linear combination of unbiased estimating equations For max-stable processes, as only the bivariate densities are known we will consider the pairwise likelihood l p (y;θ) = n log f (y (i) k,y(j) k ;θ) k=1 i<j
24 24/39 Asymptotics Instead of having n{ H(θ)} 1/2 d (ˆθ θ) N(0,Id p ), n + where H(θ) = E[ 2 l(θ;y)], (M 1/2 ) T M 1/2 = M 1 When we work under misspecification - which is the case when using composite likelihoods, we now have n{h(θ)j(θ) 1 H(θ)} 1/2 d (ˆθ θ) N(0,Id p ), n + where J(θ) = Var[ l(θ; Y)] Note that when the model is correctly specified, the so-called second Bartlett idendity holds i.e. E[ 2 l(θ;y)] + Var[ l(θ;y)] = 0 (10) In other words, we have H(θ) = J(θ) and H(θ)J(θ) 1 H(θ) = H(θ) 1 Beware I ll use this abuse of notation several times for shortness
25 25/39 Pairwise likelihood The log-pairwise likelihood is given by l p (θ;y) = n K 1 K k=1 i=1 j=i+1 log f (y (i) k,y(j) k ;θ) (11) Under the usual regularity conditions for the MLE and that θ is identifiable from the densities in (11), n{h(θ)j(θ) 1 H(θ)} 1/2 (ˆθ p θ) d N(0,Id p ), n + H(θ) and J(θ) can be consistently estimated by Ĥ(ˆθ p ) = Ĵ(ˆθ p ) = 1 n 1 n n 2 l p (ˆθ p ;y i ) i=1 n { l c (ˆθ p ;y i ) l c (ˆθ p ;y i ) T} i=1
26 From unit Fréchet margins to ordinary GEV ones Most of the max-stable processes are given with unit Fréchet margins Usually the margins will be GEV - not necessary identical Lets consider the following mapping ( t : y 1 + ξ y µ ) 1/ξ (12) σ If Y GEV(µ,σ,ξ), then Z = t(y ) Fréchet(1) Hence if Y ( ) is a max-stable process such that Y (x i ) GEV(µ i,σ i,ξ i ), the pairwise likelihood for max-stable processes defined with unit Fréchet margins becomes l p (θ;y) = where z (i) k n k=1 i<j { log f (z (i) k,z(j) k } ;θ) + log J{t(y(i) k )}J{t(y(j) k )} = t(y (i) k ) = {1 + ξ i(y (i) k µ i )/σ i } 1/ξ i and (13) J{t(y (i) k )} is the Jacobian of the transformation t 26/39
27 Thus one can define trend surfaces for the evolution of the GEV parameters in space i.e. and similarly for σ and ξ µ = X µ β µ (14) so that the parameters for our max-stable process have the form ψ = (θ,β µ,β σ,β ξ ), where θ is the parameter for the simple max-stable process - i.e. with unit Fréchet margins. The definition of trend surfaces allows 1. predictions of the GEV parameters, quantiles x K 2. correct estimates of the standard errors 2 2 the two step procedure i.e. first margins, then dependence leads to std. err. underestimations 27/39
28 28/39 Outline Introduction Max-stable processes Spatial dependence of max-stable processes Inference for max-stable processes Model Selection Application References
29 29/39 Information criteria When several models M0, M1,...are fitted to our data, one would prefer the one minimising AIC = 2l(ˆθ;y) + 2p where p is the number of parameters to be estimated. Under misspecification, one should use TIC = 2l(ˆθ;y) 2tr {Ĵ(ˆθ)Ĥ(ˆθ) 1} as the 2nd Bartlett idendity is not true anymore Note that if the model is correctly specified J(θ) = H(θ), E [ 2 l(θ;y) ] + Var [ l(θ;y)] = 0 so that TIC = 2l(ˆθ;y) + 2tr {I p } = AIC
30 30/39 Likelihood ratio When nested models are considered, it is always preferable to use the likelihood ratio statistic - cf. Neyman Pearson Let (ˆκ T, ˆφ T ) T be the MLE and (κ T 0, ˆφ T κ 0 ) T the MLE under the restriction κ = κ 0. Then { } W (κ 0 ) = 2 l(ˆκ, ˆφ) l(κ 0, ˆφ κ0 ) χ 2 p, n + (15) where p is the dimension of κ 0. Under mis-specification, this result differs { W (κ 0 ) = 2 l(ˆκ, ˆφ) l(κ 0, ˆφ } κ0 ) p λ i X i, n + i=1 where X i iid χ 2 1, λ i are the eigenvalues of (H 1 JH 1 ) κ { (H 1 ) κ } 1 and M κ is the submatrix of M corresponding to the elements of κ. (16)
31 31/39 The problem with (16) is that the limiting distribution is generally unknown. As an approximation, Rotnitski and Jewell suggest to treat pw (κ 0 )/ p i=1 λ i as a χ 2 p ] E Var [ pw (κ0 ) p i=1 λ i [ ] pw (κ0 ) p i=1 λ i p, n + 2p, n + Results based on quadratic forms for Normal random variables might be more accurate Another approach consists in adjusting the likelihood so that the usual convergence to a χ 2 p is met
32 32/39 Lets consider the following adjustment of the log-likelihood l A (θ) = l(θ ), θ = ˆθ + C(θ ˆθ) (17) The following adjustment has to interesting properties 1. ˆθ maximises l A 2. At ˆθ, l A has curvature C T H(θ 0 )C To have the usual χ 2 p convergence, we need that l A has curvature H(θ 0 )J(θ 0 ) 1 H(θ 0 ) at ˆθ If we choose C = {H(θ 0 ) 1/2 } 1 {H(θ 0 )J(θ 0 ) 1 H(θ 0 )} 1/2 (18) then C T H(θ 0 )C = H(θ 0 )J(θ 0 ) 1 H(θ 0 ) and W A (κ 0 ) = 2 Note the presence of ˆφ κ0,a! { } l A (ˆκ, ˆφ) l A (κ 0, ˆφ κ0,a) χ 2 p, n +
33 33/39 Outline Introduction Max-stable processes Spatial dependence of max-stable processes Inference for max-stable processes Model Selection Application References
34 34/39 Extreme precipitations around Zurich Elevation (meters) 35 stations (47 obs./stations) kept for validation κ = R 2 10 km Zurich Smith Covariance σ 11 σ 12 σ 22 h h + DoF l p TIC Isotropic 259 (45) 0 ( ) σ 22 = σ Anisotropic 251 (46) 64 (13) 290 (50) Schlather Covariance λ κ h h + DoF l p TIC Whittle 38.1 (10.7) 0.45 ( ) Cauchy 8.0 ( 2.2) 0.34 (0.16) Stable 34.8 (11.5) 0.95 (0.16) Exponential 34.1 ( 9.0) Bessel 0.95 (0.11) 150 ( ) Geometric Covariance σ 2 λ κ h h + DoF l p TIC Whittle 5.77 (8.64) (456.3) 0.40 ( ) Cauchy 36.0 (673) 5.19 (1.54) 9e 3 (0.16) Stable 5.74 (14.9) 237 (961) 0.79 (0.17) Exponential 2.42 (0.93) 53.2 (18.4) Bessel 1.28 (0.30) 0.61 (0.03) 150 ( ) 7.7 never h α Brown Resnick: γ(h) = λ Covariance λ α h h + DoF l p TIC (3.43) 0.74 (0.07)
35 35/39 Schlather Geometric Brown Resnick θ(h) Whittle Cauchy θ(h) h Whittle Cauchy h θ(h) Brown Resnick Smith h Figure: Comparison between the empirical extremal coefficients and the fitted extremal coefficient functions. Grey points : fitting data set. Black points : Validation data set. Although there seems to be a slight underestimation of the dependence, the fits look reasonably well The empirical estimates suggest the use of a nugget effect that has no theoretical meaning...
36 36/39 Comments on the Geometric Gaussian model θ(h) λ = 200 λ = 400 λ = 600 λ = 800 λ = h σ Figure: Impact of fixing the scale parameter λ in the Geometric gaussian model. The scale λ and variance σ 2 are far from being orthogonal - identifiability problem? Remember that σ 2 drives the upper bound of θ(h). Our application study does not cover the full range of dependence, difficulties to estimate both λ and σ 2 The safest strategy seems to fix λ and estimate σ 2 and hence the upper bound of θ(h). λ
37 Observed Observed Observed Model Model Model Observed Observed Observed Model Model Model Observed Observed Observed Model Model Model Figure: Model checking for the geometric models having a Whittle Matérn correlation function. Top row compares pairwise maxima estimated from the model and the observed ones for pairs of stations separated by 7 km (left), 45 km (middle) and 83 km. The middle row compares the observed and predicted minima (left), mean (middle) and maxima (right) for a group of 5 stations chosen randomly. The bottom row compares the observed and predicted minima (left), mean (middle) and maxima (right) for all the stations kept for model validation. The 95% confidence envelopes are also reported. The plot is transformed to the unit Gumbel scale using the empirical distribution function for the observations and the integral transform for the model. 37/39
38 Smith Schlather Geometric Precipitation (cm) Zurich Zurich Zurich Figure: One realisation of the best Smith, Schlather and geometric gaussian max-stable models. Smith s model leads unrealistic realisation while the Schlather and Geometric ones seem more plausible. Note how the Schlather model leads to larger range effects - as θ(h) 1 + 1/2 Simulations from fitted models are available (but difficult) and the modelling of spatial extremes in a subset of our region is therefore possible 38/39
39 39/39 Outline Introduction Max-stable processes Spatial dependence of max-stable processes Inference for max-stable processes Model Selection Application References
40 39/39 Chandler, R. E. and Bate, S. (2007). Inference for clustered data using the independence loglikelihood. Biometrika, 94(1): Cooley, D., Naveau, P., and Poncet, P. (2006). Variograms for spatial max-stable random fields. In Springer, editor, Dependence in Probability and Statistics, volume 187, pages Springer, New York, lecture notes in statistics edition. Davison, A. (2003). Statistical Models. Cambridge Series in Statistical and Probabilistic Mathematics. Cambridge University Press. de Haan, L. (1984). A spectral representation for max-stable processes. The Annals of Probability, 12(4): de Haan, L. and Pereira, T. (2006). Spatial extremes: Models for the stationary case. The Annals of Statistics, 34: Kent, J. (1982). Robust properties of likelihood ratio tests. Biometrika, 69: Naveau, P., Guillou, A., and Cooley, D. (2009). Modelling pairwise dependence of maxima in space. Biometrika, 96(1):1 17.
41 39/39 Padoan, S., Ribatet, M., and Sisson, S. (2009). Likelihood-based inference for max-stable processes. To appear in the Journal of the American Statistical Association (Theory & Methods). Ribatet, M. (2009). SpatialExtremes: Modelling Spatial Extremes. R package version Rotnitzky, A. and Jewell, N. (1990). Hypothesis testing of regression parameters in semiparametric generalized linear models for cluster correlated data. Biometrika, 77: Schlather, M. (2002). Models for stationary max-stable random fields. Extremes, 5(1): Schlather, M. and Tawn, J. (2003). A dependence measure for multivariate and spatial extremes: Properties and inference. Biometrika, 90(1): Smith, R. (1990). Max-stable processes and spatial extreme. Unpublished manuscript.
Extreme 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 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 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 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 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 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 informationThreshold estimation in marginal modelling of spatially-dependent non-stationary extremes
Threshold estimation in marginal modelling of spatially-dependent non-stationary extremes Philip Jonathan Shell Technology Centre Thornton, Chester philip.jonathan@shell.com Paul Northrop University College
More 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 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 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 informationOverview of Extreme Value Theory. Dr. Sawsan Hilal space
Overview of Extreme Value Theory Dr. Sawsan Hilal space Maths Department - University of Bahrain space November 2010 Outline Part-1: Univariate Extremes Motivation Threshold Exceedances Part-2: Bivariate
More 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 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 informationJoint work with Nottingham colleagues Simon Preston and Michail Tsagris.
/pgf/stepx/.initial=1cm, /pgf/stepy/.initial=1cm, /pgf/step/.code=1/pgf/stepx/.expanded=- 10.95415pt,/pgf/stepy/.expanded=- 10.95415pt, /pgf/step/.value required /pgf/images/width/.estore in= /pgf/images/height/.estore
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 informationLikelihood-based inference for max-stable processes
Likelihood-based inference for max-stable processes Simone Padoan, Mathieu Ribatet, Scott Sisson To cite this version: Simone Padoan, Mathieu Ribatet, Scott Sisson. Likelihood-based inference for max-stable
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 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 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 informationHANDBOOK OF APPLICABLE MATHEMATICS
HANDBOOK OF APPLICABLE MATHEMATICS Chief Editor: Walter Ledermann Volume VI: Statistics PART A Edited by Emlyn Lloyd University of Lancaster A Wiley-Interscience Publication JOHN WILEY & SONS Chichester
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 informationStat 5101 Lecture Notes
Stat 5101 Lecture Notes Charles J. Geyer Copyright 1998, 1999, 2000, 2001 by Charles J. Geyer May 7, 2001 ii Stat 5101 (Geyer) Course Notes Contents 1 Random Variables and Change of Variables 1 1.1 Random
More informationCovariance function estimation in Gaussian process regression
Covariance function estimation in Gaussian process regression François Bachoc Department of Statistics and Operations Research, University of Vienna WU Research Seminar - May 2015 François Bachoc Gaussian
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 informationStatistics for analyzing and modeling precipitation isotope ratios in IsoMAP
Statistics for analyzing and modeling precipitation isotope ratios in IsoMAP The IsoMAP uses the multiple linear regression and geostatistical methods to analyze isotope data Suppose the response variable
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 informationGauge Plots. Gauge Plots JAPANESE BEETLE DATA MAXIMUM LIKELIHOOD FOR SPATIALLY CORRELATED DISCRETE DATA JAPANESE BEETLE DATA
JAPANESE BEETLE DATA 6 MAXIMUM LIKELIHOOD FOR SPATIALLY CORRELATED DISCRETE DATA Gauge Plots TuscaroraLisa Central Madsen Fairways, 996 January 9, 7 Grubs Adult Activity Grub Counts 6 8 Organic Matter
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 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 informationConditional 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 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 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 informationIntroduction to Estimation Methods for Time Series models Lecture 2
Introduction to Estimation Methods for Time Series models Lecture 2 Fulvio Corsi SNS Pisa Fulvio Corsi Introduction to Estimation () Methods for Time Series models Lecture 2 SNS Pisa 1 / 21 Estimators:
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 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 informationMaximum Likelihood, Logistic Regression, and Stochastic Gradient Training
Maximum Likelihood, Logistic Regression, and Stochastic Gradient Training Charles Elkan elkan@cs.ucsd.edu January 17, 2013 1 Principle of maximum likelihood Consider a family of probability distributions
More informationMAS223 Statistical Inference and Modelling Exercises
MAS223 Statistical Inference and Modelling Exercises The exercises are grouped into sections, corresponding to chapters of the lecture notes Within each section exercises are divided into warm-up questions,
More informationOn the Estimation and Application of Max-Stable Processes
On the Estimation and Application of Max-Stable Processes Zhengjun Zhang Department of Statistics University of Wisconsin Madison, WI 53706, USA Co-author: Richard Smith EVA 2009, Fort Collins, CO Z. Zhang
More informationFall 2017 STAT 532 Homework Peter Hoff. 1. Let P be a probability measure on a collection of sets A.
1. Let P be a probability measure on a collection of sets A. (a) For each n N, let H n be a set in A such that H n H n+1. Show that P (H n ) monotonically converges to P ( k=1 H k) as n. (b) For each n
More informationInvestigation of an Automated Approach to Threshold Selection for Generalized Pareto
Investigation of an Automated Approach to Threshold Selection for Generalized Pareto Kate R. Saunders Supervisors: Peter Taylor & David Karoly University of Melbourne April 8, 2015 Outline 1 Extreme Value
More informationAssociation studies and regression
Association studies and regression CM226: Machine Learning for Bioinformatics. Fall 2016 Sriram Sankararaman Acknowledgments: Fei Sha, Ameet Talwalkar Association studies and regression 1 / 104 Administration
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 informationAccurate directional inference for vector parameters
Accurate directional inference for vector parameters Nancy Reid February 26, 2016 with Don Fraser, Nicola Sartori, Anthony Davison Nancy Reid Accurate directional inference for vector parameters York University
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 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 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 informationMA/ST 810 Mathematical-Statistical Modeling and Analysis of Complex Systems
MA/ST 810 Mathematical-Statistical Modeling and Analysis of Complex Systems Principles of Statistical Inference Recap of statistical models Statistical inference (frequentist) Parametric vs. semiparametric
More informationRecap. Vector observation: Y f (y; θ), Y Y R m, θ R d. sample of independent vectors y 1,..., y n. pairwise log-likelihood n m. weights are often 1
Recap Vector observation: Y f (y; θ), Y Y R m, θ R d sample of independent vectors y 1,..., y n pairwise log-likelihood n m i=1 r=1 s>r w rs log f 2 (y ir, y is ; θ) weights are often 1 more generally,
More informationAccurate directional inference for vector parameters
Accurate directional inference for vector parameters Nancy Reid October 28, 2016 with Don Fraser, Nicola Sartori, Anthony Davison Parametric models and likelihood model f (y; θ), θ R p data y = (y 1,...,
More informationSTATISTICAL MODELS FOR QUANTIFYING THE SPATIAL DISTRIBUTION OF SEASONALLY DERIVED OZONE STANDARDS
STATISTICAL MODELS FOR QUANTIFYING THE SPATIAL DISTRIBUTION OF SEASONALLY DERIVED OZONE STANDARDS Eric Gilleland Douglas Nychka Geophysical Statistics Project National Center for Atmospheric Research Supported
More informationExtreme Precipitation: An Application Modeling N-Year Return Levels at the Station Level
Extreme Precipitation: An Application Modeling N-Year Return Levels at the Station Level Presented by: Elizabeth Shamseldin Joint work with: Richard Smith, Doug Nychka, Steve Sain, Dan Cooley Statistics
More informationSpatial Statistics with Image Analysis. Lecture L02. Computer exercise 0 Daily Temperature. Lecture 2. Johan Lindström.
C Stochastic fields Covariance Spatial Statistics with Image Analysis Lecture 2 Johan Lindström November 4, 26 Lecture L2 Johan Lindström - johanl@maths.lth.se FMSN2/MASM2 L /2 C Stochastic fields Covariance
More informationCopula Regression RAHUL A. PARSA DRAKE UNIVERSITY & STUART A. KLUGMAN SOCIETY OF ACTUARIES CASUALTY ACTUARIAL SOCIETY MAY 18,2011
Copula Regression RAHUL A. PARSA DRAKE UNIVERSITY & STUART A. KLUGMAN SOCIETY OF ACTUARIES CASUALTY ACTUARIAL SOCIETY MAY 18,2011 Outline Ordinary Least Squares (OLS) Regression Generalized Linear Models
More informationComposite Likelihood Estimation
Composite Likelihood Estimation With application to spatial clustered data Cristiano Varin Wirtschaftsuniversität Wien April 29, 2016 Credits CV, Nancy M Reid and David Firth (2011). An overview of composite
More informationSYSM 6303: Quantitative Introduction to Risk and Uncertainty in Business Lecture 4: Fitting Data to Distributions
SYSM 6303: Quantitative Introduction to Risk and Uncertainty in Business Lecture 4: Fitting Data to Distributions M. Vidyasagar Cecil & Ida Green Chair The University of Texas at Dallas Email: M.Vidyasagar@utdallas.edu
More informationA Conditional Approach to Modeling Multivariate Extremes
A Approach to ing Multivariate Extremes By Heffernan & Tawn Department of Statistics Purdue University s April 30, 2014 Outline s s Multivariate Extremes s A central aim of multivariate extremes is trying
More informationModel Specification Testing in Nonparametric and Semiparametric Time Series Econometrics. Jiti Gao
Model Specification Testing in Nonparametric and Semiparametric Time Series Econometrics Jiti Gao Department of Statistics School of Mathematics and Statistics The University of Western Australia Crawley
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 informationMultivariate Statistics
Multivariate Statistics Chapter 2: Multivariate distributions and inference Pedro Galeano Departamento de Estadística Universidad Carlos III de Madrid pedro.galeano@uc3m.es Course 2016/2017 Master in Mathematical
More informationAn Extended BIC for Model Selection
An Extended BIC for Model Selection at the JSM meeting 2007 - Salt Lake City Surajit Ray Boston University (Dept of Mathematics and Statistics) Joint work with James Berger, Duke University; Susie Bayarri,
More informationComposite likelihood methods
1 / 20 Composite likelihood methods Nancy Reid University of Warwick, April 15, 2008 Cristiano Varin Grace Yun Yi, Zi Jin, Jean-François Plante 2 / 20 Some questions (and answers) Is there a drinks table
More informationGeneralized additive modelling of hydrological sample extremes
Generalized additive modelling of hydrological sample extremes Valérie Chavez-Demoulin 1 Joint work with A.C. Davison (EPFL) and Marius Hofert (ETHZ) 1 Faculty of Business and Economics, University of
More 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 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 informationStatistícal Methods for Spatial Data Analysis
Texts in Statistícal Science Statistícal Methods for Spatial Data Analysis V- Oliver Schabenberger Carol A. Gotway PCT CHAPMAN & K Contents Preface xv 1 Introduction 1 1.1 The Need for Spatial Analysis
More informationFractional Hot Deck Imputation for Robust Inference Under Item Nonresponse in Survey Sampling
Fractional Hot Deck Imputation for Robust Inference Under Item Nonresponse in Survey Sampling Jae-Kwang Kim 1 Iowa State University June 26, 2013 1 Joint work with Shu Yang Introduction 1 Introduction
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 informationAsymptotic Multivariate Kriging Using Estimated Parameters with Bayesian Prediction Methods for Non-linear Predictands
Asymptotic Multivariate Kriging Using Estimated Parameters with Bayesian Prediction Methods for Non-linear Predictands Elizabeth C. Mannshardt-Shamseldin Advisor: Richard L. Smith Duke University Department
More informationTopics on estimation, prediction and bounding risk for multivariate extremes
Topics on estimation, prediction and bounding risk for multivariate extremes by Robert Alohimakalani Yuen A dissertation submitted in partial fulfillment of the requirements for the degree of Doctor of
More informationLikelihood-Based Methods
Likelihood-Based Methods Handbook of Spatial Statistics, Chapter 4 Susheela Singh September 22, 2016 OVERVIEW INTRODUCTION MAXIMUM LIKELIHOOD ESTIMATION (ML) RESTRICTED MAXIMUM LIKELIHOOD ESTIMATION (REML)
More informationLecture 8: Information Theory and Statistics
Lecture 8: Information Theory and Statistics Part II: Hypothesis Testing and I-Hsiang Wang Department of Electrical Engineering National Taiwan University ihwang@ntu.edu.tw December 23, 2015 1 / 50 I-Hsiang
More informationESTIMATING THE MEAN LEVEL OF FINE PARTICULATE MATTER: AN APPLICATION OF SPATIAL STATISTICS
ESTIMATING THE MEAN LEVEL OF FINE PARTICULATE MATTER: AN APPLICATION OF SPATIAL STATISTICS Richard L. Smith Department of Statistics and Operations Research University of North Carolina Chapel Hill, N.C.,
More informationZwiers FW and Kharin VV Changes in the extremes of the climate simulated by CCC GCM2 under CO 2 doubling. J. Climate 11:
Statistical Analysis of EXTREMES in GEOPHYSICS Zwiers FW and Kharin VV. 1998. Changes in the extremes of the climate simulated by CCC GCM2 under CO 2 doubling. J. Climate 11:2200 2222. http://www.ral.ucar.edu/staff/ericg/readinggroup.html
More informationSimulating Realistic Ecological Count Data
1 / 76 Simulating Realistic Ecological Count Data Lisa Madsen Dave Birkes Oregon State University Statistics Department Seminar May 2, 2011 2 / 76 Outline 1 Motivation Example: Weed Counts 2 Pearson Correlation
More informationCHAPTER 1 INTRODUCTION
CHAPTER 1 INTRODUCTION 1.0 Discrete distributions in statistical analysis Discrete models play an extremely important role in probability theory and statistics for modeling count data. The use of discrete
More informationBayesian Model Averaging for Multivariate Extreme Values
Bayesian Model Averaging for Multivariate Extreme Values Philippe Naveau naveau@lsce.ipsl.fr Laboratoire des Sciences du Climat et l Environnement (LSCE) Gif-sur-Yvette, France joint work with A. Sabourin
More informationKriging models with Gaussian processes - covariance function estimation and impact of spatial sampling
Kriging models with Gaussian processes - covariance function estimation and impact of spatial sampling François Bachoc former PhD advisor: Josselin Garnier former CEA advisor: Jean-Marc Martinez Department
More informationComposite Likelihood
Composite Likelihood Nancy Reid January 30, 2012 with Cristiano Varin and thanks to Don Fraser, Grace Yi, Ximing Xu Background parametric model f (y; θ), y R m ; θ R d likelihood function L(θ; y) f (y;
More informationOverview of Extreme Value Analysis (EVA)
Overview of Extreme Value Analysis (EVA) Brian Reich North Carolina State University July 26, 2016 Rossbypalooza Chicago, IL Brian Reich Overview of Extreme Value Analysis (EVA) 1 / 24 Importance of extremes
More informationTesting Statistical Hypotheses
E.L. Lehmann Joseph P. Romano Testing Statistical Hypotheses Third Edition 4y Springer Preface vii I Small-Sample Theory 1 1 The General Decision Problem 3 1.1 Statistical Inference and Statistical Decisions
More informationA Very Brief Summary of Statistical Inference, and Examples
A Very Brief Summary of Statistical Inference, and Examples Trinity Term 2009 Prof. Gesine Reinert Our standard situation is that we have data x = x 1, x 2,..., x n, which we view as realisations of random
More informationHIERARCHICAL MODELS IN EXTREME VALUE THEORY
HIERARCHICAL MODELS IN EXTREME VALUE THEORY Richard L. Smith Department of Statistics and Operations Research, University of North Carolina, Chapel Hill and Statistical and Applied Mathematical Sciences
More informationSpatial statistics, addition to Part I. Parameter estimation and kriging for Gaussian random fields
Spatial statistics, addition to Part I. Parameter estimation and kriging for Gaussian random fields 1 Introduction Jo Eidsvik Department of Mathematical Sciences, NTNU, Norway. (joeid@math.ntnu.no) February
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 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 informationParametric Modelling of Over-dispersed Count Data. Part III / MMath (Applied Statistics) 1
Parametric Modelling of Over-dispersed Count Data Part III / MMath (Applied Statistics) 1 Introduction Poisson regression is the de facto approach for handling count data What happens then when Poisson
More informationOn the spatial dependence of extreme ocean storm seas
Accepted by Ocean Engineering, August 2017 On the spatial dependence of extreme ocean storm seas Emma Ross a, Monika Kereszturi b, Mirrelijn van Nee c, David Randell d, Philip Jonathan a, a Shell Projects
More informationStatistical Methods in HYDROLOGY CHARLES T. HAAN. The Iowa State University Press / Ames
Statistical Methods in HYDROLOGY CHARLES T. HAAN The Iowa State University Press / Ames Univariate BASIC Table of Contents PREFACE xiii ACKNOWLEDGEMENTS xv 1 INTRODUCTION 1 2 PROBABILITY AND PROBABILITY
More informationStat 5102 Final Exam May 14, 2015
Stat 5102 Final Exam May 14, 2015 Name Student ID The exam is closed book and closed notes. You may use three 8 1 11 2 sheets of paper with formulas, etc. You may also use the handouts on brand name distributions
More informationNotes on the Multivariate Normal and Related Topics
Version: July 10, 2013 Notes on the Multivariate Normal and Related Topics Let me refresh your memory about the distinctions between population and sample; parameters and statistics; population distributions
More informationi=1 h n (ˆθ n ) = 0. (2)
Stat 8112 Lecture Notes Unbiased Estimating Equations Charles J. Geyer April 29, 2012 1 Introduction In this handout we generalize the notion of maximum likelihood estimation to solution of unbiased estimating
More informationGMM - Generalized method of moments
GMM - Generalized method of moments GMM Intuition: Matching moments You want to estimate properties of a data set {x t } T t=1. You assume that x t has a constant mean and variance. x t (µ 0, σ 2 ) Consider
More informationAkaike Information Criterion
Akaike Information Criterion Shuhua Hu Center for Research in Scientific Computation North Carolina State University Raleigh, NC February 7, 2012-1- background Background Model statistical model: Y j =
More informationStatistical Methods for Handling Incomplete Data Chapter 2: Likelihood-based approach
Statistical Methods for Handling Incomplete Data Chapter 2: Likelihood-based approach Jae-Kwang Kim Department of Statistics, Iowa State University Outline 1 Introduction 2 Observed likelihood 3 Mean Score
More informationGeneralized Linear Models. Kurt Hornik
Generalized Linear Models Kurt Hornik Motivation Assuming normality, the linear model y = Xβ + e has y = β + ε, ε N(0, σ 2 ) such that y N(μ, σ 2 ), E(y ) = μ = β. Various generalizations, including general
More informationPart IB Statistics. Theorems with proof. Based on lectures by D. Spiegelhalter Notes taken by Dexter Chua. Lent 2015
Part IB Statistics Theorems with proof Based on lectures by D. Spiegelhalter Notes taken by Dexter Chua Lent 2015 These notes are not endorsed by the lecturers, and I have modified them (often significantly)
More informationLecture 3 September 1
STAT 383C: Statistical Modeling I Fall 2016 Lecture 3 September 1 Lecturer: Purnamrita Sarkar Scribe: Giorgio Paulon, Carlos Zanini Disclaimer: These scribe notes have been slightly proofread and may have
More informationLimit Kriging. Abstract
Limit Kriging V. Roshan Joseph School of Industrial and Systems Engineering Georgia Institute of Technology Atlanta, GA 30332-0205, USA roshan@isye.gatech.edu Abstract A new kriging predictor is proposed.
More informationMultivariate Distributions
IEOR E4602: Quantitative Risk Management Spring 2016 c 2016 by Martin Haugh Multivariate Distributions We will study multivariate distributions in these notes, focusing 1 in particular on multivariate
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 information