tests/testthat/test-moezipf.R

In zipfextR: Zipf Extended Distributions

test_that("[moezipf - preconditions] Checking preconditions", {
  expect_error(dmoezipf(1, "2.5", -1.5),
               "Incorrect alpha parameter. This parameter should be greater than one.")
  expect_error(pmoezipf(1, 0.5, 1.5),
               "Incorrect alpha parameter. This parameter should be greater than one.")
  expect_error(pmoezipf(1, 2.5, "-1.5"),
               "Incorrect beta parameter. You should provide a numeric value.")
  expect_error(dmoezipf("1", 1.5, 1.5),
               label = "The x value is not included into the support of the distribution.")
  expect_error(pmoezipf(0, 2.5, 1.5),
               "The x value is not included into the support of the distribution.")
  expect_error(dmoezipf(1.6, 2.5, 1.5),
               "The x value is not included into the support of the distribution.")
  expect_error(qmoezipf(1.3, 2.5, 1.6),
               "Wrong values for the p parameter.")
  expect_error(rmoezipf(1.6, 2.5, 1.5),
               "The x value is not included into the support of the distribution.")
})

test_that("[moezipf - pmf] The summation of all probabilities must be 1.", {
  expect_equal(1,sum(dmoezipf(1:1000, alpha = 2.5, beta = 0.10)), tolerance = 1*10^(-3))
  expect_equal(1,sum(dmoezipf(1:1000, alpha = 2.5, beta = 1.3)), tolerance = 1*10^(-3))
  expect_equal(1,sum(dmoezipf(1:1000, alpha = 2.5, beta = 3)), tolerance = 1*10^(-3))
})

test_that("[moezipf - pmf] The probability lays on the interval (0, 1).", {
  pvalue <- dmoezipf(x = 9, alpha = 2, beta = 0.10)
  expect_lt(pvalue, 1)
  expect_gt(pvalue, 0)

  pvalue <- dmoezipf(x = 9, alpha = 1.25, beta = 5.3)
  expect_lt(pvalue, 1)
  expect_gt(pvalue, 0)
})

test_that("[moezipf - cdf] The summation of all probabilities must be 1.", {
  expect_equal(1,pmoezipf(1000, alpha = 2.5, beta = 0.10), tolerance = 1*10^(-3))
  expect_equal(1,pmoezipf(1000, alpha = 2.5, beta = 1.3), tolerance = 1*10^(-3))
  expect_equal(1,pmoezipf(1000, alpha = 2.5, beta = 3), tolerance = 1*10^(-3))
})

test_that("[moezipf- cdf] The cumulative probabilities have to be in the interval (0, 1).", {
  pvalue <- pmoezipf(q = 9, alpha = 2, beta = 0.10)
  expect_lt(pvalue, 1)
  expect_gt(pvalue, 0)

  pvalue <- pmoezipf(q = 9, alpha = 1.25, beta = 5.3)
  expect_lt(pvalue, 1)
  expect_gt(pvalue, 0)
})