tests/testthat/test-POST_FATE.graphic_evolutionPixels.R
In MayaGueguen/RFate: Vegetation modelling with the FATE model
library(RFate)
context("POST_FATE.graphic_evolutionPixels() function")
## INPUTS
test_that("POST_FATE.graphic_evolutionPixels gives error with missing data", {
expect_error(POST_FATE.graphic_evolutionPixels()
, "`name.simulation` does not exist or does not contain a PARAM_SIMUL/ folder")
expect_error(POST_FATE.graphic_evolutionPixels(NA)
, "`name.simulation` does not exist or does not contain a PARAM_SIMUL/ folder")
expect_error(POST_FATE.graphic_evolutionPixels(NULL)
, "`name.simulation` does not exist or does not contain a PARAM_SIMUL/ folder")
})
## INPUTS
test_that("POST_FATE.graphic_evolutionPixels gives error with wrong data : name.simulation", {
expect_error(POST_FATE.graphic_evolutionPixels(1)
, "`name.simulation` does not exist or does not contain a PARAM_SIMUL/ folder")
expect_error(POST_FATE.graphic_evolutionPixels("a")
, "`name.simulation` does not exist or does not contain a PARAM_SIMUL/ folder")
expect_error(POST_FATE.graphic_evolutionPixels(factor(1))
, "`name.simulation` does not exist or does not contain a PARAM_SIMUL/ folder")
expect_error(POST_FATE.graphic_evolutionPixels(matrix(seq(2), ncol=2))
, "`name.simulation` does not exist or does not contain a PARAM_SIMUL/ folder")
})
## INPUTS
test_that("POST_FATE.graphic_evolutionPixels gives error with wrong data : folders", {
if (dir.exists("FATE_simulation")) unlink("FATE_simulation", recursive = TRUE)
dir.create("FATE_simulation/")
expect_error(POST_FATE.graphic_evolutionPixels(name.simulation = "FATE_simulation")
, "`name.simulation` does not exist or does not contain a PARAM_SIMUL/ folder")
dir.create("FATE_simulation/PARAM_SIMUL/")
expect_error(POST_FATE.graphic_evolutionPixels(name.simulation = "FATE_simulation")
, "`name.simulation` does not exist or does not contain a RESULTS/ folder")
dir.create("FATE_simulation/RESULTS/")
expect_error(POST_FATE.graphic_evolutionPixels(name.simulation = "FATE_simulation")
, "`name.simulation` does not exist or does not contain a DATA/ folder")
dir.create("FATE_simulation/DATA/")
})
## INPUTS
test_that("POST_FATE.graphic_evolutionPixels gives error with wrong data : file.simulParam", {
if (dir.exists("FATE_simulation")) unlink("FATE_simulation", recursive = TRUE)
PRE_FATE.skeletonDirectory()
## TEST file.simulParam : correct content
expect_error(POST_FATE.graphic_evolutionPixels(name.simulation = "FATE_simulation")
, "The folder FATE_simulation/PARAM_SIMUL/ does not contain adequate files")
expect_error(POST_FATE.graphic_evolutionPixels(name.simulation = "FATE_simulation"
, file.simulParam = NULL)
, "The folder FATE_simulation/PARAM_SIMUL/ does not contain adequate files")
expect_error(POST_FATE.graphic_evolutionPixels(name.simulation = "FATE_simulation"
, file.simulParam = NA)
, "The folder FATE_simulation/PARAM_SIMUL/ does not contain adequate files")
expect_error(POST_FATE.graphic_evolutionPixels(name.simulation = "FATE_simulation"
, file.simulParam = "")
, "The folder FATE_simulation/PARAM_SIMUL/ does not contain adequate files")
expect_error(POST_FATE.graphic_evolutionPixels(name.simulation = "FATE_simulation"
, file.simulParam = "")
, "The folder FATE_simulation/PARAM_SIMUL/ does not contain adequate files")
## TEST file.simulParam : correct value
expect_error(POST_FATE.graphic_evolutionPixels(name.simulation = "FATE_simulation"
, file.simulParam = "toto")
, "`FATE_simulation/PARAM_SIMUL/toto` does not exist")
file.create("FATE_simulation/PARAM_SIMUL/ParamSimul.txt")
})
## INPUTS
test_that("POST_FATE.graphic_evolutionPixels gives error with wrong data : csv file (abundance)", {
if (dir.exists("FATE_simulation")) unlink("FATE_simulation", recursive = TRUE)
PRE_FATE.skeletonDirectory()
library(raster)
map_mask = raster(nrows = 5, ncols = 5)
map_mask[] = 1
writeRaster(map_mask, filename = "FATE_simulation/DATA/MASK/map_mask.tif", overwrite = TRUE)
PRE_FATE.params_globalParameters(name.simulation = "FATE_simulation"
, required.no_PFG = 6
, required.no_strata = 5
, required.simul_duration = 100
, required.seeding_duration = 10
, required.seeding_timestep = 1
, required.seeding_input = 100
, required.max_abund_low = 3000
, required.max_abund_medium = 5000
, required.max_abund_high = 9000)
PRE_FATE.params_PFGsuccession(name.simulation = "FATE_simulation"
, mat.PFG.succ = data.frame(PFG = paste0("PFG",1:6)
, type = c("C", "C", "H", "H", "P", "P")
, height = c(10, 250, 36, 68, 1250, 550)
, maturity = c(5, 5, 3, 3, 8, 9)
, longevity = c(12, 200, 25, 4, 110, 70)))
PRE_FATE.params_simulParameters(name.simulation = "FATE_simulation"
, name.MASK = "map_mask.tif")
dir.create("FATE_simulation/RESULTS/SIMUL_V1")
dir.create("FATE_simulation/RESULTS/SIMUL_V1/ABUND_perPFG_allStrata")
dir.create("FATE_simulation/RESULTS/SIMUL_V1/ABUND_perPFG_perStrata")
dir.create("FATE_simulation/RESULTS/SIMUL_V1/LIGHT")
dir.create("FATE_simulation/RESULTS/SIMUL_V1/SOIL")
## TEST RESULTS folder
expect_error(POST_FATE.graphic_evolutionPixels(name.simulation = "FATE_simulation")
, "FATE_simulation/RESULTS/POST_FATE_TABLE_PIXEL_evolution_abundance_SIMUL_V1.csv` does not exist")
})
## INPUTS
test_that("POST_FATE.graphic_evolutionPixels gives error with wrong data : csv file (light)", {
if (dir.exists("FATE_simulation")) unlink("FATE_simulation", recursive = TRUE)
PRE_FATE.skeletonDirectory()
library(raster)
map_mask = raster(nrows = 5, ncols = 5)
map_mask[] = 1
writeRaster(map_mask, filename = "FATE_simulation/DATA/MASK/map_mask.tif", overwrite = TRUE)
PRE_FATE.params_globalParameters(name.simulation = "FATE_simulation"
, required.no_PFG = 6
, required.no_strata = 5
, required.simul_duration = 100
, required.seeding_duration = 10
, required.seeding_timestep = 1
, required.seeding_input = 100
, required.max_abund_low = 3000
, required.max_abund_medium = 5000
, required.max_abund_high = 9000
, doLight = TRUE
, LIGHT.thresh_medium = 4000
, LIGHT.thresh_low = 7000)
PRE_FATE.params_PFGsuccession(name.simulation = "FATE_simulation"
, mat.PFG.succ = data.frame(PFG = paste0("PFG",1:6)
, type = c("C", "C", "H", "H", "P", "P")
, height = c(10, 250, 36, 68, 1250, 550)
, maturity = c(5, 5, 3, 3, 8, 9)
, longevity = c(12, 200, 25, 4, 110, 70)))
PRE_FATE.params_PFGlight(name.simulation = "FATE_simulation"
, mat.PFG.light = data.frame(PFG = paste0("PFG",1:6)
, type = c("C", "C", "H", "H", "P", "P"))
, mat.PFG.tol = data.frame(PFG = paste0("PFG",1:6)
, strategy_tol = c("ubiquist", "undergrowth"
, "ubiquist", "semi_shade"
, "pioneer", "full_light")))
PRE_FATE.params_simulParameters(name.simulation = "FATE_simulation"
, name.MASK = "map_mask.tif")
dir.create("FATE_simulation/RESULTS/SIMUL_V1")
dir.create("FATE_simulation/RESULTS/SIMUL_V1/ABUND_perPFG_allStrata")
dir.create("FATE_simulation/RESULTS/SIMUL_V1/ABUND_perPFG_perStrata")
dir.create("FATE_simulation/RESULTS/SIMUL_V1/LIGHT")
dir.create("FATE_simulation/RESULTS/SIMUL_V1/SOIL")
cat("PFG ID.pixel X Y HAB 1\nPFG1 1 -144 72 1 10"
, file = "FATE_simulation/RESULTS/POST_FATE_TABLE_PIXEL_evolution_abundance_SIMUL_V1.csv")
## TEST RESULTS folder
expect_error(POST_FATE.graphic_evolutionPixels(name.simulation = "FATE_simulation")
, "FATE_simulation/RESULTS/POST_FATE_TABLE_PIXEL_evolution_light_SIMUL_V1.csv` does not exist")
})
## INPUTS
test_that("POST_FATE.graphic_evolutionPixels gives error with wrong data : csv file (soil)", {
if (dir.exists("FATE_simulation")) unlink("FATE_simulation", recursive = TRUE)
PRE_FATE.skeletonDirectory()
library(raster)
map_mask = raster(nrows = 5, ncols = 5)
map_mask[] = 1
writeRaster(map_mask, filename = "FATE_simulation/DATA/MASK/map_mask.tif", overwrite = TRUE)
PRE_FATE.params_globalParameters(name.simulation = "FATE_simulation"
, required.no_PFG = 6
, required.no_strata = 5
, required.simul_duration = 100
, required.seeding_duration = 10
, required.seeding_timestep = 1
, required.seeding_input = 100
, required.max_abund_low = 3000
, required.max_abund_medium = 5000
, required.max_abund_high = 9000
, doSoil = TRUE
, SOIL.init = 2.5
, SOIL.retention = 0.5)
PRE_FATE.params_PFGsuccession(name.simulation = "FATE_simulation"
, mat.PFG.succ = data.frame(PFG = paste0("PFG",1:6)
, type = c("C", "C", "H", "H", "P", "P")
, height = c(10, 250, 36, 68, 1250, 550)
, maturity = c(5, 5, 3, 3, 8, 9)
, longevity = c(12, 200, 25, 4, 110, 70)))
PRE_FATE.params_PFGsoil(name.simulation = "FATE_simulation"
, mat.PFG.soil = data.frame(PFG = paste0("PFG",1:6)
, type = c("C", "C", "H", "H", "P", "P")
, soil_contrib = c(2, 2.5, 3.5, 3, 1.5, 3)
, soil_tol_min = c(1, 2, 2, 1.5, 1, 2)
, soil_tol_max = c(3, 3, 5, 4.5, 2, 4)))
PRE_FATE.params_simulParameters(name.simulation = "FATE_simulation"
, name.MASK = "map_mask.tif")
dir.create("FATE_simulation/RESULTS/SIMUL_V1")
dir.create("FATE_simulation/RESULTS/SIMUL_V1/ABUND_perPFG_allStrata")
dir.create("FATE_simulation/RESULTS/SIMUL_V1/ABUND_perPFG_perStrata")
dir.create("FATE_simulation/RESULTS/SIMUL_V1/LIGHT")
dir.create("FATE_simulation/RESULTS/SIMUL_V1/SOIL")
cat("PFG ID.pixel X Y HAB 1\nPFG1 1 -144 72 1 10"
, file = "FATE_simulation/RESULTS/POST_FATE_TABLE_PIXEL_evolution_abundance_SIMUL_V1.csv")
cat("PFG STRATUM ID.pixel X Y HAB 1\nPFG1 1 1 -144 72 1 10"
, file = "FATE_simulation/RESULTS/POST_FATE_TABLE_PIXEL_evolution_light_SIMUL_V1.csv")
## TEST RESULTS folder
expect_error(POST_FATE.graphic_evolutionPixels(name.simulation = "FATE_simulation")
, "FATE_simulation/RESULTS/POST_FATE_TABLE_PIXEL_evolution_soil_SIMUL_V1.csv` does not exist")
})
## OUTPUTS
test_that("POST_FATE.graphic_evolutionPixels gives correct outputs :", {
{
if (dir.exists("FATE_simulation")) unlink("FATE_simulation", recursive = TRUE)
PRE_FATE.skeletonDirectory()
library(raster)
map_0 = raster(nrows = 5, ncols = 5, resolution = 50)
map_0[] = 0
map_1 = map_0
map_1[] = 1
map_hab = map_0
map_hab[] = sample(c(1, 5, 10), ncell(map_hab), replace = TRUE)
writeRaster(map_1, filename = "FATE_simulation/DATA/MASK/map_mask.tif", overwrite = TRUE)
writeRaster(map_hab, filename = "FATE_simulation/DATA/MASK/map_hab.tif", overwrite = TRUE)
PRE_FATE.params_globalParameters(name.simulation = "FATE_simulation"
, required.no_PFG = 6
, required.no_strata = 5
, required.simul_duration = 100
, required.seeding_duration = 10
, required.seeding_timestep = 1
, required.seeding_input = 100
, required.max_abund_low = 3000
, required.max_abund_medium = 5000
, required.max_abund_high = 9000
, doLight = TRUE
, LIGHT.thresh_medium = 4000
, LIGHT.thresh_low = 7000
, doSoil = TRUE
, SOIL.init = 2.5
, SOIL.retention = 0.5)
PRE_FATE.params_PFGsuccession(name.simulation = "FATE_simulation"
, mat.PFG.succ = data.frame(PFG = paste0("PFG",1:6)
, type = c("C", "C", "H", "H", "P", "P")
, height = c(10, 250, 36, 68, 1250, 550)
, maturity = c(5, 5, 3, 3, 8, 9)
, longevity = c(12, 200, 25, 4, 110, 70)))
PRE_FATE.params_PFGlight(name.simulation = "FATE_simulation"
, mat.PFG.light = data.frame(PFG = paste0("PFG",1:6)
, type = c("C", "C", "H", "H", "P", "P"))
, mat.PFG.tol = data.frame(PFG = paste0("PFG",1:6)
, strategy_tol = c("ubiquist", "undergrowth"
, "ubiquist", "semi_shade"
, "pioneer", "full_light")))
PRE_FATE.params_PFGsoil(name.simulation = "FATE_simulation"
, mat.PFG.soil = data.frame(PFG = paste0("PFG",1:6)
, type = c("C", "C", "H", "H", "P", "P")
, soil_contrib = c(2, 2.5, 3.5, 3, 1.5, 3)
, soil_tol_min = c(1, 2, 2, 1.5, 1, 2)
, soil_tol_max = c(3, 3, 5, 4.5, 2, 4)))
PRE_FATE.params_simulParameters(name.simulation = "FATE_simulation"
, name.MASK = "map_mask.tif")
dir.create("FATE_simulation/RESULTS/SIMUL_V1")
dir.create("FATE_simulation/RESULTS/SIMUL_V1/ABUND_perPFG_allStrata")
dir.create("FATE_simulation/RESULTS/SIMUL_V1/ABUND_perPFG_perStrata")
dir.create("FATE_simulation/RESULTS/SIMUL_V1/LIGHT")
dir.create("FATE_simulation/RESULTS/SIMUL_V1/SOIL")
for (i in 1:6)
{
writeRaster(map_1, filename = paste0("FATE_simulation/RESULTS/SIMUL_V1/ABUND_perPFG_allStrata/"
, "Abund_YEAR_1_PFG", i, "_STRATA_all.tif"), overwrite = TRUE)
}
writeRaster(map_1, filename = paste0("FATE_simulation/RESULTS/SIMUL_V1/LIGHT/"
, "Light_Resources_YEAR_1_STRATA_1.tif"), overwrite = TRUE)
writeRaster(map_1, filename = paste0("FATE_simulation/RESULTS/SIMUL_V1/SOIL/"
, "Soil_Resources_YEAR_1.tif"), overwrite = TRUE)
tempEvol = POST_FATE.temporalEvolution(name.simulation = "FATE_simulation"
, no_years = 10
, opt.ras_habitat = "FATE_simulation/DATA/MASK/map_hab.tif")
}
expect_message(POST_FATE.graphic_evolutionPixels(name.simulation = "FATE_simulation")
, "> POST_FATE_TABLE_PIXEL_evolution_pixels_")
expect_warning(POST_FATE.graphic_evolutionPixels(name.simulation = "FATE_simulation"
, opt.cells_ID = 42
, opt.doPlot = FALSE)
, "The values given in `opt.cells_ID` do not match with any cells of the studied area")
plotPix = POST_FATE.graphic_evolutionPixels(name.simulation = "FATE_simulation"
, opt.cells_ID = c(1,2,3)
, opt.fixedScale = FALSE)
expect_output(str(plotPix), "List")
expect_equal(length(plotPix), 1)
expect_output(str(plotPix[[1]]), "List")
expect_equal(length(plotPix[[1]]), 2)
expect_output(str(plotPix[[1]]$tab), "data.frame")
expect_output(str(plotPix[[1]]$tab), "6 variables")
expect_output(str(plotPix[[1]]$plot), "List")
})
MayaGueguen/RFate documentation built on Oct. 17, 2020, 8:06 a.m.
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
library(RFate)
context("POST_FATE.graphic_evolutionPixels() function")
## INPUTS
test_that("POST_FATE.graphic_evolutionPixels gives error with missing data", {
expect_error(POST_FATE.graphic_evolutionPixels()
, "`name.simulation` does not exist or does not contain a PARAM_SIMUL/ folder")
expect_error(POST_FATE.graphic_evolutionPixels(NA)
, "`name.simulation` does not exist or does not contain a PARAM_SIMUL/ folder")
expect_error(POST_FATE.graphic_evolutionPixels(NULL)
, "`name.simulation` does not exist or does not contain a PARAM_SIMUL/ folder")
})
## INPUTS
test_that("POST_FATE.graphic_evolutionPixels gives error with wrong data : name.simulation", {
expect_error(POST_FATE.graphic_evolutionPixels(1)
, "`name.simulation` does not exist or does not contain a PARAM_SIMUL/ folder")
expect_error(POST_FATE.graphic_evolutionPixels("a")
, "`name.simulation` does not exist or does not contain a PARAM_SIMUL/ folder")
expect_error(POST_FATE.graphic_evolutionPixels(factor(1))
, "`name.simulation` does not exist or does not contain a PARAM_SIMUL/ folder")
expect_error(POST_FATE.graphic_evolutionPixels(matrix(seq(2), ncol=2))
, "`name.simulation` does not exist or does not contain a PARAM_SIMUL/ folder")
})
## INPUTS
test_that("POST_FATE.graphic_evolutionPixels gives error with wrong data : folders", {
if (dir.exists("FATE_simulation")) unlink("FATE_simulation", recursive = TRUE)
dir.create("FATE_simulation/")
expect_error(POST_FATE.graphic_evolutionPixels(name.simulation = "FATE_simulation")
, "`name.simulation` does not exist or does not contain a PARAM_SIMUL/ folder")
dir.create("FATE_simulation/PARAM_SIMUL/")
expect_error(POST_FATE.graphic_evolutionPixels(name.simulation = "FATE_simulation")
, "`name.simulation` does not exist or does not contain a RESULTS/ folder")
dir.create("FATE_simulation/RESULTS/")
expect_error(POST_FATE.graphic_evolutionPixels(name.simulation = "FATE_simulation")
, "`name.simulation` does not exist or does not contain a DATA/ folder")
dir.create("FATE_simulation/DATA/")
})
## INPUTS
test_that("POST_FATE.graphic_evolutionPixels gives error with wrong data : file.simulParam", {
if (dir.exists("FATE_simulation")) unlink("FATE_simulation", recursive = TRUE)
PRE_FATE.skeletonDirectory()
## TEST file.simulParam : correct content
expect_error(POST_FATE.graphic_evolutionPixels(name.simulation = "FATE_simulation")
, "The folder FATE_simulation/PARAM_SIMUL/ does not contain adequate files")
expect_error(POST_FATE.graphic_evolutionPixels(name.simulation = "FATE_simulation"
, file.simulParam = NULL)
, "The folder FATE_simulation/PARAM_SIMUL/ does not contain adequate files")
expect_error(POST_FATE.graphic_evolutionPixels(name.simulation = "FATE_simulation"
, file.simulParam = NA)
, "The folder FATE_simulation/PARAM_SIMUL/ does not contain adequate files")
expect_error(POST_FATE.graphic_evolutionPixels(name.simulation = "FATE_simulation"
, file.simulParam = "")
, "The folder FATE_simulation/PARAM_SIMUL/ does not contain adequate files")
expect_error(POST_FATE.graphic_evolutionPixels(name.simulation = "FATE_simulation"
, file.simulParam = "")
, "The folder FATE_simulation/PARAM_SIMUL/ does not contain adequate files")
## TEST file.simulParam : correct value
expect_error(POST_FATE.graphic_evolutionPixels(name.simulation = "FATE_simulation"
, file.simulParam = "toto")
, "`FATE_simulation/PARAM_SIMUL/toto` does not exist")
file.create("FATE_simulation/PARAM_SIMUL/ParamSimul.txt")
})
## INPUTS
test_that("POST_FATE.graphic_evolutionPixels gives error with wrong data : csv file (abundance)", {
if (dir.exists("FATE_simulation")) unlink("FATE_simulation", recursive = TRUE)
PRE_FATE.skeletonDirectory()
library(raster)
map_mask = raster(nrows = 5, ncols = 5)
map_mask[] = 1
writeRaster(map_mask, filename = "FATE_simulation/DATA/MASK/map_mask.tif", overwrite = TRUE)
PRE_FATE.params_globalParameters(name.simulation = "FATE_simulation"
, required.no_PFG = 6
, required.no_strata = 5
, required.simul_duration = 100
, required.seeding_duration = 10
, required.seeding_timestep = 1
, required.seeding_input = 100
, required.max_abund_low = 3000
, required.max_abund_medium = 5000
, required.max_abund_high = 9000)
PRE_FATE.params_PFGsuccession(name.simulation = "FATE_simulation"
, mat.PFG.succ = data.frame(PFG = paste0("PFG",1:6)
, type = c("C", "C", "H", "H", "P", "P")
, height = c(10, 250, 36, 68, 1250, 550)
, maturity = c(5, 5, 3, 3, 8, 9)
, longevity = c(12, 200, 25, 4, 110, 70)))
PRE_FATE.params_simulParameters(name.simulation = "FATE_simulation"
, name.MASK = "map_mask.tif")
dir.create("FATE_simulation/RESULTS/SIMUL_V1")
dir.create("FATE_simulation/RESULTS/SIMUL_V1/ABUND_perPFG_allStrata")
dir.create("FATE_simulation/RESULTS/SIMUL_V1/ABUND_perPFG_perStrata")
dir.create("FATE_simulation/RESULTS/SIMUL_V1/LIGHT")
dir.create("FATE_simulation/RESULTS/SIMUL_V1/SOIL")
## TEST RESULTS folder
expect_error(POST_FATE.graphic_evolutionPixels(name.simulation = "FATE_simulation")
, "FATE_simulation/RESULTS/POST_FATE_TABLE_PIXEL_evolution_abundance_SIMUL_V1.csv` does not exist")
})
## INPUTS
test_that("POST_FATE.graphic_evolutionPixels gives error with wrong data : csv file (light)", {
if (dir.exists("FATE_simulation")) unlink("FATE_simulation", recursive = TRUE)
PRE_FATE.skeletonDirectory()
library(raster)
map_mask = raster(nrows = 5, ncols = 5)
map_mask[] = 1
writeRaster(map_mask, filename = "FATE_simulation/DATA/MASK/map_mask.tif", overwrite = TRUE)
PRE_FATE.params_globalParameters(name.simulation = "FATE_simulation"
, required.no_PFG = 6
, required.no_strata = 5
, required.simul_duration = 100
, required.seeding_duration = 10
, required.seeding_timestep = 1
, required.seeding_input = 100
, required.max_abund_low = 3000
, required.max_abund_medium = 5000
, required.max_abund_high = 9000
, doLight = TRUE
, LIGHT.thresh_medium = 4000
, LIGHT.thresh_low = 7000)
PRE_FATE.params_PFGsuccession(name.simulation = "FATE_simulation"
, mat.PFG.succ = data.frame(PFG = paste0("PFG",1:6)
, type = c("C", "C", "H", "H", "P", "P")
, height = c(10, 250, 36, 68, 1250, 550)
, maturity = c(5, 5, 3, 3, 8, 9)
, longevity = c(12, 200, 25, 4, 110, 70)))
PRE_FATE.params_PFGlight(name.simulation = "FATE_simulation"
, mat.PFG.light = data.frame(PFG = paste0("PFG",1:6)
, type = c("C", "C", "H", "H", "P", "P"))
, mat.PFG.tol = data.frame(PFG = paste0("PFG",1:6)
, strategy_tol = c("ubiquist", "undergrowth"
, "ubiquist", "semi_shade"
, "pioneer", "full_light")))
PRE_FATE.params_simulParameters(name.simulation = "FATE_simulation"
, name.MASK = "map_mask.tif")
dir.create("FATE_simulation/RESULTS/SIMUL_V1")
dir.create("FATE_simulation/RESULTS/SIMUL_V1/ABUND_perPFG_allStrata")
dir.create("FATE_simulation/RESULTS/SIMUL_V1/ABUND_perPFG_perStrata")
dir.create("FATE_simulation/RESULTS/SIMUL_V1/LIGHT")
dir.create("FATE_simulation/RESULTS/SIMUL_V1/SOIL")
cat("PFG ID.pixel X Y HAB 1\nPFG1 1 -144 72 1 10"
, file = "FATE_simulation/RESULTS/POST_FATE_TABLE_PIXEL_evolution_abundance_SIMUL_V1.csv")
## TEST RESULTS folder
expect_error(POST_FATE.graphic_evolutionPixels(name.simulation = "FATE_simulation")
, "FATE_simulation/RESULTS/POST_FATE_TABLE_PIXEL_evolution_light_SIMUL_V1.csv` does not exist")
})
## INPUTS
test_that("POST_FATE.graphic_evolutionPixels gives error with wrong data : csv file (soil)", {
if (dir.exists("FATE_simulation")) unlink("FATE_simulation", recursive = TRUE)
PRE_FATE.skeletonDirectory()
library(raster)
map_mask = raster(nrows = 5, ncols = 5)
map_mask[] = 1
writeRaster(map_mask, filename = "FATE_simulation/DATA/MASK/map_mask.tif", overwrite = TRUE)
PRE_FATE.params_globalParameters(name.simulation = "FATE_simulation"
, required.no_PFG = 6
, required.no_strata = 5
, required.simul_duration = 100
, required.seeding_duration = 10
, required.seeding_timestep = 1
, required.seeding_input = 100
, required.max_abund_low = 3000
, required.max_abund_medium = 5000
, required.max_abund_high = 9000
, doSoil = TRUE
, SOIL.init = 2.5
, SOIL.retention = 0.5)
PRE_FATE.params_PFGsuccession(name.simulation = "FATE_simulation"
, mat.PFG.succ = data.frame(PFG = paste0("PFG",1:6)
, type = c("C", "C", "H", "H", "P", "P")
, height = c(10, 250, 36, 68, 1250, 550)
, maturity = c(5, 5, 3, 3, 8, 9)
, longevity = c(12, 200, 25, 4, 110, 70)))
PRE_FATE.params_PFGsoil(name.simulation = "FATE_simulation"
, mat.PFG.soil = data.frame(PFG = paste0("PFG",1:6)
, type = c("C", "C", "H", "H", "P", "P")
, soil_contrib = c(2, 2.5, 3.5, 3, 1.5, 3)
, soil_tol_min = c(1, 2, 2, 1.5, 1, 2)
, soil_tol_max = c(3, 3, 5, 4.5, 2, 4)))
PRE_FATE.params_simulParameters(name.simulation = "FATE_simulation"
, name.MASK = "map_mask.tif")
dir.create("FATE_simulation/RESULTS/SIMUL_V1")
dir.create("FATE_simulation/RESULTS/SIMUL_V1/ABUND_perPFG_allStrata")
dir.create("FATE_simulation/RESULTS/SIMUL_V1/ABUND_perPFG_perStrata")
dir.create("FATE_simulation/RESULTS/SIMUL_V1/LIGHT")
dir.create("FATE_simulation/RESULTS/SIMUL_V1/SOIL")
cat("PFG ID.pixel X Y HAB 1\nPFG1 1 -144 72 1 10"
, file = "FATE_simulation/RESULTS/POST_FATE_TABLE_PIXEL_evolution_abundance_SIMUL_V1.csv")
cat("PFG STRATUM ID.pixel X Y HAB 1\nPFG1 1 1 -144 72 1 10"
, file = "FATE_simulation/RESULTS/POST_FATE_TABLE_PIXEL_evolution_light_SIMUL_V1.csv")
## TEST RESULTS folder
expect_error(POST_FATE.graphic_evolutionPixels(name.simulation = "FATE_simulation")
, "FATE_simulation/RESULTS/POST_FATE_TABLE_PIXEL_evolution_soil_SIMUL_V1.csv` does not exist")
})
## OUTPUTS
test_that("POST_FATE.graphic_evolutionPixels gives correct outputs :", {
{
if (dir.exists("FATE_simulation")) unlink("FATE_simulation", recursive = TRUE)
PRE_FATE.skeletonDirectory()
library(raster)
map_0 = raster(nrows = 5, ncols = 5, resolution = 50)
map_0[] = 0
map_1 = map_0
map_1[] = 1
map_hab = map_0
map_hab[] = sample(c(1, 5, 10), ncell(map_hab), replace = TRUE)
writeRaster(map_1, filename = "FATE_simulation/DATA/MASK/map_mask.tif", overwrite = TRUE)
writeRaster(map_hab, filename = "FATE_simulation/DATA/MASK/map_hab.tif", overwrite = TRUE)
PRE_FATE.params_globalParameters(name.simulation = "FATE_simulation"
, required.no_PFG = 6
, required.no_strata = 5
, required.simul_duration = 100
, required.seeding_duration = 10
, required.seeding_timestep = 1
, required.seeding_input = 100
, required.max_abund_low = 3000
, required.max_abund_medium = 5000
, required.max_abund_high = 9000
, doLight = TRUE
, LIGHT.thresh_medium = 4000
, LIGHT.thresh_low = 7000
, doSoil = TRUE
, SOIL.init = 2.5
, SOIL.retention = 0.5)
PRE_FATE.params_PFGsuccession(name.simulation = "FATE_simulation"
, mat.PFG.succ = data.frame(PFG = paste0("PFG",1:6)
, type = c("C", "C", "H", "H", "P", "P")
, height = c(10, 250, 36, 68, 1250, 550)
, maturity = c(5, 5, 3, 3, 8, 9)
, longevity = c(12, 200, 25, 4, 110, 70)))
PRE_FATE.params_PFGlight(name.simulation = "FATE_simulation"
, mat.PFG.light = data.frame(PFG = paste0("PFG",1:6)
, type = c("C", "C", "H", "H", "P", "P"))
, mat.PFG.tol = data.frame(PFG = paste0("PFG",1:6)
, strategy_tol = c("ubiquist", "undergrowth"
, "ubiquist", "semi_shade"
, "pioneer", "full_light")))
PRE_FATE.params_PFGsoil(name.simulation = "FATE_simulation"
, mat.PFG.soil = data.frame(PFG = paste0("PFG",1:6)
, type = c("C", "C", "H", "H", "P", "P")
, soil_contrib = c(2, 2.5, 3.5, 3, 1.5, 3)
, soil_tol_min = c(1, 2, 2, 1.5, 1, 2)
, soil_tol_max = c(3, 3, 5, 4.5, 2, 4)))
PRE_FATE.params_simulParameters(name.simulation = "FATE_simulation"
, name.MASK = "map_mask.tif")
dir.create("FATE_simulation/RESULTS/SIMUL_V1")
dir.create("FATE_simulation/RESULTS/SIMUL_V1/ABUND_perPFG_allStrata")
dir.create("FATE_simulation/RESULTS/SIMUL_V1/ABUND_perPFG_perStrata")
dir.create("FATE_simulation/RESULTS/SIMUL_V1/LIGHT")
dir.create("FATE_simulation/RESULTS/SIMUL_V1/SOIL")
for (i in 1:6)
{
writeRaster(map_1, filename = paste0("FATE_simulation/RESULTS/SIMUL_V1/ABUND_perPFG_allStrata/"
, "Abund_YEAR_1_PFG", i, "_STRATA_all.tif"), overwrite = TRUE)
}
writeRaster(map_1, filename = paste0("FATE_simulation/RESULTS/SIMUL_V1/LIGHT/"
, "Light_Resources_YEAR_1_STRATA_1.tif"), overwrite = TRUE)
writeRaster(map_1, filename = paste0("FATE_simulation/RESULTS/SIMUL_V1/SOIL/"
, "Soil_Resources_YEAR_1.tif"), overwrite = TRUE)
tempEvol = POST_FATE.temporalEvolution(name.simulation = "FATE_simulation"
, no_years = 10
, opt.ras_habitat = "FATE_simulation/DATA/MASK/map_hab.tif")
}
expect_message(POST_FATE.graphic_evolutionPixels(name.simulation = "FATE_simulation")
, "> POST_FATE_TABLE_PIXEL_evolution_pixels_")
expect_warning(POST_FATE.graphic_evolutionPixels(name.simulation = "FATE_simulation"
, opt.cells_ID = 42
, opt.doPlot = FALSE)
, "The values given in `opt.cells_ID` do not match with any cells of the studied area")
plotPix = POST_FATE.graphic_evolutionPixels(name.simulation = "FATE_simulation"
, opt.cells_ID = c(1,2,3)
, opt.fixedScale = FALSE)
expect_output(str(plotPix), "List")
expect_equal(length(plotPix), 1)
expect_output(str(plotPix[[1]]), "List")
expect_equal(length(plotPix[[1]]), 2)
expect_output(str(plotPix[[1]]$tab), "data.frame")
expect_output(str(plotPix[[1]]$tab), "6 variables")
expect_output(str(plotPix[[1]]$plot), "List")
})
R Package Documentation
Browse R Packages
We want your feedback!
Note that we can't provide technical support on individual packages. You should contact the package authors for that.
Embedding an R snippet on your website
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.