tests/testthat/test-logLikelihood.ESF.R

In GUILDS: Implementation of Sampling Formulas for the Unified Neutral Model of Biodiversity and Biogeography, with or without Guild Structure

context("logLikelihood.ESF")

test_that("logLikelihood.ESF: use", {
  set.seed(42)
  J <- 1000
  theta <- 100
  m <- 0.1
  I <- m * (J - 1) / (1 - m)

  abund <- generate.ESF(theta, I, J)
  LL1 <- logLikelihood.ESF(theta, m, abund)
  LL2 <- logLikelihood.ESF(theta * 2, m * 2, abund)

  testthat::expect_gt(LL1, LL2)
})

test_that("logLikelihood.ESF: abuse", {
  set.seed(42)
  J <- 1000
  theta <- 100
  m <- 0.1
  I <- m * (J - 1) / (1 - m)

  abund <- generate.ESF(theta, I, J)
  LL1 <- logLikelihood.ESF(-10, m, abund)
  expect_equal(is.infinite(LL1[[1]]), TRUE)
  LL1 <- logLikelihood.ESF(theta, 2, abund)
  expect_equal(is.infinite(LL1[[1]]), TRUE)
})