tests/testthat/test-tapered_estimator.R

In CovEsts: Nonparametric Estimators for Covariance Functions

# taper_single
test_that("taper_single(, , \"tukey\") works", {
  expect_equal(taper_single(0.5, 1, "tukey"), 1)
})

test_that("taper_single() fails for nonnumeric x", {
  expect_error(taper_single('a', 1, "tukey"))
  expect_error(taper_single(NA, 1, "tukey"))
})

test_that("taper_single() fails x of length not equal to 1", {
  expect_error(taper_single(c(), 1, "tukey"))
  expect_error(taper_single(c(1, 2), 1, "tukey"))
})

test_that("taper_single() fails for negative x", {
  expect_error(taper_single(-0.01, 1, "tukey"))
})

test_that("taper_single() fails for x great than 1", {
  expect_error(taper_single(1.01, 1, "tukey"))
})

test_that("taper_single() fails for nonnumeric rho", {
  expect_error(taper_single(0.5, 'a', "tukey"))
  expect_error(taper_single(0.5, 1i, "tukey"))
})

test_that("taper_single() fails for negative rho", {
  expect_error(taper_single(0.5, -0.01, "tukey"))
})

test_that("taper_single() fails for rho that is not a single number", {
  expect_error(taper_single(0.5, c(1, 2), "tukey"))
})

test_that("taper_single() fails for window_name that is not one of 'tukey', 'triangular', 'power_sine', 'blackman_window', 'hann_poisson', 'welch'", {
  expect_error(taper_single(0.5, 1, "window"))
})

test_that("taper_single() fails for nonboolean custom_window", {
  expect_error(taper_single(0.5, 1, "window", custom_window = 1))
  expect_error(taper_single(0.5, 1, "window", custom_window = 'TRUE'))
})


# H2n
test_that("H2n(, , \"tukey\") works", {
  expect_equal(H2n(3, 1, "tukey"), 1.5)
})

test_that("H2n(, , \"tukey\") fails for nonnumeric n", {
  expect_error(H2n('a', 1, "tukey"))
  expect_error(H2n(1i, 1, "tukey"))
})

test_that("H2n(, , \"tukey\") fails for n < 1", {
  expect_error(H2n(0, 1, "tukey"))
})

test_that("H2n(, , \"tukey\") fails for noninteger n", {
  expect_error(H2n(2.95, 1, "tukey"))
})

test_that("H2n() fails for nonnumeric rho", {
  expect_error(H2n(3, 'a', "tukey"))
  expect_error(H2n(3, 1i, "tukey"))
})

test_that("H2n() fails for negative rho", {
  expect_error(H2n(3, -0.01, "tukey"))
})

test_that("H2n() fails for rho that is not a single number", {
  expect_error(H2n(0.5, c(1, 2), "tukey"))
})

test_that("H2n() fails for window_name that is not one of 'tukey', 'triangular', 'power_sine', 'blackman_window', 'hann_poisson', 'welch'", {
  expect_error(H2n(3, 1, "window"))
})

test_that("H2n() fails for nonboolean custom_window", {
  expect_error(H2n(3, 1, "window", custom_window = 1))
  expect_error(H2n(3, 1, "window", custom_window = 'TRUE'))
})

# taper
test_that("taper(, , \"tukey\") works", {
  expect_equal(taper(c(0.25, 0.5, 0.75), 1, "tukey"), c(0.5, 1, 0.5))
})

test_that("taper(, , \"tukey\") fails for nonnumeric x", {
  expect_error(taper(c(0.25, 'a', 0.75), 1, "tukey"))
  expect_error(taper(c(0.25, 1i, 0.75), 1, "tukey"))
})

test_that("taper(, , \"tukey\") fails for empty x", {
  expect_error(taper(c(), 1, "tukey"))
})

test_that("taper(, , \"tukey\") fails for any NAs in x", {
  expect_error(taper(c(0.25, NA, 0.75), 1, "tukey"))
})

