install.packages("dismo") library(dismo) bhut<-gbif("Bhutanitis","lidderdalii") colnames(bhut) localities<-subset(bhut, select=c("lat","lon")) head(localities) b<-complete.cases(subset(bhut, select=c("lat","lon"))) localities[b,] newmap <- getMap(resolution = "low") plot(newmap) points(localities$lon, localities$lat, col = "red", cex = 2)