tests/testthat/test-clado_event.R

In joshwlambert/DAISIEmainland: Simulates Data for Phylogenetic Data on Islands for DAISIE with the Addition of Mainland Dynamics

test_that("clado_event produces correct output on non-cladogenesis species", {
  set.seed(
    1,
    kind = "Mersenne-Twister",
    normal.kind = "Inversion",
    sample.kind = "Rejection"
  )
  updated_state <- clado_event(
    timeval = 0.5,
    island_tbl = create_test_island_tbl(island_scenario = 1),
    max_spec_id = 1
  )
  expect_true(is.list(updated_state))
  expected_island_tbl <- create_test_island_tbl(island_scenario = 19)
  expect_equal(updated_state$island_tbl, expected_island_tbl)
  expect_equal(updated_state$max_spec_id, 3)
})

test_that("clado_event produces correct output on cladogenesis species", {
  set.seed(
    1,
    kind = "Mersenne-Twister",
    normal.kind = "Inversion",
    sample.kind = "Rejection"
  )
  updated_state <- clado_event(
    timeval = 0.9,
    island_tbl = create_test_island_tbl(island_scenario = 19),
    max_spec_id = 3
  )
  expect_true(is.list(updated_state))
  expected_island_tbl <- create_test_island_tbl(island_scenario = 45)
  expect_equal(updated_state$island_tbl, expected_island_tbl)
  expect_equal(updated_state$max_spec_id, 5)
})

test_that("clado_event fails with incorrect input", {
  island_tbl <- create_test_island_tbl(island_scenario = 1)
  expect_error(clado_event(
    timeval = "nonsense",
    island_tbl = island_tbl,
    max_spec_id = 1
  ))
  expect_error(clado_event(
    timeval = 0.5,
    island_tbl = "nonsense",
    max_spec_id = 1
  ))
  expect_error(clado_event(
    timeval = 0.5,
    island_tbl = island_tbl,
    max_spec_id = "nonsense"
  ))
})