# NOT RUN {
library(sf)
nuts2 <- gisco_get_nuts(nuts_level = "2", country = "ITA")
nuts3 <- gisco_get_nuts(nuts_level = "3", country = "ITA")
plot(
st_geometry(nuts3),
lwd = 1,
border = "grey40"
)
plot(st_geometry(nuts2),
add = TRUE)
title(
main = "NUTS2 and 3 Levels on Italy",
sub = gisco_attributions(),
cex.sub = 0.7,
line = 1
)
# }
Run the code above in your browser using DataLab