vis_bioclimPlot Visualize bivariate plot of BIOCLIM model
vis_bioclimPlot.RdThis functions creates a bivariate plot with two of the environmental variables used for modeling as x and y axes and occurrences as observations.
Arguments
- x
bioclim model including values for each environmental layer at each occurrence point
- a
numeric Environmental layer to be used as x axis. Default is layer 1.
- b
numeric. Environmental layer to be used as x axis. Default is layer 2.
- p
numeric. (0-1) percentile distribution to be used for plotting envelope and showing points outside of envelope. Default is 0.9
Value
A bivariate plot of environmental values for occurrences. Includes a blue rectangle representing the bioclimatic envelope given p. Occurrences that are inside the envelope for all layers (included those not plotted) are shown as green circles and those outside of the envelope for one ore more variables are plotted as orange triangles.
Details
This is a bivariate plot with x and y axes representing two of the environmental layers used for modeling (user selected although 1 and 2 as default). Occurrences used for modeling are shown with differential visualization if they are outside of the selected percentile distribution (for any variable). Plot also includes a rectangle representing the bivariate bioclimatic envelope according to a provided percentile.
Examples
if (FALSE) { # \dontrun{
envs <- envs_userEnvs(rasPath = list.files(system.file("extdata/wc",
package = "wallace"),
pattern = ".tif$", full.names = TRUE),
rasName = list.files(system.file("extdata/wc",
package = "wallace"),
pattern = ".tif$", full.names = FALSE))
occs <- read.csv(system.file("extdata/Bassaricyon_alleni.csv",
package = "wallace"))
bg <- read.csv(system.file("extdata/Bassaricyon_alleni_bgPoints.csv",
package = "wallace"))
partblock <- part_partitionOccs(occs, bg, method = 'block')
m <- model_bioclim(occs, bg, partblock, envs)
bioclimPlot <- vis_bioclimPlot(x = m@models$bioclim,
a = 1, b = 2, p = 1)
} # }