Simple Multiple Maps

library(maps)

## Warning: package 'maps' was built under R version 3.6.3

library(ggplot2)

## Warning: package 'ggplot2' was built under R version 3.6.3

library(data.table)

## Warning: package 'data.table' was built under R version 3.6.2

library(dplyr)

## Warning: package 'dplyr' was built under R version 3.6.3

## 
## Attaching package: 'dplyr'

## The following objects are masked from 'package:data.table':
## 
##     between, first, last

## The following objects are masked from 'package:stats':
## 
##     filter, lag

## The following objects are masked from 'package:base':
## 
##     intersect, setdiff, setequal, union

us_states <- map_data("state")
head(us_states)

##        long      lat group order  region subregion
## 1 -87.46201 30.38968     1     1 alabama      <NA>
## 2 -87.48493 30.37249     1     2 alabama      <NA>
## 3 -87.52503 30.37249     1     3 alabama      <NA>
## 4 -87.53076 30.33239     1     4 alabama      <NA>
## 5 -87.57087 30.32665     1     5 alabama      <NA>
## 6 -87.58806 30.32665     1     6 alabama      <NA>

load("J:/Github_Folks/socviz/data/opiates.rda")
opiates

## # A tibble: 800 x 11
##     year state  fips deaths population crude adjusted adjusted_se region abbr 
##    <int> <chr> <int>  <int>      <int> <dbl>    <dbl>       <dbl> <ord>  <chr>
##  1  1999 Alab~     1     37    4430141   0.8      0.8         0.1 South  AL   
##  2  1999 Alas~     2     27     624779   4.3      4           0.8 West   AK   
##  3  1999 Ariz~     4    229    5023823   4.6      4.7         0.3 West   AZ   
##  4  1999 Arka~     5     28    2651860   1.1      1.1         0.2 South  AR   
##  5  1999 Cali~     6   1474   33499204   4.4      4.5         0.1 West   CA   
##  6  1999 Colo~     8    164    4226018   3.9      3.7         0.3 West   CO   
##  7  1999 Conn~     9    151    3386401   4.5      4.4         0.4 North~ CT   
##  8  1999 Dela~    10     32     774990   4.1      4.1         0.7 South  DE   
##  9  1999 Dist~    11     28     570213   4.9      4.9         0.9 South  DC   
## 10  1999 Flor~    12    402   15759421   2.6      2.6         0.1 South  FL   
## # ... with 790 more rows, and 1 more variable: division_name <chr>

opiates$region <- tolower(opiates$state)
opiates_map <- left_join(us_states, opiates)

## Joining, by = "region"

library(viridis)

## Warning: package 'viridis' was built under R version 3.6.2

## Loading required package: viridisLite

## Warning: package 'viridisLite' was built under R version 3.6.2

library(ggthemes)

## Warning: package 'ggthemes' was built under R version 3.6.3

p0 <- ggplot(data = subset(opiates_map, year > 1999),
             mapping = aes(x = long, y = lat,
                           group = group,
                           fill = adjusted))

p1 <- p0 + geom_polygon(color = "gray90", size = 0.05) +
  coord_map(projection = "albers", lat0 = 39, lat1 = 45) 

p2 <- p1 + scale_fill_viridis_c(option = "plasma")

p2 + theme_map() + facet_wrap(~ year, ncol = 3) +
  theme(legend.position = "bottom",
        strip.background = element_blank()) +
  labs(fill = "Death rate per 100,000 population ",
       title = "Opiate Related Deaths by State, 2000-2014")