# NOT RUN {
### 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("vertnet", "gbif"),
limit = 50)
xx <- occ2sp(out, just_coords = TRUE)
points(xx, col = 2)
# }
Run the code above in your browser using DataLab