tests/testthat/test_spinar_boot.R

In spINAR: (Semi)Parametric Estimation and Bootstrapping of INAR Models

# spinar_boot <- function(x, p, B, setting, type = NA, distr = NA, M = 100, level = 0.05)
test_that("input", {
   tmp <- sample(1:10, 5, replace = TRUE)
   ######################## x ########################
   expect_error(spinar_boot(c(1), 1, 10, setting = "sp"), "Assertion on 'x' failed")
   expect_error(spinar_boot(c(1, 2, -3, 3, 3, 3), 1, 10, setting = "sp"))
   expect_error(spinar_boot(c(1, 2, 1.5, 3, 3, 3), 1, 10, setting = "sp"), "Assertion on 'x' failed: Must be of type 'integerish'")
   ######################## p ########################
   expect_error(spinar_boot(tmp, "1", 100, setting = "sp"), "Assertion on 'p' failed: Must be of type 'integerish', not 'character'")
   expect_error(spinar_boot(tmp, 3, 100, setting = "sp"), "Assertion on 'p' failed")
   expect_error(spinar_boot(tmp, 0, 100, setting = "sp"), "Assertion on 'p' failed")
   expect_error(spinar_boot(tmp, 1.5, 100, setting = "sp"), "Assertion on 'p' failed: Must be of type 'integerish'")
   expect_error(spinar_boot(tmp, "2", 100, setting = "sp"), "Assertion on 'p' failed: Must be of type 'integerish', not 'character'")
   ######################## B ########################
   expect_error(spinar_boot(tmp, 1, -10, setting = "sp"), "Assertion on 'B' failed")
   expect_error(spinar_boot(tmp, 1, 0, setting = "sp"), "Assertion on 'B' failed")
   expect_error(spinar_boot(tmp, 1, 10.5, setting = "sp"), "Assertion on 'B' failed: Must be of type 'integerish'")
   ######################## setting ########################
   expect_error(spinar_boot(tmp, 1, 10, setting = 0))
   ######################## type ########################
   expect_error(spinar_boot(tmp, 1, 10, setting = "p", type = 2))
   ######################## distr ########################
   expect_error(spinar_boot(tmp, 1, 10, setting = "p", type = "mom", distr = 21))
   ######################## M ########################
   expect_error(spinar_boot(tmp, 1, 10, setting = "p", type = "mom", distr = 'poi', M="text"), "Assertion on 'M' failed: Must be of type 'integerish', not 'character'.")
   expect_error(spinar_boot(tmp, 1, 10, setting = "p", type = "mom", distr = 'poi', M=-1), "Assertion on 'M' failed")
   ######################## level ########################
   expect_error(spinar_boot(tmp, 1, 10, setting = "p", type = "mom", distr = 'poi', M=2, level = 2), "Assertion on 'level' failed")
})

set.seed(123)
m <- 10
tmp = spinar_boot(x=sample(1:10, m, replace=TRUE), p=2, B=2, setting = "sp")

test_that("output", {
   ######################## size ########################
   expect_true(expect_true(dim(tmp$x_star)[1]==m))
   expect_true(expect_true(dim(tmp$x_star)[2]==2))
   expect_true(sum(tmp$parameters_star[1,3:(length(tmp$parameters_star)/2)])==1)
   expect_true(length((tmp$parameters_star)) == length(tmp$bs_ci_hall))
   expect_true(length((tmp$bs_ci_percentile)) == length(tmp$bs_ci_hall))
   ######################## type ########################
   expect_type(tmp, "list")
   ######################## values ########################
   expect_true(tmp$parameters_star[1] < 1)
   expect_true(tmp$parameters_star[1] > 0)
   tmp0 <- spinar_boot(x=sample(1:8, 10, replace=TRUE), p=1, B=2, setting = "p", type = "mom", distr="poi", progress = FALSE)
   expect_true(length(tmp0$bs_ci_hall)== 4)
   tmp1 <- spinar_boot(x=sample(1:8, 10, replace=TRUE), p=1, B=2, setting = "p", type = "mom", distr="geo")
   expect_true(length(tmp1$bs_ci_percentile)==4)
   tmp2 <- spinar_boot(x=sample(1:8, 10, replace=TRUE), p=1, B=2, setting = "p", type = "mom", distr="nb", level = 1)
   print(length(tmp2$parameters)==6)
})

test_that("print.spinar_boot", {
   expect_output(print(tmp), "spinar_boot object")
   expect_output(print(tmp), "B=2, n=10")
   expect_output(print(tmp), "Hall's Bootstrap Percentile Confidence Intervals")
   expect_output(print(tmp), "Bootstrap Percentile Confidence Intervals")
})