Proper multivariate conditional autoregressive models for spatial data analysis.

In the past decade conditional autoregressive modelling specifications have found considerable application for the analysis of spatial data. Nearly all of this work is done in the univariate case and employs an improper specification. Our contribution here is to move to multivariate conditional autoregressive models and to provide rich, flexible classes which yield proper distributions. Our approach is to introduce spatial autoregression parameters. We first clarify what classes can be developed from the family of Mardia (1988) and contrast with recent work of Kim et al. (2000). We then present a novel parametric linear transformation which provides an extension with attractive interpretation. We propose to employ these models as specifications for second-stage spatial effects in hierarchical models. Two applications are discussed; one for the two-dimensional case modelling spatial patterns of child growth, the other for a four-dimensional situation modelling spatial variation in HLA-B allele frequencies. In each case, full Bayesian inference is carried out using Markov chain Monte Carlo simulation.

[1]  D. Brook.  On the distinction between the conditional probability and the joint probability approaches in the specification of nearest-neighbour systems , 1964 .

[2]  J. Besag.  Spatial Interaction and the Statistical Analysis of Lattice Systems , 1974 .

[3]  Donald Geman,et al.  Stochastic Relaxation, Gibbs Distributions, and the Bayesian Restoration of Images , 1984, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[4]  P. Heywood,et al.  Nutritional status of young children--the 1982/83 National Nutrition Survey. , 1988, Papua and New Guinea medical journal.

[5]  K. Mardia.  Multi-dimensional multivariate Gaussian Markov random fields with application to image processing , 1988 .

[6]  J. Besag,et al.  Bayesian image restoration, with two applications in spatial statistics , 1991 .

[7]  T J Cole,et al.  Smoothing reference centile curves: the LMS method and penalized likelihood. , 1992, Statistics in medicine.

[8]  J. Waterlow,et al.  Protein Energy Malnutrition , 1993, Definitions.

[9]  C Montomoli,et al.  Spatial correlation in ecological analysis. , 1993, International journal of epidemiology.

[10]  Robert Haining,et al.  Statistics for spatial data: by Noel Cressie, 1991, John Wiley & Sons, New York, 900 p., ISBN 0-471-84336-9, US $89.95 , 1993 .

[11]  Kerrie Mengersen,et al.  [Bayesian Computation and Stochastic Systems]: Rejoinder , 1995 .

[12]  J. Besag,et al.  Bayesian Computation and Stochastic Systems , 1995 .

[13]  Alan E. Gelfand,et al.  Model choice: A minimum posterior predictive loss approach , 1998, AISTATS.

[14]  H. Goldstein,et al.  Multilevel Modelling of the Geographical Distributions of Diseases , 1999, Journal of the Royal Statistical Society. Series C, Applied statistics.

[15]  H. Goldstein,et al.  Multivariate spatial models for event data. , 2000, Statistics in medicine.

[16]  A E Gelfand,et al.  Spatial modelling of multinomial data with latent structure: an application to geographical mapping of human gene and haplotype frequencies. , 2000, Biostatistics.

[17]  Håvard Rue,et al.  On block updating in Markov random field models for disease mapping. (REVISED, May 2001) , 2000 .

[18]  P. Vounatsou,et al.  Spatial patterns of child growth in Papua New Guinea and their relation to environment, diet, socio-economic status and subsistence activities , 2001, Annals of human biology.

[19]  H. Rue,et al.  On Block Updating in Markov Random Field Models for Disease Mapping , 2002 .