Web Appendices: Hierarchical and Joint Site-Edge Methods for Medicare Hospice Service Region Boundary Analysis
|
|
- Lillian Horton
- 5 years ago
- Views:
Transcription
1 Web Appendices: Hierarchical and Joint Site-Edge Methods for Medicare Hospice Service Region Boundary Analysis Haijun Ma, Bradley P. Carlin and Sudipto Banerjee December 8, 2008 Web Appendix A: Selecting binding strength and penalty parameters Here, we address the proper selection of the binding strength parameter ν and the penalty parameter κ in the penalized SE-Ising prior (8). To do this, for a given map (say, St. Luke s) we obtain samples φ E g, g = 1,...,G = 3000 from the penalized SE-Ising distribution for different values of ν and κ. As a comparison, we also obtain i.i.d. draws φ E k,g iid Bernoulli(p), k = 1,...,K, g = 1,...,G = 3000, where K is number of edges in the map and p for each draw is the mean of the φ E g from the two priors should have the same number of on edges. draws from the penalized SE-Ising prior. Thus on average, samples Table 1: Comparison of average number of strings of on edges, M, and average string length, L, for penalized SE-Ising (Ising) and i.i.d. Bernoulli (IID) models over 3000 prior draws, St. Luke s map. Ê(φ E ) is the empirical mean for the penalized SE-Ising. ν κ M Ising M IID L Ising L IID Ê(φ E ) The simulation results using the St. Luke s map are given in Table 1. It can be seen that for ν = 1 and most values of κ, the penalized SE-Ising distribution performs better, delivering fewer and longer strings of on edges than the i.i.d. Bernoulli. For ν = 100 and most values of κ, the penalized SE-Ising instead performs worse, producing more and shorter strings of edges. For ν = 0.5 and ν = 1, as κ increases, the number of strings of on edges and the average length of the strings decrease, with the penalized SE-Ising performing better for more negative values of κ. For fixed ν, experimentation with a finer grid of κ values suggests a fairly abrupt 1
2 Table 2: Comparison of average number of strings of on edges, M, and average string length, L, for penalized SE-Ising (Ising) and i.i.d. Bernoulli (IID) models over G = 3000 prior draws, SMDC map. Ê(φ E ) is the empirical mean for the penalized SE-Ising. ν κ M Ising M IID L Ising L IID Ê(φ E) change at a particular point, as often seen with models of this sort (Winkler, 2003). In our setting, this point is roughly κ = 3 for ν = 0.5 and κ = 7 for ν = 1. Similar patterns are observed from the simulation results for SMDC, given in Table 2. In the remainder of our paper, we continue with the ν and κ values used in Table 1 (ν = 0.5 and κ = 3) of the manuscript, since they appear reasonably well justifed by the prior simulations in Tables 1 and 2, and since they deliver the appealing fitted wombled maps shown below. While we were able to obtain satisfactory boundaries in our application simply by experimenting with a few values for ν and κ, one might wonder if, like τ ψ (for CAR2), they could be estimated from the data. A hyperprior for (ν, κ) would likely need to feature some informative content, since data-based information here will be weak, and in any case we want to encourage the model to choose longer strings of boundary segments. Designing efficient updating schemes for (ν, κ) will not be easy; moreover, such schemes may need to involve the normalizing constant of the Ising model. This amounts to finding certain vertex-induced subgraphs (see, e.g., Lund and Yannakakis, 1993) for which no polynomial time algorithm exists. Since this difficult problem arises at every MCMC iteration, we considered only the case of fixed ν and κ. As an enhancement to our SE model, we might want the binding strength to increase if the areal level information supports connection of two edges. Leclerc (1989) recommends smoothing second differences across adjacent areas. Alternatively, note that our edge priors currently use only Gibbs distributions based on cliques of size two. While this framework is conceptually and computationally feasible, it does not permit unequivocal distinction of several important edge patterns. Consider for example four boundary segments arranged in a four-way edge intersection. CAR models with cliques of size two cannot differentiate between a continuation (north-south or east-west boundary only) and a cross (both boundaries crossing at the intersection), though in our context we would want the model to favor the former. The image processing literature (Dass and Nair, 2003; Winkler, 2003) uses cliques of size 4 to address this problem, but in our irregular lattice setting it is not at all clear how such cliques would be defined, nor how the problem could be managed computationally. 2
3 [0.5,0.6) [0.6,0.8) [0.8,0.86] (a)spatial Classification, p [0.5,0.6) [0.6,0.8) [0.8,0.98] (b)spatial Classification, p Figure 1: Posterior mean of p i for the spatial classification model: (a) St. Luke s; (b) SMDC. [0,0.02) [0.02,0.05) [0.05,0.13) [0.13,0.41] [0,0.02) [0.02,0.08) [0.08,0.21) [0.21,0.65] (a) Spatial Classification, p (b) Spatial Classification, p Figure 2: Posterior median of p,ij = p i p j for the spatial classification model: (a) St. Luke s; (b) SMDC. Web Appendix B: Comparisons with spatial classification Given the complex nature of our chosen boundary analysis models, one might wonder if a simple spatial classification model (say, a logistic regression model for the probability that each zip code is served) might work just as well. That is, suppose we let Z i = 1 if Y i > 0 (i.e., in there are any patients served by the hospice in zip code i), and Z i = 0 if Y i = 0. Suppose Z i Bernoulli(p i ), where p i is the probability that area i is served. We adopt the mean structure logit(p i ) = β 0 + β 1 X i + φ i, i = 1,...,n, where φ CAR(τ φ, W), and the single covariate x i is again the intercentroidal distance from the patient s zip to the nearest relevant hospice home base zip. Figure 1 plots the posterior means of p i, i = 1,..., n. It can be seen that the areas that are closer to the headquarters of the hospices have higher chances of being served, and these probabilities decrease as the distance from the headquarters increases. However, for both hospices, the totality of regions for which E(p i z) 0.8 or even 0.6 is much smaller than the self-reported service area. Figures 2 and 3 give the corresponding boundary plots based on p,ij = p i p j and log OR ij = log pi(1 pj) p j(1 p i) respectively. Note that neither of these two criteria appear suitable: the smoothness of the fitted surface in Figure 1 has apparently led to the identification of boundary segments throughout the map. Thus, this apparently sensible class of spatial classification models is not helpful for identifying service area boundaries. 3
4 [0,0.24) [0.24,0.6) [0.6,1.07) [1.07,2.12] (a) Spatial Classification, logor [0,0.35) [0.35,0.89) [0.89,1.49) [1.49,4.02] (b) Spatial Classification, logor Figure 3: Posterior median of log OR ij = log pi(1 pj) p j(1 p i) for the spatial classification model: (a) St. Luke s; (b) SMDC. Web Appendix C: Edge correction, thresholding and BLV definition A practical issue in almost any spatial analysis is the differential effect of the methods near the outer boundary of the spatial domain. As pointed out in Subsection 2.2, areas having more neighbors are subject to a greater degree of smoothing by the CAR model. Spatial units close to the outer boundary typically have fewer neighbors, and may thus be undersmoothed. In the case of areal boundary analysis, this problem may manifest as a tendency to find spurious boundary segments that connect to the outer boundary. This motivates a search for a spatial edge correction method. Various such methods have been proposed in the literature; see Lawson et al. (1999) for a review. For our random edge models, if an edge element touches the outer boundary, we extend its edge-level neighborhood structure to include those segments of the outer boundary as its neighbors. However, these additional outer boundary segments are not treated as random variables to be estimated, but mere placeholders in the algorithm. For SE-Ising, we set these φ E ij we are certain they are part of the boundary). parameters equal to 1 (since A second practical issue concerns the detection of spurious boundary segments within a homogenous region. We wish to find only segments that encompass a hospice s service area; minor boundaries separating different levels of service within the area are less important. To address this problem, we may use a form of thresholding, redefining the BLV in the µ-based case as µ,ij = (µ i µ j ) if (µ i c)(µ j c) < 0, and 0 otherwise. That is, the BLV is defined so that differences in the means of adjacent ZIP codes are recorded only when the means lie on opposite sides of c, where c is a predetermined cut-off value indicating a minimum service level. The BLV is thus 0 if the adjacent ZIP codes are both inside the service region, or both outside of it. In our analysis, we would like to list a ZIP code as served if it contains at least one served patient over our three-year observation period. As such, we somewhat arbitrarily selected c = 1.5 as the threshold value for our µ -based maps, meaning that a posterior median difference of at least one-half a patient per year will lead to the establishment of a boundary. We found this c to perform well in eliminating intra-service area clutter from our maps, and thus used it in the final maps presented in the main paper. However, at the urging of the associate editor, we also undertook a brief investigation of a few other values of c, namely 0.5, 1, 2, and 5. The results for c = 0.5 and 5 using the same shading levels as used in the main paper are shown in Figure 4; results for c = 1 and 2 were visually indistinguishable from those for c = 1.5. On the whole, the results for the St. Luke s map are very similar to before, though the smaller value of c does lead to some minor boundaries in the northwest portion of the map, where we have a few small but nonzero counts. Differences are somewhat more 4
5 [10,16.49] [10,19.81] (b) Penalized SE Ising(c=0.5), µ (b) Penalized SE Ising(c=5), µ [10,140.45] [10,142.05] (b) Penalized SE Ising(c=0.5), µ (b) Penalized SE Ising(c=5), µ Figure 4: Maps of St. Luke s and SMDC s service area boundaries given by the penalized SE-Ising model with varying thresholding cutoffs c: (a) St. Luke s service area boundaries, c = 0.5; (b) St. Luke s service area boundaries, c = 5; (c) SMDC s service area boundaries, c = 0.5; (d) SMDC s service area boundaries, c = 5. pronounced for the SMDC map, with the larger c causing the boundaries to contract somewhat around the more highly served ZIP code areas. Thus, at least within a reasonable range (here, 0.5 to 5), c acts as a boundary intensity parameter. While in our case c can be chosen from context (our half a patient per year requirement), it can also be determined in the way an ordinary CAR model smoothing parameter is sometimes set: by personal preference, perhaps based on context (in this case, how weak a definition of service one is willing to tolerate, and hence how large the final service area will be). Web Appendix D: Normalizing constants in the SE-Ising model We now return to the issue of unknown normalizing constants in (3) and (6). Several technical issues arise due to the singularity of D w W. Were the CAR distribution proper, the normalizing constants would involve the determinant of D w W and τ φ to some power. Viewing the CAR as a degenerate multivariate normal distribution defined on a hyperplane and comparing it to an n-dimensional multivariate normal distribution, we see that the analogous normalizing constant for (1) would be τ φ 2π (D w W) 1/2. But since (D w W) is rank deficient, D w W = 0. Besag and Higdon (1999) suggest keeping only D w W +, defined as the non-zero eigenvalue product of (D w W); see also Lu et al. (2007). Hodges et al. (2003) recommend using rank(d w W)/2 as the exponent for τ φ /(2π), and show that rank(d w W) = n I where I is the number of disconnected islands in the spatial structure. Since (D w W) is a function of φ E, both its determinant and I will need to be updated at every MCMC iteration. To stabilize the estimation of φ S i corresponding to singleton islands (islands consisting of just one areal unit) 5
6 that appear and disappear from the model, we assign each an independent N(0, 1/τ φ ) prior component. The rest of the φ S vector retains its joint improper CAR(τ φ, W) prior, parameterized using a sum-to-zero constraint (Besag and Kooperberg, 1995). As mentioned above, these problems can also be solved via the proper CAR that uses D ǫ = Diag(w i+ + ǫ) for some small ǫ > 0. This however adds the difficult problem of choosing the value of ǫ, to which the results are often rather sensitive (Lu et al., 2007). Additional references Besag, J. and Higdon, D. (1999). Bayesian analysis of agricultural field experiments (with discussion). J. Roy. Statist. Soc., Ser. B, 61, Besag, J. and Kooperberg, C. (1995). On conditional and intrinsic autoregressions. Biometrika, 82, Hodges, J.S., Carlin, B.P., and Fan, Q. (2003). On the precision of the conditionally autoregressive prior in spatial models. Biometrics, 59, Leclerc, Y.G. (1989). Constructing simple stable descriptions for image partitioning. International Journal of Computer Vision, 3, Lund, C. and Yannakakis, M. (1993). The approximation of maximum subgraph problems. In Proceedings of the 20th International Colloquium on Automata, Languages, and Programming, 700,
Bayesian Areal Wombling for Geographic Boundary Analysis
Bayesian Areal Wombling for Geographic Boundary Analysis Haolan Lu, Haijun Ma, and Bradley P. Carlin haolanl@biostat.umn.edu, haijunma@biostat.umn.edu, and brad@biostat.umn.edu Division of Biostatistics
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 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 informationAreal Unit Data Regular or Irregular Grids or Lattices Large Point-referenced Datasets
Areal Unit Data Regular or Irregular Grids or Lattices Large Point-referenced Datasets Is there spatial pattern? Chapter 3: Basics of Areal Data Models p. 1/18 Areal Unit Data Regular or Irregular Grids
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 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 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 informationBayesian spatial hierarchical modeling for temperature extremes
Bayesian spatial hierarchical modeling for temperature extremes Indriati Bisono Dr. Andrew Robinson Dr. Aloke Phatak Mathematics and Statistics Department The University of Melbourne Maths, Informatics
More informationPart 8: GLMs and Hierarchical LMs and GLMs
Part 8: GLMs and Hierarchical LMs and GLMs 1 Example: Song sparrow reproductive success Arcese et al., (1992) provide data on a sample from a population of 52 female song sparrows studied over the course
More informationSpatio-Temporal Threshold Models for Relating UV Exposures and Skin Cancer in the Central United States
Spatio-Temporal Threshold Models for Relating UV Exposures and Skin Cancer in the Central United States Laura A. Hatfield and Bradley P. Carlin Division of Biostatistics School of Public Health University
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 informationspbayes: An R Package for Univariate and Multivariate Hierarchical Point-referenced Spatial Models
spbayes: An R Package for Univariate and Multivariate Hierarchical Point-referenced Spatial Models Andrew O. Finley 1, Sudipto Banerjee 2, and Bradley P. Carlin 2 1 Michigan State University, Departments
More informationA spatio-temporal model for extreme precipitation simulated by a climate model
A spatio-temporal model for extreme precipitation simulated by a climate model Jonathan Jalbert Postdoctoral fellow at McGill University, Montréal Anne-Catherine Favre, Claude Bélisle and Jean-François
More informationHierarchical Modeling for non-gaussian Spatial Data
Hierarchical Modeling for non-gaussian 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 informationHierarchical Modelling for non-gaussian Spatial Data
Hierarchical Modelling for non-gaussian Spatial Data Sudipto Banerjee 1 and Andrew O. Finley 2 1 Department of Forestry & Department of Geography, Michigan State University, Lansing Michigan, U.S.A. 2
More informationProbabilistic Graphical Models
2016 Robert Nowak Probabilistic Graphical Models 1 Introduction We have focused mainly on linear models for signals, in particular the subspace model x = Uθ, where U is a n k matrix and θ R k is a vector
More informationHierarchical Modeling and Analysis for Spatial Data
Hierarchical Modeling and Analysis for Spatial Data Bradley P. Carlin, Sudipto Banerjee, and Alan E. Gelfand brad@biostat.umn.edu, sudiptob@biostat.umn.edu, and alan@stat.duke.edu University of Minnesota
More informationEffects of residual smoothing on the posterior of the fixed effects in disease-mapping models
Effects of residual smoothing on the posterior of the fixed effects in disease-mapping models By BRIAN J. REICH a, JAMES S. HODGES a, and VESNA ZADNIK b 1 a Division of Biostatistics, School of Public
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 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 informationHierarchical Modelling for non-gaussian Spatial Data
Hierarchical Modelling for non-gaussian Spatial Data Geography 890, Hierarchical Bayesian Models for Environmental Spatial Data Analysis February 15, 2011 1 Spatial Generalized Linear Models Often data
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 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 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 informationVariational Inference (11/04/13)
STA561: Probabilistic machine learning Variational Inference (11/04/13) Lecturer: Barbara Engelhardt Scribes: Matt Dickenson, Alireza Samany, Tracy Schifeling 1 Introduction In this lecture we will further
More informationHierarchical Modelling for Univariate and Multivariate Spatial Data
Hierarchical Modelling for Univariate and Multivariate Spatial Data p. 1/4 Hierarchical Modelling for Univariate and Multivariate Spatial Data Sudipto Banerjee sudiptob@biostat.umn.edu University of Minnesota
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 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 informationLearning discrete graphical models via generalized inverse covariance matrices
Learning discrete graphical models via generalized inverse covariance matrices Duzhe Wang, Yiming Lv, Yongjoon Kim, Young Lee Department of Statistics University of Wisconsin-Madison {dwang282, lv23, ykim676,
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 informationA fast sampler for data simulation from spatial, and other, Markov random fields
A fast sampler for data simulation from spatial, and other, Markov random fields Andee Kaplan Iowa State University ajkaplan@iastate.edu June 22, 2017 Slides available at http://bit.ly/kaplan-phd Joint
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 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 informationARIC Manuscript Proposal # PC Reviewed: _9/_25_/06 Status: A Priority: _2 SC Reviewed: _9/_25_/06 Status: A Priority: _2
ARIC Manuscript Proposal # 1186 PC Reviewed: _9/_25_/06 Status: A Priority: _2 SC Reviewed: _9/_25_/06 Status: A Priority: _2 1.a. Full Title: Comparing Methods of Incorporating Spatial Correlation in
More informationCompatibility of conditionally specified models
Compatibility of conditionally specified models Hua Yun Chen Division of epidemiology & Biostatistics School of Public Health University of Illinois at Chicago 1603 West Taylor Street, Chicago, IL 60612
More informationA Note on Bayesian Inference After Multiple Imputation
A Note on Bayesian Inference After Multiple Imputation Xiang Zhou and Jerome P. Reiter Abstract This article is aimed at practitioners who plan to use Bayesian inference on multiplyimputed datasets in
More informationBayesian Image Segmentation Using MRF s Combined with Hierarchical Prior Models
Bayesian Image Segmentation Using MRF s Combined with Hierarchical Prior Models Kohta Aoki 1 and Hiroshi Nagahashi 2 1 Interdisciplinary Graduate School of Science and Engineering, Tokyo Institute of Technology
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 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 informationVariable Selection in Structured High-dimensional Covariate Spaces
Variable Selection in Structured High-dimensional Covariate Spaces Fan Li 1 Nancy Zhang 2 1 Department of Health Care Policy Harvard University 2 Department of Statistics Stanford University May 14 2007
More informationNearest Neighbor Gaussian Processes for Large Spatial Data
Nearest Neighbor Gaussian Processes for Large Spatial Data Abhi Datta 1, Sudipto Banerjee 2 and Andrew O. Finley 3 July 31, 2017 1 Department of Biostatistics, Bloomberg School of Public Health, Johns
More informationIntroduction to Machine Learning. Introduction to ML - TAU 2016/7 1
Introduction to Machine Learning Introduction to ML - TAU 2016/7 1 Course Administration Lecturers: Amir Globerson (gamir@post.tau.ac.il) Yishay Mansour (Mansour@tau.ac.il) Teaching Assistance: Regev Schweiger
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 informationIntelligent Systems:
Intelligent Systems: Undirected Graphical models (Factor Graphs) (2 lectures) Carsten Rother 15/01/2015 Intelligent Systems: Probabilistic Inference in DGM and UGM Roadmap for next two lectures Definition
More informationAn EM algorithm for Gaussian Markov Random Fields
An EM algorithm for Gaussian Markov Random Fields Will Penny, Wellcome Department of Imaging Neuroscience, University College, London WC1N 3BG. wpenny@fil.ion.ucl.ac.uk October 28, 2002 Abstract Lavine
More information7 The structure of graphs excluding a topological minor
7 The structure of graphs excluding a topological minor Grohe and Marx [39] proved the following structure theorem for graphs excluding a topological minor: Theorem 7.1 ([39]). For every positive integer
More informationTAKEHOME FINAL EXAM e iω e 2iω e iω e 2iω
ECO 513 Spring 2015 TAKEHOME FINAL EXAM (1) Suppose the univariate stochastic process y is ARMA(2,2) of the following form: y t = 1.6974y t 1.9604y t 2 + ε t 1.6628ε t 1 +.9216ε t 2, (1) where ε is i.i.d.
More informationStatistical Practice
Statistical Practice A Note on Bayesian Inference After Multiple Imputation Xiang ZHOU and Jerome P. REITER This article is aimed at practitioners who plan to use Bayesian inference on multiply-imputed
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 informationAnalysis of Marked Point Patterns with Spatial and Non-spatial Covariate Information
Analysis of Marked Point Patterns with Spatial and Non-spatial Covariate Information p. 1/27 Analysis of Marked Point Patterns with Spatial and Non-spatial Covariate Information Shengde Liang, Bradley
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 information13: Variational inference II
10-708: Probabilistic Graphical Models, Spring 2015 13: Variational inference II Lecturer: Eric P. Xing Scribes: Ronghuo Zheng, Zhiting Hu, Yuntian Deng 1 Introduction We started to talk about variational
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 informationBAYESIAN HIERARCHICAL MODELS FOR MISALIGNED DATA: A SIMULATION STUDY
STATISTICA, anno LXXV, n. 1, 2015 BAYESIAN HIERARCHICAL MODELS FOR MISALIGNED DATA: A SIMULATION STUDY Giulia Roli 1 Dipartimento di Scienze Statistiche, Università di Bologna, Bologna, Italia Meri Raggi
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 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 informationSTATS 200: Introduction to Statistical Inference. Lecture 29: Course review
STATS 200: Introduction to Statistical Inference Lecture 29: Course review Course review We started in Lecture 1 with a fundamental assumption: Data is a realization of a random process. The goal throughout
More informationSpatial Smoothing in Stan: Conditional Auto-Regressive Models
Spatial Smoothing in Stan: Conditional Auto-Regressive Models Charles DiMaggio, PhD, NYU School of Medicine Stephen J. Mooney, PhD, University of Washington Mitzi Morris, Columbia University Dan Simpson,
More informationarxiv: v1 [stat.ap] 18 Nov 2015
Automatic Region-wise Spatially Varying Coefficient Regression Model: an Application to National Cardiovascular Disease Mortality and Air Pollution Association Study arxiv:1511.05924v1 [stat.ap] 18 Nov
More informationChris Fraley and Daniel Percival. August 22, 2008, revised May 14, 2010
Model-Averaged l 1 Regularization using Markov Chain Monte Carlo Model Composition Technical Report No. 541 Department of Statistics, University of Washington Chris Fraley and Daniel Percival August 22,
More informationRegularization Parameter Selection for a Bayesian Multi-Level Group Lasso Regression Model with Application to Imaging Genomics
Regularization Parameter Selection for a Bayesian Multi-Level Group Lasso Regression Model with Application to Imaging Genomics arxiv:1603.08163v1 [stat.ml] 7 Mar 016 Farouk S. Nathoo, Keelin Greenlaw,
More informationA Bayesian perspective on GMM and IV
A Bayesian perspective on GMM and IV Christopher A. Sims Princeton University sims@princeton.edu November 26, 2013 What is a Bayesian perspective? A Bayesian perspective on scientific reporting views all
More informationGaussian Process Regression Model in Spatial Logistic Regression
Journal of Physics: Conference Series PAPER OPEN ACCESS Gaussian Process Regression Model in Spatial Logistic Regression To cite this article: A Sofro and A Oktaviarina 018 J. Phys.: Conf. Ser. 947 01005
More informationFully Bayesian Spatial Analysis of Homicide Rates.
Fully Bayesian Spatial Analysis of Homicide Rates. Silvio A. da Silva, Luiz L.M. Melo and Ricardo S. Ehlers Universidade Federal do Paraná, Brazil Abstract Spatial models have been used in many fields
More informationKatsuhiro Sugita Faculty of Law and Letters, University of the Ryukyus. Abstract
Bayesian analysis of a vector autoregressive model with multiple structural breaks Katsuhiro Sugita Faculty of Law and Letters, University of the Ryukyus Abstract This paper develops a Bayesian approach
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 informationSparse Linear Models (10/7/13)
STA56: Probabilistic machine learning Sparse Linear Models (0/7/) Lecturer: Barbara Engelhardt Scribes: Jiaji Huang, Xin Jiang, Albert Oh Sparsity Sparsity has been a hot topic in statistics and machine
More informationBAYESIAN MODEL FOR SPATIAL DEPENDANCE AND PREDICTION OF TUBERCULOSIS
BAYESIAN MODEL FOR SPATIAL DEPENDANCE AND PREDICTION OF TUBERCULOSIS Srinivasan R and Venkatesan P Dept. of Statistics, National Institute for Research Tuberculosis, (Indian Council of Medical Research),
More informationHigh dimensional Ising model selection
High dimensional Ising model selection Pradeep Ravikumar UT Austin (based on work with John Lafferty, Martin Wainwright) Sparse Ising model US Senate 109th Congress Banerjee et al, 2008 Estimate a sparse
More informationA Fully Nonparametric Modeling Approach to. BNP Binary Regression
A Fully Nonparametric Modeling Approach to Binary Regression Maria Department of Applied Mathematics and Statistics University of California, Santa Cruz SBIES, April 27-28, 2012 Outline 1 2 3 Simulation
More informationLecture 5: Spatial probit models. James P. LeSage University of Toledo Department of Economics Toledo, OH
Lecture 5: Spatial probit models James P. LeSage University of Toledo Department of Economics Toledo, OH 43606 jlesage@spatial-econometrics.com March 2004 1 A Bayesian spatial probit model with individual
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 informationMarkov random fields. The Markov property
Markov random fields The Markov property Discrete time: (X k X k!1,x k!2,... = (X k X k!1 A time symmetric version: (X k! X!k = (X k X k!1,x k+1 A more general version: Let A be a set of indices >k, B
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 informationProperties of optimizations used in penalized Gaussian likelihood inverse covariance matrix estimation
Properties of optimizations used in penalized Gaussian likelihood inverse covariance matrix estimation Adam J. Rothman School of Statistics University of Minnesota October 8, 2014, joint work with Liliana
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 informationProbabilistic Graphical Models
Probabilistic Graphical Models David Sontag New York University Lecture 4, February 16, 2012 David Sontag (NYU) Graphical Models Lecture 4, February 16, 2012 1 / 27 Undirected graphical models Reminder
More informationProbabilistic Graphical Models Homework 2: Due February 24, 2014 at 4 pm
Probabilistic Graphical Models 10-708 Homework 2: Due February 24, 2014 at 4 pm Directions. This homework assignment covers the material presented in Lectures 4-8. You must complete all four problems to
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 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 Thinking in Spatial Statistics
Bayesian Thinking in Spatial Statistics Lance A. Waller Department of Biostatistics Rollins School of Public Health Emory University 1518 Clifton Road NE Atlanta, GA 30322 E-mail: lwaller@sph.emory.edu
More informationThe Ising model and Markov chain Monte Carlo
The Ising model and Markov chain Monte Carlo Ramesh Sridharan These notes give a short description of the Ising model for images and an introduction to Metropolis-Hastings and Gibbs Markov Chain Monte
More informationDownloaded from:
Camacho, A; Kucharski, AJ; Funk, S; Breman, J; Piot, P; Edmunds, WJ (2014) Potential for large outbreaks of Ebola virus disease. Epidemics, 9. pp. 70-8. ISSN 1755-4365 DOI: https://doi.org/10.1016/j.epidem.2014.09.003
More information13 : Variational Inference: Loopy Belief Propagation and Mean Field
10-708: Probabilistic Graphical Models 10-708, Spring 2012 13 : Variational Inference: Loopy Belief Propagation and Mean Field Lecturer: Eric P. Xing Scribes: Peter Schulam and William Wang 1 Introduction
More informationPrinciples of Bayesian Inference
Principles of Bayesian Inference Sudipto Banerjee University of Minnesota July 20th, 2008 1 Bayesian Principles Classical statistics: model parameters are fixed and unknown. A Bayesian thinks of parameters
More informationCOMPOSITIONAL IDEAS IN THE BAYESIAN ANALYSIS OF CATEGORICAL DATA WITH APPLICATION TO DOSE FINDING CLINICAL TRIALS
COMPOSITIONAL IDEAS IN THE BAYESIAN ANALYSIS OF CATEGORICAL DATA WITH APPLICATION TO DOSE FINDING CLINICAL TRIALS M. Gasparini and J. Eisele 2 Politecnico di Torino, Torino, Italy; mauro.gasparini@polito.it
More informationStat 5101 Lecture Notes
Stat 5101 Lecture Notes Charles J. Geyer Copyright 1998, 1999, 2000, 2001 by Charles J. Geyer May 7, 2001 ii Stat 5101 (Geyer) Course Notes Contents 1 Random Variables and Change of Variables 1 1.1 Random
More informationGaussian Multiscale Spatio-temporal Models for Areal Data
Gaussian Multiscale Spatio-temporal Models for Areal Data (University of Missouri) Scott H. Holan (University of Missouri) Adelmo I. Bertolde (UFES) Outline Motivation Multiscale factorization The multiscale
More informationStatistical Inference on Large Contingency Tables: Convergence, Testability, Stability. COMPSTAT 2010 Paris, August 23, 2010
Statistical Inference on Large Contingency Tables: Convergence, Testability, Stability Marianna Bolla Institute of Mathematics Budapest University of Technology and Economics marib@math.bme.hu COMPSTAT
More information39th Annual ISMS Marketing Science Conference University of Southern California, June 8, 2017
Permuted and IROM Department, McCombs School of Business The University of Texas at Austin 39th Annual ISMS Marketing Science Conference University of Southern California, June 8, 2017 1 / 36 Joint work
More informationIntroduction to graphical models: Lecture III
Introduction to graphical models: Lecture III Martin Wainwright UC Berkeley Departments of Statistics, and EECS Martin Wainwright (UC Berkeley) Some introductory lectures January 2013 1 / 25 Introduction
More informationScaling the Topology of Symmetric, Second-Order Planar Tensor Fields
Scaling the Topology of Symmetric, Second-Order Planar Tensor Fields Xavier Tricoche, Gerik Scheuermann, and Hans Hagen University of Kaiserslautern, P.O. Box 3049, 67653 Kaiserslautern, Germany E-mail:
More informationProbabilistic Graphical Models Lecture Notes Fall 2009
Probabilistic Graphical Models Lecture Notes Fall 2009 October 28, 2009 Byoung-Tak Zhang School of omputer Science and Engineering & ognitive Science, Brain Science, and Bioinformatics Seoul National University
More informationMinimax-Regret Sample Design in Anticipation of Missing Data, With Application to Panel Data. Jeff Dominitz RAND. and
Minimax-Regret Sample Design in Anticipation of Missing Data, With Application to Panel Data Jeff Dominitz RAND and Charles F. Manski Department of Economics and Institute for Policy Research, Northwestern
More informationBayesian analysis of conditional autoregressive models
Ann Inst Stat Math 202) 64:07 33 DOI 0.007/s0463-00-0298- Bayesian analysis of conditional autoregressive models Victor De Oliveira Received: 8 December 2008 / Revised: 4 January 200 / Published online:
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 informationAdding Spatially-Correlated Errors Can Mess Up The Fixed Effect You Love
Adding Spatially-Correlated Errors Can Mess Up The Fixed Effect You Love James S. Hodges 1, Brian J. Reich 2 1 Division of Biostatistics, University of Minnesota, Minneapolis, Minnesota USA 55414 2 Department
More information5. Discriminant analysis
5. Discriminant analysis We continue from Bayes s rule presented in Section 3 on p. 85 (5.1) where c i is a class, x isap-dimensional vector (data case) and we use class conditional probability (density
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 informationKneib, Fahrmeir: Supplement to "Structured additive regression for categorical space-time data: A mixed model approach"
Kneib, Fahrmeir: Supplement to "Structured additive regression for categorical space-time data: A mixed model approach" Sonderforschungsbereich 386, Paper 43 (25) Online unter: http://epub.ub.uni-muenchen.de/
More information