trees_3.R

zoom <- subset(trees, zipcode %in% c(11239, 11206, 11212, 11224, 11221, 
                                     11201, 11215, 11217, 11231, 11234))
zoom$hood <- as.factor(ifelse(zoom$zipcode %in% c(11201, 11215, 11217, 11231, 11234), 1,0))

map <- get_map(location = c(lon = -73.95, lat = 40.64), zoom = 12,
               maptype = "satellite", source = "google")
ggmap(map) + geom_point(data=zoom, aes(x = longitude, y = latitude, col = hood),
                        size = 0.5, shape = 16, alpha = 0.1, show.legend = F)