A full scale, non stationary approach for the kriging of large spatio(-temporal) datasets
|
|
- Brianna Fowler
- 6 years ago
- Views:
Transcription
1 A full scale, non stationary approach for the kriging of large spatio(-temporal) datasets Thomas Romary, Nicolas Desassis & Francky Fouedjio Mines ParisTech Centre de Géosciences, Equipe Géostatistique Réseau RESSTE 11 mai 2015 T. Romary et al. (Mines ParisTech) RESSTE 1 / 25
2 Introduction Introduction Z(x, t) = µ(x, t) + S(x, t) second order spatio-(temporal) random field Inference Possible non-stationarity in the covariance structure Prediction complexity O(n 3 ), storage O(n 2 ) T. Romary et al. (Mines ParisTech) RESSTE 2 / 25
3 Introduction Notations n size of the dataset z = {z(x 1 ),..., z(x n )} : data C covariance fonction C covariance matrix associated to the data f deterministic covariates T. Romary et al. (Mines ParisTech) RESSTE 3 / 25
4 Introduction Universal kriging system Matrix form We have f 00 f 0 =. f 0L, F = Kriging system and variance f f 1L..... f n0... f nl Z = Λ t z, µ = µ 1. µ L ( ) ( ) C F Λ F t 0 µ = ( C0 f 0 ) and σ 2 UK = C(0) Λ t C 0 µ t f 0 T. Romary et al. (Mines ParisTech) RESSTE 4 / 25
5 Introduction Universal kriging system Z = C 0C 1 z + (f 0 C 0C 1 F)(F C 1 F) 1 F C 1 z Idea C = A + P BP t where A is sparse P is a n p matrix with p n Sherman-Woodbury-Morrison formula C 1 = A 1 A 1 P (B 1 + P t A 1 P ) 1 P t A 1 T. Romary et al. (Mines ParisTech) RESSTE 5 / 25
6 Introduction Universal kriging system Z = C 0C 1 z + (f 0 C 0C 1 F)(F C 1 F) 1 F C 1 z Idea C = A + P BP t where A is sparse P is a n p matrix with p n Sherman-Woodbury-Morrison formula C 1 = A 1 A 1 P (B 1 + P t A 1 P ) 1 P t A 1 T. Romary et al. (Mines ParisTech) RESSTE 5 / 25
7 Introduction Spatio-temporal framework Satellite / Remote sensing data {Z(x, t) : x X R d, t = 1,..., T }, d = 2, 3 Spatio-temporal mixed effect model (Katzfuss and Cressie [2011]) where µ mean function Z(x, t) = µ(x, t) + S(x, t) + S second order random function p ξ i (t)p i (x) Ξ t = (ξ t 1,..., ξt p) centered, square-integrable random variables P i (x) basis functions i=1 T. Romary et al. (Mines ParisTech) RESSTE 6 / 25
8 Introduction Spatio-temporal framework Temporal Structure µ(x, t) = f t (x) β t Ξ t = H t Ξ t 1 + u t, such that var(ξ t Ξ t 1 ) = B t S(x, t) = σ(x, t)ε xt, where ε is a Gaussian white noise Inference by EM algorithm, see Katzfuss and Cressie [2011] Prediction by kriging, filtering, smoothing T. Romary et al. (Mines ParisTech) RESSTE 7 / 25
9 Introduction Spatial-only framework Assume E(Z(x)) = 0 where Z(x) = S(x) + p ξ i P i (x) + ε(x) i=1 cov(z(x), Z(y)) = C(x, y) = D(x, y) + P (x) t BP (y) + σ 2 xδ xy cov(s(x), S(y)) = D(x, y) = C NS (x, y)c CS (x, y) C NS is a non stationary covariance function (Fouedjio et al. [2014]) C CS is a compactly supported covariance function (Furrer et al. [2006], Gneiting [2002]), with fixed range parameter T. Romary et al. (Mines ParisTech) RESSTE 8 / 25
10 Inference Inference Several steps 1 Adjust D 2 Select p basis function from a dictionary by minimizing 3 Estimate B z P β λ β 1 T. Romary et al. (Mines ParisTech) RESSTE 9 / 25
11 Inference D term Compute and fit a local kernel variogram estimator at a set of knots Interpolate the parameters over the domain to get D(x, y) = C NS (x, y)c CS (x, y) see Fouedjio et al. [2014] for more details T. Romary et al. (Mines ParisTech) RESSTE 10 / 25
12 Inference Selection of basis functions Dictionary Multi resolution basis made of columns of compactly supported covariance matrices (including anisotropic ones) Selection by LASSO (glmnet) λ is set by cross-validation z P β λ β 1 T. Romary et al. (Mines ParisTech) RESSTE 11 / 25
13 Inference Estimation of B 1 Compute a non-parametric estimator of the covariance at a set of knots Ĉ K (x, y) = 1 K ij (z(x i ))(z(x j )) Σ Kxy where K ij = K h (x x i, y y j ) and Σ Kxy = i,j K ij 2 Compute P = QR 3 Compute B = R 1 Q T (ĈK A θ)q(r 1 ) T i,j T. Romary et al. (Mines ParisTech) RESSTE 12 / 25
14 Example Example Synthetic data on [0, 1] 2, Gaussian RF, centered with covariance where C(x, y) = 0.3sph( x y, 0.1) + P (x) t BP (y) + 0.1δ xy P 50 basis functions taken at random among 736 B drawn from a Wishart distribution Comparison between full model, small scale only, basis functions only, true model range of the small scale term = 0.05 T. Romary et al. (Mines ParisTech) RESSTE 13 / 25
15 Example Reference realization T. Romary et al. (Mines ParisTech) RESSTE 14 / 25
16 Example Reference realization sample (5000) T. Romary et al. (Mines ParisTech) RESSTE 15 / 25
17 Example Building the dictionary Ingredients models: cubic, Wendland 2 ranges: 0.5, 0.3 angles: 0, π/6, π/3, π/2, 2π/3, 5π/6 736 basis functions in total T. Romary et al. (Mines ParisTech) RESSTE 16 / 25
18 Example Selection T. Romary et al. (Mines ParisTech) RESSTE 17 / 25
19 Example Kriging map: full model T. Romary et al. (Mines ParisTech) RESSTE 18 / 25
20 Example Kriging map: full model T. Romary et al. (Mines ParisTech) RESSTE 19 / 25
21 Example Kriging map: short scale only T. Romary et al. (Mines ParisTech) RESSTE 20 / 25
22 Example Kriging map: basis only T. Romary et al. (Mines ParisTech) RESSTE 21 / 25
23 Example Kriging map: true model T. Romary et al. (Mines ParisTech) RESSTE 22 / 25
24 Example Results Model MSE MSE (true) Full Taper only Basis only (36) True model T. Romary et al. (Mines ParisTech) RESSTE 23 / 25
25 Example Discussion Work in progress Lots of tuning: knots (basis functions, variogram estimation, covariance estimation) parameters of the basis functions bandwidth parameters Prediction variances T. Romary et al. (Mines ParisTech) RESSTE 24 / 25
26 Example References S. Banerjee, A. E. Gelfand, A. O. Finley, and H. Sang. Gaussian predictive process models for large spatial data sets. Journal of the Royal Statistical Society. Series B, 70(4): , N. Cressie and G. Johannesson. Fixed rank kriging for very large spatial data sets. Journal of the Royal Statistical Society. Series B, 70: , Francky Fouedjio, Nicolas Desassis, and Jacques Rivoirard. A generalized convolution model and estimation for non-stationary random functions. arxiv preprint arxiv: , R. Furrer, M. G. Genton, and D. W. Nychka. Covariance tapering for interpolation of large spatial datasets. Journal of Computational and Graphical Statistics, 15(3): , T. Gneiting. Compactly supported correlation functions. Journal of Multivariate Analysis, 83(2): , Matthias Katzfuss and Noel Cressie. Spatio-temporal smoothing and em estimation for massive remote-sensing data sets. Journal of Time Series Analysis, 32(4): , H. Sang and J.Z. Huang. A full scale approximation of covariance functions for large spatial data sets. Journal of the Royal Statistical Society. Series B, 74(1): , M. L. Stein. A modeling approach for large spatial datasets. Journal of the Korean Statistical Society, 37:3 10, T. Romary et al. (Mines ParisTech) RESSTE 25 / 25
27 Sparse approximation: Covariance Tapering Principle C tap (x, y) = C(x, y)c θ (x, y) where C θ (x, y) is a compactly supported kernel = C tap (x, y) generates a sparse matrix Complexity O(n), storage O(n) R. Furrer, M. G. Genton, and D. W. Nychka. Covariance tapering for interpolation of large spatial datasets. Journal of Computational and Graphical Statistics, 15(3): , 2006 T. Romary et al. (Mines ParisTech) RESSTE 26 / 25
28 Low rank approximation Principle C LR (x, y) = σ 2 δ xy + P (x) B P (y) where σ 2 is the variance of the nugget, P is a vector of p basis functions, B R p p with p n Complexity O(n p 2 ), storage O(n p) N. Cressie and G. Johannesson. Fixed rank kriging for very large spatial data sets. Journal of the Royal Statistical Society. Series B, 70: , 2008 S. Banerjee, A. E. Gelfand, A. O. Finley, and H. Sang. Gaussian predictive process models for large spatial data sets. Journal of the Royal Statistical Society. Series B, 70(4): , 2008 T. Romary et al. (Mines ParisTech) RESSTE 27 / 25
29 Combination Principle C(x, y) C tap (x, y) + P (x) B P (y) Inversion by the Sherman-Morrison-Woodbury formula Complexity O(n p 2 ), storage O(n p) M. L. Stein. A modeling approach for large spatial datasets. Journal of the Korean Statistical Society, 37:3 10, 2008 H. Sang and J.Z. Huang. A full scale approximation of covariance functions for large spatial data sets. Journal of the Royal Statistical Society. Series B, 74(1): , 2012 T. Romary et al. (Mines ParisTech) RESSTE 28 / 25
Geostatistical Modeling for Large Data Sets: Low-rank methods
Geostatistical Modeling for Large Data Sets: Low-rank methods Whitney Huang, Kelly-Ann Dixon Hamil, and Zizhuang Wu Department of Statistics Purdue University February 22, 2016 Outline Motivation Low-rank
More informationNonparametric Bayesian Methods
Nonparametric Bayesian Methods Debdeep Pati Florida State University October 2, 2014 Large spatial datasets (Problem of big n) Large observational and computer-generated datasets: Often have spatial and
More informationBayesian Modeling and Inference for High-Dimensional Spatiotemporal Datasets
Bayesian Modeling and Inference for High-Dimensional Spatiotemporal Datasets Sudipto Banerjee University of California, Los Angeles, USA Based upon projects involving: Abhirup Datta (Johns Hopkins University)
More informationOn Gaussian Process Models for High-Dimensional Geostatistical Datasets
On Gaussian Process Models for High-Dimensional Geostatistical Datasets Sudipto Banerjee Joint work with Abhirup Datta, Andrew O. Finley and Alan E. Gelfand University of California, Los Angeles, USA May
More informationParameter Estimation in the Spatio-Temporal Mixed Effects Model Analysis of Massive Spatio-Temporal Data Sets
Parameter Estimation in the Spatio-Temporal Mixed Effects Model Analysis of Massive Spatio-Temporal Data Sets Matthias Katzfuß Advisor: Dr. Noel Cressie Department of Statistics The Ohio State University
More informationExperimental variogram of the residuals in the universal kriging (UK) model
Experimental variogram of the residuals in the universal kriging (UK) model Nicolas DESASSIS Didier RENARD Technical report R141210NDES MINES ParisTech Centre de Géosciences Equipe Géostatistique 35, rue
More informationGeostatistics for Gaussian processes
Introduction Geostatistical Model Covariance structure Cokriging Conclusion Geostatistics for Gaussian processes Hans Wackernagel Geostatistics group MINES ParisTech http://hans.wackernagel.free.fr Kernels
More informationAn Additive Gaussian Process Approximation for Large Spatio-Temporal Data
An Additive Gaussian Process Approximation for Large Spatio-Temporal Data arxiv:1801.00319v2 [stat.me] 31 Oct 2018 Pulong Ma Statistical and Applied Mathematical Sciences Institute and Duke University
More informationBayesian dynamic modeling for large space-time weather datasets using Gaussian predictive processes
Bayesian dynamic modeling for large space-time weather datasets using Gaussian predictive processes Alan Gelfand 1 and Andrew O. Finley 2 1 Department of Statistical Science, Duke University, Durham, North
More informationarxiv: v1 [stat.me] 3 Dec 2014
A Generalized Convolution Model and Estimation for Non-stationary Random Functions F. Fouedjio, N. Desassis, J. Rivoirard Équipe Géostatistique - Centre de Géosciences, MINES ParisTech 35, rue Saint-Honoré,
More informationBayesian dynamic modeling for large space-time weather datasets using Gaussian predictive processes
Bayesian dynamic modeling for large space-time weather datasets using Gaussian predictive processes Andrew O. Finley 1 and Sudipto Banerjee 2 1 Department of Forestry & Department of Geography, Michigan
More informationBayesian dynamic modeling for large space-time weather datasets using Gaussian predictive processes
Bayesian dynamic modeling for large space-time weather datasets using Gaussian predictive processes Sudipto Banerjee 1 and Andrew O. Finley 2 1 Biostatistics, School of Public Health, University of Minnesota,
More informationA Fused Lasso Approach to Nonstationary Spatial Covariance Estimation
Supplementary materials for this article are available at 10.1007/s13253-016-0251-8. A Fused Lasso Approach to Nonstationary Spatial Covariance Estimation Ryan J. Parker, Brian J. Reich,andJoEidsvik Spatial
More informationKarhunen-Loeve Expansion and Optimal Low-Rank Model for Spatial Processes
TTU, October 26, 2012 p. 1/3 Karhunen-Loeve Expansion and Optimal Low-Rank Model for Spatial Processes Hao Zhang Department of Statistics Department of Forestry and Natural Resources Purdue University
More informationA Note on the comparison of Nearest Neighbor Gaussian Process (NNGP) based models
A Note on the comparison of Nearest Neighbor Gaussian Process (NNGP) based models arxiv:1811.03735v1 [math.st] 9 Nov 2018 Lu Zhang UCLA Department of Biostatistics Lu.Zhang@ucla.edu Sudipto Banerjee UCLA
More informationBayesian dynamic modeling for large space-time weather datasets using Gaussian predictive processes
Bayesian dynamic modeling for large space-time weather datasets using Gaussian predictive processes Andrew O. Finley Department of Forestry & Department of Geography, Michigan State University, Lansing
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 informationLow-rank methods and predictive processes for spatial models
Low-rank methods and predictive processes for spatial models Sam Bussman, Linchao Chen, John Lewis, Mark Risser with Sebastian Kurtek, Vince Vu, Ying Sun February 27, 2014 Outline Introduction and general
More informationMultivariate modelling and efficient estimation of Gaussian random fields with application to roller data
Multivariate modelling and efficient estimation of Gaussian random fields with application to roller data Reinhard Furrer, UZH PASI, Búzios, 14-06-25 NZZ.ch Motivation Microarray data: construct alternative
More informationA Framework for Daily Spatio-Temporal Stochastic Weather Simulation
A Framework for Daily Spatio-Temporal Stochastic Weather Simulation, Rick Katz, Balaji Rajagopalan Geophysical Statistics Project Institute for Mathematics Applied to Geosciences National Center for Atmospheric
More informationA Multi-resolution Gaussian process model for the analysis of large spatial data sets.
1 2 3 4 A Multi-resolution Gaussian process model for the analysis of large spatial data sets. Douglas Nychka, Soutir Bandyopadhyay, Dorit Hammerling, Finn Lindgren, and Stephan Sain August 13, 2013 5
More informationA multi-resolution Gaussian process model for the analysis of large spatial data sets.
National Science Foundation A multi-resolution Gaussian process model for the analysis of large spatial data sets. Doug Nychka Soutir Bandyopadhyay Dorit Hammerling Finn Lindgren Stephen Sain NCAR/TN-504+STR
More informationThe Matrix Reloaded: Computations for large spatial data sets
The Matrix Reloaded: Computations for large spatial data sets The spatial model Solving linear systems Matrix multiplication Creating sparsity Doug Nychka National Center for Atmospheric Research Sparsity,
More informationHierarchical Nearest-Neighbor Gaussian Process Models for Large Geo-statistical Datasets
Hierarchical Nearest-Neighbor Gaussian Process Models for Large Geo-statistical Datasets Abhirup Datta 1 Sudipto Banerjee 1 Andrew O. Finley 2 Alan E. Gelfand 3 1 University of Minnesota, Minneapolis,
More informationGaussian predictive process models for large spatial data sets.
Gaussian predictive process models for large spatial data sets. Sudipto Banerjee, Alan E. Gelfand, Andrew O. Finley, and Huiyan Sang Presenters: Halley Brantley and Chris Krut September 28, 2015 Overview
More informationDouglas Nychka, Soutir Bandyopadhyay, Dorit Hammerling, Finn Lindgren, and Stephan Sain. October 10, 2012
A multi-resolution Gaussian process model for the analysis of large spatial data sets. Douglas Nychka, Soutir Bandyopadhyay, Dorit Hammerling, Finn Lindgren, and Stephan Sain October 10, 2012 Abstract
More informationAn Introduction to Spatial Statistics. Chunfeng Huang Department of Statistics, Indiana University
An Introduction to Spatial Statistics Chunfeng Huang Department of Statistics, Indiana University Microwave Sounding Unit (MSU) Anomalies (Monthly): 1979-2006. Iron Ore (Cressie, 1986) Raw percent data
More informationThe Matrix Reloaded: Computations for large spatial data sets
The Matrix Reloaded: Computations for large spatial data sets Doug Nychka National Center for Atmospheric Research The spatial model Solving linear systems Matrix multiplication Creating sparsity Sparsity,
More informationarxiv: v4 [stat.me] 14 Sep 2015
Does non-stationary spatial data always require non-stationary random fields? Geir-Arne Fuglstad 1, Daniel Simpson 1, Finn Lindgren 2, and Håvard Rue 1 1 Department of Mathematical Sciences, NTNU, Norway
More informationA Bayesian Spatio-Temporal Geostatistical Model with an Auxiliary Lattice for Large Datasets
Statistica Sinica (2013): Preprint 1 A Bayesian Spatio-Temporal Geostatistical Model with an Auxiliary Lattice for Large Datasets Ganggang Xu 1, Faming Liang 1 and Marc G. Genton 2 1 Texas A&M University
More informationLatent Gaussian Processes and Stochastic Partial Differential Equations
Latent Gaussian Processes and Stochastic Partial Differential Equations Johan Lindström 1 1 Centre for Mathematical Sciences Lund University Lund 2016-02-04 Johan Lindström - johanl@maths.lth.se Gaussian
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 informationBasics in Geostatistics 2 Geostatistical interpolation/estimation: Kriging methods. Hans Wackernagel. MINES ParisTech.
Basics in Geostatistics 2 Geostatistical interpolation/estimation: Kriging methods Hans Wackernagel MINES ParisTech NERSC April 2013 http://hans.wackernagel.free.fr Basic concepts Geostatistics Hans Wackernagel
More informationNouveaux développements en géostatistique minière. New developments in mining geostatistics
Ressources Minérales : la vision du Mineur 1 1-3 Février 2012 Nouveaux développements en géostatistique minière New developments in mining geostatistics Jacques Rivoirard Responsable de l Équipe Géostatistique
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 informationPositive Definite Functions on Spheres
Positive Definite Functions on Spheres Tilmann Gneiting Institute for Applied Mathematics Heidelberg University, Germany t.gneiting@uni-heidelberg.de www.math.uni-heidelberg.de/spatial/tilmann/ International
More informationJournal of Multivariate Analysis. Nonstationary modeling for multivariate spatial processes
Journal of Multivariate Analysis () 76 9 Contents lists available at SciVerse ScienceDirect Journal of Multivariate Analysis journal homepage: wwwelseviercom/locate/jmva Nonstationary modeling for multivariate
More informationHierarchical Low Rank Approximation of Likelihoods for Large Spatial Datasets
Hierarchical Low Rank Approximation of Likelihoods for Large Spatial Datasets Huang Huang and Ying Sun CEMSE Division, King Abdullah University of Science and Technology, Thuwal, 23955, Saudi Arabia. July
More informationHandbook of Spatial Statistics Chapter 2: Continuous Parameter Stochastic Process Theory by Gneiting and Guttorp
Handbook of Spatial Statistics Chapter 2: Continuous Parameter Stochastic Process Theory by Gneiting and Guttorp Marcela Alfaro Córdoba August 25, 2016 NCSU Department of Statistics Continuous Parameter
More informationProceedings of the 2015 Winter Simulation Conference L. Yilmaz, W. K. V. Chan, I. Moon, T. M. K. Roeder, C. Macal, and M. D. Rossetti, eds.
Proceedings of the 2015 Winter Simulation Conference L. Yilmaz, W. K. V. Chan, I. Moon, T. M. K. Roeder, C. Macal, and M. D. Rossetti, eds. AN ADDITIVE GLOBAL AND LOCAL GAUSSIAN PROCESS MODEL FOR LARGE
More informationSummary STK 4150/9150
STK4150 - Intro 1 Summary STK 4150/9150 Odd Kolbjørnsen May 22 2017 Scope You are expected to know and be able to use basic concepts introduced in the book. You knowledge is expected to be larger than
More informationSpatial Backfitting of Roller Measurement Values from a Florida Test Bed
Spatial Backfitting of Roller Measurement Values from a Florida Test Bed Daniel K. Heersink 1, Reinhard Furrer 1, and Mike A. Mooney 2 1 Institute of Mathematics, University of Zurich, CH-8057 Zurich 2
More informationModelling Multivariate Spatial Data
Modelling Multivariate Spatial Data Sudipto Banerjee 1 Biostatistics, School of Public Health, University of Minnesota, Minneapolis, Minnesota, U.S.A. June 20th, 2014 1 Point-referenced spatial data often
More informationA full-scale approximation of covariance functions for large spatial data sets
A full-scale approximation of covariance functions for large spatial data sets Huiyan Sang Department of Statistics, Texas A&M University, College Station, USA. Jianhua Z. Huang Department of Statistics,
More informationA Divide-and-Conquer Bayesian Approach to Large-Scale Kriging
A Divide-and-Conquer Bayesian Approach to Large-Scale Kriging Cheng Li DSAP, National University of Singapore Joint work with Rajarshi Guhaniyogi (UC Santa Cruz), Terrance D. Savitsky (US Bureau of Labor
More informationIsotropic Covariance Functions on Spheres: Some Properties and Modeling Considerations
Isotropic Covariance Functions on Spheres: Some Properties and Modeling Considerations Joseph Guinness a, Montserrat Fuentes a a NC State University, Department of Statistics Abstract Introducing flexible
More informationCOVARIANCE APPROXIMATION FOR LARGE MULTIVARIATE SPATIAL DATA SETS WITH AN APPLICATION TO MULTIPLE CLIMATE MODEL ERRORS 1
The Annals of Applied Statistics 2011, Vol. 5, No. 4, 2519 2548 DOI: 10.1214/11-AOAS478 Institute of Mathematical Statistics, 2011 COVARIANCE APPROXIMATION FOR LARGE MULTIVARIATE SPATIAL DATA SETS WITH
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 informationNearest Neighbor Gaussian Processes for Large Spatial Data
Nearest Neighbor Gaussian Processes for Large Spatial Data Abhi Datta 1, Sudipto Banerjee 2 and Andrew O. Finley 3 July 31, 2017 1 Department of Biostatistics, Bloomberg School of Public Health, Johns
More informationCovariance Tapering for Interpolation of Large Spatial Datasets
Covariance Tapering for Interpolation of Large Spatial Datasets Reinhard Furrer, Marc G. Genton and Douglas Nychka Interpolation of a spatially correlated random process is used in many areas. The best
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 informationENVIRONMENTAL DATA ANALYSIS WILLIAM MENKE JOSHUA MENKE WITH MATLAB COPYRIGHT 2011 BY ELSEVIER, INC. ALL RIGHTS RESERVED.
ENVIRONMENTAL DATA ANALYSIS WITH MATLAB WILLIAM MENKE PROFESSOR OF EARTH AND ENVIRONMENTAL SCIENCE COLUMBIA UNIVERSITY JOSHUA MENKE SOFTWARE ENGINEER JOM ASSOCIATES COPYRIGHT 2011 BY ELSEVIER, INC. ALL
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 informationWrapped Gaussian processes: a short review and some new results
Wrapped Gaussian processes: a short review and some new results Giovanna Jona Lasinio 1, Gianluca Mastrantonio 2 and Alan Gelfand 3 1-Università Sapienza di Roma 2- Università RomaTRE 3- Duke University
More informationStatistical modeling of MODIS cloud data using the spatial random effects model
University of Wollongong Research Online Centre for Statistical & Survey Methodology Working Paper Series Faculty of Engineering and Information Sciences 2013 Statistical modeling of MODIS cloud data using
More informationPredictive spatio-temporal models for spatially sparse environmental data. Umeå University
Seminar p.1/28 Predictive spatio-temporal models for spatially sparse environmental data Xavier de Luna and Marc G. Genton xavier.deluna@stat.umu.se and genton@stat.ncsu.edu http://www.stat.umu.se/egna/xdl/index.html
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 informationarxiv: v1 [stat.me] 24 May 2010
The role of the nugget term in the Gaussian process method Andrey Pepelyshev arxiv:1005.4385v1 [stat.me] 24 May 2010 Abstract The maximum likelihood estimate of the correlation parameter of a Gaussian
More informationKriging and Alternatives in Computer Experiments
Kriging and Alternatives in Computer Experiments C. F. Jeff Wu ISyE, Georgia Institute of Technology Use kriging to build meta models in computer experiments, a brief review Numerical problems with kriging
More informationA BAYESIAN SPATIO-TEMPORAL GEOSTATISTICAL MODEL WITH AN AUXILIARY LATTICE FOR LARGE DATASETS
Statistica Sinica 25 (2015), 61-79 doi:http://dx.doi.org/10.5705/ss.2013.085w A BAYESIAN SPATIO-TEMPORAL GEOSTATISTICAL MODEL WITH AN AUXILIARY LATTICE FOR LARGE DATASETS Ganggang Xu 1, Faming Liang 1
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 informationGaussian processes for spatial modelling in environmental health: parameterizing for flexibility vs. computational efficiency
Gaussian processes for spatial modelling in environmental health: parameterizing for flexibility vs. computational efficiency Chris Paciorek March 11, 2005 Department of Biostatistics Harvard School of
More informationTheory and Computation for Gaussian Processes
University of Chicago IPAM, February 2015 Funders & Collaborators US Department of Energy, US National Science Foundation (STATMOS) Mihai Anitescu, Jie Chen, Ying Sun Gaussian processes A process Z on
More informationNew Flexible Compact Covariance Model on a Sphere
New Flexible Compact Covariance Model on a Sphere Alexander Gribov, Konstantin Krivoruchko Esri 380 New York Street Redlands, CA 9373 Emails: agribov@esri.com; kkrivoruchko@esri.com arxiv:1701.03405v1
More informationSpatial Statistics with Image Analysis. Outline. A Statistical Approach. Johan Lindström 1. Lund October 6, 2016
Spatial Statistics Spatial Examples More Spatial Statistics with Image Analysis Johan Lindström 1 1 Mathematical Statistics Centre for Mathematical Sciences Lund University Lund October 6, 2016 Johan Lindström
More informationStatistical Models for Monitoring and Regulating Ground-level Ozone. Abstract
Statistical Models for Monitoring and Regulating Ground-level Ozone Eric Gilleland 1 and Douglas Nychka 2 Abstract The application of statistical techniques to environmental problems often involves a tradeoff
More informationCovariance Tapering for Interpolation of Large Spatial Datasets
Covariance Tapering for Interpolation of Large Spatial Datasets Reinhard FURRER, Marc G. GENTON, and Douglas NYCHKA Interpolation of a spatially correlated random process is used in many scientific areas.
More informationSeparable approximations of space-time covariance matrices
ENVIRONMETRICS Environmetrics 2007; 18: 681 695 Published online 8 May 2007 in Wiley InterScience (www.interscience.wiley.com).854 Separable approximations of space-time covariance matrices Marc G. Genton
More informationSTATISTICAL INTERPOLATION METHODS RICHARD SMITH AND NOEL CRESSIE Statistical methods of interpolation are all based on assuming that the process
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 STATISTICAL INTERPOLATION METHODS
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 informationX random; interested in impact of X on Y. Time series analogue of regression.
Multiple time series Given: two series Y and X. Relationship between series? Possible approaches: X deterministic: regress Y on X via generalized least squares: arima.mle in SPlus or arima in R. We have
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 information7 Geostatistics. Figure 7.1 Focus of geostatistics
7 Geostatistics 7.1 Introduction Geostatistics is the part of statistics that is concerned with geo-referenced data, i.e. data that are linked to spatial coordinates. To describe the spatial variation
More informationFitting Large-Scale Spatial Models with Applications to Microarray Data Analysis
Fitting Large-Scale Spatial Models with Applications to Microarray Data Analysis Stephan R Sain Department of Mathematics University of Colorado at Denver Denver, Colorado ssain@mathcudenveredu Reinhard
More informationSpatial and Environmental Statistics
Spatial and Environmental Statistics Dale Zimmerman Department of Statistics and Actuarial Science University of Iowa January 17, 2019 Dale Zimmerman (UIOWA) Spatial and Environmental Statistics January
More informationPractical Bayesian Optimization of Machine Learning. Learning Algorithms
Practical Bayesian Optimization of Machine Learning Algorithms CS 294 University of California, Berkeley Tuesday, April 20, 2016 Motivation Machine Learning Algorithms (MLA s) have hyperparameters that
More informationMultivariate Gaussian Random Fields with SPDEs
Multivariate Gaussian Random Fields with SPDEs Xiangping Hu Daniel Simpson, Finn Lindgren and Håvard Rue Department of Mathematics, University of Oslo PASI, 214 Outline The Matérn covariance function and
More informationAn Inversion-Free Estimating Equations Approach for. Gaussian Process Models
An Inversion-Free Estimating Equations Approach for Gaussian Process Models Mihai Anitescu Jie Chen Michael L. Stein November 29, 2015 Abstract One of the scalability bottlenecks for the large-scale usage
More informationHilbert Space Methods for Reduced-Rank Gaussian Process Regression
Hilbert Space Methods for Reduced-Rank Gaussian Process Regression Arno Solin and Simo Särkkä Aalto University, Finland Workshop on Gaussian Process Approximation Copenhagen, Denmark, May 2015 Solin &
More informationDecision-Oriented Environmental Mapping with Radial Basis Function Neural Networks
Decision-Oriented Environmental Mapping with Radial Basis Function Neural Networks V. Demyanov (1), N. Gilardi (2), M. Kanevski (1,2), M. Maignan (3), V. Polishchuk (1) (1) Institute of Nuclear Safety
More informationProbabilistic Graphical Models
Probabilistic Graphical Models Brown University CSCI 2950-P, Spring 2013 Prof. Erik Sudderth Lecture 12: Gaussian Belief Propagation, State Space Models and Kalman Filters Guest Kalman Filter Lecture by
More informationSpace-time data. Simple space-time analyses. PM10 in space. PM10 in time
Space-time data Observations taken over space and over time Z(s, t): indexed by space, s, and time, t Here, consider geostatistical/time data Z(s, t) exists for all locations and all times May consider
More informationENVIRONMENTAL ADAPTIVE SAMPLING FOR MOBILE SENSOR NETWORKS USING GAUSSIAN PROCESSES. Yunfei Xu
ENVIRONMENTAL ADAPTIVE SAMPLING FOR MOBILE SENSOR NETWORKS USING GAUSSIAN PROCESSES By Yunfei Xu A DISSERTATION Submitted to Michigan State University in partial fulfillment of the requirements for the
More informationBayesian Transgaussian Kriging
1 Bayesian Transgaussian Kriging Hannes Müller Institut für Statistik University of Klagenfurt 9020 Austria Keywords: Kriging, Bayesian statistics AMS: 62H11,60G90 Abstract In geostatistics a widely used
More informationImproving Spatial Data Interoperability
Improving Spatial Data Interoperability A Framework for Geostatistical Support-To To-Support Interpolation Michael F. Goodchild, Phaedon C. Kyriakidis, Philipp Schneider, Matt Rice, Qingfeng Guan, Jordan
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 informationChapter 4 - Fundamentals of spatial processes Lecture notes
TK4150 - Intro 1 Chapter 4 - Fundamentals of spatial processes Lecture notes Odd Kolbjørnsen and Geir Storvik January 30, 2017 STK4150 - Intro 2 Spatial processes Typically correlation between nearby sites
More informationSimulation of a Gaussian random vector: A propagative approach to the Gibbs sampler
Simulation of a Gaussian random vector: A propagative approach to the Gibbs sampler C. Lantuéjoul, N. Desassis and F. Fouedjio MinesParisTech Kriging and Gaussian processes for computer experiments CHORUS
More informationBruno Sansó. Department of Applied Mathematics and Statistics University of California Santa Cruz bruno
Bruno Sansó Department of Applied Mathematics and Statistics University of California Santa Cruz http://www.ams.ucsc.edu/ bruno Climate Models Climate Models use the equations of motion to simulate changes
More informationLecture 9: Introduction to Kriging
Lecture 9: Introduction to Kriging Math 586 Beginning remarks Kriging is a commonly used method of interpolation (prediction) for spatial data. The data are a set of observations of some variable(s) of
More informationMassive MIMO: Signal Structure, Efficient Processing, and Open Problems II
Massive MIMO: Signal Structure, Efficient Processing, and Open Problems II Mahdi Barzegar Communications and Information Theory Group (CommIT) Technische Universität Berlin Heisenberg Communications and
More informationGaussian processes. Chuong B. Do (updated by Honglak Lee) November 22, 2008
Gaussian processes Chuong B Do (updated by Honglak Lee) November 22, 2008 Many of the classical machine learning algorithms that we talked about during the first half of this course fit the following pattern:
More informationGeostatistics in Hydrology: Kriging interpolation
Chapter Geostatistics in Hydrology: Kriging interpolation Hydrologic properties, such as rainfall, aquifer characteristics (porosity, hydraulic conductivity, transmissivity, storage coefficient, etc.),
More informationA STATISTICAL TECHNIQUE FOR MODELLING NON-STATIONARY SPATIAL PROCESSES
A STATISTICAL TECHNIQUE FOR MODELLING NON-STATIONARY SPATIAL PROCESSES JOHN STEPHENSON 1, CHRIS HOLMES, KERRY GALLAGHER 1 and ALEXANDRE PINTORE 1 Dept. Earth Science and Engineering, Imperial College,
More informationSome Topics in Convolution-Based Spatial Modeling
Some Topics in Convolution-Based Spatial Modeling Catherine A. Calder 1 and Noel Cressie Department of Statistics The Ohio State University 1958 Neil Ave. Columbus, OH 43210 USA Appears in the Proceedings
More informationGaussian Processes. Le Song. Machine Learning II: Advanced Topics CSE 8803ML, Spring 2012
Gaussian Processes Le Song Machine Learning II: Advanced Topics CSE 8803ML, Spring 01 Pictorial view of embedding distribution Transform the entire distribution to expected features Feature space Feature
More informationModified Linear Projection for Large Spatial Data Sets
Modified Linear Projection for Large Spatial Data Sets Toshihiro Hirano July 15, 2014 arxiv:1402.5847v2 [stat.me] 13 Jul 2014 Abstract Recent developments in engineering techniques for spatial data collection
More informationParallel inference for massive distributed spatial data using low-rank models
Parallel inference for massive distributed spatial data using low-rank models Matthias Katfuss Dorit Hammerling arxiv:1402.1472v4 [stat.co] 5 Feb 2016 Published in Statistics and Computing (DOI: 10.1007/s11222-016-9627-4)
More informationExtremes and Atmospheric Data
Extremes and Atmospheric Data Eric Gilleland Research Applications Laboratory National Center for Atmospheric Research 2007-08 Program on Risk Analysis, Extreme Events and Decision Theory, opening workshop
More informationZudi Lu, Dag Tjøstheim and Qiwei Yao Spatial smoothing, Nugget effect and infill asymptotics
Zudi Lu, Dag Tjøstheim and Qiwei Yao Spatial smoothing, Nugget effect and infill asymptotics Article (Submitted version) (Pre-refereed) Original citation: Lu, Zudi and Tjostheim, Dag and Yao, Qiwei (2008)
More information