Fusing space-time data under measurement error for computer model output
|
|
- Gertrude Hutchinson
- 5 years ago
- Views:
Transcription
1 for computer model output SAMSI joint work with Alan E. Gelfand and David M. Holland
2 Introduction In many environmental disciplines data come from two sources: monitoring networks and numerical models Numerical models are deterministic mathematical models used to predict environmental spatio-temporal processes Describe the underlying physical and chemical processes via partial differential equations Equations solved via numerical methods by discretizing space and time Predictions are given in terms of averages over grid cells
3 Introduction (a) Observed ozone concentration (b) Predicted ozone concentration Sparse locations Missing data Essentially, true value Large spatial domains No missingness Calibration concerns
4 Our goal Fuse the two sources of data Address the following issues Spatial scale of outputs from numerical models Outputs from numerical models are given in terms of predictions over grid cells, but predictions at points are more useful Calibration of numerical models Correct outputs from numerical models Problem: Comparing averages over grid cells with points
5 Downscaler: main idea (a) Observed ozone concentration (b) Predicted ozone concentration
6 Downscaler: main idea (a) Observed ozone concentration (b) Predicted ozone concentration
7 Downscaler: main idea (a) Observed ozone concentration (b) Predicted ozone concentration To each point s in the domain S with observation Y (s) we associate the numerical model output at grid cell B, x(b), where B is such that s B: Y (s) x(b)
8 Downscaler Time t is fixed. Y (s) observation at point s, x(b) numerical model output at grid cell B. For s in B: Y (s) = β 0 (s) + β 1 (s)x(b) + ε(s) with β i (s) = β i + β i (s), i=0,1. β 0 (s) and β 1 (s) correlated mean-zero GP. ε(s) ind N(0,τ 2 )
9 Downscaler Time t is fixed. Y (s) observation at point s, x(b) numerical model output at grid cell B. For s in B: Y (s) = β 0 (s) + β 1 (s)x(b) + ε(s) with β i (s) = β i + β i (s), i=0,1. β 0 (s) and β 1 (s) correlated mean-zero GP. ε(s) ind N(0,τ 2 ) Extension to space-time: For s in B and for each t Y (s,t) = β 0 (s,t)+ β 1 (s,t)x(b,t)+ε(s,t) with β i (s,t) = β i,t + β i (s,t), i=0,1. Temporal dependence in β i,t and β i (s,t): dynamic or independent ε(s,t) ind N(0,τ 2 )
10 Downscaler Driven by true station data rather than uncalibrated model output Computationally feasible also for large spatial domains Allows local calibration of the numerical model output Endows the spatial process Y (s) with a non-stationary covariance structure Straightforward prediction at an unmonitored sites Better predictive performance than other methods (geostatistical and model-based)
11 Extending the downscaler Improve the input to the downscaler provided by the numerical model output accounting for uncertainty in the numerical model output accounting for uncertainty in the association x(b) Y (s) with s in B
12 Extending the downscaler Improve the input to the downscaler provided by the numerical model output accounting for uncertainty in the numerical model output accounting for uncertainty in the association x(b) Y (s) with s in B Will consider two possibilities: a Measurement Error Model (MEM) a Smoother with spatially-varying weights
13 Static setting
14 A MEM for outputs from computer model Downscaler: Y (s) = β 0 (s) + β 1 (s)x(b) + ε(s)
15 A MEM for outputs from computer model Downscaler: Y (s) = β 0 (s) + β 1 (s)x(b) + ε(s) Model the numerical model output as a stochastic process x(b) = Ṽ (B) + η(b) η(b) ind N(0,σ 2 ) with Ṽ (B) = µ + V (B) and V (B) ICAR(ξ2,W ).
16 A MEM for outputs from computer model Downscaler: Y (s) = β 0 (s) + β 1 (s)x(b) + ε(s) Model the numerical model output as a stochastic process x(b) = Ṽ (B) + η(b) η(b) ind N(0,σ 2 ) with Ṽ (B) = µ + V (B) and V (B) ICAR(ξ2,W ). For s B: Y (s) = β 0 (s) + β 1 Ṽ (B) + ε(s) ε(s) ind N(0,τ 2 ) with β 0 (s) = β 0 + β 0 (s) and β 0 (s) mean-zero GP with exponential covariance function.
17 A MEM for outputs from computer model Downscaler: Y (s) = β 0 (s) + β 1 (s)x(b) + ε(s) Model the numerical model output as a stochastic process x(b) = Ṽ (B) + η(b) η(b) ind N(0,σ 2 ) with Ṽ (B) = µ + V (B) and V (B) ICAR(ξ2,W ). For s B: Y (s) = β 0 (s) + β 1 Ṽ (B) + ε(s) ε(s) ind N(0,τ 2 ) with β 0 (s) = β 0 + β 0 (s) and β 0 (s) mean-zero GP with exponential covariance function. Then: X MEM Ṽ downscaler Y
18 A MEM for outputs from computer model (a) Observed ozone concentration (b) Predicted ozone concentration
19 A MEM for outputs from computer model (a) Observed ozone concentration (b) Predicted ozone concentration
20 A MEM for outputs from computer model (a) Observed ozone concentration (b) Predicted ozone concentration
21 A MEM for outputs from computer model (a) Observed ozone concentration (b) Predicted ozone concentration
22 A MEM for outputs from computer model /5 1/5 1/5 1/5 1/ (a) Observed ozone concentration (b) Predicted ozone concentration To each point s in the domain S with observation Y (s) we associate the latent Gaussian Markov Random Field Ṽ (B), where B is such that s B and Ṽ (B) also associated with x(b) Y (s) x(b ),B δb
23 Smoother with spatially-varying weights Downscaler: Y (s) = β 0 (s) + β 1 (s)x(b) + ε(s)
24 Smoother with spatially-varying weights Downscaler: Y (s) = β 0 (s) + β 1 (s)x(b) + ε(s) The numerical model output, {x(b)}, is NOT modeled as a stochastic process. For s B : Y (s) = β 0 (s) + β 1 x(s) + ε(s) ε(s) ind N(0,τ 2 ) with x(s) = g k=1 w k(s)x(b k ), and w k (s) random and spatially-varying.
25 Smoother with spatially-varying weights Downscaler: Y (s) = β 0 (s) + β 1 (s)x(b) + ε(s) The numerical model output, {x(b)}, is NOT modeled as a stochastic process. For s B : Y (s) = β 0 (s) + β 1 x(s) + ε(s) ε(s) ind N(0,τ 2 ) with x(s) = g k=1 w k(s)x(b k ), and w k (s) random and spatially-varying. Let r k, k = 1,...,g be the centroids of the numerical model grid cells B k, then w k (s) := K (s r k;ψ) exp(q(r k )) g l=1 K (s r l;ψ) exp(q(r l ))
26 Smoother with spatially-varying weights Downscaler: Y (s) = β 0 (s) + β 1 (s)x(b) + ε(s) The numerical model output, {x(b)}, is NOT modeled as a stochastic process. For s B : Y (s) = β 0 (s) + β 1 x(s) + ε(s) ε(s) ind N(0,τ 2 ) with x(s) = g k=1 w k(s)x(b k ), and w k (s) random and spatially-varying. Let r k, k = 1,...,g be the centroids of the numerical model grid cells B k, then w k (s) := K (s r k;ψ) exp(q(r k )) g l=1 K (s r l;ψ) exp(q(r l )) Q( ) mean-zero GP with variance σ 2 Q, exponential covariance function. K (s r k ;ψ) = exp( ψ s r k ).
27 Smoother with spatially-varying weights (a) Observed ozone concentration (b) Predicted ozone concentration
28 Smoother with spatially-varying weights (a) Observed ozone concentration (b) Predicted ozone concentration
29 Smoother with spatially-varying weights (a) Observed ozone concentration (b) Predicted ozone concentration
30 w_6676 w_6675 w_6674 w_6673 w_6672 w_6671 w_6670 w_6669 w_6688 w_6687 w_6686 w_6685 w_6684 w_6683 w_6682 w_6681 w_6680 w_6679 w_6678 w_6677 w_6864 w_6863 w_6862 w_6861 w_6860 w_6859 w_6858 w_6857 w_6876 w_6875 w_6874 w_6873 w_6872 w_6871 w_6870 w_6869 w_6868 w_6867 w_6866 w_6865 w_7053 w_7052 w_7051 w_7050 w_7049 w_7048 w_7047 w_7046 w_7045 w_7064 w_7063 w_7062 w_7061 w_7060 w_7059 w_7058 w_7057 w_7056 w_7055 w_7054 w_7242 w_7241 w_7240 w_7239 w_7238 w_7237 w_7236 w_7235 w_7234 w_7233 w_7252 w_7251 w_7250 w_7249 w_7248 w_7247 w_7246 w_7245 w_7244 w_7243 w_7430 w_7429 w_7428 w_7427 w_7426 w_7425 w_7424 w_7423 w_7422 w_7421 w_7440 w_7439 w_7438 w_7437 w_7436 w_7435 w_7434 w_7433 w_7432 w_7431 w_7609 w_7619 w_7618 w_7617 w_7616 w_7615 w_7614 w_7613 w_7612 w_7611 w_7610 w_7628 w_7627 w_7626 w_7625 w_7624 w_7623 w_7622 w_7621 w_7620 w_7798 w_7797 w_7807 w_7806 w_7805 w_7804 w_7803 w_7802 w_7801 w_7800 w_7799 w_7816 w_7815 w_7814 w_7813 w_7812 w_7811 w_7810 w_7809 w_7808 w_7986 w_7985 w_7995 w_7994 w_7993 w_7992 w_7991 w_7990 w_7989 w_7988 w_7987 w_8004 w_8003 w_8002 w_8001 w_8000 w_7999 w_7998 w_7997 w_7996 w_8176 w_8175 w_8174 w_8173 w_8184 w_8183 w_8182 w_8181 w_8180 w_8179 w_8178 w_8177 w_8192 w_8191 w_8190 w_8189 w_8188 w_8187 w_8186 w_8185 w_8364 w_8363 w_8362 w_8361 w_8372 w_8371 w_8370 w_8369 w_8368 w_8367 w_8366 w_8365 w_8380 w_8379 w_8378 w_8377 w_8376 w_8375 w_8374 w_8373 w_8552 w_8551 w_8550 w_8549 w_8560 w_8559 w_8558 w_8557 w_8556 w_8555 w_8554 w_8553 w_8568 w_8567 w_8566 w_8565 w_8564 w_8563 w_8562 w_8561 w_8742 w_8741 w_8740 w_8739 w_8738 w_8737 w_8749 w_8748 w_8747 w_8746 w_8745 w_8744 w_8743 w_8756 w_8755 w_8754 w_8753 w_8752 w_8751 w_8750 w_8930 w_8929 w_8928 w_8927 w_8926 w_8925 w_8937 w_8936 w_8935 w_8934 w_8933 w_8932 w_8931 w_8944 w_8943 w_8942 w_8941 w_8940 w_8939 w_8938 w_9118 w_9117 w_9116 w_9115 w_9114 w_9113 w_9125 w_9124 w_9123 w_9122 w_9121 w_9120 w_9119 w_9132 w_9131 w_9130 w_9129 w_9128 w_9127 w_9126 w_9306 w_9305 w_9304 w_9303 w_9302 w 9301 w_9313 w_9312 w_9311 w_9310 w_9309 w_9308 w_9307 w_9320 w_9319 w_9318 w_9317 w_9316 w_9315 w_9314 w_ w_94 w 948 w w w_ w_ w_9 w_94 w_949 Smoother with spatially-varying weights (a) Observed ozone concentration (b) Predicted ozone concentration To each point s in the domain S with observation Y (s) we associate x(s), where x(s)= g k=1 w k(s)x(b k ) Y (s) x(s)
31 Spatio-temporal setting
32 A MEM for outputs from computer model Y (s,t) observations at s and time t, x(b,t) numerical model output at grid cell B and time t x(b,t) = Ṽ (B,t) + η(b,t) with Ṽ (B,t) = µ t + V (B,t) and V t = {V (B,t)}. ind η(b,t) N(0,σ 2 )
33 A MEM for outputs from computer model Y (s,t) observations at s and time t, x(b,t) numerical model output at grid cell B and time t x(b,t) = Ṽ (B,t) + η(b,t) with Ṽ (B,t) = µ t + V (B,t) and V t = {V (B,t)}. For s B: ind η(b,t) N(0,σ 2 ) Y (s,t) = β 0 (s,t) + β 1,t Ṽ (B,t) + ε(s,t) ε(s,t) ind N(0,τ 2 ) with β 0 (s,t) = β 0,t + β 0 (s,t).
34 A MEM for outputs from computer model Y (s,t) observations at s and time t, x(b,t) numerical model output at grid cell B and time t x(b,t) = Ṽ (B,t) + η(b,t) with Ṽ (B,t) = µ t + V (B,t) and V t = {V (B,t)}. For s B: ind η(b,t) N(0,σ 2 ) Y (s,t) = β 0 (s,t) + β 1,t Ṽ (B,t) + ε(s,t) ε(s,t) ind N(0,τ 2 ) with β 0 (s,t) = β 0,t + β 0 (s,t). Temporal dependence in β i,t, i = 0,1, β 0 (s,t) and V t : either (i) independent in time, or (ii) dynamic.
35 Smoother with spatially-varying weights Y (s,t) observations at s and time t, x(b,t) numerical model output at grid cell B and time t For s B : Y (s,t) = β 0 (s,t) + β 1,t x(s,t) + ε(s,t) ε(s,t) ind N(0,τ 2 ) with β 0 (s,t) = β 0,t + β 0 (s,t), x(s,t) = g k=1 w k,t(s)x(b k,t), and w k,t (s) random and varying both in space and time.
36 Smoother with spatially-varying weights Y (s,t) observations at s and time t, x(b,t) numerical model output at grid cell B and time t For s B : Y (s,t) = β 0 (s,t) + β 1,t x(s,t) + ε(s,t) ε(s,t) ind N(0,τ 2 ) with β 0 (s,t) = β 0,t + β 0 (s,t), x(s,t) = g k=1 w k,t(s)x(b k,t), and w k,t (s) random and varying both in space and time. w k,t (s) := K (s r k;ψ) exp(q(r k,t)) g l=1 K (s r l;ψ) exp(q(r l,t))
37 Smoother with spatially-varying weights Y (s,t) observations at s and time t, x(b,t) numerical model output at grid cell B and time t For s B : Y (s,t) = β 0 (s,t) + β 1,t x(s,t) + ε(s,t) ε(s,t) ind N(0,τ 2 ) with β 0 (s,t) = β 0,t + β 0 (s,t), x(s,t) = g k=1 w k,t(s)x(b k,t), and w k,t (s) random and varying both in space and time. w k,t (s) := K (s r k;ψ) exp(q(r k,t)) g l=1 K (s r l;ψ) exp(q(r l,t)) Temporal dependence in β i,t, i = 0,1, β 0 (s,t) and Q(s,t): either (i) independent in time, or (ii) dynamic.
38 Application Training sites Validation sites Daily maximum 8-hour ozone concentration (ppb): observations and CMAQ model output. Model output on 12-km grid cells (g =40,044). Period: June 1- August 31, monitoring sites: 400 training sites and 400 validation sites Transformation to square root scale; predictions were back-transformed Time-varying parameters: independent in time
39 Results Assessed predictive performance at the 400 validation sites for the period June 1 - August 31, Coverage Avg. length Method MSPE MAPE 95% PI 95% PI Downscaler % 31.6 MEM % 29.1 Smoother % 28.8
40 Results Downscaler: x(b,t) Y (s,t) MEM: Ṽ (B,t) Y (s,t) Smoother: x(s,t) Y (s,t)
41 Results Downscaler: x(b,t) Y (s,t) MEM: Ṽ (B,t) Y (s,t) Smoother: x(s,t) Y (s,t) Computed the correlation between Y (s,t) and: (i) x(b,t) and the posterior predictive mean of (ii) Ṽ (B,t) and (iii) x(s,t).
42 Results Downscaler: x(b,t) Y (s,t) MEM: Ṽ (B,t) Y (s,t) Smoother: x(s,t) Y (s,t) Computed the correlation between Y (s,t) and: (i) x(b,t) and the posterior predictive mean of (ii) Ṽ (B,t) and (iii) x(s,t). Posterior mean of Posterior mean Day x(b,t) Ṽ (B,t) of x(s,t) July 4, July 20, August 9,
43 Results Posterior mean of the random and spatially and temporally- varying weights w k,t (s) for July 4, 2001 at four sites s3 s4 s2 s
44 Results Posterior mean of the random and spatially and temporally- varying weights w k,t (s) for July 4, 2001 at four sites
45 Results Posterior mean of the random and spatially and temporally- varying weights w k,t (s) for July 4, 2001, July 20, and August 9, 2001 at two sites
46 Conclusions Presented two methods to extend the downscaler model to fuse numerical model ouput and monitoring data Both are computationally feasibile even for large computer model output (smoother + predictive processes (Banerjee et al., JRSS B, 2008)) Smoother with spatially-varying weights yielded better predictions than other methods Clear directionality in the weights of the smoother model
Space-time downscaling under error in computer model output
Space-time downscaling under error in computer model output University of Michigan Department of Biostatistics joint work with Alan E. Gelfand, David M. Holland, Peter Guttorp and Peter Craigmile Introduction
More informationSpace-Time Data fusion Under Error in Computer Model Output: An Application to Modeling Air Quality
Biometrics 68, 837 848 September 2012 DOI: 10.1111/j.1541-0420.2011.01725.x Space-Time Data fusion Under Error in Computer Model Output: An Application to Modeling Air Quality Veronica J. Berrocal, 1,
More informationA Spatio-Temporal Downscaler for Output From Numerical Models
Supplementary materials for this article are available at 10.1007/s13253-009-0004-z. A Spatio-Temporal Downscaler for Output From Numerical Models Veronica J. BERROCAL,AlanE.GELFAND, and David M. HOLLAND
More informationFusing point and areal level space-time data. data with application to wet deposition
Fusing point and areal level space-time data with application to wet deposition Alan Gelfand Duke University Joint work with Sujit Sahu and David Holland Chemical Deposition Combustion of fossil fuel produces
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 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 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 informationBayesian hierarchical modelling for data assimilation of past observations and numerical model forecasts
Bayesian hierarchical modelling for data assimilation of past observations and numerical model forecasts Stan Yip Exeter Climate Systems, University of Exeter c.y.yip@ex.ac.uk Joint work with Sujit Sahu
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 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 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 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 informationState Space Representation of Gaussian Processes
State Space Representation of Gaussian Processes Simo Särkkä Department of Biomedical Engineering and Computational Science (BECS) Aalto University, Espoo, Finland June 12th, 2013 Simo Särkkä (Aalto University)
More informationNRCSE. Misalignment and use of deterministic models
NRCSE Misalignment and use of deterministic models Work with Veronica Berrocal Peter Craigmile Wendy Meiring Paul Sampson Gregory Nikulin The choice of spatial scale some questions 1. Which spatial scale
More informationBayesian spatial quantile regression
Brian J. Reich and Montserrat Fuentes North Carolina State University and David B. Dunson Duke University E-mail:reich@stat.ncsu.edu Tropospheric ozone Tropospheric ozone has been linked with several adverse
More informationStatistics for extreme & sparse data
Statistics for extreme & sparse data University of Bath December 6, 2018 Plan 1 2 3 4 5 6 The Problem Climate Change = Bad! 4 key problems Volcanic eruptions/catastrophic event prediction. Windstorms
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 informationGroundwater permeability
Groundwater permeability Easy to solve the forward problem: flow of groundwater given permeability of aquifer Inverse problem: determine permeability from flow (usually of tracers With some models enough
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 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 informationSTA414/2104 Statistical Methods for Machine Learning II
STA414/2104 Statistical Methods for Machine Learning II Murat A. Erdogdu & David Duvenaud Department of Computer Science Department of Statistical Sciences Lecture 3 Slide credits: Russ Salakhutdinov Announcements
More informationHierarchical Bayesian auto-regressive models for large space time data with applications to ozone concentration modelling
Hierarchical Bayesian auto-regressive models for large space time data with applications to ozone concentration modelling Sujit Kumar Sahu School of Mathematics, University of Southampton, Southampton,
More informationDisease mapping with Gaussian processes
EUROHEIS2 Kuopio, Finland 17-18 August 2010 Aki Vehtari (former Helsinki University of Technology) Department of Biomedical Engineering and Computational Science (BECS) Acknowledgments Researchers - Jarno
More informationCombining Monitoring Data and Computer model Output in Assessing Environmental Exposure
Combining Monitoring Data and Computer model Output in Assessing Environmental Exposure Alan E. GELFAND and Sujit K. SAHU December 21, 28 Abstract In the environmental exposure community numerical models
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 informationFully Bayesian Deep Gaussian Processes for Uncertainty Quantification
Fully Bayesian Deep Gaussian Processes for Uncertainty Quantification N. Zabaras 1 S. Atkinson 1 Center for Informatics and Computational Science Department of Aerospace and Mechanical Engineering University
More informationHierarchical Bayesian models for space-time air pollution data
Hierarchical Bayesian models for space-time air pollution data Sujit K. Sahu June 14, 2011 Abstract Recent successful developments in Bayesian modeling and computation for large and complex data sets have
More informationOnline Updating of Computer Model Output Using Real-Time Sensor Data
Supplementary materials for this article are available online. Please go to http://www.tandfonline.com/r/tech Online Updating of Computer Model Output Using Real-Time Sensor Data Huijing JIANG IBM Watson
More informationLecture 16 Deep Neural Generative Models
Lecture 16 Deep Neural Generative Models CMSC 35246: Deep Learning Shubhendu Trivedi & Risi Kondor University of Chicago May 22, 2017 Approach so far: We have considered simple models and then constructed
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 informationModel Selection for Geostatistical Models
Model Selection for Geostatistical Models Richard A. Davis Colorado State University http://www.stat.colostate.edu/~rdavis/lectures Joint work with: Jennifer A. Hoeting, Colorado State University Andrew
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 informationFusing point and areal level space-time data with application
Fusing point and areal level space-time data with application to wet deposition Sujit K. Sahu University of Southampton, UK Alan E. Gelfand Duke University, USA David M. Holland Environmental Protection
More informationGaussian Processes (10/16/13)
STA561: Probabilistic machine learning Gaussian Processes (10/16/13) Lecturer: Barbara Engelhardt Scribes: Changwei Hu, Di Jin, Mengdi Wang 1 Introduction In supervised learning, we observe some inputs
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 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 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 informationHierarchical Modeling for Multivariate Spatial Data
Hierarchical Modeling for Multivariate Spatial Data Sudipto Banerjee 1 and Andrew O. Finley 2 1 Biostatistics, School of Public Health, University of Minnesota, Minneapolis, Minnesota, U.S.A. 2 Department
More informationData are collected along transect lines, with dense data along. Spatial modelling using GMRFs. 200m. Today? Spatial modelling using GMRFs
Centre for Mathematical Sciences Lund University Engineering geology Lund University Results A non-stationary extension The model Estimation Gaussian Markov random fields Basics Approximating Mate rn covariances
More informationCalibration of numerical model output using nonparametric spatial density functions
Calibration of numerical model output using nonparametric spatial density functions Jingwen Zhou 1, Montserrat Fuentes 1, and Jerry Davis 2 1 North Carolina State University, Department of Statistics,
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 informationNonparametric Bayesian Methods (Gaussian Processes)
[70240413 Statistical Machine Learning, Spring, 2015] Nonparametric Bayesian Methods (Gaussian Processes) Jun Zhu dcszj@mail.tsinghua.edu.cn http://bigml.cs.tsinghua.edu.cn/~jun State Key Lab of Intelligent
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 SAE using Complex Survey Data Lecture 4A: Hierarchical Spatial Bayes Modeling
Bayesian SAE using Complex Survey Data Lecture 4A: Hierarchical Spatial Bayes Modeling Jon Wakefield Departments of Statistics and Biostatistics University of Washington 1 / 37 Lecture Content Motivation
More informationSpatial Lasso with Application to GIS Model Selection. F. Jay Breidt Colorado State University
Spatial Lasso with Application to GIS Model Selection F. Jay Breidt Colorado State University with Hsin-Cheng Huang, Nan-Jung Hsu, and Dave Theobald September 25 The work reported here was developed under
More informationIntroduction. Chapter 1
Chapter 1 Introduction In this book we will be concerned with supervised learning, which is the problem of learning input-output mappings from empirical data (the training dataset). Depending on the characteristics
More informationADVANCED MACHINE LEARNING ADVANCED MACHINE LEARNING. Non-linear regression techniques Part - II
1 Non-linear regression techniques Part - II Regression Algorithms in this Course Support Vector Machine Relevance Vector Machine Support vector regression Boosting random projections Relevance vector
More informationGaussian Processes. 1 What problems can be solved by Gaussian Processes?
Statistical Techniques in Robotics (16-831, F1) Lecture#19 (Wednesday November 16) Gaussian Processes Lecturer: Drew Bagnell Scribe:Yamuna Krishnamurthy 1 1 What problems can be solved by Gaussian Processes?
More informationStatistical downscaling daily rainfall statistics from seasonal forecasts using canonical correlation analysis or a hidden Markov model?
Statistical downscaling daily rainfall statistics from seasonal forecasts using canonical correlation analysis or a hidden Markov model? Andrew W. Robertson International Research Institute for Climate
More informationNon-gaussian spatiotemporal modeling
Dec, 2008 1/ 37 Non-gaussian spatiotemporal modeling Thais C O da Fonseca Joint work with Prof Mark F J Steel Department of Statistics University of Warwick Dec, 2008 Dec, 2008 2/ 37 1 Introduction Motivation
More informationA Method for Estimating Urban Background Concentrations in Support of Hybrid Air Pollution Modeling for Environmental Health Studies
Supplementary Information A Method for Estimating Urban Background Concentrations in Support of Hybrid Air Pollution Modeling for Environmental Health Studies 1. The 64 Unique Monitoring Locations Used
More informationModels for models. Douglas Nychka Geophysical Statistics Project National Center for Atmospheric Research
Models for models Douglas Nychka Geophysical Statistics Project National Center for Atmospheric Research Outline Statistical models and tools Spatial fields (Wavelets) Climate regimes (Regression and clustering)
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 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 informationSpatial bias modeling with application to assessing remotely-sensed aerosol as a proxy for particulate matter
Spatial bias modeling with application to assessing remotely-sensed aerosol as a proxy for particulate matter Chris Paciorek Department of Biostatistics Harvard School of Public Health application joint
More informationWhat is (certain) Spatio-Temporal Data?
What is (certain) Spatio-Temporal Data? A spatio-temporal database stores triples (oid, time, loc) In the best case, this allows to look up the location of an object at any time 2 What is (certain) Spatio-Temporal
More informationSTOCHASTIC MODELING OF ENVIRONMENTAL TIME SERIES. Richard W. Katz LECTURE 5
STOCHASTIC MODELING OF ENVIRONMENTAL TIME SERIES Richard W Katz LECTURE 5 (1) Hidden Markov Models: Applications (2) Hidden Markov Models: Viterbi Algorithm (3) Non-Homogeneous Hidden Markov Model (1)
More informationBayesian data analysis in practice: Three simple examples
Bayesian data analysis in practice: Three simple examples Martin P. Tingley Introduction These notes cover three examples I presented at Climatea on 5 October 0. Matlab code is available by request to
More informationSequential Monte Carlo Methods for Bayesian Computation
Sequential Monte Carlo Methods for Bayesian Computation A. Doucet Kyoto Sept. 2012 A. Doucet (MLSS Sept. 2012) Sept. 2012 1 / 136 Motivating Example 1: Generic Bayesian Model Let X be a vector parameter
More informationGeostatistical 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 informationLecture 6: Bayesian Inference in SDE Models
Lecture 6: Bayesian Inference in SDE Models Bayesian Filtering and Smoothing Point of View Simo Särkkä Aalto University Simo Särkkä (Aalto) Lecture 6: Bayesian Inference in SDEs 1 / 45 Contents 1 SDEs
More informationGaussian Process Approximations of Stochastic Differential Equations
Gaussian Process Approximations of Stochastic Differential Equations Cédric Archambeau Dan Cawford Manfred Opper John Shawe-Taylor May, 2006 1 Introduction Some of the most complex models routinely run
More informationof the 7 stations. In case the number of daily ozone maxima in a month is less than 15, the corresponding monthly mean was not computed, being treated
Spatial Trends and Spatial Extremes in South Korean Ozone Seokhoon Yun University of Suwon, Department of Applied Statistics Suwon, Kyonggi-do 445-74 South Korea syun@mail.suwon.ac.kr Richard L. Smith
More informationSpatio-temporal prediction of site index based on forest inventories and climate change scenarios
Forest Research Institute Spatio-temporal prediction of site index based on forest inventories and climate change scenarios Arne Nothdurft 1, Thilo Wolf 1, Andre Ringeler 2, Jürgen Böhner 2, Joachim Saborowski
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 informationPILCO: A Model-Based and Data-Efficient Approach to Policy Search
PILCO: A Model-Based and Data-Efficient Approach to Policy Search (M.P. Deisenroth and C.E. Rasmussen) CSC2541 November 4, 2016 PILCO Graphical Model PILCO Probabilistic Inference for Learning COntrol
More informationGeneralized Spatial Dirichlet Process Models
Generalized Spatial Dirichlet Process Models By JASON A. DUAN Institute of Statistics and Decision Sciences at Duke University, Durham, North Carolina, 27708-0251, U.S.A. jason@stat.duke.edu MICHELE GUINDANI
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 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 informationStudy Notes on the Latent Dirichlet Allocation
Study Notes on the Latent Dirichlet Allocation Xugang Ye 1. Model Framework A word is an element of dictionary {1,,}. A document is represented by a sequence of words: =(,, ), {1,,}. A corpus is a collection
More informationBayesian Dynamic Linear Modelling for. Complex Computer Models
Bayesian Dynamic Linear Modelling for Complex Computer Models Fei Liu, Liang Zhang, Mike West Abstract Computer models may have functional outputs. With no loss of generality, we assume that a single computer
More informationBayesian Forecasting Using Spatio-temporal Models with Applications to Ozone Concentration Levels in the Eastern United States
1 Bayesian Forecasting Using Spatio-temporal Models with Applications to Ozone Concentration Levels in the Eastern United States Sujit Kumar Sahu Southampton Statistical Sciences Research Institute, University
More informationA short introduction to INLA and R-INLA
A short introduction to INLA and R-INLA Integrated Nested Laplace Approximation Thomas Opitz, BioSP, INRA Avignon Workshop: Theory and practice of INLA and SPDE November 7, 2018 2/21 Plan for this talk
More informationStatistical analysis of regional climate models. Douglas Nychka, National Center for Atmospheric Research
Statistical analysis of regional climate models. Douglas Nychka, National Center for Atmospheric Research National Science Foundation Olso workshop, February 2010 Outline Regional models and the NARCCAP
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 informationSTAT 518 Intro Student Presentation
STAT 518 Intro Student Presentation Wen Wei Loh April 11, 2013 Title of paper Radford M. Neal [1999] Bayesian Statistics, 6: 475-501, 1999 What the paper is about Regression and Classification Flexible
More informationBayesian spatial hierarchical modeling for temperature extremes
Bayesian spatial hierarchical modeling for temperature extremes Indriati Bisono Dr. Andrew Robinson Dr. Aloke Phatak Mathematics and Statistics Department The University of Melbourne Maths, Informatics
More informationGaussian Processes 1. Schedule
1 Schedule 17 Jan: Gaussian processes (Jo Eidsvik) 24 Jan: Hands-on project on Gaussian processes (Team effort, work in groups) 31 Jan: Latent Gaussian models and INLA (Jo Eidsvik) 7 Feb: Hands-on project
More informationFactor Analysis and Kalman Filtering (11/2/04)
CS281A/Stat241A: Statistical Learning Theory Factor Analysis and Kalman Filtering (11/2/04) Lecturer: Michael I. Jordan Scribes: Byung-Gon Chun and Sunghoon Kim 1 Factor Analysis Factor analysis is used
More informationProbabilistic Graphical Models
Probabilistic Graphical Models Brown University CSCI 2950-P, Spring 2013 Prof. Erik Sudderth Lecture 13: Learning in Gaussian Graphical Models, Non-Gaussian Inference, Monte Carlo Methods Some figures
More informationLecture 6: Multiple Model Filtering, Particle Filtering and Other Approximations
Lecture 6: Multiple Model Filtering, Particle Filtering and Other Approximations Department of Biomedical Engineering and Computational Science Aalto University April 28, 2010 Contents 1 Multiple Model
More informationGaussian Process Approximations of Stochastic Differential Equations
Gaussian Process Approximations of Stochastic Differential Equations Cédric Archambeau Centre for Computational Statistics and Machine Learning University College London c.archambeau@cs.ucl.ac.uk CSML
More informationLaplace Distribution Based Stochastic Models
Laplace Distribution Based Stochastic Models Anastassia Baxevani Department of Mathematics and Statistics University of Cyprus results based on work of/with Podgorski, Wegener, Aberg, Bolin, Walin July
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 informationParametric Techniques Lecture 3
Parametric Techniques Lecture 3 Jason Corso SUNY at Buffalo 22 January 2009 J. Corso (SUNY at Buffalo) Parametric Techniques Lecture 3 22 January 2009 1 / 39 Introduction In Lecture 2, we learned how to
More informationSpatiotemporal Analysis of Solar Radiation for Sustainable Research in the Presence of Uncertain Measurements
Spatiotemporal Analysis of Solar Radiation for Sustainable Research in the Presence of Uncertain Measurements Alexander Kolovos SAS Institute, Inc. alexander.kolovos@sas.com Abstract. The study of incoming
More informationA graph contains a set of nodes (vertices) connected by links (edges or arcs)
BOLTZMANN MACHINES Generative Models Graphical Models A graph contains a set of nodes (vertices) connected by links (edges or arcs) In a probabilistic graphical model, each node represents a random variable,
More informationNew Insights into History Matching via Sequential Monte Carlo
New Insights into History Matching via Sequential Monte Carlo Associate Professor Chris Drovandi School of Mathematical Sciences ARC Centre of Excellence for Mathematical and Statistical Frontiers (ACEMS)
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 informationLinear Dynamical Systems
Linear Dynamical Systems Sargur N. srihari@cedar.buffalo.edu Machine Learning Course: http://www.cedar.buffalo.edu/~srihari/cse574/index.html Two Models Described by Same Graph Latent variables Observations
More informationCBMS Lecture 1. Alan E. Gelfand Duke University
CBMS Lecture 1 Alan E. Gelfand Duke University Introduction to spatial data and models Researchers in diverse areas such as climatology, ecology, environmental exposure, public health, and real estate
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 informationA statistical framework to combine multivariate spatial data and physical models for hurricane surface wind prediction
A statistical framework to combine multivariate spatial data and physical models for hurricane surface wind prediction Kristen M. Foley 1, Montserrat Fuentes 1 1 Department of Statistics, North Carolina
More informationOn Statistical Approaches to Generate Level 3 Products from Satellite Remote Sensing Retrievals. and Clint Shumack
remote sensing Article On Statistical Approaches to Generate Level 3 Products from Satellite Remote Sensing Retrievals Andrew Zammit-Mangion * ID, Noel Cressie ID and Clint Shumack National Institute for
More informationParametric Techniques
Parametric Techniques Jason J. Corso SUNY at Buffalo J. Corso (SUNY at Buffalo) Parametric Techniques 1 / 39 Introduction When covering Bayesian Decision Theory, we assumed the full probabilistic structure
More informationspbayes: An R Package for Univariate and Multivariate Hierarchical Point-referenced Spatial Models
spbayes: An R Package for Univariate and Multivariate Hierarchical Point-referenced Spatial Models Andrew O. Finley 1, Sudipto Banerjee 2, and Bradley P. Carlin 2 1 Michigan State University, Departments
More informationThree data analysis problems
Three data analysis problems Andreas Zezas University of Crete CfA Two types of problems: Fitting Source Classification Fitting: complex datasets Fitting: complex datasets Maragoudakis et al. in prep.
More informationEmpirical Bayes methods for the transformed Gaussian random field model with additive measurement errors
1 Empirical Bayes methods for the transformed Gaussian random field model with additive measurement errors Vivekananda Roy Evangelos Evangelou Zhengyuan Zhu CONTENTS 1.1 Introduction......................................................
More informationLearning to Control an Octopus Arm with Gaussian Process Temporal Difference Methods
Learning to Control an Octopus Arm with Gaussian Process Temporal Difference Methods Yaakov Engel Joint work with Peter Szabo and Dmitry Volkinshtein (ex. Technion) Why use GPs in RL? A Bayesian approach
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 information