Description Usage Arguments Details Value Author(s) See Also Examples
Plot the values of LISAs at species occurrence locations, which can be used to identify the locations that need positional uncertainty treatment.
1 2 3 4 5 6 |
x |
an object of class |
y |
optional. Boundary map of the study area, an object of class |
... |
additional argument (see details) |
This function generates a map (i.e. a bubble plot) in which the species points present the magnitude of LISA in predictors at the location as open or filled circles with different sizes.
Below are additional arguments:
cex
the maximum symbol size (circle) in the plot.
levels
specifies the number of LISA levels at which the points are presented .
xyLegend
a vector with two numbers, specifying the coordinates of the legend. If missing, the function tries to find the appropriate location for it.
...
xlab
, ylab
and main
same as the base plot
function.
plots the bubble plot.
Babak Naimi naimi.b@gmail.com
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 | file <- system.file("external/predictors.grd", package="usdm")
r <- brick(file) # reading a RasterBrick including 5 rasters (predictor variables)
r
plot(r) # visualize the raster layers
#files.path <- system.file("external", package="usdm") # path to location of example files
#require(rgdal)
#sp <- readOGR(dsn=files.path,layer="species_nl") # reading species data (shapefile)
#splisa <- speciesLisa(x=r,y=sp,uncertainty=15000,weights=c(0.22,0.2,0.38,0.2))
#splisa
#plot(splisa)
#bnd <- readOGR(dsn=files.path,layer="boundary") # reading the boundary map
#plot(splisa,bnd)
#plot(splisa,bnd,levels=c(2,4,6,8))
#plot(splisa,bnd,levels=c(-5,-3,0,3,5))
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.