Creating higher order neighbours lists

Plotting the second order neighbours list:

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