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