tests/testthat/test-load_param_space.R

In Neves-P/DAISIErobustness: Test the Robustness of DAISIE to Geodynamics and Traits

test_that("test load_param_space loads oceanic_ontogeny
          parameter space", {
  oceanic_ontogeny <- load_param_space(
    param_space_name = "oceanic_ontogeny_cs"
    )
  expect_equal(nrow(oceanic_ontogeny), 256)
  expect_equal(ncol(oceanic_ontogeny), 22)
  expect_equal(names(oceanic_ontogeny),
               c("time", "M", "lac", "mu", "K", "gam", "laa", "x", "d",
                 "island_ontogeny", "sea_level", "max_area", "current_area",
                 "peak_time", "total_island_age", "sea_level_amplitude",
                 "sea_level_frequency", "island_gradient_angle",
                 "extcutoff", "x_s", "x_nonend", "divdepmodel"))
  expect_true(is.numeric(oceanic_ontogeny[[1]]))
  expect_true(is.numeric(oceanic_ontogeny[[2]]))
  expect_true(is.numeric(oceanic_ontogeny[[3]]))
  expect_true(is.numeric(oceanic_ontogeny[[4]]))
  expect_true(is.numeric(oceanic_ontogeny[[5]]))
  expect_true(is.numeric(oceanic_ontogeny[[6]]))
  expect_true(is.numeric(oceanic_ontogeny[[7]]))
  expect_true(is.numeric(oceanic_ontogeny[[8]]))
  expect_true(is.numeric(oceanic_ontogeny[[9]]))
  expect_true(is.character(oceanic_ontogeny[[10]]))
  expect_true(is.character(oceanic_ontogeny[[11]]))
  expect_true(is.numeric(oceanic_ontogeny[[12]]))
  expect_true(is.numeric(oceanic_ontogeny[[13]]))
  expect_true(is.numeric(oceanic_ontogeny[[14]]))
  expect_true(is.numeric(oceanic_ontogeny[[15]]))
  expect_true(is.numeric(oceanic_ontogeny[[16]]))
  expect_true(is.numeric(oceanic_ontogeny[[17]]))
  expect_true(is.numeric(oceanic_ontogeny[[18]]))
  expect_true(is.numeric(oceanic_ontogeny[[19]]))
  expect_true(is.numeric(oceanic_ontogeny[[20]]))
  expect_true(is.numeric(oceanic_ontogeny[[21]]))
  expect_true(is.character(oceanic_ontogeny[[22]]))
})

test_that("test load_param_space loads oceanic_sea_level
          parameter space", {
  oceanic_sea_level <- load_param_space(
    param_space_name = "oceanic_sea_level_cs"
  )
  expect_equal(nrow(oceanic_sea_level), 512)
  expect_equal(ncol(oceanic_sea_level), 22)
  expect_equal(names(oceanic_sea_level),
               c("time", "M", "lac", "mu", "K", "gam", "laa", "x", "d",
                 "island_ontogeny", "sea_level", "max_area", "current_area",
                 "peak_time", "total_island_age", "sea_level_amplitude",
                 "sea_level_frequency", "island_gradient_angle",
                 "extcutoff", "x_s", "x_nonend", "divdepmodel"))
  expect_true(is.numeric(oceanic_sea_level[[1]]))
  expect_true(is.numeric(oceanic_sea_level[[2]]))
  expect_true(is.numeric(oceanic_sea_level[[3]]))
  expect_true(is.numeric(oceanic_sea_level[[4]]))
  expect_true(is.numeric(oceanic_sea_level[[5]]))
  expect_true(is.numeric(oceanic_sea_level[[6]]))
  expect_true(is.numeric(oceanic_sea_level[[7]]))
  expect_true(is.numeric(oceanic_sea_level[[8]]))
  expect_true(is.numeric(oceanic_sea_level[[9]]))
  expect_true(is.character(oceanic_sea_level[[10]]))
  expect_true(is.character(oceanic_sea_level[[11]]))
  expect_true(is.numeric(oceanic_sea_level[[12]]))
  expect_true(is.numeric(oceanic_sea_level[[13]]))
  expect_true(is.numeric(oceanic_sea_level[[14]]))
  expect_true(is.numeric(oceanic_sea_level[[15]]))
  expect_true(is.numeric(oceanic_sea_level[[16]]))
  expect_true(is.numeric(oceanic_sea_level[[17]]))
  expect_true(is.numeric(oceanic_sea_level[[18]]))
  expect_true(is.numeric(oceanic_sea_level[[19]]))
  expect_true(is.numeric(oceanic_sea_level[[20]]))
  expect_true(is.numeric(oceanic_sea_level[[21]]))
  expect_true(is.character(oceanic_sea_level[[22]]))
})

