DEM 7093 - Lab 1/Homework 2: Using R for Basic Map making

Using tmap to make Maps with Quantile, Jenks, Pretty (and Fixed) Breaks

Reference:

• More styles here: https://geocompr.github.io/post/2019/tmap-color-scales/
• Used Color Brewer’s codes for the palettes from: https://colorbrewer2.org/
• Data is 2017 5-Year(?) ACS for Bexar County Census Tracts

Maps for Median Household Income (MHI)

library(tmap)
library(tmaptools)

#Pretty Breaks - MHI
tm_shape(sa_acs2)+
  tm_polygons("MHI", 
              title="MHI", 
              palette="YlGnBu", 
              border.col="#FFFFFF", 
              style="pretty",
              legend.hist=T )+
  tm_format("World", title="San Antonio MHI Estimates - Pretty Breaks", legend.outside=T)+
  tm_scale_bar(position="LEFT", breaks=c(0,2.5,5))+
  tm_compass()

#Fixed Breaks - MHI
tm_shape(sa_acs2, unit="mi")+
  tm_polygons("MHI", 
              title="MHI", 
              palette="YlGnBu", 
              border.col="#FFFFFF", 
              style="fixed", 
              breaks= c(0,25000,50000,75000,100000,250000),
              labels=c("$0 to 25K","$25 to 50K","$50 to 75K","$75 to 100K","$100K +"), 
              legend.hist=T )+
  tm_format("World", title="San Antonio MHI Estimates - Fixed Breaks", legend.outside=T)+
  tm_scale_bar(position="LEFT", breaks=c(0,2.5,5))+
  tm_compass(type="arrow")

Maps for Poverty Rate

#Quantile Breaks - Poverty
tm_shape(sa_acs2)+
  tm_polygons("ppov", 
              title="% in Poverty", 
              palette="BuPu", 
              border.col="#FFFFFF", 
              style="quantile", 
              n=5 ,
              legend.hist=T)+
  tm_format("World", title="San Antonio Poverty Estimates - Quantile Breaks", legend.outside=T)+
  tm_scale_bar()+
  tm_compass()

#Jenks Breaks - Poverty
tm_shape(sa_acs2)+
  tm_polygons("ppov", 
              title="% in Poverty", 
              palette="BuPu", 
              border.col="#FFFFFF", 
              style="jenks", 
              n=5,
              legend.hist=T )+
  tm_format("World", title="San Antonio Poverty Estimates - Jenks Breaks", legend.outside=T)+
  tm_scale_bar()+
  tm_compass()

#Pretty Breaks - Poverty
tm_shape(sa_acs2)+
  tm_polygons("ppov", 
              title="% in Poverty", 
              palette="BuPu", 
              border.col="#FFFFFF", 
              style="pretty", 
              n=5,
              legend.hist=T )+
  tm_format("World", title="San Antonio Poverty Estimates - Pretty Breaks", legend.outside=T)+
  tm_scale_bar()+
  tm_compass()