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.

An important issue in many cities is related to crime events, and spatio–temporal Bayesian approach leads to identifying crime patterns and hotspots. In Bayesian analysis of spatio–temporal crime data, there is no closed form for posterior distribution because of its non-Gaussian distribution and existence of latent variables. In this case, we face different challenges such as high dimensional parameters, extensive simulation and time-consuming computation in applying MCMC methods. In this paper, we use INLA to analyze crime data in Colombia. The advantages of this method can be the estimation of criminal events at a specific time and location and exploring unusual patterns in places.