test_that("test load_param_space loads oceanic_ontogeny_sea_level
          parameter space", {
  oceanic_ontogeny_sea_level <- load_param_space(
    param_space_name = "oceanic_ontogeny_sea_level_cs"
  )
  expect_equal(nrow(oceanic_ontogeny_sea_level), 512)
  expect_equal(ncol(oceanic_ontogeny_sea_level), 22)
  expect_equal(names(oceanic_ontogeny_sea_level),
               c("time", "M", "lac", "mu", "K", "gam", "laa", "x", "d",
                 "island_ontogeny", "sea_level", "max_area", "current_area",
                 "peak_time", "total_island_age", "sea_level_amplitude",
                 "sea_level_frequency", "island_gradient_angle",
                 "extcutoff", "x_s", "x_nonend", "divdepmodel"))
  expect_true(is.numeric(oceanic_ontogeny_sea_level[[1]]))
  expect_true(is.numeric(oceanic_ontogeny_sea_level[[2]]))
  expect_true(is.numeric(oceanic_ontogeny_sea_level[[3]]))
  expect_true(is.numeric(oceanic_ontogeny_sea_level[[4]]))
  expect_true(is.numeric(oceanic_ontogeny_sea_level[[5]]))
  expect_true(is.numeric(oceanic_ontogeny_sea_level[[6]]))
  expect_true(is.numeric(oceanic_ontogeny_sea_level[[7]]))
  expect_true(is.numeric(oceanic_ontogeny_sea_level[[8]]))
  expect_true(is.numeric(oceanic_ontogeny_sea_level[[9]]))
  expect_true(is.character(oceanic_ontogeny_sea_level[[10]]))
  expect_true(is.character(oceanic_ontogeny_sea_level[[11]]))
  expect_true(is.numeric(oceanic_ontogeny_sea_level[[12]]))
  expect_true(is.numeric(oceanic_ontogeny_sea_level[[13]]))
  expect_true(is.numeric(oceanic_ontogeny_sea_level[[14]]))
  expect_true(is.numeric(oceanic_ontogeny_sea_level[[15]]))
  expect_true(is.numeric(oceanic_ontogeny_sea_level[[16]]))
  expect_true(is.numeric(oceanic_ontogeny_sea_level[[17]]))
  expect_true(is.numeric(oceanic_ontogeny_sea_level[[18]]))
  expect_true(is.numeric(oceanic_ontogeny_sea_level[[19]]))
  expect_true(is.numeric(oceanic_ontogeny_sea_level[[20]]))
  expect_true(is.numeric(oceanic_ontogeny_sea_level[[21]]))
  expect_true(is.character(oceanic_ontogeny_sea_level[[22]]))
})

test_that("test load_param_space loads continental parameter space", {
  continental <- load_param_space(
    param_space_name = "continental_cs"
  )
  expect_equal(nrow(continental), 384)
  expect_equal(ncol(continental), 10)
  expect_equal(names(continental),
               c("time", "M", "lac", "mu", "K", "gam", "laa", "x_s",
                 "x_nonend", "divdepmodel"))
  expect_true(is.numeric(continental[[1]]))
  expect_true(is.numeric(continental[[2]]))
  expect_true(is.numeric(continental[[3]]))
  expect_true(is.numeric(continental[[4]]))
  expect_true(is.numeric(continental[[5]]))
  expect_true(is.numeric(continental[[6]]))
  expect_true(is.numeric(continental[[7]]))
  expect_true(is.numeric(continental[[8]]))
  expect_true(is.numeric(continental[[9]]))
  expect_true(is.character(continental[[10]]))
})

test_that("test load_param_space loads continental_land_bridge parameter
          space", {
  continental_land_bridge <- load_param_space(
    param_space_name = "continental_land_bridge_cs"
  )
  expect_equal(nrow(continental_land_bridge), 128)
  expect_equal(ncol(continental_land_bridge), 16)
  expect_equal(names(continental_land_bridge),
               c("time", "M", "lac_1", "mu_1", "K_1", "gam_1",
                 "laa_1", "x_s", "x_nonend", "shift_times", "divdepmodel",
                 "lac_2", "mu_2", "K_2", "gam_2", "laa_2"))
  expect_true(is.numeric(continental_land_bridge[[1]]))
  expect_true(is.numeric(continental_land_bridge[[2]]))
  expect_true(is.numeric(continental_land_bridge[[3]]))
  expect_true(is.numeric(continental_land_bridge[[4]]))
  expect_true(is.numeric(continental_land_bridge[[5]]))
  expect_true(is.numeric(continental_land_bridge[[6]]))
  expect_true(is.numeric(continental_land_bridge[[7]]))
  expect_true(is.numeric(continental_land_bridge[[8]]))
  expect_true(is.numeric(continental_land_bridge[[9]]))
  expect_true(is.list(continental_land_bridge[[10]]))
  expect_true(is.character(continental_land_bridge[[11]]))
  expect_true(is.numeric(continental_land_bridge[[12]]))
  expect_true(is.numeric(continental_land_bridge[[13]]))
  expect_true(is.numeric(continental_land_bridge[[14]]))
  expect_true(is.numeric(continental_land_bridge[[15]]))
})

test_that("test load_param_space loads trait parameter space", {
  trait <- load_param_space(
    param_space_name = "trait_CES"
  )
  expect_equal(nrow(trait), 1664)
  expect_equal(ncol(trait), 14)
  expect_equal(names(trait),
               c("time", "M", "M2", "lac", "mu", "gam", "laa", "trans",
                 "trans2", "K",  "lac2", "mu2", "gam2", "laa2"))
  expect_true(is.numeric(trait[[1]]))
  expect_true(is.numeric(trait[[2]]))
  expect_true(is.numeric(trait[[3]]))
  expect_true(is.numeric(trait[[4]]))
  expect_true(is.numeric(trait[[5]]))
  expect_true(is.numeric(trait[[6]]))
  expect_true(is.numeric(trait[[7]]))
  expect_true(is.numeric(trait[[8]]))
  expect_true(is.numeric(trait[[9]]))
  expect_true(is.numeric(trait[[10]]))
  expect_true(is.numeric(trait[[11]]))
  expect_true(is.numeric(trait[[12]]))
  expect_true(is.numeric(trait[[13]]))
  expect_true(is.numeric(trait[[14]]))
})