tests/testthat/test-LSA.R

#   MIT License
#
#   Copyright(c) 2017-2018 Ilia Kosenkov [ilia.kosenkov.at.gm@gmail.com]
#
#   Permission is hereby granted, free of charge, to any person obtaining a copy
#   of this software and associated documentation files(the "Software"), to deal
#   in the Software without restriction, including without limitation the rights
#   to use, copy, modify, merge, publish, distribute, sublicense, and / or sell
#   copies of the Software, and to permit persons to whom the Software is
#   furnished to do so, subject to the following conditions:
#
#   The above copyright notice and this permission
#   notice shall be included in all
#   copies or substantial portions of the Software.
#
#   THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
#   IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
#   FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT.
#   IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM,
#   DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION OF CONTRACT,
#   TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION WITH
#   THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.

context("Lomb-Scaragle analysis tests")

test_that("`LSAPeriodogramEx` works as `LSAPeriodogram`", {

    t <- runif(200) %>% sort
    x <- rnorm(200)

    w <- seq(0, 5 * 2 * pi, by = 0.1)

    p1 <- lsa_periodogram(w, t, x - mean(x))

    p2 <- lsa_periodogram_ex(tibble(T = t, X = x), T, X, w)
    expect_equal(p1, p2$PSD, tolerance = 1e-16)

})
Ilia-Kosenkov/RLibs documentation built on Jan. 26, 2020, 2:21 p.m.