“…Their methods allowed for spatial autocorrelation but have not been able to benefit from covariates. Still, making inferences for autologistic models is notoriously difficult, requiring numerical approximaton methods for estimation ( Griffith, 2004;Hughes et al, 2011;Hardouin and Cressie, 2018 ). The main disadvantage is that the extension to multiple scale models increases the dimension of the parameter spaces while the calibration data available do not increase in most cases.…”