---

Gaussian random field is usually used to model Gaussian spatial data. In practice, we may encounter non-Gaussian data that are skewed. One solution to model skew spatial data is to use a skew random field. Recently, many skew random fields have been proposed to model this type of data, some of which have problems such as complexity, non-identifiability, and non-stationarity. In this article, a flexible class of closed skew-normal distribution is introduced to construct valid stationary random fields, and some important properties of this class such as identifiability and closedness under marginalization and conditioning are examined. The reasons for developing valid spatial models based on these skew random fields are also explained. Additionally, the identifiability of the spatial correlation model based on empirical variogram is investigated in a simulation study with the stationary skew random field as a competing model. Furthermore, spatial predictions using a likelihood approach are presented on these skew random fields and a simulation study is performed to evaluate the likelihood estimation of their parameters. 

---

In some instances, the occurrence of an event can be influenced by its spatial location, giving rise to spatial survival data. The accurate and precise estimation of parameters in a spatial survival model poses a challenge due to the complexity of the likelihood function, highlighting the significance of employing a Bayesian approach in survival analysis. In a Bayesian spatial survival model, the spatial correlation between event times is elucidated using a geostatistical model. This article presents a simulation study to estimate the parameters of classical and spatial survival models, evaluating the performance of each model in fitting simulated survival data. Ultimately, it is demonstrated that the spatial survival model exhibits superior efficacy in analyzing blood cancer data compared to conventional models.

---

The spatial generalized linear mixed models are often used, where the latent variables representing spatial correlations are modeled through a Gaussian random field to model the categorical spatial data. The violation of the Gaussian assumption affects the accuracy of predictions and parameter estimates in these models. In this paper, the spatial generalized linear mixed models are fitted and analyzed by utilizing a stationary skew Gaussian random field and employing an approximate Bayesian approach. The performance of the model and the approximate Bayesian approach is examined through a simulation example, and implementation on an actual data set is presented.

---

Spatial regression models are used to analyze quantitative spatial responses based on linear and non-linear relationships with explanatory variables. Usually, the spatial correlation of responses is modeled with a Gaussian random field based on a multivariate normal distribution. However, in practice, we encounter skewed responses, which are analyzed using skew-normal distributions. Closed skew-normal distribution is one of the extended families of skew-normal distributions, which has similar properties to normal distributions. This article presents a hierarchical Bayesian analysis based on a flexible subclass of closed skew-normal distributions. Given the time-consuming nature of Monte Carlo methods in hierarchical Bayes analysis, we have opted to use the variational Bayes approach to approximate the posterior distribution. This decision was made to expedite the analysis process without compromising the accuracy of our results. Then, the proposed model is implemented and analyzed based on the real earthquake data of Iran.

---

In this article, autoregressive spatial regression and second-order moving average will be presented to model the outputs of a heavy-tailed skewed spatial random field resulting from the developed multivariate generalized Skew-Laplace distribution. The model parameters are estimated by the maximum likelihood method using the Kolbeck-Leibler divergence criterion. Also, the best spatial predictor will be provided. Then, a simulation study is conducted to validate and evaluate the performance of the proposed model. The method is applied to analyze a real data.

---

Financial and economic indicators, such as housing prices, often show spatial correlation and heterogeneity. While spatial econometric models effectively address spatial dependency, they face challenges in capturing heterogeneity. Geographically weighted regression is naturally used to model this heterogeneity, but it can become too complex when data show homogeneity across subregions. In this paper, spatially homogeneous subareas are identified through spatial clustering, and Bayesian spatial econometric models are then fitted to each subregion. The integrated nested Laplace approximation method is applied to overcome the computational complexity of posterior inference and the difficulties of MCMC algorithms. The proposed methodology is assessed through a simulation study and applied to analyze housing prices in Mashhad City.