Nothing
context("maxLikelihood.ESF")
test_that("maxLikelihood.ESF: use", {
set.seed(42)
J <- 10000
theta <- 100
m <- 0.1
I <- m * (J - 1) / (1 - m)
v <- generate.ESF(theta, I, J)
LL <- maxLikelihood.ESF(init_vals = c(100, 0.1),
abund = c(v))
m_est <- LL$par[2]
expect_equal(m,
m_est,
tolerance = 0.05)
expect_equal(theta,
LL$par[1],
tolerance = 5)
})
test_that("maxLikelihood.ESF: abuse", {
set.seed(42)
J <- 100
theta <- 100
m <- 0.1
I <- m * (J - 1) / (1 - m)
v <- generate.ESF(theta, I, J)
expect_error(
LL <- maxLikelihood.ESF(init_vals = c(-100, 0.1),
abund = v , verbose = FALSE),
"initial theta can not be below one"
)
expect_error(
LL <- maxLikelihood.ESF(init_vals = c(100, -0.1),
abund = v, verbose = FALSE),
"initial m can not be below zero"
)
expect_error(
LL <- maxLikelihood.ESF(init_vals = c(100, 100),
abund = v, verbose = FALSE),
"initial m can not be above 1"
)
expect_error(
LL <- maxLikelihood.ESF(init_vals = c(100, 0.1),
abund = c(1), verbose = FALSE),
"Need more than 1 species in the dataset"
)
})
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.