Creating higher order neighbours lists

Plotting the third order neighbours list:

plot(map_geom, border = "lightgray")
plot.nb(nblags[[3]], map_geom, add = TRUE)