test_that("taper() fails for nonnumeric rho", {
  expect_error(taper(3, 'a', "tukey"))
  expect_error(taper(3, 1i, "tukey"))
})

test_that("taper() fails for negative rho", {
  expect_error(taper(3, -0.01, "tukey"))
})

test_that("taper() fails for rho that is not a single number", {
  expect_error(taper(0.5, c(1, 2), "tukey"))
})

test_that("taper() fails for nonboolean custom_window", {
  expect_error(taper(3, 1, "window", custom_window = 1))
  expect_error(taper(3, 1, "window", custom_window = 'TRUE'))
})

# tapered_single
test_that("tapered_single(, , \"tukey\") works", {
  expect_equal(tapered_single(test_X, test_meanX, 1, test_h2n, test_taperVals_t, test_taperVals_h), 0)
})

test_that("tapered_single() fails for nonnumeric X", {
  expect_error(tapered_single(c(1, 'a', 3), test_meanX, 1, test_h2n, test_taperVals_t, test_taperVals_h))
  expect_error(tapered_single(c(1, 1i, 3), test_meanX, 1, test_h2n, test_taperVals_t, test_taperVals_h))
})

test_that("tapered_single() fails for empty X", {
  expect_error(tapered_single(c(), test_meanX, 1, test_h2n, test_taperVals_t, test_taperVals_h))
})

test_that("tapered_single() fails for at least one NA in X", {
  expect_error(tapered_single(c(1, NA, 3), test_meanX, 1, test_h2n, test_taperVals_t, test_taperVals_h))
})

test_that("tapered_single() fails for nonnumeric meanX", {
  expect_error(tapered_single(test_X, 'a', 1, test_h2n, test_taperVals_t, test_taperVals_h))
  expect_error(tapered_single(test_X, 1i, 1, test_h2n, test_taperVals_t, test_taperVals_h))
})

test_that("tapered_single() fails for meanX not of length 1", {
  expect_error(tapered_single(test_X, c(), 1, test_h2n, test_taperVals_t, test_taperVals_h))
  expect_error(tapered_single(test_X, c(1, 2), 1, test_h2n, test_taperVals_t, test_taperVals_h))
})

test_that("tapered_single() fails for nonnumeric h", {
  expect_error(tapered_single(test_X, test_meanX, 'a', test_h2n, test_taperVals_t, test_taperVals_h))
  expect_error(tapered_single(test_X, test_meanX, 1i, test_h2n, test_taperVals_t, test_taperVals_h))
})

test_that("tapered_single() fails for h not of length 1", {
  expect_error(tapered_single(test_X, test_meanX, c(), test_h2n, test_taperVals_t, test_taperVals_h))
  expect_error(tapered_single(test_X, test_meanX, c(1, 2), test_h2n, test_taperVals_t, test_taperVals_h))
})

test_that("tapered_single() fails for noninteger h", {
  expect_error(tapered_single(test_X, test_meanX, 0.1, test_h2n, test_taperVals_t, test_taperVals_h))
})

test_that("tapered_single() fails for nonnumeric h2n", {
  expect_error(tapered_single(test_X, test_meanX, 1, 'a', test_taperVals_t, test_taperVals_h))
  expect_error(tapered_single(test_X, test_meanX, 1, 1i, test_taperVals_t, test_taperVals_h))
})

test_that("tapered_single() fails for h2n not of length 1", {
  expect_error(tapered_single(test_X, test_meanX, 1, c(), test_taperVals_t, test_taperVals_h))
  expect_error(tapered_single(test_X, test_meanX, 1, c(1, 2), test_taperVals_t, test_taperVals_h))
})

test_that("tapered_single() fails for nonnumeric taperVals_t", {
  expect_error(tapered_single(test_X, test_meanX, 1, test_h2n, c(0.1, 'a', 0.3), test_taperVals_h))
  expect_error(tapered_single(test_X, test_meanX, 1, test_h2n, c(0.1, 1i, 0.3), test_taperVals_h))
})

