tests/test-all.R

In willpearse/phest: Calculate PHEnological ESTimates

library(testthat)
library(phest)

test_that("Weibull estimates are sensible",{
    # These are the Roberts & Solow data (2003)
    expect_equal(
        weib.limit(c(1662, 1638, 1631, 1628, 1628, 1611, 1607, 1602, 1601, 1598), upper=TRUE),
        setNames(c(1690, 1669, 1799), c("estimate","lower-ci","upper-ci")),
        tolerance=0.0005
    )
    expect_equal(
        weib.limit(c(1662, 1638, 1631, 1628, 1628, 1611, 1607, 1602, 1601, 1598)),
        setNames(c(1595, 1575, 1597), c("estimate","lower-ci","upper-ci")),
        tolerance=0.0005
    )    
})

test_that("Weibull CIs are checked",{
    model <- weib.limit(c(1, 2, 3, rep(4,50)))
    expect_true(all(is.na(model[-1])))
})