Using mapview

inat <- read.csv("data/inat-observations.csv")

map_locations <- st_as_sf(inat, coords = c("longitude", "latitude"), crs = 4326)
ggplot() +
  geom_sf(data = map_locations)

map_data <- inat %>% 
  st_as_sf(coords = c("longitude", "latitude"),   crs = 4326, remove=FALSE) 


mapview(map_data)

This is end of the map.