test_that("tapered_single() fails for empty taperVals_t", {
  expect_error(tapered_single(test_X, test_meanX, 1, test_h2n, c(), test_taperVals_h))
})

test_that("tapered_single() fails for any taperVals_t between [-1, 1]", {
  expect_error(tapered_single(test_X, test_meanX, 1, test_h2n, c(-1.1, 0.2, 0.3), test_taperVals_h))
  expect_error(tapered_single(test_X, test_meanX, 1, test_h2n, c(0.1, 0.2, 1.1), test_taperVals_h))
})

test_that("tapered_single() fails for nonnumeric taperVals_t", {
  expect_error(tapered_single(test_X, test_meanX, 1, test_h2n, test_taperVals_t, c(0.1, 'a', 0.3)))
  expect_error(tapered_single(test_X, test_meanX, 1, test_h2n, test_taperVals_t, c(0.1, 1i, 0.3)))
})

test_that("tapered_single() fails for empty taperVals_t", {
  expect_error(tapered_single(test_X, test_meanX, 1, test_h2n, test_taperVals_h, c()))
})

test_that("tapered_single() fails for any taperVals_t between [-1, 1]", {
  expect_error(tapered_single(test_X, test_meanX, 1, test_h2n, test_taperVals_t, c(-1.1, 0.2, 0.3)))
  expect_error(tapered_single(test_X, test_meanX, 1, test_h2n, test_taperVals_t, c(0.1, 0.2, 1.1)))
})

# tapered_est
test_that("tapered_est(, , , \"tukey\") works", {
  expect_equal(tapered_est(test_X, 1, "tukey", maxLag = 2), c(1/12, 0, -1/24))
})

test_that("tapered_est() fails for nonnumeric X", {
  expect_error(tapered_est(c(1, 'a', 3), 1, "tukey", maxLag = 2))
  expect_error(tapered_est(c(1, 1i, 3), 1, "tukey", maxLag = 2))
})

test_that("tapered_est() fails for empty X", {
  expect_error(tapered_est(c(), 1, "tukey", maxLag = 2))
})

test_that("tapered_est() fails for at least one NA in X", {
  expect_error(tapered_est(c(1, NA, 3), 1, "tukey", maxLag = 2))
})

test_that("tapered_est() fails for at nonnumeric maxLag", {
  expect_error(tapered_est(test_X, 1, "tukey", maxLag = 'a'))
  expect_error(tapered_est(test_X, 1, "tukey", maxLag = 1i))
})

test_that("tapered_est() fails for maxLag whose length is not 1", {
  expect_error(tapered_est(test_X, 1, "tukey", maxLag = c()))
  expect_error(tapered_est(test_X, 1, "tukey", maxLag = c(1, 2)))
})

test_that("tapered_est() fails for maxLag less than 1", {
  expect_error(tapered_est(test_X, 1, "tukey", , maxLag = 0))
})

test_that("tapered_est() fails for maxLag greater than or equal to length(X)", {
  expect_error(tapered_est(test_X, 1, "tukey", maxLag = 3))
})

test_that("tapered_est() fails for noninteger maxLag", {
  expect_error(tapered_est(test_X, 1, "tukey", maxLag = 1.2))
})

test_that("tapered_est() fails for nonnumeric rho", {
  expect_error(tapered_est(test_X, 'a', "tukey", maxLag = 2))
  expect_error(tapered_est(test_X, 1i, "tukey", maxLag = 2))
})

test_that("tapered_est() fails for rho not of length 1", {
  expect_error(tapered_est(test_X, c(), "tukey", maxLag = 2))
  expect_error(tapered_est(test_X, c(1, 2), "tukey", maxLag = 2))
})

test_that("tapered_est() fails for nonboolean custom_window", {
  expect_error(tapered_est(test_X, 1, "tukey", custom_window = 1, , maxLag = 2))
  expect_error(tapered_est(test_X, 1, "tukey", custom_window = 'TRUE', maxLag = 2))
})