Test: jack

Ivan Jacob Agaloos Pesigan

2020-07-16

library(testthat)
library(jeksterslabRboot)
context("Test jack.")

Parameters

data <- c(
  11,
  22,
  33
)
thetahat <- mean(data)
jack_samples <- jack(data = data)
thetahatstarjack <- sapply(
  X = jack_samples,
  FUN = mean
)

Results

results <- jack_hat(
  thetahatstarjack = thetahatstarjack,
  thetahat = thetahat
)

testthat

test_that("mean", {
  expect_equivalent(
    results[["hat"]][["mean"]],
    22
  )
})

test_that("bias", {
  expect_equivalent(
    results[["hat"]][["bias"]],
    0
  )
})

test_that("se", {
  expect_equivalent(
    results[["hat"]][["se"]],
    0
  )
})

test_that("thetahatjack", {
  expect_equivalent(
    results[["hat"]][["thetahatjack"]],
    22
  )
})

Example

data <- c(
  17.23,
  13.93,
  15.78,
  14.91,
  18.21,
  14.28,
  18.83,
  13.45,
  18.71,
  18.81,
  11.29,
  13.39,
  11.57,
  10.94,
  15.52,
  15.25
)
ps <- c(
  1.605,
  1.151,
  0.998,
  0.942,
  2.416,
  1.043,
  3.122,
  1.362,
  2.972,
  3.097,
  3.308,
  1.393,
  2.951,
  3.806,
  0.958,
  0.937
)
mean_ps <- 2.00389
ci_ll <- 1.45
ci_ul <- 2.56
example <- jack_hat(
  thetahatstarjack = sapply(X = jack(data), FUN = function(x) log(var(x))),
  thetahat = log(var(data))
)

test_that("mean_ps", {
  expect_equivalent(
    round(
      x = mean_ps,
      digits = 3
    ),
    round(
      x = example[["hat"]][["mean_ps"]],
      digits = 3
    )
  )
})

test_that("ps", {
  expect_equivalent(
    round(
      x = ps,
      digits = 3
    ),
    round(
      x = example[["ps"]],
      digits = 3
    )
  )
})

test_that("ll", {
  expect_equivalent(
    ci_ll,
    round(
      x = example[["ci"]][["ci_2.5"]],
      digits = 2
    )
  )
})

test_that("ul", {
  expect_equivalent(
    ci_ul,
    round(
      x = example[["ci"]][["ci_97.5"]],
      digits = 2
    )
  )
})