Map of Europe at NUTS 2 Level
Miluji Sb <milujisb <at> gmail.com> writes:
Dear all. I would like to draw a map of France, Italy, Spain, and Portugal at NUTS 2 level. I used the following code:
...
# Subset NUTS 2 level data map_nuts2 <- subset(EU_NUTS, STAT_LEVL_ == 2) # Draw basic plot plot(map_nuts2)
country <- substring(as.character(map_nuts2$NUTS_ID), 1, 2)
map <- c("ES", "FR", "IT", "PT")
map_nuts2a <- map_nuts2[country %in% map,]
plot(map_nuts2)
but this includes all the overseas territories and islands. Use xlim=, ylim=
to exclude these, or select based on knowing their NUTS2 codes:
plot(nuts2a, axes=TRUE, xlim=c(-11.62699, 19.20174), ylim=c(34.35470, 52.04182))
Roger
PS. One gets more, and more rapid replies to this kind of question on R-sig-geo.
This does produce a plot but its rather 'ugle'. Is there any way I can subset the data further and draw a map for France, Italy, Spain, and Portugal only? Thank you very much! Sincerely, Milu