Skip to content

Latest commit

 

History

History
99 lines (65 loc) · 3.17 KB

README.md

File metadata and controls

99 lines (65 loc) · 3.17 KB

ggcountymx

Generate ggplot2 geom_map Mexico county maps

This is a simple package with one purpose:

  • Make it easier to generate Mexico Municipality Choroplet maps

For installation, do:

library(devtools)
# Download Library (read Notes)
install_github("nicholasbucher/ggcountymx")

library(ggcountymx)

# Load dataset to Enviroment (read Notes)
data(geomx)

BCS <- ggcountymx("03")  # Entity code as Character 
BCS$gg

To get:

map1

The BCS object in the above code contains

  • the gg ggplot2 object
  • a map object which is the polygon points data frame
  • a mun.names object which is a list of all municipality names with FIPS codes in that municipality
  • a geom object (geom) for the municipality map

This lets you add further map layers (e.g. for a choropleth):

# built-in MEX socioeconomic by FIPS code data set from ENIGH 2016
data(socioeconomic)

# define appropriate (& nicely labeled) socioeconomic breaks

socioeconomic$brk<-cut(socioeconomic$Level,
                       breaks=c(0, 0.5 , 1 , 1.5 , 2 , 2.5 , 3 ,3.5 , 4),
                       labels=c("Outcast",
                                "Low",
                                "Lower Middle",
                                "Middle",
                                "Upper Middle",
                                "Lower High",
                                "Middle High",
                                "Upper High"),
                      include.lowest=TRUE)
                 
# get the Mexico municipality map
BCS <- ggcountymx("03")

# start the plot with our base map
gg <- BCS$gg

# add a new geom with our population (choropleth)

choro_geom <- geom_map(data=socioeconomic,
                       map=BCS$map,
                       aes(map_id=FIPS, fill=brk),
                       color="white", size=0.125)

gg<-gg + choro_geom

# define nice colors

choro_color<-scale_fill_manual(values=c("#ff0000", "#ff4000", "#ff8000",
                                        "#ffbf00", "#ffff00", "#bfff00",
                                        "#80ff00","#40ff00"),
                               name="Socioeconomic Class")

gg<-gg + choro_color

# plot the map
gg

map2

And, combining individual maps is pretty straightforward:

BCN <- ggcountymx("02", fill="#c7e9b4", color="white")
BCS <- ggcountymx("03", fill="#41b6c4", color="white")
SON <- ggcountymx("26", fill="#253494", color="white")

BCS$gg + BCN$geom + SON$geom

map2

Notes:

  • Mexico's territory is defined by more than three million geografic points, file weights 46mb, hence, download might take longer than usual.
  • File is to big to be read on every function call, after loading library always load dataset to Enviroment data(geomx).
  • resolution parameter reduces the percetage of points used.