tests/testthat/test-P_hnorm.R

In PeteyCoco/ppdeconv: Poisson Process Deconvolution Estimation

test_that("the dimension of P is correct", {
  # Test using the length.out argument
  l_breaks = seq(0,100, length.out = 201)
  r_breaks = seq(0,100, length.out = 101)

  P <- P_hnorm(l_breaks = l_breaks,
                 r_breaks = r_breaks,
                 sd = 10)

  expect_equal(dim(P), c( 100, 200))

  # Test using the by argument
  l_breaks <- seq(0, 100, by = 0.5)
  r_breaks <- seq(0, 100, by = 1)

  P <- P_hnorm(l_breaks = l_breaks,
                 r_breaks = r_breaks,
                 sd = 10)

  expect_equal(dim(P), c( 100, 200))
})

test_that("the entries of P are positive", {
  expect(all(P_hnorm(l_breaks = seq(0,10, length.out = 11),
                       r_breaks = seq(0,20, length.out = 21),
                       sd = 1) >= 0),
         failure_message = "Some entries of P are negative")
})