Hierarchical spatio-temporal models are used for modeling space-time responses and temporally and spatially correlations of the data is considered via Gaussian latent random field with Matérn covariance function. The most important interest in these models is estimation of the model parameters and the latent variables, and is predict of the response variables at new locations and times. In this paper, to analyze these models, the Bayesian approach is presented. Because of the complexity of the posterior distributions and the full conditional distributions of these models and the use of Monte Carlo samples in a Bayesian analysis, the computation time is too long. For solving this problem, Gaussian latent random field with Matern covariance function are represented as a Gaussian Markov Random Field (GMRF) through the Stochastic Partial Differential Equations (SPDE) approach. Approximatin Baysian method and Integrated Nested Laplace Approximation (INLA) are used to obtain an approximation of the posterior distributions and to inference about the model. Finally, the presented methods are applied to a case study on rainfall data observed in the weather stations of Semnan in 2013.

Geostatistical spatial count data in finite populations can be seen in many applications, such as urban management and medicine. The traditional model for analyzing these data is the spatial logit-binomial model. In the most applied situations, these data have overdispersion alongside the spatial variability. The binomial model is not the appropriate candidate to account for the overdispersion. The proper alternative is a beta-binomial model that has sufficient flexibility to account for the extra variability due to the possible overdispersion of counts. In this paper, we describe a Bayesian spatial beta-binomial for geostatistical count data by using a combination of the integrated nested Laplace approximation and the stochastic partial differential equations methods. We apply the methodology for analyzing the number of people injured/killed in car crashes in Mashhad, Iran. We further evaluate the performance of the model using a simulation study.