tests/testthat/test-gbdcd.R

In leandromineti/gbdcd: Implementation of the Bayesian Detection of Clusters and Discontinuities

context("Gaussian BDCD")

context("RcppArmadillo functions")

test_that("Frequency matrix", {
  t <- c(1:64)
  dt <- rcpp_freqmatrix(t)
  expect_is(dt, "matrix")
  expect_equal(nrow(dt), 64)
  expect_equal(ncol(dt), 64)
})

data("aneeldata", package = "gbdcd")
data("aneelshape", package = "gbdcd")

test_that("Partitioning", {
  dt <- aneeldata[[1]]
  part <- rcpp_partition(dt, c(3, 40)) 
  expect_equal(length(part), 57)
})

context("Gbdcd")

test_that("gbdcd function", {
  res <- gbdcd(y = aneelshape$z_Precipitation, neigh = aneeldata$connections,
                      n_iterations = 100, burn_in = 50, c = 0.35, 0, sigma_0 = sqrt(2))
  expect_is(res, "list")
})