library(tidyverse) library(ggplot2) my_map <- read.csv(“EUvaccine.csv”) my_map

mapdata <- map_data(“world”) mapdata

mapdata <- left_join(mapdata, my_map, by = “region”)

mapdata

mapdata1 <- mapdata %>% filter(!is.na(mapdata$Perc_vaccinated)) mapdata1

map1 <- ggplot(mapdata1, aes(x = long, y = lat, group = group))+ geom_polygon(aes(fill= Perc_vaccinated), color = “black”)

map1

map2 <- map1 + scale_fill_gradient(name = “% vaccinated”, low = “yellow”, high = “red”, na.value = “grey50”)+ theme(axis.text.x = element_blank(), axis.text.y = element_blank(), axis.ticks = element_blank(), axis.title.y = element_blank(), axis.title.x = element_blank(), rect = element_blank()) map2