st_coordinates()
sf
inla()
Step 2:
nb4 <- knn2nb(knn_centroids) head(nb4)
## [[1]] ## [1] 2 3 4 ## ## [[2]] ## [1] 1 4 8 ## ## [[3]] ## [1] 1 4 5 ## ## [[4]] ## [1] 2 7 8 ## ## [[5]] ## [1] 3 8 11 ## ## [[6]] ## [1] 5 9 10