I want to fit a GLM and control for spatial relationships within the data. To this end, I wanted to include the term spaMM::Matern(1 | x + y) where x and y are the column and row index from the raster stack where the other data points were sampled.

I have run a number of hurdle models and these models have the lowest AIC value, within 2 AIC of eachother. I am trying to see if there is a way to visualise model fit for both of these models.