Case: housing prices in Boston

map <- st_read(system.file("shapes/boston_tracts.shp",
                           package = "spData"), quiet = TRUE)
map$vble <- log(map$MEDV)
dim(map)

## [1] 506  38