tests/testthat/test-set_posterior_by_index.R

In richelbilderbeek/wiritttes: WIRITTE simulation

context("set_posterior_by_index")

test_that("set_posterior_by_index: use", {
  file <- read_file(find_path("toy_example_1.RDa"))
  posterior <- get_posterior_by_index(file = file, i = 1)
  testthat::expect_true(tracerer::is_posterior(posterior))
  posterior_new <- tracerer::parse_beast_posterior(
    trees_filename = find_path(filename = "beast2_example_output.trees"),
    log_filename   = find_path(filename = "beast2_example_output.log")
  )
  names(posterior_new) <- c("trees", "estimates")
  testthat::expect_true(tracerer::is_posterior(posterior_new))
  testthat::expect_false(are_identical_posteriors(posterior, posterior_new))
  file <- set_posterior_by_index(
    file = file,
    i = 1,
    posterior = posterior_new
  )
  posterior_new_again <- get_posterior_by_index(file = file, i = 1)
  names(posterior_new_again) <- c("trees", "estimates")
  testthat::expect_true(
    are_identical_posteriors(
      posterior_new,
      posterior_new_again
    )
  )
})


test_that("set_posterior_by_index: abuse", {

  expect_error(
    set_posterior_by_index(
      file = read_file(find_path("toy_example_1.RDa")),
      i = 0,
      posterior = ape::rcoal(10)
    ),
    "index must be at least 1"
  )

  expect_error(
    set_posterior_by_index(
      file = read_file(find_path("toy_example_1.RDa")),
      i = 3,
      posterior = ape::rcoal(10)
    ),
    "index must be less than number of posteriors"
  )

})