tests/testthat/test-bwLscvg.R

In circularKDE: Recent Methods for Kernel Density Estimation of Circular Data

test_that("bwLscvg returns a numeric value for valid input", {
  set.seed(60)
  x <- rvonmises(50, circular(pi / 2), 1, control.circular = list(units = "radians"))
  result <- bwLscvg(x, g = 4)$minimum
  expect_equal(result, 6.6326176)
  expect_type(result, "double")
  expect_length(result, 1)
})

test_that("bwLscvg returns a numeric value for valid input", {
  set.seed(123)
  x <- rvonmises(50, circular(pi / 2), 1)
  result <- bwLscvg(x)$minimum
  expect_equal(result, 5.09811023)
  expect_type(result, "double")
  expect_length(result, 1)
})

test_that("bwLscvg throws error on empty input", {
  expect_error(bwLscvg(numeric(0)), "`x` must be a non-empty object.")
})

test_that("bwLscvg throws error if input is not numeric", {
  expect_error(bwLscvg(c("a", "b")), "must be a numeric vector")
})

test_that("bwLscvg throws error if x contains only NAs", {
  x <- c(NA, NA)
  expect_error(bwLscvg(x), "`x` contains all missing values.")
})

test_that("bwLscvg removes NA values and returns result", {
  x <- circular(c(0, pi / 2, NA, pi))
  result <- bwLscvg(x, g = 4)$minimum
  expect_type(result, "double")
  expect_cli_warning(bwLscvg(x, g = 4),
                     1,
                     "! `x` contains missing values, which will be removed.")
})

test_that("bwLscvg handles non-numeric g", {
  x <- circular(seq(0, 2 * pi, length.out = 5))
  expect_cli_warning(
    result <- bwLscvg(x, g = "wrong")$minimum,
    1,
    "! Argument `g` must be positive numeric number and not equal to 2. Value 4 for coefficient was used."
  )
  expect_type(result, "double")
})

test_that("bwLscvg handles non-numeric lower", {
  x <- circular(seq(0, 2 * pi, length.out = 5))
  expect_cli_warning(
    result <- bwLscvg(x, lower = "zero", g = 4)$minimum,
    1,
    "! Argument `lower` must be numeric. Default value 0 for lower boundary was used."
  )
  expect_type(result, "double")
})

test_that("bwLscvg handles non-numeric upper", {
  x <- circular(seq(0, 2 * pi, length.out = 5))
  expect_cli_warning(
    result <- bwLscvg(x, upper = "sixty", g = 4)$minimum,
    1,
    "! Argument `upper` must be numeric. Default value 60 for upper boundary was used."
  )
  expect_type(result, "double")
})

test_that("bwLscvg warns and resets invalid boundary values", {
  x <- circular(seq(0, 2 * pi, length.out = 5))
  expect_cli_warning(
    result <- bwLscvg(x, lower = -5, upper = 5, g = 4)$minimum,
    1,
    "! The boundaries must be positive numbers and 'lower' must be smaller than 'upper'. Default boundaries lower=0, upper=60 were used."
  )
  expect_type(result, "double")

  expect_cli_warning(
    result <- bwLscvg(x, lower = 10, upper = 5, g = 4)$minimum,
    1,
    "! The boundaries must be positive numbers and 'lower' must be smaller than 'upper'. Default boundaries lower=0, upper=60 were used."
  )
  expect_type(result, "double")
})

test_that("bwLscvg warns when minimum is at edge of the range", {
  x <- circular(rep(0, 10))
  expect_cli_warning(bwLscvg(x, g = 4)$minimum,
                     1,
                     "! Minimum/maximum occurred at one end of the range.")
})