st_coordinates()
sf
inla()
BYM2 = different parametrisation of BYM
bi=1√τ⋅(√1−ϕ⋅vi+√ϕ⋅ui)
inla.doc("^besag$", section = "latent") inla.doc("^bym$", section = "latent") inla.doc("^bym2$", section = "latent")