st_coordinates()
sf
inla()
map <- st_read(system.file("shapes/boston_tracts.shp", package = "spData"), quiet = TRUE) map$vble <- log(map$MEDV) dim(map)
## [1] 506 38