Logo Questions Linux Laravel Mysql Ubuntu Git Menu
 

ggplot2 mapping county boundries in one color and state boundries in another on the same map

Tags:

r

ggplot2

I am creating a choropleth county map with grey borders, and I also want to include the state boundries in black. Does anyone know how I go about adding a second layer of state mapping to an existing county map?

Here's the data set and code I ended up using:

#load libraries
    library(ggplot2)
    library(ggmap)
    library(maps)
    library(plyr)

#get wif file
wip <- read.csv("wip.csv") 

#get map data for US counties and states
county_map <- map_data("county")
state_map <- map_data("state")

#merge wip and county_map
wip_map <- merge(county_map, wip, by.x=c("region", "subregion"), 
    by.y=c("region","subregion"), all.x=TRUE)

#resort merged data
wip_map <- arrange(wip_map, group, order)

#relpace NA with 0's
wip_map[is.na(wip_map)] <- 0

#generate a disctrete color pallette    
pal <- c("#F7FCF5","#74C476","#41AB5D","#238B45","#006D2C","#00441B")


theme_clean <- function(base_size = 12) {
    require(grid)
    theme_grey(base_size) %+replace%
    theme(
        axis.title      =   element_blank(),
        axis.text       =   element_blank(),
        panel.background    =   element_blank(),
        panel.grid      =   element_blank(),
        axis.ticks.length   =   unit(0,"cm"),
        axis.ticks.margin   =   unit(0,"cm"),
        panel.margin    =   unit(0,"lines"),
        plot.margin     =   unit(c(0,0,0,0),"lines"),
        complete = TRUE
        )
    }

final_map <- ggplot(wip_map, aes(x=long, y=lat, group=group, fill=factor(CATEGORY))) +
            geom_polygon(colour="grey", aes(fill=factor(CATEGORY))) +
            scale_fill_manual(values=pal) +
            expand_limits(x = wip_map$long, y = wip_map$lat) +
            coord_map("polyconic") + 
            labs(fill="Number Per\nCounty") + 
                theme_clean()

 final_map + geom_path( data = state_map , colour = "red")

Thanks!

like image 729
huxley Avatar asked May 19 '13 04:05

huxley


1 Answers

Just add a geom_path to your code...

I used red to highlight the boundaries but you can easily just set it to black.

    ggplot( wip_map, aes( x = long , y = lat , group=group ) ) +
        geom_polygon( colour = "grey" , aes( fill = factor( CATEGORY ) ) ) +
        scale_fill_manual( values = pal ) +
        expand_limits( x = wip_map$long, y = wip_map$lat ) +
        coord_map( "polyconic" ) + 
        labs(fill="Number Per\nCounty") + 
        theme_clean( ) +
        geom_path( data = state_map , colour = "red")

enter image description here

like image 176
Simon O'Hanlon Avatar answered Oct 06 '22 01:10

Simon O'Hanlon