We did 100 simulations with **data** size = 500 and = 6. We computed SAR probit parameter estimation for each method; MLE, linearized GMM, Bayes, and conditional approximate likelihood. We constructed a cumulative mean plot for every SAR probit parameter estimation with a variable value of spatial lag as the following:  Figure 1-3, we can see that parameter estimation from linearized GMM is the worst method because the estimation is so far away from the real value of SAR probit parameter (beta). We also can check this by a cumulative mean plot of w x y − ; z, = 0,1,2 in figure 4-6.    