Mapping

Data Source

we got US population by race from https://worldpopulationreview.com/states/states-by-race and we got the resident population from www.census.gov

Those both datasets are open source that anyone can legally download it.

Loading data

df <- read.csv("Mapping.csv")
# df$State <- df$State %>%
#   abbr2state() %>%
#   tolower()

df$Resident.Population <- as.numeric(gsub(",", "", df$Resident.Population))

for (i in  6:56){
  df[i, 1] = sub('^.(.*)', '\\1' , df[i, 1])
}
df <- df[6:56, ]

df$State <- df$State %>% tolower()

# a = "39,512,223"
# as.numeric(gsub(",", "", a))

US population

us_pop <- ggplot(data = map1,
                 mapping = aes(x = long, y = lat,
                               group = group,
                               fill = Resident.Population)) +
  geom_polygon(color = "white", size = .1) +
  theme_bw() +
  coord_map() +
  scale_fill_continuous(type = "viridis")+
  ggtitle("US resident population") +
    theme(axis.title.x = element_blank(),
          axis.title.y = element_blank(),
          legend.title = element_blank())


us_pop

Asian population in US

us_pop <- ggplot(data = map1,
                 mapping = aes(x = long, y = lat,
                               group = group,
                               fill = Asian)) +
  geom_polygon(color = "white", size = .1) +
  theme_bw() +
  coord_map() +
  scale_fill_continuous()+
  ggtitle("US resident population") +
    theme(axis.title.x = element_blank(),
          axis.title.y = element_blank(),
          legend.title = element_blank())


us_pop

Those plot show the population in United State. We couls see that there are so many Asian race in CA, which make others state end up having the save colors

us_pop <- ggplot(data = map1,
                 mapping = aes(x = long, y = lat,
                               group = group,
                               fill = AsianTotalPerc)) +
  geom_polygon(color = "white", size = .1) +
  theme_bw() +
  coord_map() +
  scale_fill_continuous(type = "viridis")+
  ggtitle("US resident population") +
    theme(axis.title.x = element_blank(),
          axis.title.y = element_blank(),
          legend.title = element_blank())

us_pop

After re-scale to see only the percentage, we can see that it is much easier to define Asian population in each State