I'm trying to plot the locations of restaurants in NYC onto the plot of NYC boroughs that I made from the shapefile provided by the NYC Taxi Data website using Rstudio. I don't understand why it is plotting the points of the restaurant locations wrong.
aoi_boundary = the shapefile provided by the NYC Taxi data website. (https://www1.nyc.gov/site/tlc/about/tlc-trip-record-data.page) rest_list = the list of restaurants with columns with the longitude and latitude coordinates where the restaurant is located
This is what I've attempted along with the attached result:
aoi_boundary_HARV <- st_read("/Users/rodrigotrujillo/Documents/College/Fall 2019/STAT 405/Project/taxi_zones/taxi_zones.shp") ggplot(data = aoi_boundary_HARV) + theme_bw() + geom_sf(size = .2, color = "black", fill = NA) + ggtitle("Yearly Pickup Location Density Map") + coord_sf() + geom_point(data = rest_list, aes(x = Longitude, y = Latitude, color = GRADE, size = 2))