if (FALSE) {
# simulate data
sim_pus <- sim.pus(225L)
sim_spp <- terra::rast(
lapply(c("uniform", "normal", "bimodal"),
sim.species, n = 1, res = 1, x = sim_pus)
)
# calculate average for 1 species
puvspr1.dat <- calcSpeciesAverageInPus(sim_pus, sim_spp[[1]])
# calculate average for multiple species
puvspr2.dat <- calcSpeciesAverageInPus(sim_pus, sim_spp)
}
Run the code above in your browser using DataLab