plot.speciesLisa: Plot positional uncertainty based on LISA

Description Usage Arguments Details Value Author(s) See Also Examples

Description

Plot the values of LISAs at species occurrence locations, which can be used to identify the locations that need positional uncertainty treatment.

Usage

1
2
3
4
5
6
## S4 method for signature 'speciesLISA,missing'
plot(x, y, ...)
## S4 method for signature 'speciesLISA,SpatialPolygons'
plot(x, y, ...)
## S4 method for signature 'speciesLISA,SpatialPolygonsDataFrame'
plot(x, y, ...)

Arguments

x

an object of class speciesLISA, which is the output of speciesLisa function.

y

optional. Boundary map of the study area, an object of class SpatialPolygons.

...

additional argument (see details)

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.

Value

plots the bubble plot.

Author(s)

Babak Naimi naimi.b@gmail.com

http://r-gis.net

See Also

speciesLisa, lisa

Examples

 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))

usdm documentation built on May 2, 2019, 5:50 p.m.