Abstract
This study applies computationally intensive methods for Bayesian analysis of spatially distributed data. It is assumed that the space is divided in contiguous and disjoint regions or areas. The neighboring structure in a given problem may indicate a wide range of number of neighbors per area, ranging from very few neighbors to cases where all areas neighbor each other. The main aim of this work is to evaluate the influence of neighborhood on results of Markov Chain Monte Carlo (MCMC) methods. Proper and improper prior specifications for state parameters are compared. Three schemes, proposed in the literature, for sampling from the joint posterior distribution are also compared. The comparison criterion is based on the autocorrelation structure of the chains. Two classes of models are studied: the first one is characterized by a simple model without any explanatory variables and the second one is an extension with multiple regression components. Initially, sensitivity of the analysis to different prior distributions is addressed. Finally, extensive empirical analyses confront the outcomes obtained with different neighboring arrangements of the units. Results are shown to generalize those obtained with dynamic or state space models.
Similar content being viewed by others
References
Assunção JJ, Gamerman D, Assunção RM (1999) Regional differences in factor productivities of Brazilian agriculture: a Bayesian space varying parameter approach. In: Proceedings of the XVII Latin American meeting of the Econometric Society, Technical Report 127, Statistical Laboratory, Universidade Federal do Rio de Janeiro, Rio de Janeiro
Bernardinelli L, Clayton D, Pascutto C, Montomoli C, Ghislandi M, Songini M (1995) Bayesian analysis of space-time variation in disease risk. Stat Med 14: 2433–2443
Besag J, Kooperberg CL (1995) On conditional and intrinsic autoregressions. Biometrika 82: 733–746
Besag J, York J, Mollié A (1991) Bayesian image restoration, with two applications in spatial statistics. (with discussion). Ann Inst Stat Math 43: 1–59
Besag J, Green P, Higdon D, Mengersen K (1995) Bayesian computation and stochastic system. Stat Sci 10(1): 3–66
Carter C, Kohn R (1994) On Gibbs sampling for state space models. Biometrika 81: 541–553
Doornik JA, Ooms M (2006) Introduction to Ox: an object-oriented matrix language. Timberlake Consultants Press, London
Doucet A, Briers M, Senecal S (2006) Efficient block sampling strategies for sequential Monte Carlo methods. J Comput Graph Stat 15(3): 693–711
Fernandez C, Green P (2002) Modelling spatially correlated data via mixtures: a Bayesian approach. J R Stat Soc Ser B 64(4): 805–826
Ferreira MAR, Oliveira V (2007) Bayesian reference analysis for Gaussian Markov random fields. J Multivar Anal 98(4): 789–812
Fruhwirth-Schnatter S (1994) Data augmentation and dynamic linear models. J Time Ser Anal 15: 183–202
Gamerman D (1998) Markov Chain Monte Carlo for dynamic generalised linear models. Biometrika 85: 215–227
Gamerman D, Lopes HF (2006) Markov Chain Monte Carlo: stochastic simulation for Bayesian inference, vol 68, 2nd edn. Chapman & Hall/CRC, London
Gamerman D, Moreira ARB, Rue H (2003) Space varying regression models: specifications and simulation. Comput Stat Data Anal 42: 513–533
Gelfand AE, Vounatsou P (2002) Proper multivariate conditional autoregressive models for spatial data analysis. Biostatistics 4: 11–25
Knorr-Held L (1999) Conditional prior proposals in dynamic models. Scand J Stat 26: 129–144
Knorr-Held L, Rue H (2002) On block updating in Markov random field models for disease mapping. Scand J Stat 29: 597–614
Knorr-Held L, Raber G, Becker N (2002) Disease mapping of stage-specific cancer incidence data. Biometrics 58: 492–501
Liu JS, Wong WH, Kong A (1994) Correlation structure and convergence rate of the Gibbs sampler: applications to the comparison of estimators and augmentation schemes. Biometrika 81: 27–40
Moreira ARB, Migon HS (1999) Spatial heterogeneity of the agricultural productivity: Brazil 1970–1996. In: Proceedings of the XXI Brazilian meeting on econometrics, Sociedade Brasileira de Econometria, Rio de Janeiro
Reis EA, Salazar E, Gamerman D (2006) Comparison of sampling schemes for dynamic linear models. Int Stat Rev 74(2): 203–214
Rue H (2001) Fast sampling of Gaussian Markov random fields. J R Stat Soc Ser B 63(2): 325–338
Rue H, Tjelmeland H (2002) Fitting Gaussian Markov random fields to Gaussian random fields. Scand J Stat 29: 31–49
Rue H, Steinsland I, Sveinung E (2004) Approximating hidden Gaussian Markov random fields. J R Stat Soc Ser B 66(4): 877–892
Shephard N (1994) Partial non-Gaussian state space. Biometrika 81: 115–131
van Dyk DA, Park T (2008) Partially collapsed Gibbs samplers: theory and methods. J Am Stat Assoc 103: 790–796
West M, Harrison PJ (1997) Bayesian forecasting and dynamic models, 2nd edn. Springer-Verlag, New York
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
About this article
Cite this article
Mayrink, V.D., Gamerman, D. On computational aspects of Bayesian spatial models: influence of the neighboring structure in the efficiency of MCMC algorithms. Comput Stat 24, 641–669 (2009). https://doi.org/10.1007/s00180-009-0153-0
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00180-009-0153-0