tests/testthat/test_cline_summary.R

In tjthurman/BAHZ: Bayesian Analysis of Hybrid Zones

context("cline summary")

load("ref_geno_stanfit.Rda")
load("ref_pheno_stanfit.Rda")

# Test of argument checks
test_that("cline summary argument checks work", {
  expect_error(cline_summary(stanfit = "a"), "stanfit")
  expect_error(cline_summary(ref_geno_stanfit, prob = "x"), "numeric")
  expect_error(cline_summary(ref_geno_stanfit, prob = c(0.6, 0.9)), "length")
  expect_error(cline_summary(ref_geno_stanfit, prob = -.5), "between")
  expect_error(cline_summary(ref_geno_stanfit, prob = 1.1), "between")
  expect_error(cline_summary(ref_geno_stanfit, method = "xxx"), "method")
  expect_error(cline_summary(ref_geno_stanfit, show.all = "xxx"), "show.all")
})


# Test of functionality
test_that("cline summary provides correct summary", {
  expect_equal_to_reference(cline_summary(ref_geno_stanfit,
                                          prob = 0.5), file = "ref_summary1.Rda")
  expect_equal_to_reference(cline_summary(ref_geno_stanfit, method = "ET",
                                          prob = 0.875), file = "ref_summary2.Rda")
  expect_equal_to_reference(cline_summary(ref_geno_stanfit, show.all = T,
                                          prob = 0.65), file = "ref_summary3.Rda")
  expect_equal_to_reference(cline_summary(ref_pheno_stanfit, show.all = F,
                                          prob = 0.65), file = "ref_summary4.Rda")
})