st_coordinates()
sf
inla()
nblag(<neighbours list>, maxlag =): to produce maxlag higher order neighbours lists
nblag(<neighbours list>, maxlag =)
maxlag