tests/testthat/test-sig_noise.R

In wincowgerDEV/OpenSpecy: Analyze, Process, Identify, and Share Raman and (FT)IR Spectra

data("raman_hdpe")

test_that("sig_noise() handles input errors corretly", {
  sig_noise(1:1000) |> expect_error()
  restrict_range(raman_hdpe, 300, 320) |> sig_noise() |>
    expect_warning()
})

test_that("sig_noise() can restrict the range correctly", {
  big_peak <- sig_noise(raman_hdpe, metric = "sig", sig_max = 3000,
                        sig_min = 2500)
  low_region <- sig_noise(raman_hdpe, metric = "sig", sig_max = 1000,
                          sig_min = 500)
  expect_true(big_peak > low_region)
  restric_peak <- restrict_range(raman_hdpe, min = 2500, max = 3000,
                                 make_rel = F) |>
    sig_noise(metric = "sig")
  expect_identical(big_peak, restric_peak)
})

test_that("sig_noise() returns correct values", {
  sig_noise(raman_hdpe, metric = "sig") |> round(2) |> unname() |>
    expect_equal(101.17)

  sig_noise(raman_hdpe, metric = "noise") |> round(2) |> unname() |>
    expect_equal(61.01)

  sig_noise(raman_hdpe, metric = "sig_times_noise") |> round(2) |> unname() |>
    expect_equal(6172.1)

  sig_noise(raman_hdpe, metric = "sig_over_noise") |> round(2) |> unname() |>
    expect_equal(1.66)

  sig_noise(raman_hdpe, metric = "tot_sig") |> round(2) |> unname() |>
    expect_equal(97527)
})

test_that("entropy results in accurate info", {
    sig_noise(raman_hdpe, 
              metric = "entropy",
              breaks = 10) |> 
        round(2) |>
        unname() |>
        expect_equal(1.17)
    
    noise <- raman_hdpe
    set.seed(10)
    noise$spectra[[1]] <- runif(n = length(noise$spectra[[1]]))
    
    sig_noise(noise, metric = "entropy",
              breaks = 10) |> 
        round(2) |>
        unname() |>
        expect_equal(3.32)
    
    expect_true(sig_noise(raman_hdpe, 
                          metric = "entropy",
                          breaks = 10) |> 
                    unname() < 
                  sig_noise(noise, metric = "entropy",
                                       breaks = 10) |> unname())
})