Spatial generalized linear mixed models with multivariate CAR models for areal data

Abstract Disease mapping studies have been widely performed with considering only one disease in the estimated models. Simultaneous modeling of different diseases can also be a valuable tool both from the epidemiological and from the statistical point of view. In particular, when we have several measurements recorded at each spatial location, we need to consider multivariate models in order to handle the dependence among the multivariate components as well as the spatial dependence between locations. These models can be studied in the class of spatial generalized linear mixed models (SGLMMs). It is well known that the frequentist analysis of SGLMMs is computationally difficult. Recently, there are a few papers which explored multivariate spatial models for areal data adopting the Bayesian framework as the natural inferential approach. We propose to use an approach, which yields to maximum likelihood estimation, to conduct frequentist analysis of SGLMMs with multivariate conditional autoregressive models for areal data. The performance of the proposed approach is evaluated through simulation studies and also by a real dataset.

[1]  Dongchu Sun,et al.  A Bivariate Bayes Method for Improving the Estimates of Mortality Rates With a Twofold Conditional Autoregressive Model , 2001 .

[2]  Bradley P Carlin,et al.  Order‐free co‐regionalized areal data models with application to multiple‐disease mapping , 2007, Journal of the Royal Statistical Society. Series B, Statistical methodology.

[3]  J. Andrew Royle,et al.  A hierarchical approach to multivariate spatial modeling and prediction , 1999 .

[4]  Mahmoud Torabi,et al.  Spatio‐temporal modelling of disease mapping of rates , 2010 .

[5]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[6]  Mahmoud Torabi,et al.  Hierarchical Bayes estimation of spatial statistics for rates , 2012 .

[7]  Mahmoud Torabi,et al.  Hierarchical Bayesian Spatiotemporal Analysis of Childhood Cancer Trends , 2012 .

[8]  Subhash R. Lele,et al.  Estimability and Likelihood Inference for Generalized Linear Mixed Models Using Data Cloning , 2010 .

[9]  A. M. Walker On the Asymptotic Behaviour of Posterior Distributions , 1969 .

[10]  A. Gelfand,et al.  Proper multivariate conditional autoregressive models for spatial data analysis. , 2003, Biostatistics.

[11]  Mahmoud Torabi,et al.  Spatio-temporal modelling using B-spline for disease mapping: analysis of childhood cancer trends , 2011 .

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

[13]  Natural variability of benthic species composition in the Delaware Bay , 2004, Environmental and Ecological Statistics.

[14]  James D. Hamilton A standard error for the estimated state vector of a state-space model , 1986 .

[15]  Sudipto Banerjee,et al.  SMOOTHED ANOVA WITH SPATIAL EFFECTS AS A COMPETITOR TO MCAR IN MULTIVARIATE SPATIAL SMOOTHING. , 2010, The annals of applied statistics.

[16]  Bradley P Carlin,et al.  Generalized Hierarchical Multivariate CAR Models for Areal Data , 2005, Biometrics.

[17]  Brian Dennis,et al.  Data cloning: easy maximum likelihood estimation for complex ecological models using Bayesian Markov chain Monte Carlo methods. , 2007, Ecology letters.

[18]  Péter Sólymos,et al.  dclone: Data Cloning in R , 2010, R J..