test_jfst_r.R
In neodistr: Neo-Normal Distribution

library(neodistr)
library(testthat)

# ----------------------------A. Test R Functions------------------------------------------------
# ----------------------------Test Inapproriate parameters-------------------------------------------

test_that("1. Wrong parameter values in R PDF and PMF functions", {
  
  expect_error(expect_true(is.nan(djfst(1, 1, -1, 1, 1))))
  expect_error(expect_true(is.nan(djfst(1, 1, 1, -1, 1))))
  expect_error(expect_true(is.nan(djfst(1, 1, 1, 1, -1))))
  
})

test_that("2. Wrong parameter values in R CDF functions", {
  
  expect_error(expect_true(is.nan(pjfst(1, 1, -1, 1, 1))))
  expect_error(expect_true(is.nan(pjfst(1, 1, 1, -1, 1))))
  expect_error(expect_true(is.nan(pjfst(1, 1, 1, 1, -1))))
  
})


test_that("3. Wrong parameter values in R quantile functions", {
  
  expect_error(expect_true(is.nan(qjfst(0.5, 1, -1, 1, 1))))
  expect_error(expect_true(is.nan(qjfst(0.5, 1, 1, -1, 1))))
  expect_error(expect_true(is.nan(qjfst(0.5, 1, 1, 1, -1))))
  
})


test_that("4. Wrong parameter values in R RNG functions", {
  
  expect_error(expect_true(is.nan(rjfst(1, 1, -1, 1, 1))))
  expect_error(expect_true(is.nan(rjfst(1, 1, 1, -1, 1))))
  expect_error(expect_true(is.nan(rjfst(1, 1, 1, 1, -1))))
  
  expect_error(expect_true(is.nan(rng_jfst(1, 1, -1, 1, 1))))
  expect_error(expect_true(is.nan(rng_jfst(1, 1, 1, -1, 1))))
  expect_error(expect_true(is.nan(rng_jfst(1, 1, 1, 1, -1))))
  
})

# ------------------------------------Infinity Test-----------------------------------------------
test_that("5. Testing PDFs  against infinite values", {
  
  expect_true(!is.nan(djfst(Inf, 0, 1, 1, 1)) && is.finite(djfst(Inf, 0, 1, 1, 1)))
  
})

test_that("6. Testing CDFs against infinite values", {
  
  expect_true(!is.nan(pjfst(Inf, 0, 1, 1, 1)) && is.finite(pjfst(Inf, 0, 1, 1, 1)))
  
})

test_that("7. Testing PDFs  against negatively infinite values", {
  
  expect_true(!is.nan(djfst(-Inf, 0, 1, 1, 1)) && is.finite(djfst(-Inf, 0, 1, 1, 1)))
  
})

test_that("8. Testing CDFs against negatively infinite values", {
  
  expect_true(!is.nan(pjfst(-Inf, 0, 1, 1, 1)) && is.finite(pjfst(-Inf, 0, 1, 1, 1)))
  
})

# -------------------------------Test Log Probability--------------------------------------------------

test_that("9. Check if log-probabilities are logs of probabilities (PDF's)", {
  
  x <- c(-Inf, -15, -10, -5, -1, -0.5, 0, 0.5, 1, 5, 10, 15, Inf)
  
  expect_equal(suppressWarnings(djfst(x, 0, 1, 1, 1, log = TRUE)),
               log(suppressWarnings(djfst(x, 0, 1, 1, 1))))
  
})

test_that("10. Check if log-probabilities are logs of probabilities (CDF's)", {
  
  x <- c(-Inf, -15, -10, -5, -1, 0, 1, 5, 10, 15, Inf)
  
  expect_equal(suppressWarnings(pjfst(x, 0, 1, 1, 1, log.p = TRUE)),
               log(suppressWarnings(pjfst(x, 0, 1, 1, 1))))
  
})

# -------------------------------------Test NAs-------------------------------------

test_that("11. Missing values in PDF and PMF functions", {
  
  expect_true(is.na(djfst(NA, 0, 1, 1, 1)))
  expect_error(expect_true(is.na(djfst(1, NA, 1, 1, 1))))
  expect_error(expect_true(is.na(djfst(1, 0, NA, 1, 1))))
  expect_error(expect_true(is.na(djfst(1, 0, 1, NA, 1))))
  expect_error(expect_true(is.na(djfst(1, 0, 1, 1, NA))))
  
})

