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])))
})
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.