tests/testthat/test_analyse_SAR.NCF.R

In Luminescence: Comprehensive Luminescence Dating Data Analysis

## load data
file <- system.file("extdata/NCF.binx", package = "Luminescence")
ncf <- read_BIN2R(file, fastForward = TRUE, verbose = FALSE)

test_that("input validation", {
  testthat::skip_on_cran()

  expect_error(analyse_SAR.NCF(iris),
               "'object' should be of class 'RLum.Analysis' or a 'list' of such")
  expect_error(analyse_SAR.NCF(set_RLum("RLum.Analysis")),
               "'object' cannot be an empty RLum.Analysis")

  data(ExampleData.BINfileData, envir = environment())
  object <- Risoe.BINfileData2RLum.Analysis(CWOSL.SAR.Data, pos = 1:2)
  expect_error(analyse_SAR.NCF(object, signal_integral = 1:2,
                               background_integral = 100:250),
               "No additional dose point found, check that the NCF-SAR protocol")

  expect_message(expect_message(expect_null(
      analyse_SAR.NCF(set_RLum("RLum.Analysis",
                               records = list(set_RLum("RLum.Data.Curve",
                                                       recordType = "false",
                                                       info = list(IRR_TIME = 10)))),
                      signal_integral = 1:2, background_integral = 100:200)),
      "Error: No record of type 'OSL', 'IRSL', 'POSL' detected, NULL returned"),
      "Error: CW-OSL analysis skipped: check your sequence, NULL returned")
})

test_that("check functionality", {
  testthat::skip_on_cran()

  SW({
  res <- analyse_SAR.NCF(ncf, signal_integral = 1:2,
                         background_integral = 100:250,
                         dose_rate_source = 0.1)
  })
  expect_type(res, "list")
  expect_length(res, 1)
  expect_s4_class(res[[1]], "RLum.Results")

  expect_silent(analyse_SAR.NCF(ncf, signal_integral = 1:2,
                                background_integral = 100:250,
                                plot = FALSE, verbose = FALSE))
})

test_that("snapshot tests", {
  testthat::skip_on_cran()

  snapshot.tolerance <- 1.5e-5

  SW({
  expect_snapshot_RLum(
      analyse_SAR.NCF(
          ncf,
          signal_integral = 1:2,
          background_integral = 100:250,
          dose_rate_source = 0.1),
      tolerance = snapshot.tolerance)
  })
})

test_that("graphical snapshot tests", {
  testthat::skip_on_cran()
  testthat::skip_if_not_installed("vdiffr")

  set.seed(1)

  SW({
  vdiffr::expect_doppelganger("default",
                              analyse_SAR.NCF(ncf, signal_integral = 1:2,
                                              background_integral = 100:250,
                                              dose_rate_source = 0.1))
  })
})