R help - Jun 2021 - Plot NUTS regions with color

Dear R-help ML,

I have downloaded the NUTS shapefiles from here: 
https://ec.europa.eu/eurostat/web/gisco/geodata/reference-data/administrative-units-statistical-units/nuts

I can load them with the following code:

    library(rgdal)
    shp_bdir <- [PATH TO SHAPE FILE]
    layername <- "NUTS_RG_01M_2021_4326"
    shp_folder <- file.path(shp_bdir, paste0(layername,".shp"))
    EU_NUTS <- readOGR(dsn = shp_folder, layer = layername)

Then I can plot the regions with:

    plot(EU_NUTS)

Now I have a list of NUTS - 3 Level, for example:

    l <- c("AT223", "AT212", "AT212",
"AT121")

I would like a plot where the NUTS regions are colored, and the more a 
particular NUTS is present is the list, the darker its color.

For example "AT212" should be darker than "AT223", because
the former is
present two times in the list.

How can I achieve that?

Thank you and best regards.


	[[alternative HTML version deleted]]

I suggest you post this in the r-sig-geo list rather than here. The
expertise you seek is more likely to be there.

Bert Gunter

"The trouble with having an open mind is that people keep coming along and
sticking things into it."
-- Opus (aka Berkeley Breathed in his "Bloom County" comic strip )


On Wed, Jun 23, 2021 at 10:36 AM Phillips Rogfield <thebudget72 at
gmail.com>
wrote:
> Dear R-help ML,
>
> I have downloaded the NUTS shapefiles from here:
>
>
https://ec.europa.eu/eurostat/web/gisco/geodata/reference-data/administrative-units-statistical-units/nuts
>
> I can load them with the following code:
>
>     library(rgdal)
>     shp_bdir <- [PATH TO SHAPE FILE]
>     layername <- "NUTS_RG_01M_2021_4326"
>     shp_folder <- file.path(shp_bdir,
paste0(layername,".shp"))
>     EU_NUTS <- readOGR(dsn = shp_folder, layer = layername)
>
> Then I can plot the regions with:
>
>     plot(EU_NUTS)
>
> Now I have a list of NUTS - 3 Level, for example:
>
>     l <- c("AT223", "AT212", "AT212",
"AT121")
>
> I would like a plot where the NUTS regions are colored, and the more a
> particular NUTS is present is the list, the darker its color.
>
> For example "AT212" should be darker than "AT223",
because the former is
> present two times in the list.
>
> How can I achieve that?
>
> Thank you and best regards.
>
>
>         [[alternative HTML version deleted]]
>
> ______________________________________________
> R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide
> http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>
	[[alternative HTML version deleted]]