In spatial generalized linear mixed models, spatial correlation is assumed by adding normal latent variables to the model. In these models because of the non-Gaussian spatial response and the presence of latent variables, the likelihood function cannot usually be given in a closed form, thus the maximum likelihood approach is very challenging. The main purpose of this paper is to introduce two new algorithms for the maximum likelihood estimations of parameters and to compare them in terms of speed and accuracy with existing algorithms. The presented algorithms are applied to a simulation study and their performances are compared.