Bayesian Transformed Gaussian Random Field: A Review
|
|
- Ethel Caldwell
- 6 years ago
- Views:
Transcription
1 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) 1
2 Berger, De Oliveira, Sansó (2001). Objective Bayesian Analysis of Spatially Correlated Data. JASA, 96, ( ) Bindel, De Oliveira, Kedem (1997). An Implementation of the Bayesian Transformed Gaussian Spatial Prediction Model. bnk/btg page.html. Box, Cox (1964). An Analysis of Transformations (with discussion). JRSS, B 26, De Oliveira (2000). Bayesian Prediction of Clipped Gaussian Random Fields. Comp. Data Analysis, 34, ( ) De Oliveira, Kedem, Short (1997). Bayesian Prediction of Transformed Gaussian Random Fields. JASA, 92,
3 Handcock, Stein (1993). A Bayesian Analysis of Kriging. Technometrics, 35, ( ) De Oliveira, Ecker (2002). Bayesian Hot Spot Detection in the Presence of a Spatial Trend: Application to Total Nitrogen Concentration in the Chesapeake Bay. Environmetrics, 13, ( ) Kedem, Fokianos (2002). Regression Models for Time Series Analysis, New York: Wiley. Kozintsev, Kedem (2000). Generation of Similar Images From a Given Discrete Image. J. Comp. Graphical Stat., 2000, Vol. 9, No. 2, Kozintseva (1999). Comparison of Three Methods of Spatial Prediction. M.A. Thesis, Department of Mathematics, University of Maryland, College Park. 3
4 Spatial/temporal geostatistical data often display: Non-Gaussian skewed sampling distributions. Positive continuous data. Heavy right tails. Bounded support. Small data sets observed irregularly (gaps). 4
5 Possible remedies: Bayesian Transformed Gaussian (BTG): A Bayesian approach combined with parametric families of nonlinear transformations to Gaussian data. BTG provides a unified framework for inference and prediction/interpolation in a wide variety of models, Gaussian and non- Gaussian. Will describe BTG and illustrate it using spatial and temporal data. 5
6 Stationary isotropic Gaussian random field. Let {Z(s)}, s D R d, be a spatial process or a random field. A random field {Z(s)} is Gaussian if for all s 1,...,s n D, the vector (Z(s 1 ),..., Z(s n )) has a multivariate normal distribution. {Z(s)} is (second order) stationary when for s,s + h D we have ( ) E(Z(s)) = µ, ( ) Cov(Z(s + h), Z(s)) C(h). The function C( ) is called the covariogram or covariance function. 6
7 We shall assume that C(h) depends only on the distance h between the locations s + h and s but not on the direction of h. In this case the covariance function as well as the process are called isotropic. The corresponding isotropic correlation function is given by K(l) = C(l)/C(0), where l is the distance between points. Useful special case: Matérn correlation K θ (l) = ( ) 1 l θ2 ( ) 2 θ 2 1 Γ(θ 2 ) θ κ l 1 θ2 θ if l if l = 0 where θ 1 > 0, θ 2 > 0, and κ θ2 is a modified Bessel function of the third kind of order θ 2. 7
8 Z K θ (l) = ( ) 1 l θ2 ( ) 2 θ 2 1 Γ(θ 2 ) θ κ l 1 θ2 θ if l if l = 0 Matérn (θ 1 = 8, θ 2 = 3) Y X
9 Z Spherical correlation: K θ (l) = ( lθ ) ( lθ ) 3 if l θ 0 if l > θ where θ > 0 controls the correlation range. Spherical (θ = 120) Y X 9
10 Z Exponential correlation: θ 1 (0,1), θ 2 (0,2] K θ (l) = exp(l θ 2 log θ 1 ) Exponential (θ 1 = 0.5, θ 2 = 1) Y X
11 Z Rational quadratic: θ 1 > 0, θ 2 > 0. K θ (l) = ( ) θ2 1 + l2 θ1 2 Rational quadratic (θ 1 = 12, θ 2 = 8) Y X 11
12 Clipped, at 3 levels, realizations from Gaussian random fields. Top left: Matérn (8,3). Top right: spherical (120). Bottom left: exponential (0.5,1). Bottom right: rational quadratic (12,8). generate.cgi?4 Kozintsev(1999), Kozintsev and Kedem(2000). 12
13 Ordinary Kriging. Given the data Z (Z(s 1 ),..., Z(s n )) observed at locations {s 1,...,s n } in D, the problem is to predict (or estimate) Z(s 0 ) at location s 0 using the best linear unbiased predictor (BLUP) obtained by minimizing E(Z(s 0 ) n i=1 λ i Z(s i )) 2 subject to n i=1 λ i = 1 13
14 Define Then 1 = (1,1,...,1), 1 n vector c = (C(s 0 s 1 ),..., C(s 0 s n )) C = (C(s i s j )), i, j = 1,..., n λ = (λ 1, λ 2,..., λ n ) ( ˆλ = C 1 c C 1 ) c 1 C The ordinary kriging predictor is then Ẑ(s 0 ) = ˆλ Z. 14
15 Define, m = 1 1 C 1 c 1 C 1 1 and the kriging variance σ 2 k (s 0) = E(Z(s 0 ) Ẑ(s 0 )) 2 = C(0) ˆλ c + m. Under the Gaussian assumption, Ẑ(s 0 ) ± 1.96σ k (s 0 ) is a 95% prediction interval for Z(s 0 ). For non- Gaussian fields this may not hold. 15
16 Bayesian Spatial Prediction: The BTG Model. RF {Z(s), s D} observed at s 1,...,s n D Parametric family of monotone transformations G = {g λ (.) : λ Λ}. Assumption: Z(.) can be transformed into a Gaussian random field by a member of G. A useful parametric family of transformations often used in applications to normalize positive data is the Box-Cox (1964) family of power transformations, g λ (x) = { x λ 1 λ if λ 0 log(x) if λ = 0. 16
17 For some unknown transformation parameter λ Λ, {g λ (Z(s)), s D} is a Gaussian random field with E{g λ (Z(s))} = p j=1 β j f j (s), cov{g λ (Z(s)), g λ (Z(u))} = τ 1 K θ (s,u), Regression parameters: β = (β 1,..., β p ) Covariates: f(s) = (f 1 (s),..., f p (s)) Variance: τ 1 = var{g λ (Z(s))} Simplifying assumption: Isotropy, K θ (s,u) = K θ ( s u ) θ = (θ 1,..., θ q ) Θ R q. 17
18 Data: Z obs = (Z 1,obs,..., Z n,obs ) g λ (Z i,obs ) = g λ (Z(s i )) + ǫ i ; i = 1,..., n, ǫ 1,..., ǫ n are i.i.d. N(0, ξ τ ). Parameters: η = (β, τ, ξ, θ, λ). Prediction problem: Predict Z 0 = (Z(s 01 ),..., Z(s 0k )) from the predictive density function, defined by p(z o z obs ) = = Ω p(z o, η z obs )dη Ω p(z o η,z obs )p(η z obs )dη, where Ω = R p (0, ) 2 Θ Λ. 18
19 Notation: For a = (a 1,..., a n ), we write g λ (a) (g λ (a 1 ),..., g λ (a n )). The Likelihood: L(η;z obs ) = ( τ )n 2 Ψξ,θ 12 exp { τ2 } 2π Q J λ, Q = ( g λ (z obs ) Xβ ) Ψ 1 ξ,θ ( gλ (z obs ) Xβ ). X n p design matrix, X ij = f j (s i ). Ψ ξ,θ = Σ θ + ξi, n n matrix. Σ θ;ij = K θ (s i,s j ). I identity matrix. J λ = n i=1 g (z i,obs ), the Jacobian. 19
20 The Prior. Insightful arguments in Box and Cox(1964), De Oliveira, Kedem, Short (1997), as well as practical experience lead us to the prior p(η) p(ξ)p(θ)p(λ) p, τj n where p(ξ), p(θ) and p(λ) are the prior marginals of ξ, θ and λ, respectively, which are assumed to be proper. λ Unusual prior: it depends on the data through the Jacobian. For more on prior selection see Berger, De Oliveira, Sansó (2001). 20
21 Simplifying assumption: No measurement noise (ξ = 0). g λ (Z i,obs ) = g λ (Z(s i )), i = 1,..., n, p(β, τ, θ, λ) p(θ)p(λ) τj p/n λ (1) η = (β, τ, θ, λ). Also, write z = z obs 21
22 The Posterior. p(η z) = p(β, τ, θ, λ z) = p(β, τ θ, λ,z)p(θ, λ z). To get the first factor: (β τ, θ, λ,z) N p (ˆβ θ,λ, 1 τ (X Σ 1 θ X) 1 ) where (τ θ, λ,z) Ga( n p 2, 2 q θ,λ ) ˆβ θ,λ = (X Σ 1 θ X) 1 X Σ 1 θ g λ(z) q θ,λ = (g λ (z) Xˆβ θ,λ ) Σ 1 θ (g λ(z) Xˆβ θ,λ ). and p(β, τ θ, λ,z) = p(β τ, θ, λ,z)p(τ θ, λ,z) is Normal-Gamma. 22
23 To get the second factor: p(θ, λ z) Σ θ 1/2 X Σ 1 θ X 1/2 q n p 2 θ,λ J1 p n λ p(θ)p(λ) In addition to the joint posterior distribution p(η z) derived above, the predictive density p(z o z) also requires p(z o η,z). We have p(z o η,z) = ( τ 2π )k/2 D θ 1/2 g λ (z oj) j=1 { exp τ 2 (g λ(z o ) M β,θ,λ ) D 1 θ (g λ(z) M β,θ,λ ) where M β,θ,λ,d θ are known. k } 23
24 We now have the integrand p(z o η,z)p(η z) needed for p(z o z). By integrating out β and τ we obtain the simplified form of the predictive density: p(z o z) = where = Λ Λ Θ p(z o θ, λ,z)p(θ, λ z)dθdλ Θ p(z o θ, λ,z)p(z θ, λ)p(θ)p(λ)dθdλ p(z θ, λ)p(θ)p(λ)dθdλ Λ Θ p(z o θ, λ,z) = Γ(n p+k 2 ) k j=1 g λ (z oj) Γ( n p 2 )πk/2 q θ,λ C θ 1/2 and from Bayes theorem, [1 + (g λ (z o ) m θ,λ ) ( q θ,λ C θ ) 1 (g λ (z o ) m θ,λ )] n p+k 2 p(z θ, λ) Σ θ 1/2 X Σ 1 θ where m θ,λ,c θ are known. n p X 1/2 q 2 θ,λ J1 p n λ 24
25 BTG Algorithm: Predictive Density Approximation 1. Let S = {z o (j) : j = 1,..., r} be the set of values obtained by discretizing the effective range of Z Generate independently θ 1,..., θ m i.i.d. p(θ) and λ 1,..., λ m i.i.d. p(λ). 3. For z o S, the approximation to p(z o z) is given by mi=1 p(z o θ ˆp m (z o z) = i, λ i,z)p(z θ i, λ i ) mi=1 p(z θ i, λ i ) p(z o θ, λ,z) and p(z θ, λ) given above. ( ) Ẑ 0 = Median of (Z 0 Z) 25
26 Software: tkbtg application. Hybrid of C++, Tcl/Tk, and FORTRAN 77 (Bindel et al (1997)). The tkbtg Interface Layout 26
27 Example: Spatial Rainfall Prediction Rain gauge positions and weekly rainfall totals in mm, Darwin, Australia, y x 27
28 1. Use the Box-Cox transformation family. 2. λ Unif( 3,3). 3. m = Correlation: Matérn and exponential. 5. No covariate information. Assume constant regression: E{g λ (Z(s))} = β Data apparently not normal Mean=33.94 Median=29.73 SD= Total Rainfall in mm 28
29 No. z ẑ 95% PI (10.70, 56.78) (15.29, 53.35) (20.93, 52.49) (18.49, 48.29) (9.53, 56.45) (16.46, 61.79) (15.40, 33.59) (11.52, 34.66) (28.25, 100) (18.01, 52.30) (15.62, 33.40) (15.35, 37.54) (44.14, 100) (18.58, 62.63) (14.00, 30.63) (13.43, 29.81) (19.25, 73.84) (16.57, 42.46) (10.79, 27.21) (20.33, 54.39) (14.71, 31.22) (11.64, 33.74) (12.21, 40.91) (10.85, 32.81) 29
30 Z Spatial prediction and contour maps from the Darwin data using Matérn correlation Y X y x 30
31 Z Spatial prediction and contour maps from the Darwin data using exponential correlation Y X y x 31
32 Predictive densities, 95% prediction intervals, and cross-validation: Predicting a true value from the remaining 23 observations using Matérn correlation. The vertical line marks the location of the true value. p(x) True: Median: (15.60, 53.02) p(x) True: Median: (41.52, 100) x x p(x) True: Median: 22.1 (13.26, 30.94) p(x) True: Median: (14.84, 31.26) x x 32
33 Comparison of BTG With Kriging and Trans- Gaussian kriging (Kozintseva (1999)). Cross validation results using artificial data on grid. Data obtained by transforming a Gaussian (0,1) RF using inverse Box-Cox transformation. In Kriging and TG kriging λ, θ, were known. Not in BTG (!) λ = 0: Log-Normal. λ = 1: Normal. λ = 0.5: Between Normal and Log-Normal. In most cases BTG has more reliable but larger prediction intervals. BTG predicts at the original scale. TG kriging does not. 33
34 Matérn(1,10) λ KRG MSE TGK MSE BTG MSE KRG AvePI TGK AvePI BTG AvePI KRG % out 100% 48% 6% TGK % out 18% 8% 6% BTG % out 12% 6% 6% Exponential(e 0.03,1) λ KRG MSE TGK MSE BTG MSE KRG AvePI TGK AvePI BTG AvePI KRG % out 98% 64% 2% TGK % out 20% 4% 2% BTG % out 6% 2% 2% 34
35 Application of BTG to Time Series Prediction. Short time series observed irregularly. Set: s = (x, y) = (t,0). Can predict/interpolate as in state space prediction: k step prediction forward, backward, and in the middle. Example 1: Monthly data of unemployed women 20 years of age and older, Data source: Bureau of Labor Statistics. N = 48. Example 2: Monthly airline passenger data, Data source: Box-Jenkins (1976). Use only N = 36 out of 144 observations, t = 51,...,86. 35
36 Example: Prediction of Monthly number of unemployed women (Age 20), Data in hundreds of thousands. Cross validation and 95% prediction intervals. Observations at times t = 12,13,36 are outside their 95% PI s. N = 48 1 = 47. Unemployed Women No. Unemployed Women True Predicted Lower Upper t 36
37 Forward and backward one and two step prediction in the unemployed women series. In 2 step higher dispersion. One step forward Two step forward p(x) True: Median: (17.24, 26.33) p(x) True: Median: (16.36, 28.05) x x One step backward Two step backward p(x) True: Median: (21.94, 30.40) p(x) True: Median: (19.77, 29.61) x x 37
38 Example: Time series of monthly international airline passengers in thousands, January December N=144. Seasonal time series. No. Passengers t 38
39 BTG cross validation and prediction intervals for the monthly airline passengers series, t = 51,..., 86, using Matérn correlation. Observations at t = 62,63 are outside the PI. N = 36 1 = 35. No. Passengers True Predicted Lower Upper t 39
40 Application to Rainfall: Heuristic Argument. Let X n represent the area average rain rate over a region such that X n = mx n 1 + λ + ǫ n, n = 1,2,3,, where the noise {ǫ n } is a martingale difference. It can be argued that necessary conditions for X n LogNormal, n, are that m 1 and λ
41 The monotone increase in ˆm (Curve a) and the monotone decrease in ˆλ (Curve b) as a function of the square root of the area. Source: Kedem and Chiu(1987). This suggests that the lognormal distribution as a model for averages or rainfall amounts over large areas or long periods. 41
42 It is interesting to obtain the posterior p(λ z) of λ, the transformation parameter in the Box- Cox family, given the data, where the data are weekly rainfall totals from Darwin, Australia. With a uniform prior for λ, the medians of p(λ z) in 5 different weeks are: Week 1 median = 0.45 (to the left of 0). Week 2 median = 0.45 (to the right of 0). Week 3 median = 0.15 (Not far from 0). Week 4 median = 0.20 (Not far from 0). Week 5 median = 0.95 (Close to 1). 42
43 Weekly posterior p(λ z) of λ given rainfall totals from Darwin, Australia, for 5 different weeks. (a) (b) density MAE = density MAE = lambda lambda (c) (d) density MAE = 0.15 density MAE = lambda lambda (e) density MAE = lambda 43
Comparing 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 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 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 informationComparing Non-informative Priors for Estimation and. Prediction in Spatial Models
Comparing Non-informative Priors for Estimation and Prediction in Spatial Models Vigre Semester Report by: Regina Wu Advisor: Cari Kaufman January 31, 2010 1 Introduction Gaussian random fields with specified
More informationModeling and Interpolation of Non-Gaussian Spatial Data: A Comparative Study
Modeling and Interpolation of Non-Gaussian Spatial Data: A Comparative Study Gunter Spöck, Hannes Kazianka, Jürgen Pilz Department of Statistics, University of Klagenfurt, Austria hannes.kazianka@uni-klu.ac.at
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 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 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 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 informationIndex. Regression Models for Time Series Analysis. Benjamin Kedem, Konstantinos Fokianos Copyright John Wiley & Sons, Inc. ISBN.
Regression Models for Time Series Analysis. Benjamin Kedem, Konstantinos Fokianos Copyright 0 2002 John Wiley & Sons, Inc. ISBN. 0-471-36355-3 Index Adaptive rejection sampling, 233 Adjacent categories
More informationPoint-Referenced Data Models
Point-Referenced Data Models Jamie Monogan University of Georgia Spring 2013 Jamie Monogan (UGA) Point-Referenced Data Models Spring 2013 1 / 19 Objectives By the end of these meetings, participants should
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 informationFoundations of Statistical Inference
Foundations of Statistical Inference Julien Berestycki Department of Statistics University of Oxford MT 2015 Julien Berestycki (University of Oxford) SB2a MT 2015 1 / 16 Lecture 16 : Bayesian analysis
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 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 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 informationPh.D. Qualifying Exam Friday Saturday, January 6 7, 2017
Ph.D. Qualifying Exam Friday Saturday, January 6 7, 2017 Put your solution to each problem on a separate sheet of paper. Problem 1. (5106) Let X 1, X 2,, X n be a sequence of i.i.d. observations from a
More informationReliability Monitoring Using Log Gaussian Process Regression
COPYRIGHT 013, M. Modarres Reliability Monitoring Using Log Gaussian Process Regression Martin Wayne Mohammad Modarres PSA 013 Center for Risk and Reliability University of Maryland Department of Mechanical
More informationA Bayesian Treatment of Linear Gaussian Regression
A Bayesian Treatment of Linear Gaussian Regression Frank Wood December 3, 2009 Bayesian Approach to Classical Linear Regression In classical linear regression we have the following model y β, σ 2, X N(Xβ,
More informationAnalysing geoadditive regression data: a mixed model approach
Analysing geoadditive regression data: a mixed model approach Institut für Statistik, Ludwig-Maximilians-Universität München Joint work with Ludwig Fahrmeir & Stefan Lang 25.11.2005 Spatio-temporal regression
More informationMotivation Scale Mixutres of Normals Finite Gaussian Mixtures Skew-Normal Models. Mixture Models. Econ 690. Purdue University
Econ 690 Purdue University In virtually all of the previous lectures, our models have made use of normality assumptions. From a computational point of view, the reason for this assumption is clear: combined
More informationSpatial Modeling and Prediction of County-Level Employment Growth Data
Spatial Modeling and Prediction of County-Level Employment Growth Data N. Ganesh Abstract For correlated sample survey estimates, a linear model with covariance matrix in which small areas are grouped
More informationParametric Models. Dr. Shuang LIANG. School of Software Engineering TongJi University Fall, 2012
Parametric Models Dr. Shuang LIANG School of Software Engineering TongJi University Fall, 2012 Today s Topics Maximum Likelihood Estimation Bayesian Density Estimation Today s Topics Maximum Likelihood
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 informationBasics of Point-Referenced Data Models
Basics of Point-Referenced Data Models Basic tool is a spatial process, {Y (s), s D}, where D R r Chapter 2: Basics of Point-Referenced Data Models p. 1/45 Basics of Point-Referenced Data Models Basic
More informationPracticum : Spatial Regression
: Alexandra M. Schmidt Instituto de Matemática UFRJ - www.dme.ufrj.br/ alex 2014 Búzios, RJ, www.dme.ufrj.br Exploratory (Spatial) Data Analysis 1. Non-spatial summaries Numerical summaries: Mean, median,
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 informationAn Introduction to Bayesian Linear Regression
An Introduction to Bayesian Linear Regression APPM 5720: Bayesian Computation Fall 2018 A SIMPLE LINEAR MODEL Suppose that we observe explanatory variables x 1, x 2,..., x n and dependent variables y 1,
More informationRiemann Manifold Methods in Bayesian Statistics
Ricardo Ehlers ehlers@icmc.usp.br Applied Maths and Stats University of São Paulo, Brazil Working Group in Statistical Learning University College Dublin September 2015 Bayesian inference is based on Bayes
More informationBAYESIAN KRIGING AND BAYESIAN NETWORK DESIGN
BAYESIAN KRIGING AND BAYESIAN NETWORK DESIGN Richard L. Smith Department of Statistics and Operations Research University of North Carolina Chapel Hill, N.C., U.S.A. J. Stuart Hunter Lecture TIES 2004
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 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 Linear Regression
Bayesian Linear Regression Sudipto Banerjee 1 Biostatistics, School of Public Health, University of Minnesota, Minneapolis, Minnesota, U.S.A. September 15, 2010 1 Linear regression models: a Bayesian perspective
More informationTo Estimate or Not to Estimate?
To Estimate or Not to Estimate? Benjamin Kedem and Shihua Wen In linear regression there are examples where some of the coefficients are known but are estimated anyway for various reasons not least of
More information1 Isotropic Covariance Functions
1 Isotropic Covariance Functions Let {Z(s)} be a Gaussian process on, ie, a collection of jointly normal random variables Z(s) associated with n-dimensional locations s The joint distribution of {Z(s)}
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 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 informationHierarchical Modelling for Univariate Spatial Data
Spatial omain Hierarchical Modelling for Univariate Spatial ata Sudipto Banerjee 1 and Andrew O. Finley 2 1 Biostatistics, School of Public Health, University of Minnesota, Minneapolis, Minnesota, U.S.A.
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 informationREGRESSION WITH SPATIALLY MISALIGNED DATA. Lisa Madsen Oregon State University David Ruppert Cornell University
REGRESSION ITH SPATIALL MISALIGNED DATA Lisa Madsen Oregon State University David Ruppert Cornell University SPATIALL MISALIGNED DATA 10 X X X X X X X X 5 X X X X X 0 X 0 5 10 OUTLINE 1. Introduction 2.
More informationModelling geoadditive survival data
Modelling geoadditive survival data Thomas Kneib & Ludwig Fahrmeir Department of Statistics, Ludwig-Maximilians-University Munich 1. Leukemia survival data 2. Structured hazard regression 3. Mixed model
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 informationClassical and Bayesian inference
Classical and Bayesian inference AMS 132 Claudia Wehrhahn (UCSC) Classical and Bayesian inference January 8 1 / 11 The Prior Distribution Definition Suppose that one has a statistical model with parameter
More information1 Mixed effect models and longitudinal data analysis
1 Mixed effect models and longitudinal data analysis Mixed effects models provide a flexible approach to any situation where data have a grouping structure which introduces some kind of correlation between
More informationClassical and Bayesian inference
Classical and Bayesian inference AMS 132 January 18, 2018 Claudia Wehrhahn (UCSC) Classical and Bayesian inference January 18, 2018 1 / 9 Sampling from a Bernoulli Distribution Theorem (Beta-Bernoulli
More informationBayesian Machine Learning
Bayesian Machine Learning Andrew Gordon Wilson ORIE 6741 Lecture 2: Bayesian Basics https://people.orie.cornell.edu/andrew/orie6741 Cornell University August 25, 2016 1 / 17 Canonical Machine Learning
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 informationDefault Priors and Effcient Posterior Computation in Bayesian
Default Priors and Effcient Posterior Computation in Bayesian Factor Analysis January 16, 2010 Presented by Eric Wang, Duke University Background and Motivation A Brief Review of Parameter Expansion Literature
More informationModeling daily precipitation in Space and Time
Space and Time SWGen - Hydro Berlin 20 September 2017 temporal - dependence Outline temporal - dependence temporal - dependence Stochastic Weather Generator Stochastic Weather Generator (SWG) is a stochastic
More informationThe linear model is the most fundamental of all serious statistical models encompassing:
Linear Regression Models: A Bayesian perspective Ingredients of a linear model include an n 1 response vector y = (y 1,..., y n ) T and an n p design matrix (e.g. including regressors) X = [x 1,..., x
More informationKernels for Automatic Pattern Discovery and Extrapolation
Kernels for Automatic Pattern Discovery and Extrapolation Andrew Gordon Wilson agw38@cam.ac.uk mlg.eng.cam.ac.uk/andrew University of Cambridge Joint work with Ryan Adams (Harvard) 1 / 21 Pattern Recognition
More informationLECTURE 5 NOTES. n t. t Γ(a)Γ(b) pt+a 1 (1 p) n t+b 1. The marginal density of t is. Γ(t + a)γ(n t + b) Γ(n + a + b)
LECTURE 5 NOTES 1. Bayesian point estimators. In the conventional (frequentist) approach to statistical inference, the parameter θ Θ is considered a fixed quantity. In the Bayesian approach, it is considered
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 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 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 informationBayesian Regression Linear and Logistic Regression
When we want more than point estimates Bayesian Regression Linear and Logistic Regression Nicole Beckage Ordinary Least Squares Regression and Lasso Regression return only point estimates But what if we
More informationGauge Plots. Gauge Plots JAPANESE BEETLE DATA MAXIMUM LIKELIHOOD FOR SPATIALLY CORRELATED DISCRETE DATA JAPANESE BEETLE DATA
JAPANESE BEETLE DATA 6 MAXIMUM LIKELIHOOD FOR SPATIALLY CORRELATED DISCRETE DATA Gauge Plots TuscaroraLisa Central Madsen Fairways, 996 January 9, 7 Grubs Adult Activity Grub Counts 6 8 Organic Matter
More informationHuman Pose Tracking I: Basics. David Fleet University of Toronto
Human Pose Tracking I: Basics David Fleet University of Toronto CIFAR Summer School, 2009 Looking at People Challenges: Complex pose / motion People have many degrees of freedom, comprising an articulated
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 informationThe Bayesian approach to inverse problems
The Bayesian approach to inverse problems Youssef Marzouk Department of Aeronautics and Astronautics Center for Computational Engineering Massachusetts Institute of Technology ymarz@mit.edu, http://uqgroup.mit.edu
More informationIntroduction to Spatial Data and Models
Introduction to Spatial Data and Models Researchers in diverse areas such as climatology, ecology, environmental health, and real estate marketing are increasingly faced with the task of analyzing data
More informationSpatial Dynamic Factor Analysis
Spatial Dynamic Factor Analysis Esther Salazar Federal University of Rio de Janeiro Department of Statistical Methods Sixth Workshop on BAYESIAN INFERENCE IN STOCHASTIC PROCESSES Bressanone/Brixen, Italy
More informationMCMC algorithms for fitting Bayesian models
MCMC algorithms for fitting Bayesian models p. 1/1 MCMC algorithms for fitting Bayesian models Sudipto Banerjee sudiptob@biostat.umn.edu University of Minnesota MCMC algorithms for fitting Bayesian models
More informationIntroduction to Systems Analysis and Decision Making Prepared by: Jakub Tomczak
Introduction to Systems Analysis and Decision Making Prepared by: Jakub Tomczak 1 Introduction. Random variables During the course we are interested in reasoning about considered phenomenon. In other words,
More informationConsistent high-dimensional Bayesian variable selection via penalized credible regions
Consistent high-dimensional Bayesian variable selection via penalized credible regions Howard Bondell bondell@stat.ncsu.edu Joint work with Brian Reich Howard Bondell p. 1 Outline High-Dimensional Variable
More informationCopula-basierte räumliche Interpolation
Copula-basierte räumliche Interpolation Jürgen Pilz 1 und Hannes Kazianka 1 1 Alpen-Adria-University Klagenfurt, Department of Statistics, Austria Herbsttagung R in Ausbildung, Forschung und Anwendung
More informationBasic math for biology
Basic math for biology Lei Li Florida State University, Feb 6, 2002 The EM algorithm: setup Parametric models: {P θ }. Data: full data (Y, X); partial data Y. Missing data: X. Likelihood and maximum likelihood
More information1 Data Arrays and Decompositions
1 Data Arrays and Decompositions 1.1 Variance Matrices and Eigenstructure Consider a p p positive definite and symmetric matrix V - a model parameter or a sample variance matrix. The eigenstructure is
More informationEco517 Fall 2004 C. Sims MIDTERM EXAM
Eco517 Fall 2004 C. Sims MIDTERM EXAM Answer all four questions. Each is worth 23 points. Do not devote disproportionate time to any one question unless you have answered all the others. (1) We are considering
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 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 informationLinear Models A linear model is defined by the expression
Linear Models A linear model is defined by the expression x = F β + ɛ. where x = (x 1, x 2,..., x n ) is vector of size n usually known as the response vector. β = (β 1, β 2,..., β p ) is the transpose
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 informationHypothesis Testing. Econ 690. Purdue University. Justin L. Tobias (Purdue) Testing 1 / 33
Hypothesis Testing Econ 690 Purdue University Justin L. Tobias (Purdue) Testing 1 / 33 Outline 1 Basic Testing Framework 2 Testing with HPD intervals 3 Example 4 Savage Dickey Density Ratio 5 Bartlett
More informationUniversity of Oxford. Statistical Methods Autocorrelation. Identification and Estimation
University of Oxford Statistical Methods Autocorrelation Identification and Estimation Dr. Órlaith Burke Michaelmas Term, 2011 Department of Statistics, 1 South Parks Road, Oxford OX1 3TG Contents 1 Model
More informationBayesian inference. Rasmus Waagepetersen Department of Mathematics Aalborg University Denmark. April 10, 2017
Bayesian inference Rasmus Waagepetersen Department of Mathematics Aalborg University Denmark April 10, 2017 1 / 22 Outline for today A genetic example Bayes theorem Examples Priors Posterior summaries
More informationModeling Real Estate Data using Quantile Regression
Modeling Real Estate Data using Semiparametric Quantile Regression Department of Statistics University of Innsbruck September 9th, 2011 Overview 1 Application: 2 3 4 Hedonic regression data for house prices
More informationBayesian linear regression
Bayesian linear regression Linear regression is the basis of most statistical modeling. The model is Y i = X T i β + ε i, where Y i is the continuous response X i = (X i1,..., X ip ) T is the corresponding
More informationASSESSING A VECTOR PARAMETER
SUMMARY ASSESSING A VECTOR PARAMETER By D.A.S. Fraser and N. Reid Department of Statistics, University of Toronto St. George Street, Toronto, Canada M5S 3G3 dfraser@utstat.toronto.edu Some key words. Ancillary;
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 informationChris Bishop s PRML Ch. 8: Graphical Models
Chris Bishop s PRML Ch. 8: Graphical Models January 24, 2008 Introduction Visualize the structure of a probabilistic model Design and motivate new models Insights into the model s properties, in particular
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 informationStatistics & Data Sciences: First Year Prelim Exam May 2018
Statistics & Data Sciences: First Year Prelim Exam May 2018 Instructions: 1. Do not turn this page until instructed to do so. 2. Start each new question on a new sheet of paper. 3. This is a closed book
More 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 informationProblem Selected Scores
Statistics Ph.D. Qualifying Exam: Part II November 20, 2010 Student Name: 1. Answer 8 out of 12 problems. Mark the problems you selected in the following table. Problem 1 2 3 4 5 6 7 8 9 10 11 12 Selected
More informationHow to build an automatic statistician
How to build an automatic statistician James Robert Lloyd 1, David Duvenaud 1, Roger Grosse 2, Joshua Tenenbaum 2, Zoubin Ghahramani 1 1: Department of Engineering, University of Cambridge, UK 2: Massachusetts
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 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 informationMore on nuisance parameters
BS2 Statistical Inference, Lecture 3, Hilary Term 2009 January 30, 2009 Suppose that there is a minimal sufficient statistic T = t(x ) partitioned as T = (S, C) = (s(x ), c(x )) where: C1: the distribution
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 informationPhysics 403. Segev BenZvi. Parameter Estimation, Correlations, and Error Bars. Department of Physics and Astronomy University of Rochester
Physics 403 Parameter Estimation, Correlations, and Error Bars Segev BenZvi Department of Physics and Astronomy University of Rochester Table of Contents 1 Review of Last Class Best Estimates and Reliability
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 informationPRODUCING PROBABILITY MAPS TO ASSESS RISK OF EXCEEDING CRITICAL THRESHOLD VALUE OF SOIL EC USING GEOSTATISTICAL APPROACH
PRODUCING PROBABILITY MAPS TO ASSESS RISK OF EXCEEDING CRITICAL THRESHOLD VALUE OF SOIL EC USING GEOSTATISTICAL APPROACH SURESH TRIPATHI Geostatistical Society of India Assumptions and Geostatistical Variogram
More informationOutline lecture 2 2(30)
Outline lecture 2 2(3), Lecture 2 Linear Regression it is our firm belief that an understanding of linear models is essential for understanding nonlinear ones Thomas Schön Division of Automatic Control
More informationStatistics for analyzing and modeling precipitation isotope ratios in IsoMAP
Statistics for analyzing and modeling precipitation isotope ratios in IsoMAP The IsoMAP uses the multiple linear regression and geostatistical methods to analyze isotope data Suppose the response variable
More informationSTAT 425: Introduction to Bayesian Analysis
STAT 425: Introduction to Bayesian Analysis Marina Vannucci Rice University, USA Fall 2017 Marina Vannucci (Rice University, USA) Bayesian Analysis (Part 1) Fall 2017 1 / 10 Lecture 7: Prior Types Subjective
More informationSTA 4273H: Statistical Machine Learning
STA 4273H: Statistical Machine Learning Russ Salakhutdinov Department of Statistics! rsalakhu@utstat.toronto.edu! http://www.utstat.utoronto.ca/~rsalakhu/ Sidney Smith Hall, Room 6002 Lecture 3 Linear
More informationStochastic Spectral Approaches to Bayesian Inference
Stochastic Spectral Approaches to Bayesian Inference Prof. Nathan L. Gibson Department of Mathematics Applied Mathematics and Computation Seminar March 4, 2011 Prof. Gibson (OSU) Spectral Approaches to
More information9 Multi-Model State Estimation
Technion Israel Institute of Technology, Department of Electrical Engineering Estimation and Identification in Dynamical Systems (048825) Lecture Notes, Fall 2009, Prof. N. Shimkin 9 Multi-Model State
More informationBayesian Inference. Chapter 4: Regression and Hierarchical Models
Bayesian Inference Chapter 4: Regression and Hierarchical Models Conchi Ausín and Mike Wiper Department of Statistics Universidad Carlos III de Madrid Master in Business Administration and Quantitative
More information