---

In this paper parameters of spatial covariance functions have been estimated using block composite likelihood method. In this method, the block composite likelihood is constructed from the joint densities of paired spatial blocks. For this purpose, after differencing data, large data sets are splited into many smaller data sets. Then each separated blocks evaluated separately and finally combined through a simple summation. The advantage of this method is that there is no need to inverse and to find determination of high dimensional matrices. The simulation shows that the block composite likelihood estimates as well as the pair composite likelihood. Finally a real data is analysed.

---

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.