tests/testthat/test-immig_list.R

In Yangshen0325/specmutual: Evolution of Mutualism Based on Speciation, Immigration and Extinction

# without mutualism, immigration rates for plant species should be equaling to
# that of animal species (nrow(M0) = ncol(M0)).
M0 <- matrix(1, 2500, ncol = 50, nrow = 50)
M0[lower.tri(M0)] <- 0
M0 <- M0[nrow(M0):1, ]
testit::assert(colSums(M0) == rowSums(M0))
gam_pars <- c(1.0, 1.0)

test_that("Immigration rates for plant and animal species are identical
          without mutualism", {
  K_pars <- c(500, 500, Inf, Inf)
  immig_list <- get_immig_rate(
    M0 = M0,
    Mt = M0,
    K_pars = K_pars,
    status_a = matrix(1, ncol = 1, nrow = nrow(M0)),
    status_p = matrix(1, ncol = 1, nrow = ncol(M0)),
    gam_pars = gam_pars
  )

  expect_equal(
    immig_list$immig_p,
    immig_list$immig_a
  )
})

# Keep plants and animals have the same number of muutalistic partners
test_that("Immigration rates for plant and animal species are identical
          with mutualism", {
  K_pars <- c(500, 500, 1.0, 1.0)
  immig_list <- get_immig_rate(
    M0 = M0,
    Mt = M0,
    K_pars = K_pars,
    status_a = matrix(1, ncol = 1, nrow = nrow(M0)),
    status_p = matrix(1, ncol = 1, nrow = ncol(M0)),
    gam_pars = gam_pars
  )

  expect_equal(
    immig_list$immig_p,
    immig_list$immig_a
  )
})