Parameter Estimation in Spatial Generalized Linear Mixed Models with Skew Gaussian Random Effects using Laplace Approximation

 Spatial generalized linear mixed models are used commonly for modelling non-Gaussian discrete spatial responses. We present an algorithm for parameter estimation of the models using Laplace approximation of likelihood function. In these models, the spatial correlation structure of data is carried out by random effects or latent variables. In most spatial analysis, it is assumed that random effects have Gaussian distribution, but the assumption is questionable. This assumption is replaced in the present work, using a skew Gaussian distribution for the latent variables, which is more flexible and includes Gaussian distribution. We examine the proposed method using a real discrete data set.