tests/testthat/test-treeselection.R

In VincyaneBadouard/LoggingLab: An R Package to Simulate Forest Logging

test_that("treeselection", {

  # Data loading
  data(Paracou6_2016)
  data(DTMParacou)
  data(SpeciesCriteria)
  data(MainTrails)
  data(HarvestableAreaOutputsCable)

  Paracou6_2016 <- dplyr::slice(Paracou6_2016, 1:2000)
  MatrixInventory <- as.matrix(Paracou6_2016)
  Matrixspeciescriteria <- as.matrix(SpeciesCriteria)

  # Check the function arguments
  expect_error(treeselection(MatrixInventory,
                             scenario ="manual", fuel = "0", objective = 20,
                             diversification = TRUE, specieslax = FALSE, objectivelax = FALSE,
                             speciescriteria = Matrixspeciescriteria),
               regexp = "The 'inventory' and 'speciescriteria' arguments of the 'treeselection' function must be data.frame")

  expect_error(treeselection(Paracou6_2016, speciescriteria = SpeciesCriteria,
                             scenario ="manual", fuel = "0", diversification = TRUE,
                             objective = "RIL"),
               regexp = "The 'objective' argument of the 'treeselection' function must be numeric")

  expect_error(treeselection(Paracou6_2016, speciescriteria = SpeciesCriteria,
                             scenario ="manual", fuel = "0", objective = 20,
                             diversification = "1"),
               regexp = "The 'diversification' argument of the 'treeselection' function must be logical or NULL")

  expect_error(treeselection(Paracou6_2016, speciescriteria = SpeciesCriteria,
                             scenario ="manual", fuel = "0", objective = 20,
                             diversification = FALSE, specieslax = 2, objectivelax = "a"),
               regexp = "The 'specieslax' and 'objectivelax' arguments of the 'treeselection' function must be logicals")

  expect_error(treeselection(Paracou6_2016, speciescriteria = SpeciesCriteria,
                             fuel = "0", diversification = TRUE, objective = 20,
                             scenario = "CL"),
               regexp = "The 'scenario' argument of the 'treeselection' function must be
         'RIL1', 'RIL2broken', 'RIL2', 'RIL3', 'RIL3fuel', 'RIL3fuelhollow' or 'manual'")

  expect_error(treeselection(Paracou6_2016, speciescriteria = SpeciesCriteria,
                             scenario ="manual", objective = 20, diversification = TRUE,
                             fuel = TRUE),
               regexp = "The 'fuel' argument of the 'treeselection' function must be '0', '1', or '2'")

  expect_error(treeselection(Paracou6_2016, speciescriteria = SpeciesCriteria,
                             objective = 20, scenario ="manual", fuel = "2",
                             diversification = FALSE, specieslax = FALSE, objectivelax = FALSE,
                             topography = NULL, plotslope = NULL,maintrails = NULL))


  expect_error(treeselection(Paracou6_2016, speciescriteria = SpeciesCriteria,
                             objective = 20, scenario ="manual", fuel = "0", diversification = TRUE,
                             topography = DTMParacou, plotslope = HarvestableAreaOutputsCable$PlotSlope,
                             maintrails = MainTrails,harvestablearea = HarvestableAreaOutputsCable$HarvestableArea,
                             harvestablepolygons = HarvestableAreaOutputsCable$HarvestablePolygons,
                             advancedloggingparameters = as.matrix(loggingparameters())),
               regexp = "The 'advancedloggingparameters' argument of the 'treeselection' function must be a list")

  expect_error(treeselection(Paracou6_2016, speciescriteria = SpeciesCriteria,
                             scenario = "manual",
                             objective = 20, fuel = NULL, diversification = T,
                             topography = DTMParacou, plotslope = HarvestableAreaOutputsCable$PlotSlope,
                             maintrails = MainTrails,harvestablepolygons = HarvestableAreaOutputsCable$HarvestablePolygons),
               regexp = "If you choose the 'manual' mode, you must fill in the arguments
         'objective', 'fuel' , 'winching' and 'diversification'")

  # Test data preparation
  inventory = addtreedim(cleaninventory(Paracou6_2016, PlotMask), volumeparameters = ForestZoneVolumeParametersTable)

  testinventory <- suppressMessages(treeselection(inventory,
                                                  topography = DTMParacou,
                                                  speciescriteria = SpeciesCriteria,
                                                  scenario = "manual", objective = 20, fuel = "2", diversification = TRUE,
                                                  winching = "0", specieslax = FALSE, objectivelax = TRUE,
                                                  harvestablearea = HarvestableAreaOutputsCable$HarvestableArea,
                                                  plotslope = HarvestableAreaOutputsCable$PlotSlope,
                                                  maintrails = MainTrails,
                                                  harvestablepolygons = HarvestableAreaOutputsCable$HarvestablePolygons,
                                                  advancedloggingparameters = loggingparameters())$inventory)

  # All alived
  if ("DeathCause" %in% names(inventory)) expect_true(all(is.null(DeathCause)))

  # Objective Volume:
  objective <- 20
  expected <- objective * HarvestableAreaOutputsCable$HarvestableArea
  VO <- suppressMessages(treeselection(inventory,
                                       topography = DTMParacou,
                                       speciescriteria = SpeciesCriteria,
                                       scenario = "manual", objective = objective, fuel = "2", winching = "0",
                                       diversification = TRUE, specieslax = FALSE, objectivelax = TRUE,
                                       harvestablearea = HarvestableAreaOutputsCable$HarvestableArea,
                                       plotslope = HarvestableAreaOutputsCable$PlotSlope,
                                       maintrails = MainTrails,
                                       harvestablepolygons = HarvestableAreaOutputsCable$HarvestablePolygons
                                       ))$VO


  expect_true(VO == expected)

  expect_true(nrow(inventory) == nrow(testinventory)) # check that the rows number don't change

})

# check args
# treeselection(inventory,
#               scenario ="manual", fuel = "0",objective = 20, diversification = TRUE, specieslax = FALSE,
#               objectivelax = TRUE, SpeciesCriteria, advancedloggingparameters = loggingparameters())
# inventory, speciescriteria (data.frame)
# objective (numeric)
# scenario, fuel  (character)
# diversification, specieslax, objectivelax (logical)
# advancedloggingparameters (liste)

# All alived: if ("DeathCause" %in% names(inventory)) expect_true(all(is.null(DeathCause)))
# Objective Volume: if (scenario == "RIL3fuelhollow"| (scenario == "manual"& fuel =="2")) {expect_true(VO = objective)
# }else{
# expect_true(VO = objective + advancedloggingparameters$ObjectiveBonus)