tests/testthat/test-as-mcmcarray.R

In poissonconsulting/mcmcr: Manipulate MCMC Samples

test_that("as.mcmcarray.mcmc with multiple parameters", {
  expect_error(as.mcmcarray(as.mcmc(subset(mcmcr::mcmcr_example, chains = 1L))),
    "^`x` must only have 1 parameter[.]$",
    class = "chk_error"
  )
})

test_that("as.mcmcarray.mcmc", {
  mcmcarray <- subset(mcmcr::mcmcr_example, chains = 1L)$beta
  expect_identical(as.mcmcarray(coda::as.mcmc(mcmcarray)), mcmcarray)

  mcmcarray <- subset(mcmcr::mcmcr_example, chains = 1L)$sigma

  x <- coda::as.mcmc(mcmcarray)
  expect_identical(as.mcmcarray(x), mcmcarray)
  colnames(x) <- "1"
  expect_warning(expect_identical(as.mcmcarray(x), structure(numeric(0), .Dim = c(1L, 400L, 0L), class = "mcmcarray")))
})

test_that("as.mcmcarray.mcmc with missing values", {
  mcmcarray <- subset(mcmcr::mcmcr_example, chains = 1L)$beta
  mcmc <- coda::as.mcmc(mcmcarray)
  mcmc <- mcmc[, c("parameter[1,1]", "parameter[2,1]", "parameter[2,2]")]
  mcmc2 <- coda::as.mcmc(as.mcmcarray(mcmc))
  mcmc2 <- mcmc2[, c("parameter[1,1]", "parameter[2,1]", "parameter[2,2]")]
  expect_identical(mcmc2, mcmc)
})