## Not run:
# gb <- keepmainfields(gbifdat,ID='',Species='species',x='decimallongitude',y='decimallatitude',
# others=c('gbifid','elevation')) # Convert example data to biogeo format
# gb <- gb[ gb$Species=='Heterotheca villosa', ] # Keep data for only one species
# dem<-raster(dem,xmn=-180, xmx=180, ymn=-60, ymx=90)
# gba<-elevcheck(gb,dem,elevc="elevation",diff=50)
# ## End(Not run)
Run the code above in your browser using DataLab