tests/testthat/test_pcomp.R

In CompGLM: Conway-Maxwell-Poisson GLM and Distribution Functions

context("pcomp")

test_that("pcomp matches ppois with nu = 1", {
    expect_equal(ppois(-5:5, 2.5, TRUE, FALSE), pcomp(-5:5, 2.5, 1, 100L, TRUE, FALSE))
    expect_equal(ppois(1:2, 2:5, TRUE, FALSE), pcomp(1:2, 2:5, 1, 100L, TRUE, FALSE))
    expect_equal(ppois(2:5, 1:2, TRUE, FALSE), pcomp(2:5, 1:2, 1, 100L, TRUE, FALSE))

    expect_equal(ppois(-5:5, 2.5, FALSE, TRUE), pcomp(-5:5, 2.5, 1, 100L, FALSE, TRUE))
    expect_equal(ppois(1:2, 2:5, FALSE, FALSE), pcomp(1:2, 2:5, 1, 100L, FALSE, FALSE))
    expect_equal(ppois(1:2, 2:5, TRUE, TRUE), pcomp(1:2, 2:5, 1, 100L, TRUE, TRUE))
})

test_that("pcomp throws", {
    expect_error(pcomp(1, -1, 2))
    expect_error(pcomp(1, 1, -2))
})