st_coordinates()
sf
inla()
Plotting the second order neighbours list:
plot(map_geom, border = "lightgray") plot.nb(nblags[[2]], map_geom, add = TRUE)