In previous research, a logistic regression of landslide occurrence on several explanatory variables was fitted and used to map landslide susceptibility for a small area in the central Apennines, Italy. Here, the spatial dependence or spatial correlation in the residuals from the fitted regression model was accounted for by inserting an autocovariate into the model. The autocovariate was estimated by applying a Gibbs sampler to the susceptibilities for neighbouring pixels. As in any landslide susceptibility analysis, accuracy was difficult to assess because of the requirement for data on future landslides. However, by comparing the ordinary logistic model to the autologistic model obtained on the same set of data, it was possible to assess the influence of the autocovariate. The autocovariate rendered the model simpler because several variables lost their significance and were, therefore, omitted from the model. Further, areas of high landslide susceptibility estimated from the autologistic model were geographically clustered, as one would expect, and this may be advantageous in terms of (i) interpreting the model and (ii) displaying the results to non-experts.