data-raw/virtualSp.R

library(virtualspecies)

# Load predictors from dismo pkg
files <- list.files(path = paste0(system.file(package = "dismo"), "/ex"),
                    pattern = "grd", full.names = TRUE)
predictors <- terra::rast(files)
virtualSp <- list()

set.seed(516516)
virtual_sp <- generateRandomSp(predictors)
virtual_sp <- convertToPA(virtual_sp, beta = 0.8)

locs <- sampleOccurrences(virtual_sp, n = 3000, type = "presence-absence")
virtualSp$presence = locs$sample.points[locs$sample.points$Real == 1, 1:2][1:400, ]
virtualSp$absence = locs$sample.points[locs$sample.points$Real == 0, 1:2][1:300, ]
bg <- dismo::randomPoints(predictors[[9]], 5000)
virtualSp$background <- as.data.frame(bg)
virtualSp$pa_map <- virtual_sp$pa.raster
virtualSp$prob_map <- virtual_sp$probability.of.occurrence

# Save object in data folder
usethis::use_data(virtualSp, overwrite = TRUE)
sgvignali/SDMtune documentation built on July 20, 2023, 1:45 a.m.