tests/testthat/test-spec-utils.R

In EarthSystemDiagnostics/proxysnr: Separate Signal and Noise in Climate Proxy Records

test_that("spectral estimation works", {

  m <- "missing data"
  expect_error(SpecMTM(stats::ts(c(rnorm(41), rep(NA, 9), rnorm(50)))), m)

  # only simple testing here on correct output structure

  data <- rnorm(1000)

  actual <- SpecMTM(stats::ts(data))
  spec.mtm <- multitaper::spec.mtm(stats::ts(data), plot = FALSE)

  expect_true(utils::hasName(actual, "dof"))

  actual$dof <- NULL

  expect_equal(names(actual), names(spec.mtm))

})

test_that("log-smoothing works", {

  # only simple testing here on correct output structure

  data <- rnorm(1000)

  spec <- SpecMTM(stats::ts(data))
  actual <- LogSmooth(spec)

  expect_type(actual, "list")
  expect_s3_class(actual, "spec")

  expect_true(all(utils::hasName(actual, c("freq", "spec", "dof"))))

  expect_equal(spec$freq, actual$freq)

})

test_that("averaging spectra works", {

  m <- "MeanSpectrum: Spectra are of different lengths."
  s1 <- SpecMTM(stats::ts(rnorm(100)))
  s2 <- SpecMTM(stats::ts(rnorm(1000)))
  s3 <- SpecMTM(stats::ts(rnorm(1000)))

  expect_error(MeanSpectrum(list(s1, s2, s3)), m, fixed = TRUE)

  s1 <- SpecMTM(stats::ts(rnorm(1000)))
  s2 <- SpecMTM(stats::ts(rnorm(1000)))
  s3 <- SpecMTM(stats::ts(rnorm(1000)))

  m <- MeanSpectrum(list(s1, s2, s3))

  expect_type(m, "list")
  expect_s3_class(m, "spec")

  expect_true(all(utils::hasName(m, c("freq", "spec", "dof"))))

  expect_equal(m$freq, s1$freq)
  expect_equal(m$spec, (s1$spec + s2$spec + s3$spec) / 3)
  expect_equal(m$dof, s1$dof + s2$dof + s3$dof)

})

test_that("running SimPLS works", {

  N <- 100

  expect_length(SimPLS(N), N)
  expect_length(SimPLS(N, alpha = 1), N)
  expect_length(SimPLS(N, beta = 1, alpha = 1), N)
  expect_length(SimPLS(N, beta = -1, alpha = 1), N)

})