tests/testthat/test-peak.R

In pmcharrison/har18: Implementation of Harrison & Pearce's (2018) Harmonicity Model

context("test-peak")

library(hrep)
library(magrittr)

test_that("example", {
  c(0, 0, 3, 4, 1) %>%
    .milne_pc_spectrum() %>%
    peak %>%
    expect_equal(4)
})

test_that("invariance to doubling the number of bins", {
  for (i in 1:10) {
    spec1 <- rnorm(10) %>% .milne_pc_spectrum()
    spec2 <- rep(as.numeric(spec1), each = 2) %>% .milne_pc_spectrum()
    expect_equal(
      peak(spec1),
      peak(spec2)
    )
  }
})

test_that("more pointy distributions get higher peaks", {
  x <- seq(from = -100, to = 100)
  pointy <- dnorm(x, sd = 1)
  less_pointy <- dnorm(x, sd = 4)
  expect_gt(
    peak(.milne_pc_spectrum(pointy)),
    peak(.milne_pc_spectrum(less_pointy))
  )
})