ensemble_range_map {S4DM} | R Documentation |
Generate ensemble predictions from S4DM range maps
Description
This function evaluates model quality and creates an ensemble of the model outputs. This function uses 5-fold, spatially stratified, cross-validation to evaluate distribution model quality.
Usage
ensemble_range_map(
occurrences,
env,
method = NULL,
presence_method = NULL,
background_method = NULL,
bootstrap = "none",
bootstrap_reps = 100,
quantile = 0.05,
constraint_regions = NULL,
background_buffer_width = NULL,
...
)
Arguments
occurrences |
Presence coordinates in long,lat format. |
env |
Environmental SpatRaster(s) |
method |
Optional. If supplied, both presence and background density estimation will use this method. |
presence_method |
Optional. Method for estimation of presence density. |
background_method |
Optional. Method for estimation of background density. |
bootstrap |
Character. One of "none" (the default, no bootstrapping), "numbag" (presence function is bootstrapped), or "doublebag" (presence and background functions are bootstrapped). |
bootstrap_reps |
Integer. Number of bootstrap replicates to use (default is 100) |
quantile |
Quantile to use for thresholding. Default is 0.05 (5 pct training presence). Set to 0 for minimum training presence (MTP). |
constraint_regions |
See get_env_bg documentation |
background_buffer_width |
Numeric or NULL. Width (meters or map units) of buffer to use to select background environment. If NULL, uses max dist between nearest occurrences. |
... |
Additional parameters passed to internal functions. |
Details
Current plug-and-play methods include: "gaussian", "kde","vine","rangebagging", "lobagoc", and "none". Current density ratio methods include: "ulsif", "rulsif".
Value
List object containing elements (1) spatRaster ensemble layer showing the proportion of maps that are included in the range across the ensemble, (2) spatRasters for individual models, and (3) model quality information.
Note
Either method
or both presence_method
and background_method
must be supplied.
Examples
# load in sample data
library(S4DM)
library(terra)
# occurrence points
data("sample_points")
occurrences <- sample_points
# environmental data
env <- rast(system.file('ex/sample_env.tif', package="S4DM"))
# rescale the environmental data
env <- scale(env)
ensemble <- ensemble_range_map(occurrences = occurrences,
env = env,
method = NULL,
presence_method = c("gaussian", "kde"),
background_method = "gaussian",
quantile = 0.05,
background_buffer_width = 100000 )