tests/testthat/test-thr_L2.R

In ccid: Cross-Covariance Isolate Detect: a New Change-Point Method for Estimating Dynamic Functional Connectivity

context("The change-points in the mean structure for a multivariate
        time series of multiplicative structure via taking the square
        root of the average of squares of contrast function values")

# In this file we test the thr_L2 function

test_that("The correct change-points are given when the function is used correctly", {
  set.seed(1)
  num.nodes <- 40 # number of nodes
  etaA.1    <- 0.95
  etaA.2    <- 0.05
  pcor1     <- GeneNet::ggm.simulate.pcor(num.nodes, etaA = etaA.1)
  pcor2     <- GeneNet::ggm.simulate.pcor(num.nodes, etaA = etaA.2)

  n <- 100
  data1 <- GeneNet::ggm.simulate.data(n, pcor1)
  data2 <- GeneNet::ggm.simulate.data(n, pcor2)

  X <- rbind(data1, data2, data1, data2) ## change-points at 100, 200, 300
  sgn <- sign(stats::cor(X))

  expect_equal(thr_L2(t(hdbinseg::gen.input(t(X), -1, TRUE, diag = FALSE, sgn))), c(100, 200, 300))
})

test_that("Error and warning messages are given correctly", {
  set.seed(100)
  test_data <- matrix(data = rchisq(1000, df = 5), nrow = 500, ncol = 2)

  # A string vector is given as the input for X.
  expect_error(thr_L2(X = "I am not a numeric matrix"))

  # The threshold constant is set equal to zero.
  expect_error(thr_L2(X = test_data, thr_const = 0))

  # The value for lambda is negative.
  expect_error(thr_L2(X = test_data, points = -2))

  # The value for lambda is a positive real number instead of a positive integer.
  expect_warning(thr_L2(X = test_data, points = 3.4))
})