Gaussian Multiscale Spatio-temporal Models for Areal Data
|
|
- Emery Greer
- 5 years ago
- Views:
Transcription
1 Gaussian Multiscale Spatio-temporal Models for Areal Data (University of Missouri) Scott H. Holan (University of Missouri) Adelmo I. Bertolde (UFES)
2 Outline Motivation Multiscale factorization The multiscale spatio-temporal model Bayesian analysis Application: Agricultural Production in Espírito Santo Unknown multiscale structure Conclusions
3 Outline Motivation Multiscale factorization The multiscale spatio-temporal model Bayesian analysis Application: Agricultural Production in Espírito Santo Unknown multiscale structure Conclusions
4 Missouri: square root of standardized mortality ratio per 100, 000 inhabitants 1990 under over 33
5 Missouri: square root of standardized mortality ratio per 100, 000 inhabitants 1991 under over 33
6 Missouri: square root of standardized mortality ratio per 100, 000 inhabitants 1992 under over 33
7 Missouri: square root of standardized mortality ratio per 100, 000 inhabitants 1993 under over 33
8 Missouri: square root of standardized mortality ratio per 100, 000 inhabitants 1994 under over 33
9 Missouri: square root of standardized mortality ratio per 100, 000 inhabitants 1995 under over 33
10 Missouri: square root of standardized mortality ratio per 100, 000 inhabitants 1996 under over 33
11 Missouri: square root of standardized mortality ratio per 100, 000 inhabitants 1997 under over 33
12 Missouri: square root of standardized mortality ratio per 100, 000 inhabitants 1998 under over 33
13 Missouri: square root of standardized mortality ratio per 100, 000 inhabitants 1999 under over 33
14 Missouri: square root of standardized mortality ratio per 100, 000 inhabitants 2000 under over 33
15 Missouri: square root of standardized mortality ratio per 100, 000 inhabitants 2001 under over 33
16 Missouri: square root of standardized mortality ratio per 100, 000 inhabitants 2002 under over 33
17 Missouri: square root of standardized mortality ratio per 100, 000 inhabitants 2003 under over 33
18 Missouri: square root of standardized mortality ratio per 100, 000 inhabitants 2004 under over 33
19 Missouri: square root of standardized mortality ratio per 100, 000 inhabitants 2005 under over 33
20 Missouri: square root of standardized mortality ratio per 100, 000 inhabitants 2006 under over 33
21 Some background Many processes of interest are naturally spatio-temporal. Frequently, quantities related to these processes are available as areal data. These processes may often be considered at several different levels of spatial resolution. Related work on dynamic spatio-temporal multiscale modeling: Berliner, Wikle and Milliff (1999), Johannesson, Cressie and Huang (2007).
22 Data Structure Here, the region of interest is divided in geographic subregions or blocks, and the data may be averages or sums over these subregions. Moreover, we assume the existence of a hierarchical multiscale structure. For example, each state in Brazil is divided into counties, microregions and macroregions.
23 Geopolitical organization (a) (b) (c) Figure: Geopolitical organization of Espírito Santo State by (a) counties, (b) microregions, and (c) macroregions.
24 Outline Motivation Multiscale factorization The multiscale spatio-temporal model Bayesian analysis Application: Agricultural Production in Espírito Santo Unknown multiscale structure Conclusions
25 Multiscale factorization At each time point we decompose the data into empirical multiscale coefficients using the spatial multiscale modeling framework of Kolaczyk and Huang (2001). See also Chapter 9 of Ferreira and Lee (2007). Interest lies in agricultural production observed at the county level, which we assume is the L th level of resolution (i.e. the finest level of resolution), on a partition of a domain S R 2. For the j th county, let y Lj, µ Lj = E(y Lj ), and σ 2 Lj = V (y Lj) respectively denote agricultural production, its latent expected value and variance.
26 Let D lj be the set of descendants of subregion (l, j). The aggregated measurements at the l th level of resolution are recursively defined by y lj = p l+1,j y l+1,j. (l+1,j ) D lj Analogously, the aggregated mean process is defined by µ lj = p l+1,j µ l+1,j. (l+1,j ) D lj Assuming conditional independence, σlj 2 = pl+1,j 2 σ2 l+1,j. (l+1,j ) D lj
27 Define σ 2 l = (σ 2 l1,..., σ2 ln l ), ρ l = p l σ 2 l, and Σ l = diag(σ 2 l ), where denotes the Hadamard product. Then with and y Dlj ylj, µ L, σ 2 L N(ν ljy lj + θ lj, Ω lj ), ν lj = ρ Dlj /σ 2 lj, θ lj = µ Dlj ν lj µ lj, Ω lj = Σ Dlj σ 2 lj ρ Dlj ρ D lj.
28 Consider θ e lj = y D lj ν lj y lj, which is an empirical estimate of θ lj. Then (Kolaczyk and Huang, 2001) θ e lj y lj, µ L, σ 2 L N(θ lj, Ω lj ), which is a singular Gaussian distribution (Anderson, 1984).
29 Exploratory Multiscale Data Analysis Macroregion 1 total Disaggregated by microregion Microregion 1 Microregion 2 Microregion 3 Microregion 4 Microregion year year
30 Empirical multiscale coefficient for Macroregion year year year θ e t111 θ e t112 θ e t year year θ e t114 θ e t115
31 Outline Motivation Multiscale factorization The multiscale spatio-temporal model Bayesian analysis Application: Agricultural Production in Espírito Santo Unknown multiscale structure Conclusions
32 The multiscale spatio-temporal model Observation equation: y tl = µ tl + v tl, v tl N(0, Σ L ) where Σ L = diag(σl1 2,..., σ2 Ln L ). Multiscale decomposition of the observation equation: y t1k µ t1k N(µ t1k, σ1k 2 ) θ e tlj θ tlj N(θ tlj, Ω lj )
33 System equations: µ t1k = µ t 1,1k + w t1k, w t1k N(0, ξ k σ 2 1k ) θ tlj = θ t 1,lj + ω tlj, ω tlj N(0, ψ lj Ω lj )
34 Priors µ 01k D 0 N(m 01k, c 01k σ 2 1k ), θ 0lj D 0 N(m 0lj, C 0lj Ω lj ), ξ k IG(0.5τ k, 0.5κ k ), ψ lj IG(0.5ϱ lj, 0.5ς lj ),
35 Outline Motivation Multiscale factorization The multiscale spatio-temporal model Bayesian analysis Application: Agricultural Production in Espírito Santo Unknown multiscale structure Conclusions
36 Empirical Bayes estimation of ν lj and Ω lj ν lj : vector of relative volatilities of the descendants of (l, j), Ω lj : singular covariance matrix of the empirical multiscale coefficient of subregion (l, j) In order to obtain an initial estimate of σlj 2, we perform a univariate time series analysis for each county using first-order dynamic linear models (West and Harrison, 1997). These analyses yield estimates σ Lj 2. Let ρ l = p l σ 2 l. We estimate ν lj and Ω lj by ν lj = ρ Dlj / σ 2 lj, Ω lj = Σ Dlj σ 2 lj ρ Dlj ρ D lj.
37 Posterior exploration Let θ lj = (θ 0lj,..., θ Tlj ), θ t j = (θ t1j,..., θ tlj ), θ = (θ 11,..., θ 1n 1, θ 21,..., θ 2n 2,..., θ L1,..., θ Ln L ), with analogous definitions for the other quantities in the model. It can be shown that, given σ 2, ξ, and ψ, the vectors µ 11,..., µ 1n1, θ 11,..., θ 1n1,..., θ L1,..., θ LnL, are conditionally independent a posteriori.
38 Gibbs sampler µ 1k : Forward Filter Backward Sampler (FFBS) (Carter and Kohn, 1994; Fruhwirth-Schnatter,1994). ξ k µ 1k, σ 2 1k, D T IG (0.5τ k, 0.5κ k ), where τ k = τ k + T and κ k = κ k + σ 2 1k T t=1 (µ t1k µ t 1,1k ) 2. ψ lj θ lj, D T IG(0.5ϱ lj, 0.5ς lj ), where ϱ lj = ϱ lj + T (m lj 1) and ςlj = ς lj + T t=1 (θ tlj θ t 1,lj ) Ω lj (θ tlj θ t 1,lj ), where Ω lj is a generalized inverse of Ω lj. θ lj : Singular FFBS.
39 Singular FFBS 1. Use the singular forward filter to obtain the mean and covariance matrix of f (θ 1lj σ 2, ψ lj, D 1 ),..., f (θ Tlj σ 2, ψ lj, D T ): posterior at t 1: θ t 1,lj D t 1 N (m t 1,lj, C t 1,lj Ω lj ) ; prior at t: θtlj D t 1 N (a tlj, R tlj Ω lj ), where a tlj = m t 1,lj and R tlj = C t 1,lj + ψ lj ; posterior at t: θtlj D t N (m tlj, C tlj ( Ω lj ), where ) C tlj = (1 + R 1 tlj ) 1 and m tlj = C tlj θ e tlj + R 1 tlj a tlj. 2. Simulate θ Tlj from θ Tlj σ 2, ψ lj, D T N(m Tlj, C Tlj Ω lj ). 3. Recursively simulate θ tlj, t = T 1,..., 0, from θ tlj θ t+1,lj,..., θ Tlj, D T θ tlj θ t+1,lj, D t N(h tlj, H tlj Ω lj ), ( ) 1 where H tlj = C 1 tlj + ψ 1 lj and ( ) h tlj = H tlj C 1 tlj m tlj + ψ 1 lj θ t+1,lj.
40 Reconstruction of the latent mean process One of the main interests of any multiscale analysis is the estimation of the latent mean process at each scale of resolution. From the g th draw from the posterior distribution, we can recursively compute the corresponding latent mean process at each level of resolution using the equation µ (g) t,d lj = θ (g) tlj + ν tlj µ (g) tlj, proceeding from the coarsest to the finest resolution level. With these draws, we can then compute the posterior mean, standard deviation and credible intervals for the latent mean process.
41 Outline Motivation Multiscale factorization The multiscale spatio-temporal model Bayesian analysis Application: Agricultural Production in Espírito Santo Unknown multiscale structure Conclusions
42 Posterior densities Density Density ξ 1 ξ 2 ξ 3 ξ 4 Density ψ 11 ψ 12 ψ 13 ψ σ 2 ξ 1,..., ξ 4 ψ 11,..., ψ 14
43 Multiscale coefficient for Macroregion year year year θ t111 θ t112 θ t year year θ t114 θ t115
44 Outline Motivation Multiscale factorization The multiscale spatio-temporal model Bayesian analysis Application: Agricultural Production in Espírito Santo Unknown multiscale structure Conclusions
45 Discrepancy measurement model Consider two subregions at the same level of resolution with observations y t1 and y t2 at time t, observational variances σ 2 1 and σ 2 2, and aggregation weights p 1 and p 2. We define the Discrepancy Measurement Model (DMM) by the dynamic linear model ( yt1 y t2 ( µt t ) ) = = ( p1 σ 2 1 /(p2 1 σ2 1 + p2 2 σ2 2 ) p 2 p 2 σ2 2/(p2 1 σ2 1 + p2 2 σ2 2 ) p 1 ) + w t, ( µt 1 t 1 ) ( µt t ) + v t, where v t N(0, V), w t N(0, W), V = diag(σ1 2, σ2 2 ), and W = diag(w µ, W ).
46 Relative discrepancy Consider the DMM model and two subregions at the same resolution level. We define the relative discrepancy between these two subregions as W /W µ.
47 Multiscale clustering algorithm 1. For each county, create a link to one of its neighbors which minimizes the estimated relative discrepancy. 2. Create intermediate-level clusters of counties. Counties are in the same cluster if and only if there is a path connecting them. 3. Obtain the data at the intermediate-level through aggregation. 4. Apply steps 1 and 2 to the intermediate-level clusters in order to obtain the coarse-level clusters. We illustrate the algorithm with a dataset on mortality in Missouri.
48 Estimated multiscale structure for the State of Missouri Kansas City Columbia Jefferson City Saint Louis Springfield
49 Missouri: square root of standardized mortality ratio per 100, 000 inhabitants 1990 Observed Fitted under over 33
50 Missouri: square root of standardized mortality ratio per 100, 000 inhabitants 1991 Observed Fitted under over 33
51 Missouri: square root of standardized mortality ratio per 100, 000 inhabitants 1992 Observed Fitted under over 33
52 Missouri: square root of standardized mortality ratio per 100, 000 inhabitants 1993 Observed Fitted under over 33
53 Missouri: square root of standardized mortality ratio per 100, 000 inhabitants 1994 Observed Fitted under over 33
54 Missouri: square root of standardized mortality ratio per 100, 000 inhabitants 1995 Observed Fitted under over 33
55 Missouri: square root of standardized mortality ratio per 100, 000 inhabitants 1996 Observed Fitted under over 33
56 Missouri: square root of standardized mortality ratio per 100, 000 inhabitants 1997 Observed Fitted under over 33
57 Missouri: square root of standardized mortality ratio per 100, 000 inhabitants 1998 Observed Fitted under over 33
58 Missouri: square root of standardized mortality ratio per 100, 000 inhabitants 1999 Observed Fitted under over 33
59 Missouri: square root of standardized mortality ratio per 100, 000 inhabitants 2000 Observed Fitted under over 33
60 Missouri: square root of standardized mortality ratio per 100, 000 inhabitants 2001 Observed Fitted under over 33
61 Missouri: square root of standardized mortality ratio per 100, 000 inhabitants 2002 Observed Fitted under over 33
62 Missouri: square root of standardized mortality ratio per 100, 000 inhabitants 2003 Observed Fitted under over 33
63 Missouri: square root of standardized mortality ratio per 100, 000 inhabitants 2004 Observed Fitted under over 33
64 Missouri: square root of standardized mortality ratio per 100, 000 inhabitants 2005 Observed Fitted under over 33
65 Missouri: square root of standardized mortality ratio per 100, 000 inhabitants 2006 Observed Fitted under over 33
66 State of Missouri application. Residuals check. Fitted Sample Quantiles Observed Theoretical Quantiles
67 Outline Motivation Multiscale factorization The multiscale spatio-temporal model Bayesian analysis Application: Agricultural Production in Espírito Santo Unknown multiscale structure Conclusions
68 Conclusions New multiscale spatio-temporal model for areal data. Estimated multiscale coefficients shed light on similarities and differences between regions within each scale of resolution. Modeling strategy naturally respects nonsmooth transitions between geographic subregions. Our divide-and-conquer modeling strategy leads to computational procedures that are scalable and fast.
Lecture 4: Dynamic models
linear s Lecture 4: s Hedibert Freitas Lopes The University of Chicago Booth School of Business 5807 South Woodlawn Avenue, Chicago, IL 60637 http://faculty.chicagobooth.edu/hedibert.lopes hlopes@chicagobooth.edu
More informationSpatio-Temporal Models for Areal Data
Spatio-Temporal Models for Areal Data Juan C. Vivar (jcvivar@dme.ufrj.br) and Marco A. R. Ferreira (marco@im.ufrj.br) Departamento de Métodos Estatísticos - IM Universidade Federal do Rio de Janeiro (UFRJ)
More informationLecture 14 Bayesian Models for Spatio-Temporal Data
Lecture 14 Bayesian Models for Spatio-Temporal Data Dennis Sun Stats 253 August 13, 2014 Outline of Lecture 1 Recap of Bayesian Models 2 Empirical Bayes 3 Case 1: Long-Lead Forecasting of Sea Surface Temperatures
More informationOnline appendix to On the stability of the excess sensitivity of aggregate consumption growth in the US
Online appendix to On the stability of the excess sensitivity of aggregate consumption growth in the US Gerdie Everaert 1, Lorenzo Pozzi 2, and Ruben Schoonackers 3 1 Ghent University & SHERPPA 2 Erasmus
More informationBayesian Model Comparison:
Bayesian Model Comparison: Modeling Petrobrás log-returns Hedibert Freitas Lopes February 2014 Log price: y t = log p t Time span: 12/29/2000-12/31/2013 (n = 3268 days) LOG PRICE 1 2 3 4 0 500 1000 1500
More informationBusiness Cycle Comovements in Industrial Subsectors
Business Cycle Comovements in Industrial Subsectors Michael T. Owyang 1 Daniel Soques 2 1 Federal Reserve Bank of St. Louis 2 University of North Carolina Wilmington The views expressed here are those
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 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 informationA Complete Spatial Downscaler
A Complete Spatial Downscaler Yen-Ning Huang, Brian J Reich, Montserrat Fuentes 1 Sankar Arumugam 2 1 Department of Statistics, NC State University 2 Department of Civil, Construction, and Environmental
More informationRegionalization of Multiscale Spatial Processes Using a Criterion for Spatial Aggregation Error. Christopher K. Wikle
Regionalization of Multiscale Spatial Processes Using a Criterion for Spatial Aggregation Error Christopher K. Wikle University of Missouri TAMU Spatial Workshop, January 2015 Collaborators: Jonathan Bradley
More informationBayesian Monte Carlo Filtering for Stochastic Volatility Models
Bayesian Monte Carlo Filtering for Stochastic Volatility Models Roberto Casarin CEREMADE University Paris IX (Dauphine) and Dept. of Economics University Ca Foscari, Venice Abstract Modelling of the financial
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 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 informationEstimating marginal likelihoods from the posterior draws through a geometric identity
Estimating marginal likelihoods from the posterior draws through a geometric identity Johannes Reichl Energy Institute at the Johannes Kepler University Linz E-mail for correspondence: reichl@energieinstitut-linz.at
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 informationLecture 13 Fundamentals of Bayesian Inference
Lecture 13 Fundamentals of Bayesian Inference Dennis Sun Stats 253 August 11, 2014 Outline of Lecture 1 Bayesian Models 2 Modeling Correlations Using Bayes 3 The Universal Algorithm 4 BUGS 5 Wrapping Up
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 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 informationThe Wishart distribution Scaled Wishart. Wishart Priors. Patrick Breheny. March 28. Patrick Breheny BST 701: Bayesian Modeling in Biostatistics 1/11
Wishart Priors Patrick Breheny March 28 Patrick Breheny BST 701: Bayesian Modeling in Biostatistics 1/11 Introduction When more than two coefficients vary, it becomes difficult to directly model each element
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 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 informationThe Effects of Monetary Policy on Stock Market Bubbles: Some Evidence
The Effects of Monetary Policy on Stock Market Bubbles: Some Evidence Jordi Gali Luca Gambetti ONLINE APPENDIX The appendix describes the estimation of the time-varying coefficients VAR model. The model
More informationBayesian Linear Models
Bayesian Linear 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 & Department
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 informationAreal data models. Spatial smoothers. Brook s Lemma and Gibbs distribution. CAR models Gaussian case Non-Gaussian case
Areal data models Spatial smoothers Brook s Lemma and Gibbs distribution CAR models Gaussian case Non-Gaussian case SAR models Gaussian case Non-Gaussian case CAR vs. SAR STAR models Inference for areal
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 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 informationOverall Objective Priors
Overall Objective Priors Jim Berger, Jose Bernardo and Dongchu Sun Duke University, University of Valencia and University of Missouri Recent advances in statistical inference: theory and case studies University
More informationRate Maps and Smoothing
Rate Maps and Smoothing Luc Anselin Spatial Analysis Laboratory Dept. Agricultural and Consumer Economics University of Illinois, Urbana-Champaign http://sal.agecon.uiuc.edu Outline Mapping Rates Risk
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 informationPartial factor modeling: predictor-dependent shrinkage for linear regression
modeling: predictor-dependent shrinkage for linear Richard Hahn, Carlos Carvalho and Sayan Mukherjee JASA 2013 Review by Esther Salazar Duke University December, 2013 Factor framework The factor framework
More informationParticle Learning and Smoothing
Particle Learning and Smoothing Hedibert Freitas Lopes The University of Chicago Booth School of Business September 18th 2009 Instituto Nacional de Pesquisas Espaciais São José dos Campos, Brazil Joint
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 informationSmall Area Estimation via Multivariate Fay-Herriot Models with Latent Spatial Dependence
Small Area Estimation via Multivariate Fay-Herriot Models with Latent Spatial Dependence Aaron T. Porter 1, Christopher K. Wikle 2, Scott H. Holan 2 arxiv:1310.7211v1 [stat.me] 27 Oct 2013 Abstract The
More informationA Divide-and-Conquer Bayesian Approach to Large-Scale Kriging
A Divide-and-Conquer Bayesian Approach to Large-Scale Kriging Cheng Li DSAP, National University of Singapore Joint work with Rajarshi Guhaniyogi (UC Santa Cruz), Terrance D. Savitsky (US Bureau of Labor
More informationModeling conditional distributions with mixture models: Theory and Inference
Modeling conditional distributions with mixture models: Theory and Inference John Geweke University of Iowa, USA Journal of Applied Econometrics Invited Lecture Università di Venezia Italia June 2, 2005
More informationResearch Division Federal Reserve Bank of St. Louis Working Paper Series
Research Division Federal Reserve Bank of St Louis Working Paper Series Kalman Filtering with Truncated Normal State Variables for Bayesian Estimation of Macroeconomic Models Michael Dueker Working Paper
More informationSteven L. Scott. Presented by Ahmet Engin Ural
Steven L. Scott Presented by Ahmet Engin Ural Overview of HMM Evaluating likelihoods The Likelihood Recursion The Forward-Backward Recursion Sampling HMM DG and FB samplers Autocovariance of samplers Some
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 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 informationHierarchical Modeling for Spatio-temporal Data
Hierarchical Modeling for Spatio-temporal 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 of
More informationScaling Neighbourhood Methods
Quick Recap Scaling Neighbourhood Methods Collaborative Filtering m = #items n = #users Complexity : m * m * n Comparative Scale of Signals ~50 M users ~25 M items Explicit Ratings ~ O(1M) (1 per billion)
More informationCore Inflation and Trend Inflation. Appendix
Core Inflation and Trend Inflation Appendix June 2015 (Revised November 2015) James H. Stock Department of Economics, Harvard University and the National Bureau of Economic Research and Mark W. Watson
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 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 informationNORGES TEKNISK-NATURVITENSKAPELIGE UNIVERSITET. Fitting dynamic models using integrated nested Laplace approximations INLA
NORGES TEKNISK-NATURVITENSKAPELIGE UNIVERSITET Fitting dynamic models using integrated nested Laplace approximations INLA by Ramiro Ruiz-Cárdenas, Elias T. Krainski and Håvard Rue PREPRINT STATISTICS NO.
More informationApproaches for Multiple Disease Mapping: MCAR and SANOVA
Approaches for Multiple Disease Mapping: MCAR and SANOVA Dipankar Bandyopadhyay Division of Biostatistics, University of Minnesota SPH April 22, 2015 1 Adapted from Sudipto Banerjee s notes SANOVA vs MCAR
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 informationAlexandra M. Schmidt and Hedibert F. Lopes. Dynamic models
Alexandra M. Schmidt and Hedibert F. Lopes Dynamic models ii Contents 1 Dynamic models 1 Alexandra M. Schmidt and Hedibert F. Lopes 1.1 Introduction........................... 1 1.2 Univariate Normal Dynamic
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 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 informationBayesian Modeling of Conditional Distributions
Bayesian Modeling of Conditional Distributions John Geweke University of Iowa Indiana University Department of Economics February 27, 2007 Outline Motivation Model description Methods of inference Earnings
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 informationApproximate Bayesian computation for spatial extremes via open-faced sandwich adjustment
Approximate Bayesian computation for spatial extremes via open-faced sandwich adjustment Ben Shaby SAMSI August 3, 2010 Ben Shaby (SAMSI) OFS adjustment August 3, 2010 1 / 29 Outline 1 Introduction 2 Spatial
More 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 Advanced Statistics and Data Mining Summer School
More informationBayesian Hierarchical Models
Bayesian Hierarchical Models Gavin Shaddick, Millie Green, Matthew Thomas University of Bath 6 th - 9 th December 2016 1/ 34 APPLICATIONS OF BAYESIAN HIERARCHICAL MODELS 2/ 34 OUTLINE Spatial epidemiology
More informationBayesian Estimation of DSGE Models
Bayesian Estimation of DSGE Models Stéphane Adjemian Université du Maine, GAINS & CEPREMAP stephane.adjemian@univ-lemans.fr http://www.dynare.org/stepan June 28, 2011 June 28, 2011 Université du Maine,
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 informationGibbs Sampling in Linear Models #2
Gibbs Sampling in Linear Models #2 Econ 690 Purdue University Outline 1 Linear Regression Model with a Changepoint Example with Temperature Data 2 The Seemingly Unrelated Regressions Model 3 Gibbs sampling
More informationAppendix for Monetary Policy Shifts and the Term Structure
Appendix for Monetary Policy Shifts and the Term Structure Andrew Ang Columbia University and NBER Jean Boivin Bank of Canada, CIRPÉE, CIRANO and NBER Sen Dong Ortus Capital Management Rudy Loo-Kung Columbia
More informationState-space Model. Eduardo Rossi University of Pavia. November Rossi State-space Model Fin. Econometrics / 53
State-space Model Eduardo Rossi University of Pavia November 2014 Rossi State-space Model Fin. Econometrics - 2014 1 / 53 Outline 1 Motivation 2 Introduction 3 The Kalman filter 4 Forecast errors 5 State
More informationWavelet-Based Nonparametric Modeling of Hierarchical Functions in Colon Carcinogenesis
Wavelet-Based Nonparametric Modeling of Hierarchical Functions in Colon Carcinogenesis Jeffrey S. Morris University of Texas, MD Anderson Cancer Center Joint wor with Marina Vannucci, Philip J. Brown,
More informationPart I State space models
Part I State space models 1 Introduction to state space time series analysis James Durbin Department of Statistics, London School of Economics and Political Science Abstract The paper presents a broad
More informationChapter 3: Maximum-Likelihood & Bayesian Parameter Estimation (part 1)
HW 1 due today Parameter Estimation Biometrics CSE 190 Lecture 7 Today s lecture was on the blackboard. These slides are an alternative presentation of the material. CSE190, Winter10 CSE190, Winter10 Chapter
More informationSpatially Varying Autoregressive Processes
Spatially Varying Autoregressive Processes Aline A. Nobre, Bruno Sansó and Alexandra M. Schmidt Abstract We develop a class of models for processes indexed in time and space that are based on autoregressive
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 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 informationBayesian Nonparametric Autoregressive Models via Latent Variable Representation
Bayesian Nonparametric Autoregressive Models via Latent Variable Representation Maria De Iorio Yale-NUS College Dept of Statistical Science, University College London Collaborators: Lifeng Ye (UCL, London,
More informationBayesian Linear Models
Bayesian Linear 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, School of Public
More informationA Spatio-Temporal Point Process Model for Ambulance Demand
A Spatio-Temporal Point Process Model for Ambulance Demand David S. Matteson Department of Statistical Science Department of Social Statistics Cornell University matteson@cornell.edu http://www.stat.cornell.edu/~matteson/
More informationThe Kalman filter, Nonlinear filtering, and Markov Chain Monte Carlo
NBER Summer Institute Minicourse What s New in Econometrics: Time Series Lecture 5 July 5, 2008 The Kalman filter, Nonlinear filtering, and Markov Chain Monte Carlo Lecture 5, July 2, 2008 Outline. Models
More informationVCMC: Variational Consensus Monte Carlo
VCMC: Variational Consensus Monte Carlo Maxim Rabinovich, Elaine Angelino, Michael I. Jordan Berkeley Vision and Learning Center September 22, 2015 probabilistic models! sky fog bridge water grass object
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 informationParticle Learning and Smoothing
Particle Learning and Smoothing Carlos Carvalho, Michael Johannes, Hedibert Lopes and Nicholas Polson This version: September 2009 First draft: December 2007 Abstract In this paper we develop particle
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 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 informationLocal Likelihood Bayesian Cluster Modeling for small area health data. Andrew Lawson Arnold School of Public Health University of South Carolina
Local Likelihood Bayesian Cluster Modeling for small area health data Andrew Lawson Arnold School of Public Health University of South Carolina Local Likelihood Bayesian Cluster Modelling for Small Area
More informationSpatio-Temporal Modelling of Credit Default Data
1/20 Spatio-Temporal Modelling of Credit Default Data Sathyanarayan Anand Advisor: Prof. Robert Stine The Wharton School, University of Pennsylvania April 29, 2011 2/20 Outline 1 Background 2 Conditional
More informationSequential Monte Carlo Methods (for DSGE Models)
Sequential Monte Carlo Methods (for DSGE Models) Frank Schorfheide University of Pennsylvania, PIER, CEPR, and NBER October 23, 2017 Some References These lectures use material from our joint work: Tempered
More informationFrailty Modeling for Spatially Correlated Survival Data, with Application to Infant Mortality in Minnesota By: Sudipto Banerjee, Mela. P.
Frailty Modeling for Spatially Correlated Survival Data, with Application to Infant Mortality in Minnesota By: Sudipto Banerjee, Melanie M. Wall, Bradley P. Carlin November 24, 2014 Outlines of the talk
More informationSpatio-temporal precipitation modeling based on time-varying regressions
Spatio-temporal precipitation modeling based on time-varying regressions Oleg Makhnin Department of Mathematics New Mexico Tech Socorro, NM 87801 January 19, 2007 1 Abstract: A time-varying regression
More informationA Spatio-temporal Model for Agricultural Yield Prediction. Oleksiy Tokovenko 1 Jeffrey H. Dorfman Lewell F. Gunter
A Spatio-temporal Model for Agricultural Yield Prediction Oleksiy Tokovenko 1 Jeffrey H. Dorfman Lewell F. Gunter Selected Paper prepared for presentation at the Agricultural and Applied Economics Association
More informationSPATIAL ANALYSIS & MORE
SPATIAL ANALYSIS & MORE Thomas A. Louis, PhD Department of Biostatistics Johns Hopins Bloomberg School of Public Health www.biostat.jhsph.edu/ tlouis/ tlouis@jhu.edu 1/28 Outline Prediction vs Inference
More informationMultiscale detection of localized anomalous structure in aggregate disease incidence data
Multiscale detection of localized anomalous structure in aggregate disease incidence data Mary M. Louie 1, Eric D. Kolaczyk 2 1 National Center for Health Statistics, 3311 Toledo Road, Room 3215, Hyattsville,
More informationBayesian non-parametric model to longitudinally predict churn
Bayesian non-parametric model to longitudinally predict churn Bruno Scarpa Università di Padova Conference of European Statistics Stakeholders Methodologists, Producers and Users of European Statistics
More informationBayesian Semiparametric GARCH Models
Bayesian Semiparametric GARCH Models Xibin (Bill) Zhang and Maxwell L. King Department of Econometrics and Business Statistics Faculty of Business and Economics xibin.zhang@monash.edu Quantitative Methods
More informationClimate Change: the Uncertainty of Certainty
Climate Change: the Uncertainty of Certainty Reinhard Furrer, UZH JSS, Geneva Oct. 30, 2009 Collaboration with: Stephan Sain - NCAR Reto Knutti - ETHZ Claudia Tebaldi - Climate Central Ryan Ford, Doug
More informationComputer Intensive Methods in Mathematical Statistics
Computer Intensive Methods in Mathematical Statistics Department of mathematics johawes@kth.se Lecture 16 Advanced topics in computational statistics 18 May 2017 Computer Intensive Methods (1) Plan of
More informationECE521 Tutorial 11. Topic Review. ECE521 Winter Credits to Alireza Makhzani, Alex Schwing, Rich Zemel and TAs for slides. ECE521 Tutorial 11 / 4
ECE52 Tutorial Topic Review ECE52 Winter 206 Credits to Alireza Makhzani, Alex Schwing, Rich Zemel and TAs for slides ECE52 Tutorial ECE52 Winter 206 Credits to Alireza / 4 Outline K-means, PCA 2 Bayesian
More informationBayesian Semiparametric GARCH Models
Bayesian Semiparametric GARCH Models Xibin (Bill) Zhang and Maxwell L. King Department of Econometrics and Business Statistics Faculty of Business and Economics xibin.zhang@monash.edu Quantitative Methods
More informationCross-sectional space-time modeling using ARNN(p, n) processes
Cross-sectional space-time modeling using ARNN(p, n) processes W. Polasek K. Kakamu September, 006 Abstract We suggest a new class of cross-sectional space-time models based on local AR models and nearest
More informationWeb Appendix for Hierarchical Adaptive Regression Kernels for Regression with Functional Predictors by D. B. Woodard, C. Crainiceanu, and D.
Web Appendix for Hierarchical Adaptive Regression Kernels for Regression with Functional Predictors by D. B. Woodard, C. Crainiceanu, and D. Ruppert A. EMPIRICAL ESTIMATE OF THE KERNEL MIXTURE Here we
More informationTimevarying VARs. Wouter J. Den Haan London School of Economics. c Wouter J. Den Haan
Timevarying VARs Wouter J. Den Haan London School of Economics c Wouter J. Den Haan Time-Varying VARs Gibbs-Sampler general idea probit regression application (Inverted Wishart distribution Drawing from
More informationOutline. Clustering. Capturing Unobserved Heterogeneity in the Austrian Labor Market Using Finite Mixtures of Markov Chain Models
Capturing Unobserved Heterogeneity in the Austrian Labor Market Using Finite Mixtures of Markov Chain Models Collaboration with Rudolf Winter-Ebmer, Department of Economics, Johannes Kepler University
More informationExpression Data Exploration: Association, Patterns, Factors & Regression Modelling
Expression Data Exploration: Association, Patterns, Factors & Regression Modelling Exploring gene expression data Scale factors, median chip correlation on gene subsets for crude data quality investigation
More informationImperfect Data in an Uncertain World
Imperfect Data in an Uncertain World James B. Elsner Department of Geography, Florida State University Tallahassee, Florida Corresponding author address: Dept. of Geography, Florida State University Tallahassee,
More informationTechnical Vignette 5: Understanding intrinsic Gaussian Markov random field spatial models, including intrinsic conditional autoregressive models
Technical Vignette 5: Understanding intrinsic Gaussian Markov random field spatial models, including intrinsic conditional autoregressive models Christopher Paciorek, Department of Statistics, University
More informationState-space Model. Eduardo Rossi University of Pavia. November Rossi State-space Model Financial Econometrics / 49
State-space Model Eduardo Rossi University of Pavia November 2013 Rossi State-space Model Financial Econometrics - 2013 1 / 49 Outline 1 Introduction 2 The Kalman filter 3 Forecast errors 4 State smoothing
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 Probabilistic Graphical Models: Exercises
Introduction to Probabilistic Graphical Models: Exercises Cédric Archambeau Xerox Research Centre Europe cedric.archambeau@xrce.xerox.com Pascal Bootcamp Marseille, France, July 2010 Exercise 1: basics
More information