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