region <-
geolevel(name = "region",
layer = layer_us_region,
key = c("geoid"))
division <-
geolevel(name = "division",
layer = layer_us_division,
key = c("geoid"))
gd <-
geodimension(name = "gd_us",
level = region) |>
add_level(division)
gd <- gd |>
relate_levels(lower_level_name = "division",
upper_level_name = "region",
by_geography = TRUE)
Run the code above in your browser using DataLab