test_that("12. Missing values  in CDF functions", {
  
  expect_true(is.na(pjfst(NA, 0, 1, 1, 1)))
  expect_error(expect_true(is.na(pjfst(1, NA, 1, 1, 1))))
  expect_error(expect_true(is.na(pjfst(1, 0, NA, 1, 1))))
  expect_error(expect_true(is.na(pjfst(1, 0, 1, NA, 1))))
  expect_error(expect_true(is.na(pjfst(1, 0, 1, 1, NA))))
  
})

test_that("13.  Missing values in inverse CDF functions", {
  
  expect_error(expect_true(is.na(qjfst(NA, 0, 1, 1, 1))))
  expect_error(expect_true(is.na(qjfst(1, NA, 1, 1, 1))))
  expect_error(expect_true(is.na(qjfst(1, 0, NA, 1, 1))))
  expect_error(expect_true(is.na(qjfst(1, 0, 1, NA, 1))))
  expect_error(expect_true(is.na(qjfst(1, 0, 1, 1, NA))))
  
})

test_that("14. Missing values in RNG functions", {
  
  expect_error(expect_true(is.na(rjfst(1, NA, 1, 1, 1))))
  expect_error(expect_true(is.na(rjfst(1, 0, NA, 1, 1))))
  expect_error(expect_true(is.na(rjfst(1, 0, 1, NA, 1))))
  expect_error(expect_true(is.na(rjfst(1, 0, 1, 1, NA))))
  
})

# --------------------- Test Probabillity -------------------------------------------

test_that("15. All probabilities/densities >= 0", {
  
  y <- c( -100, -10, -5, -1, -0.5, 0, 0.5, 1, 5, 10, 100)
  
  expect_true(suppressWarnings(all(djfst(y, 0, 1, 1, 1) >= 0)))
  
})

test_that("16. All cumulative probabilities >= 0 and <= 1", {
  
  y <- c( -100, -10, -5, -1, -0.5, 0, 0.5, 1, 5, 10, 100)
  
  expect_true(all(pjfst(y, 0, 1, 1, 1) >= 0 & pjfst(y, 0, 1, 1, 1) <= 1))
  
})


# -------------------------------------Test Quantile Function ------------------------------

test_that("17. Zeros in quantile functions", {
  
  expect_error(expect_true(!is.nan(qjfst(0, 0, 1, 1, 1))))
  
})

test_that("18. Ones in quantile functions", {
  
  expect_error(expect_true(!is.nan(qjfst(1, 0, 1, 1, 1))))
  
})

test_that("19. Checking p = F(F^-1(p))", {
  
  pp <- seq(0.00001, 1, by = 0.001)
  
  expect_equal(pp, pjfst(qjfst(pp, 0, 1, 1, 1), 0, 1, 1, 1))
  
})

# ----------------------------- Test RNG Converge -------------------------------------

probCoverage <- function(stub, ..., n = 5000L) {
  rfoo <- eval(as.name(paste0("r", stub)))
  pfoo <- eval(as.name(paste0("p", stub)))
  diff(pfoo(range(rfoo(n, ...)), ...))
}



test_that("20. Coverage of RNG's", {
  
  skip_on_cran()
  
  expect_gte(probCoverage("jfst", 0, 1, 1, 1), 0.99)
  
})


# -------------------------- PDF and CDF Test -----------------------------------------
library(neodistr)
library (testthat)

test_that("21. integrate PDF from -inf to inf == 1", {
  
  expect_equal(integrate(djfst,lower=-Inf, upper = Inf,mu=0,sigma=1,alpha=1,beta=1 )$value,1)
  
})

test_that("22. integrate PDF from - Inf until x equal to cdf x", {
  
  expect_equal((integrate(djfst, lower= -Inf, upper=4, mu=0, sigma=1,alpha=1, beta=1)$value), (pjfst(4, mu=0, sigma=1, alpha=1, beta=1)))
  
})