What s for today. Continue to discuss about nonstationary models Moving windows Convolution model Weighted stationary model
|
|
- Josephine Miles
- 5 years ago
- Views:
Transcription
1 What s for today Continue to discuss about nonstationary models Moving windows Convolution model Weighted stationary model c Mikyoung Jun (Texas A&M) Stat647 Lecture 11 October 2, / 23
2 Nonstationary covariance models If you suspect the process is nonstationary, you need to have certain ideas about the types of nonstationarity For example, as in the example from the previous lecture, you may see if the covariance has dependence on latitudes (or certain coordinates) or whether the fact that the domain is land or sea Even if you know the types of nonstationarity, it is often hard to come up with valid classes of nonstationary covariance functions This is why I started from a process in the previous example, not covariance functions directly c Mikyoung Jun (Texas A&M) Stat647 Lecture 11 October 2, / 23
3 An example of nonstationary covariance function An interesting application problem is the study of point source problem for air pollution The model in Hughes-Oliver et al. (1998a), Statistics & Probability Letters: Cor(Z(s i,s j )) = exp( θ 1 s i s j exp (θ 2 c i c j + θ 3 min[c i,c j ])) where c i = s i c and c j = s j c are the distances between sites and the point source c Note that if θ 2 = θ 3 = 0, then it reduces to an isotropic model What is the range of this model in general? c Mikyoung Jun (Texas A&M) Stat647 Lecture 11 October 2, / 23
4 An example of nonstationary covariance function It is difficult to show that this is a valid model What should you do then? c Mikyoung Jun (Texas A&M) Stat647 Lecture 11 October 2, / 23
5 Another example of nonstationary covariance model Let s consider another simple nonstationary model that was considered in the previous lecture For example, suppose Z 0 (s) is an isotropic process and we consider Z(s) = a(s)z 0 (s) which has nonstationary covariance function What is the condition on a required to have a valid nonstationary covariance model? Note that the covariance function of Z is nonstationary but the correlation function of Z is not a(s 1 )a(s 2 )Cov(Z 0 (s 1 ), Z 0 (s 2 )) Cor(Z(s 1 ), Z(s 2 )) = a(s 1 )a(s 2 ) Var(Z 0 (s 1 ))Var(Z 0 (s 1 )) It is in general harder to develop nonstationary correlation model, since nonstationary covariance does not necessarily mean nonstationary correlation c Mikyoung Jun (Texas A&M) Stat647 Lecture 11 October 2, / 23
6 A note on nonstationarity You may have nonstationary process purely due to the mean part if the mean changes over space In this class, we focus on nonstationarity in the covariance structure since dealing with nonstationarity in the mean is fairly easier However, you need to have good model for the mean so that you can focus on correctly modeling the covariance (practically with the residuals) c Mikyoung Jun (Texas A&M) Stat647 Lecture 11 October 2, / 23
7 Local stationarity Even if the process is nonstationary in general, in some situations it may be reasonable to assume that 1 the process is stationary in smaller subregions 2 the process can be decomposed into simpler, stationary processes There are generally three approaches for local stationarity: 1 moving windows 2 convolutions 3 weighted stationary processes c Mikyoung Jun (Texas A&M) Stat647 Lecture 11 October 2, / 23
8 Moving windows The basic idea is to split the spatial region into smaller subregions and fit (isotropic) covariance models in the subregions separately Then you have sort of nonstationary model since your parameter values are different in different subregions The number of the windows as well as the size and shape of the windows should be determined This approach can save some computational time in the sense that you only fit small number of data c Mikyoung Jun (Texas A&M) Stat647 Lecture 11 October 2, / 23
9 Moving windows Suppose our model is Z(s) = µ1 + e(s) and e (0,Σ(θ)) The ordinary kriging predictor at a new location s 0 given the observations at the locations s i, i = 1,,n is given by where Z = (Z(s 1 ),,Z(s n )) T Ẑ(s 0 ) = ˆµ + c T Σ 1 (Z µ1) It may be the case that Z at some location s j has little impact on the kriging predictor at the location s 0 if s j is far away from s 0 Using the moving windows idea, you can get the kriging predictor on the local neighborhood c Mikyoung Jun (Texas A&M) Stat647 Lecture 11 October 2, / 23
10 Convolution method We start from a mean zero white noise X(s) such that E(X(s)) = 0, Var(X(s)) = σ 2, and Cov(X(s 1 ),X(s 2 )) = 0 if i j Then we can construct a stationary random field Z(s) by convolving X(s) with a kernel function K s (u) centered at s such that Z(s) = K s (u)x(u)du all u What are the mean and the covariance of Z? Another way of looking at it in terms of summation: c Mikyoung Jun (Texas A&M) Stat647 Lecture 11 October 2, / 23
11 Convolution method Note that the covariance function depends on the choice of the kernel Here are some ways of choosing the kernel function Isotropic convolution: K s (u) = (2πσ 2 ) 1 exp ( 0.5((u x s x ) 2 + (u y s y ) 2 )) Note that the result of the convolution is a stationary, isotropic random field with Gaussian covariance structure Anisotropic convolution: K s (u) = A 1/2 (2π) 1 exp ( 0.5(u s) T A 1 (u s)) where A is a positive definite (but not identity) matrix This creates a geometrically anisotropic Gaussian covariance structure Nonstationary convolution: K s (u) = A(s) 1/2 (2π) 1 exp ( 0.5(u s) T A(s) 1 (u s)) This creates nonstationary covariance function. However, you have so many parameters to estimate c Mikyoung Jun (Texas A&M) Stat647 Lecture 11 October 2, / 23
12 More on convolution method This idea mainly started from the paper, Higdon, Swall, and Kern (1999), Bayesian statistics There is also a nice paper that introduces the method and has an R code (Higdon, 2002, Quantitative methods for current environmental issues c Mikyoung Jun (Texas A&M) Stat647 Lecture 11 October 2, / 23
13 More on convolution method As we ve seen, the convolution method also applies for isotropic or stationary model For stationary model, we can write Z(s) = K(u s)x(u)du all u Then the resulting covariance function is given by Cov(Z(s),Z(t)) = K(u (s t))k(u)du c Mikyoung Jun (Texas A&M) Stat647 Lecture 11 October 2, / 23
14 Kernels and resulting covariance function (picture taken from Higdon (2002)) c Mikyoung Jun (Texas A&M) Stat647 Lecture 11 October 2, / 23
15 Nonstationary kernels (picture taken from Chris Paciorek s website) c Mikyoung Jun (Texas A&M) Stat647 Lecture 11 October 2, / 23
16 Higdon (1998), Environmental and Ecological Statistics The paper is on A process-convolution approach to modelling temperatures in the North Atlantic Ocean (note this paper deals with space-time data) First define a grid process x = (x 1,,x m ) on (w 1,τ 1 ),,(w m,τ m ) For x j independent N(µ x,σx 2 ), we let m Z(s,t) = K s (s w j,t τ j ) x j j=1 Here K s ( s, t) = C S ( s) R( t) For geometrically anisotropic Gaussian kernels on (w 1,τ 1 ),,(w m,τ m ), C lk, m we let C s ( s) = w k (s)c lk ( S) and w k (s) exp ( 1 2 s l k 2 ) k=1 c Mikyoung Jun (Texas A&M) Stat647 Lecture 11 October 2, / 23
17 Data, Higdon (1998) c Mikyoung Jun (Texas A&M) Stat647 Lecture 11 October 2, / 23
18 Grid locations, Higdon (1998) c Mikyoung Jun (Texas A&M) Stat647 Lecture 11 October 2, / 23
19 Estimated mean, Higdon (1998) c Mikyoung Jun (Texas A&M) Stat647 Lecture 11 October 2, / 23
20 Weighted stationary processes We write a nonstationary process Z(s) as the weighted sum of stationary processes Z 1 (s),,z k (s): Z(s) = k w i (s)z i (s) i=1 We assume Z i s are uncorrelated In practice, we assume Z i s have same covariance function with different parameter values You can check the covariance of Z(s) is a weighed sum of the covariances of Z i c Mikyoung Jun (Texas A&M) Stat647 Lecture 11 October 2, / 23
21 Fuentes (2001), Environmetrics c Mikyoung Jun (Texas A&M) Stat647 Lecture 11 October 2, / 23
22 Fuentes (2001), Environmetrics c Mikyoung Jun (Texas A&M) Stat647 Lecture 11 October 2, / 23
What s for today. Random Fields Autocovariance Stationarity, Isotropy. c Mikyoung Jun (Texas A&M) stat647 Lecture 2 August 30, / 13
What s for today Random Fields Autocovariance Stationarity, Isotropy c Mikyoung Jun (Texas A&M) stat647 Lecture 2 August 30, 2012 1 / 13 Stochastic Process and Random Fields A stochastic process is a family
More informationWhat s for today. All about Variogram Nugget effect. Mikyoung Jun (Texas A&M) stat647 lecture 4 September 6, / 17
What s for today All about Variogram Nugget effect Mikyoung Jun (Texas A&M) stat647 lecture 4 September 6, 2012 1 / 17 What is the variogram? Let us consider a stationary (or isotropic) random field Z
More informationWhat s for today. Introduction to Space-time models. c Mikyoung Jun (Texas A&M) Stat647 Lecture 14 October 16, / 19
What s for today Introduction to Space-time models c Mikyoung Jun (Texas A&M) Stat647 Lecture 14 October 16, 2012 1 / 19 Space-time Data So far we looked at the data that vary over space Now we add another
More informationSimple example of analysis on spatial-temporal data set
Simple example of analysis on spatial-temporal data set I used the ground level ozone data in North Carolina (from Suhasini Subba Rao s website) The original data consists of 920 days of data over 72 locations
More informationAsymptotic standard errors of MLE
Asymptotic standard errors of MLE Suppose, in the previous example of Carbon and Nitrogen in soil data, that we get the parameter estimates For maximum likelihood estimation, we can use Hessian matrix
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 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 informationChapter 4 - Fundamentals of spatial processes Lecture notes
Chapter 4 - Fundamentals of spatial processes Lecture notes Geir Storvik January 21, 2013 STK4150 - Intro 2 Spatial processes Typically correlation between nearby sites Mostly positive correlation Negative
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 informationEstimation of non-stationary spatial covariance structure
Estimation of non-stationary spatial covariance structure DAVID J NOTT Department of Statistics, University of New South Wales, Sydney 2052, Australia djn@mathsunsweduau WILLIAM T M DUNSMUIR Division of
More informationLesson 4: Stationary stochastic processes
Dipartimento di Ingegneria e Scienze dell Informazione e Matematica Università dell Aquila, umberto.triacca@univaq.it Stationary stochastic processes Stationarity is a rather intuitive concept, it means
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 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 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 informationECEN 5612, Fall 2007 Noise and Random Processes Prof. Timothy X Brown NAME: CUID:
Midterm ECE ECEN 562, Fall 2007 Noise and Random Processes Prof. Timothy X Brown October 23 CU Boulder NAME: CUID: You have 20 minutes to complete this test. Closed books and notes. No calculators. If
More informationNon-stationary Cross-Covariance Models for Multivariate Processes on a Globe
Scandinavian Journal of Statistics, Vol. 38: 726 747, 2011 doi: 10.1111/j.1467-9469.2011.00751.x Published by Blackwell Publishing Ltd. Non-stationary Cross-Covariance Models for Multivariate Processes
More informationMean square continuity
Mean square continuity Suppose Z is a random field on R d We say Z is mean square continuous at s if lim E{Z(x) x s Z(s)}2 = 0 If Z is stationary, Z is mean square continuous at s if and only if K is continuous
More informationStatistics 910, #15 1. Kalman Filter
Statistics 910, #15 1 Overview 1. Summary of Kalman filter 2. Derivations 3. ARMA likelihoods 4. Recursions for the variance Kalman Filter Summary of Kalman filter Simplifications To make the derivations
More informationCross-covariance Functions for Tangent Vector Fields on the Sphere
Cross-covariance Functions for Tangent Vector Fields on the Sphere Minjie Fan 1 Tomoko Matsuo 2 1 Department of Statistics University of California, Davis 2 Cooperative Institute for Research in Environmental
More informationStatistics STAT:5100 (22S:193), Fall Sample Final Exam B
Statistics STAT:5 (22S:93), Fall 25 Sample Final Exam B Please write your answers in the exam books provided.. Let X, Y, and Y 2 be independent random variables with X N(µ X, σ 2 X ) and Y i N(µ Y, σ 2
More informationNonstationary cross-covariance models for multivariate processes on a globe
Nonstationary cross-covariance models for multivariate processes on a globe Mikyoung Jun 1 April 15, 2011 Abstract: In geophysical and environmental problems, it is common to have multiple variables of
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 informationA test for stationarity of spatio-temporal random fields on planar and spherical domains
A test for stationarity of spatio-temporal random fields on planar and spherical domains Mikyoung Jun and Marc G. Genton 1 June 13, 2010 ABSTRACT: A formal test for weak stationarity of spatial and spatio-temporal
More informationFlexible Gaussian Processes via Convolution
Flexible Gaussian Processes via Convolution Herbert K. H. Lee, Christopher H. Holloman, Catherine A. Calder, and Dave M. Higdon Duke University (and LANL) June 21, 2002 Abstract Spatial and spatio-temporal
More informationA non-stationary spatial weather generator for statistical downscaling of precipitation. Pradeebane VAITTINADA AYAR and Mathieu VRAC.
A non-stationary spatial weather generator for statistical downscaling of precipitation Pradeebane VAITTINADA AYAR and Mathieu VRAC StaRMIP Climate change impacts in the Mediterranean region Montpellier,
More informationHierarchical Modeling for Univariate Spatial Data
Hierarchical Modeling for Univariate Spatial Data Geography 890, Hierarchical Bayesian Models for Environmental Spatial Data Analysis February 15, 2011 1 Spatial Domain 2 Geography 890 Spatial Domain This
More informationModule 9: Stationary Processes
Module 9: Stationary Processes Lecture 1 Stationary Processes 1 Introduction A stationary process is a stochastic process whose joint probability distribution does not change when shifted in time or space.
More informationComparing Non-informative Priors for Estimation and Prediction in Spatial Models
Environmentrics 00, 1 12 DOI: 10.1002/env.XXXX Comparing Non-informative Priors for Estimation and Prediction in Spatial Models Regina Wu a and Cari G. Kaufman a Summary: Fitting a Bayesian model to spatial
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 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 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 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 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 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 informationEconometría 2: Análisis de series de Tiempo
Econometría 2: Análisis de series de Tiempo Karoll GOMEZ kgomezp@unal.edu.co http://karollgomez.wordpress.com Segundo semestre 2016 II. Basic definitions A time series is a set of observations X t, each
More informationPaper Review: NONSTATIONARY COVARIANCE MODELS FOR GLOBAL DATA
Paper Review: NONSTATIONARY COVARIANCE MODELS FOR GLOBAL DATA BY MIKYOUNG JUN AND MICHAEL L. STEIN Presented by Sungkyu Jung April, 2009 Outline 1 Introduction 2 Covariance Models 3 Application: Level
More informationSAMPLE PATH AND ASYMPTOTIC PROPERTIES OF SPACE-TIME MODELS. Yun Xue
SAMPLE PATH AND ASYMPTOTIC PROPERTIES OF SPACE-TIME MODELS By Yun Xue A DISSERTATION Submitted to Michigan State University in partial fulfillment of the requirements for the degree of DOCTOR OF PHILOSOPHY
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 informationOverview of Spatial Statistics with Applications to fmri
with Applications to fmri School of Mathematics & Statistics Newcastle University April 8 th, 2016 Outline Why spatial statistics? Basic results Nonstationary models Inference for large data sets An example
More informationNotes on Random Processes
otes on Random Processes Brian Borchers and Rick Aster October 27, 2008 A Brief Review of Probability In this section of the course, we will work with random variables which are denoted by capital letters,
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 informationKriging Luc Anselin, All Rights Reserved
Kriging Luc Anselin Spatial Analysis Laboratory Dept. Agricultural and Consumer Economics University of Illinois, Urbana-Champaign http://sal.agecon.uiuc.edu Outline Principles Kriging Models Spatial Interpolation
More informationEfficient Posterior Inference and Prediction of Space-Time Processes Using Dynamic Process Convolutions
Efficient Posterior Inference and Prediction of Space-Time Processes Using Dynamic Process Convolutions Catherine A. Calder Department of Statistics The Ohio State University 1958 Neil Avenue Columbus,
More informationLinear models. Rasmus Waagepetersen Department of Mathematics Aalborg University Denmark. October 5, 2016
Linear models Rasmus Waagepetersen Department of Mathematics Aalborg University Denmark October 5, 2016 1 / 16 Outline for today linear models least squares estimation orthogonal projections estimation
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 informationLecture 13: Simple Linear Regression in Matrix Format. 1 Expectations and Variances with Vectors and Matrices
Lecture 3: Simple Linear Regression in Matrix Format To move beyond simple regression we need to use matrix algebra We ll start by re-expressing simple linear regression in matrix form Linear algebra is
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 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 informationMultivariate Bayesian Linear Regression MLAI Lecture 11
Multivariate Bayesian Linear Regression MLAI Lecture 11 Neil D. Lawrence Department of Computer Science Sheffield University 21st October 2012 Outline Univariate Bayesian Linear Regression Multivariate
More informationHierarchical Modelling for Univariate Spatial Data
Hierarchical Modelling for Univariate 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 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 informationPaul Sampson Doris Damian Peter Guttorp NRCSE. T e c h n i c a l R e p o r t S e r i e s. NRCSE-TRS No January 18, 2001
Advances in Modeling and Inference for Environmental Processes with Nonstationary Spatial Covariance Paul Sampson Doris Damian Peter Guttorp NRCSE T e c h n i c a l R e p o r t S e r i e s NRCSE-TRS No.
More informationx. Figure 1: Examples of univariate Gaussian pdfs N (x; µ, σ 2 ).
.8.6 µ =, σ = 1 µ = 1, σ = 1 / µ =, σ =.. 3 1 1 3 x Figure 1: Examples of univariate Gaussian pdfs N (x; µ, σ ). The Gaussian distribution Probably the most-important distribution in all of statistics
More informationA TEST FOR STATIONARITY OF SPATIO-TEMPORAL RANDOM FIELDS ON PLANAR AND SPHERICAL DOMAINS
Statistica Sinica 22 (2012), 1737-1764 doi:http://dx.doi.org/10.5705/ss.2010.251 A TEST FOR STATIONARITY OF SPATIO-TEMPORAL RANDOM FIELDS ON PLANAR AND SPHERICAL DOMAINS Mikyoung Jun and Marc G. Genton
More informationTomoko Matsuo s collaborators
Tomoko Matsuo s collaborators Thermosphere-Ionosphere GCM Modeling A. Richmond (NCAR-HAO) T. Fuller-Rowell and M. Codrescu (NOAA) Polar Ionosphere Data Assimilation A. Richmond, G. Lu, and B. Emery (NCAR-HAO)
More information[POLS 8500] Review of Linear Algebra, Probability and Information Theory
[POLS 8500] Review of Linear Algebra, Probability and Information Theory Professor Jason Anastasopoulos ljanastas@uga.edu January 12, 2017 For today... Basic linear algebra. Basic probability. Programming
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 informationRESEARCH REPORT. Estimation of sample spacing in stochastic processes. Anders Rønn-Nielsen, Jon Sporring and Eva B.
CENTRE FOR STOCHASTIC GEOMETRY AND ADVANCED BIOIMAGING www.csgb.dk RESEARCH REPORT 6 Anders Rønn-Nielsen, Jon Sporring and Eva B. Vedel Jensen Estimation of sample spacing in stochastic processes No. 7,
More informationSpace-time analysis using a general product-sum model
Space-time analysis using a general product-sum model De Iaco S., Myers D. E. 2 and Posa D. 3,4 Università di Chieti, Pescara - ITALY; sdeiaco@tiscalinet.it 2 University of Arizona, Tucson AZ - USA; myers@math.arizona.edu
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 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 informationNew Classes of Asymmetric Spatial-Temporal Covariance Models. Man Sik Park and Montserrat Fuentes 1
New Classes of Asymmetric Spatial-Temporal Covariance Models Man Sik Park and Montserrat Fuentes 1 Institute of Statistics Mimeo Series# 2584 SUMMARY Environmental spatial data often show complex spatial-temporal
More informationReview Session: Econometrics - CLEFIN (20192)
Review Session: Econometrics - CLEFIN (20192) Part II: Univariate time series analysis Daniele Bianchi March 20, 2013 Fundamentals Stationarity A time series is a sequence of random variables x t, t =
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 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 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 informationFinal Exam. Economics 835: Econometrics. Fall 2010
Final Exam Economics 835: Econometrics Fall 2010 Please answer the question I ask - no more and no less - and remember that the correct answer is often short and simple. 1 Some short questions a) For each
More informationA Generalized Convolution Model for Multivariate Nonstationary Spatial Processes
A Generalized Convolution Model for Multivariate Nonstationary Spatial Processes Anandamayee Majumdar, Debashis Paul and Dianne Bautista Department of Mathematics and Statistics, Arizona State University,
More informationSpatial Modelling Using a New Class of Nonstationary Covariance Functions
Spatial Modelling Using a New Class of Nonstationary Covariance Functions Christopher J. Paciorek and Mark J. Schervish 13th June 2005 Abstract We introduce a new class of nonstationary covariance functions
More informationRandom vectors X 1 X 2. Recall that a random vector X = is made up of, say, k. X k. random variables.
Random vectors Recall that a random vector X = X X 2 is made up of, say, k random variables X k A random vector has a joint distribution, eg a density f(x), that gives probabilities P(X A) = f(x)dx Just
More informationECE 673-Random signal analysis I Final
ECE 673-Random signal analysis I Final Q ( point) Comment on the following statement "If cov(x ; X 2 ) 0; then the best predictor (without limitation on the type of predictor) of the random variable X
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 informationBrandon C. Kelly (Harvard Smithsonian Center for Astrophysics)
Brandon C. Kelly (Harvard Smithsonian Center for Astrophysics) Probability quantifies randomness and uncertainty How do I estimate the normalization and logarithmic slope of a X ray continuum, assuming
More informationHW5 Solutions. (a) (8 pts.) Show that if two random variables X and Y are independent, then E[XY ] = E[X]E[Y ] xy p X,Y (x, y)
HW5 Solutions 1. (50 pts.) Random homeworks again (a) (8 pts.) Show that if two random variables X and Y are independent, then E[XY ] = E[X]E[Y ] Answer: Applying the definition of expectation we have
More informationChapter 2: Fundamentals of Statistics Lecture 15: Models and statistics
Chapter 2: Fundamentals of Statistics Lecture 15: Models and statistics Data from one or a series of random experiments are collected. Planning experiments and collecting data (not discussed here). Analysis:
More informationStatistical Inference and Visualization in Scale-Space for Spatially Dependent Images
Statistical Inference and Visualization in Scale-Space for Spatially Dependent Images Amy Vaughan College of Business and Public Administration, Drake University, Des Moines, IA 0311, USA Mikyoung Jun
More informationIf we want to analyze experimental or simulated data we might encounter the following tasks:
Chapter 1 Introduction If we want to analyze experimental or simulated data we might encounter the following tasks: Characterization of the source of the signal and diagnosis Studying dependencies Prediction
More informationUnivariate Normal Distribution; GLM with the Univariate Normal; Least Squares Estimation
Univariate Normal Distribution; GLM with the Univariate Normal; Least Squares Estimation PRE 905: Multivariate Analysis Spring 2014 Lecture 4 Today s Class The building blocks: The basics of mathematical
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 information5.2 Expounding on the Admissibility of Shrinkage Estimators
STAT 383C: Statistical Modeling I Fall 2015 Lecture 5 September 15 Lecturer: Purnamrita Sarkar Scribe: Ryan O Donnell Disclaimer: These scribe notes have been slightly proofread and may have typos etc
More informationEconometría 2: Análisis de series de Tiempo
Econometría 2: Análisis de series de Tiempo Karoll GOMEZ kgomezp@unal.edu.co http://karollgomez.wordpress.com Segundo semestre 2016 IX. Vector Time Series Models VARMA Models A. 1. Motivation: The vector
More informationTwo hours. Statistical Tables to be provided THE UNIVERSITY OF MANCHESTER. 14 January :45 11:45
Two hours Statistical Tables to be provided THE UNIVERSITY OF MANCHESTER PROBABILITY 2 14 January 2015 09:45 11:45 Answer ALL four questions in Section A (40 marks in total) and TWO of the THREE questions
More informationMore Spectral Clustering and an Introduction to Conjugacy
CS8B/Stat4B: Advanced Topics in Learning & Decision Making More Spectral Clustering and an Introduction to Conjugacy Lecturer: Michael I. Jordan Scribe: Marco Barreno Monday, April 5, 004. Back to spectral
More informationLecture 4a: ARMA Model
Lecture 4a: ARMA Model 1 2 Big Picture Most often our goal is to find a statistical model to describe real time series (estimation), and then predict the future (forecasting) One particularly popular model
More informationSymmetry and Separability In Spatial-Temporal Processes
Symmetry and Separability In Spatial-Temporal Processes Man Sik Park, Montserrat Fuentes Symmetry and Separability In Spatial-Temporal Processes 1 Motivation In general, environmental data have very complex
More informationBayesian Transformed Gaussian Random Field: A Review
Bayesian Transformed Gaussian Random Field: A Review Benjamin Kedem Department of Mathematics & ISR University of Maryland College Park, MD (Victor De Oliveira, David Bindel, Boris and Sandra Kozintsev)
More informationA test for stationarity for irregularly spaced spatial data
A test for stationarity for irregularly spaced spatial data Soutir Bandyopadhyay and Suhasini Subba Rao Department of Mathematics, Lehigh University Bethlehem, PA, U.S.A. Department of Statistics, Texas
More informationITK Filters. Thresholding Edge Detection Gradients Second Order Derivatives Neighborhood Filters Smoothing Filters Distance Map Image Transforms
ITK Filters Thresholding Edge Detection Gradients Second Order Derivatives Neighborhood Filters Smoothing Filters Distance Map Image Transforms ITCS 6010:Biomedical Imaging and Visualization 1 ITK Filters:
More informationChapter 2. Some basic tools. 2.1 Time series: Theory Stochastic processes
Chapter 2 Some basic tools 2.1 Time series: Theory 2.1.1 Stochastic processes A stochastic process is a sequence of random variables..., x 0, x 1, x 2,.... In this class, the subscript always means time.
More information18 Bivariate normal distribution I
8 Bivariate normal distribution I 8 Example Imagine firing arrows at a target Hopefully they will fall close to the target centre As we fire more arrows we find a high density near the centre and fewer
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 informationSTAT 248: EDA & Stationarity Handout 3
STAT 248: EDA & Stationarity Handout 3 GSI: Gido van de Ven September 17th, 2010 1 Introduction Today s section we will deal with the following topics: the mean function, the auto- and crosscovariance
More informationMS&E 226: Small Data. Lecture 11: Maximum likelihood (v2) Ramesh Johari
MS&E 226: Small Data Lecture 11: Maximum likelihood (v2) Ramesh Johari ramesh.johari@stanford.edu 1 / 18 The likelihood function 2 / 18 Estimating the parameter This lecture develops the methodology behind
More informationHierarchical Modelling for Multivariate Spatial Data
Hierarchical Modelling for Multivariate Spatial Data Geography 890, Hierarchical Bayesian Models for Environmental Spatial Data Analysis February 15, 2011 1 Point-referenced spatial data often come as
More informationOn fixed effects estimation in spline-based semiparametric regression for spatial data
Libraries Conference on Applied Statistics in Agriculture 015-7th Annual Conference Proceedings On fixed effects estimation in spline-based semiparametric regression for spatial data Guilherme Ludwig University
More informationA Generalized Convolution Model for Multivariate Nonstationary Spatial Processes
A Generalized Convolution Model for Multivariate Nonstationary Spatial Processes Anandamayee Majumdar, Debashis Paul and Dianne Bautista Department of Mathematics and Statistics, Arizona State University,
More informationPeter Hoff Linear and multilinear models April 3, GLS for multivariate regression 5. 3 Covariance estimation for the GLM 8
Contents 1 Linear model 1 2 GLS for multivariate regression 5 3 Covariance estimation for the GLM 8 4 Testing the GLH 11 A reference for some of this material can be found somewhere. 1 Linear model Recall
More informationChapter 4 - Spatial processes R packages and software Lecture notes
TK4150 - Intro 1 Chapter 4 - Spatial processes R packages and software Lecture notes Odd Kolbjørnsen and Geir Storvik February 13, 2017 STK4150 - Intro 2 Last time General set up for Kriging type problems
More informationGeostatistics for Seismic Data Integration in Earth Models
2003 Distinguished Instructor Short Course Distinguished Instructor Series, No. 6 sponsored by the Society of Exploration Geophysicists European Association of Geoscientists & Engineers SUB Gottingen 7
More informationX t = a t + r t, (7.1)
Chapter 7 State Space Models 71 Introduction State Space models, developed over the past 10 20 years, are alternative models for time series They include both the ARIMA models of Chapters 3 6 and the Classical
More information