Title

library("ggplot2")
library("ggmap")
library("maps")

Rural <- read.csv("~/Desktop/Rural.csv")
county_map <- map_data("county")
state_map <- map_data("state")

wip_map <- merge(county_map, Rural, by = c("region", "subregion"))

wip_map[is.na(wip_map)] <- 0

theme_clean <- function(base_size = 12) {
    require(grid)
    theme_grey(base_size)
    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)
}


ggplot(wip_map, aes(x = long, y = lat, group = group)) + geom_polygon(aes(fill = factor(population))) + 
    expand_limits(x = wip_map$long, y = wip_map$lat) + coord_map("polyconic") + 
    theme_clean() + geom_path(data = state_map, colour = "red")

## Loading required package: grid