Bayesian analysis of conditional autoregressive models

Conditional autoregressive (CAR) models have been extensively used for the analysis of spatial data in diverse areas, such as demography, economy, epidemiology and geography, as models for both latent and observed variables. In the latter case, the most common inferential method has been maximum likelihood, and the Bayesian approach has not been used much. This work proposes default (automatic) Bayesian analyses of CAR models. Two versions of Jeffreys prior, the independence Jeffreys and Jeffreys-rule priors, are derived for the parameters of CAR models and properties of the priors and resulting posterior distributions are obtained. The two priors and their respective posteriors are compared based on simulated data. Also, frequentist properties of inferences based on maximum likelihood are compared with those based on the Jeffreys priors and the uniform prior. Finally, the proposed Bayesian analysis is illustrated by fitting a CAR model to a phosphate dataset from an archaeological region.

[1]  Noel A. C. Cressie,et al.  Statistics for Spatial Data: Cressie/Statistics , 1993 .

[2]  T. Raghavan,et al.  Nonnegative Matrices and Applications , 1997 .

[3]  L D Broemeling,et al.  A bayesian analysis for spatial processes with application to disease mapping. , 2000, Statistics in medicine.

[4]  Anthony N. Pettitt,et al.  A Conditional Autoregressive Gaussian Process for Irregularly Spaced Multivariate Data with Application to Modelling Large Sets of Binary Data , 2002, Stat. Comput..

[5]  Leonhard Held,et al.  Gaussian Markov Random Fields: Theory and Applications , 2005 .

[6]  D. Clayton,et al.  Empirical Bayes estimates of age-standardized relative risks for use in disease mapping. , 1987, Biometrics.

[7]  Sylvia Richardson,et al.  Spatial Linear Models with Autocorrelated Error Structure , 1992 .

[8]  Harvey S. Rosen,et al.  Budget spillovers and fiscal policy interdependence: Evidence from the states , 1993 .

[9]  N. Cressie,et al.  Spatial Modeling of Regional Variables , 1993 .

[10]  James R. Schott,et al.  Matrix Analysis for Statistics , 2005 .

[11]  Sw. Banerjee,et al.  Hierarchical Modeling and Analysis for Spatial Data , 2003 .

[12]  J. Berger,et al.  Objective Bayesian Analysis of Spatially Correlated Data , 2001 .

[13]  W. Gilks,et al.  Adaptive Rejection Metropolis Sampling Within Gibbs Sampling , 1995 .

[14]  M. D. Ugarte,et al.  Alternative Models for Describing Spatial Dependence among Dwelling Selling Prices , 2004 .

[15]  Mike Rees,et al.  5. Statistics for Spatial Data , 1993 .

[16]  N. Cressie,et al.  Some Diagnostics for Markov Random Fields , 2008 .

[17]  K. Mardia,et al.  Maximum likelihood estimation of models for residual covariance in spatial regression , 1984 .

[18]  C. D. Litton,et al.  The spatial analysis of site phos-phate data , 1988 .

[19]  J. Bernardo Reference Posterior Distributions for Bayesian Inference , 1979 .

[20]  Noel A Cressie,et al.  Likelihood-based estimation for Gaussian MRFs , 2005 .

[21]  J. Besag,et al.  On conditional and intrinsic autoregressions , 1995 .

[22]  Noel A Cressie,et al.  Statistics for Spatial Data. , 1992 .

[23]  P. Speckman,et al.  Posterior distribution of hierarchical models using CAR(1) distributions , 1999 .

[24]  Andreas Karlsson,et al.  Matrix Analysis for Statistics (2nd ed.), by James R. Schott , 2007 .