Propensity score can be described as a probability of certain treatments conditional to the given observed covariates. Propensity score is one of the known methods to allows an observational study emulating certain characteristics from that of a randomized trial. The most common method used to estimate this score is the logistic regression model. Logistic regression can be used to model the probability of a certain event. With the advancement that is happening to spatial statistics, one can also build a logistic regression model that takes into consideration to that of spatial dependence. Thus, accommodate the spatial effect that is likely happening on observation data that came from different places. Problem arises from this model, that is the estimation of the parameters on the spatial logistic model. EM algorithm which is needed for this problem, still requires another adjustment since the expectation in the E-step is not available in closed form. Variational method modification is then proposed as an alternative for this problem. This paper reviews the propensity score estimation using spatial logistic regression and discusses the variational method as an alternative method to tackle the problem in estimating the parameters on the spatial logistic regression model in a theoretical study.