tests/testthat/test-metric_midpointCE.R

In LightLogR: Process Data from Wearable Light Loggers and Optical Radiation Dosimeters

test_that("Calculation works", {
  MEDI = c(rep(0, 6), rep(250, 12), rep(0, 6))
  Datetime = lubridate::as_datetime(lubridate::dhours(0:23), tz = "UTC")
  expect_equal(midpointCE(MEDI, Datetime), lubridate::as_datetime(lubridate::dhours(11), tz = "UTC"))
})

test_that("Missing values", {
  MEDI = c(rep(0, 6), rep(250, 10), NA, NA, rep(0, 6))
  Datetime = lubridate::as_datetime(lubridate::dhours(0:23), tz = "UTC")
  expect_equal(midpointCE(MEDI, Datetime), lubridate::as_datetime(NA))
  expect_equal(midpointCE(MEDI, Datetime, na.rm = TRUE), lubridate::as_datetime(lubridate::dhours(10), tz = "UTC"))
})

test_that("Works with different date time units", {
  MEDI = c(rep(0, 6), rep(250, 12), rep(0, 6))
  Datetime = lubridate::as_datetime(lubridate::dhours(0:23), tz = "UTC")
  HMS = hms::as_hms(Datetime)
  Duration = lubridate::as.duration(HMS)
  Difftime = Datetime-Datetime[1]
  expect_equal(midpointCE(MEDI, HMS), hms::hms(0,0,11))
  expect_equal(midpointCE(MEDI, Duration), lubridate::dhours(11))
  expect_equal(midpointCE(MEDI, Difftime), lubridate::as.difftime(lubridate::dhours(11)))
})

test_that("Return data frame works", {
  MEDI = c(rep(0, 6), rep(250, 12), rep(0, 6))
  Datetime = lubridate::as_datetime(lubridate::dhours(0:23), tz = "UTC")
  expect_equal(
    midpointCE(MEDI, Datetime, as.df = TRUE),
    tibble::tibble(
      "midpointCE" = lubridate::as_datetime(lubridate::dhours(11), tz = "UTC"),
    )
  )
})