### See points on a map
library("maptools")
library("spocc")
data(wrld_simpl)
plot(wrld_simpl[wrld_simpl$NAME == "United States", ], xlim = c(-70, -60))
out <- occ(query = "Accipiter striatus", from = c("ecoengine", "gbif"), limit = 50)
sp_points <- occ2sp(out, just_coords = TRUE)
points(sp_points, col = 2)
Run the code above in your browser using DataLab