tests/testthat/test-new_seeds.R

In bain: Bayes Factors for Informative Hypotheses

# Due to a change in CRAN policy, we are no longer allowed to generate random
# seeds in the way we previously did in bain.
# We have implemented a new way to generate seeds. This test checks whether those
# provide the same results as previously, and whether they provide consistent
# results upon repeat evaluation.
estimate <- c(1,1,1)
names(estimate)<-c("a", "b", "c")
sampN<- 100
cov <- matrix(c(1,0,0,0,1,0,0,0,1),nrow=3,ncol=3)

set.seed(2)
y<-bain(estimate,"a>-.96 & a < 2.96 & b>-.96 & b < 2.96 & c>-.96 & c < 2.96",n=sampN,Sigma=cov,group_parameters=0,joint_parameters = 3)

# This test intentionally fails, because we allow bain 0.2.9's random seeds to
# be not compatible with < 0.2.9:
test_that("New seed same result as old bain", {
  expect_false(abs(y$fit$Fit_in[1] - 0.866728935817668) < .0000000000000001)
  })

set.seed(2)
y2 <- bain(estimate,"a>-.96 & a < 2.96 & b>-.96 & b < 2.96 & c>-.96 & c < 2.96",n=sampN,Sigma=cov,group_parameters=0,joint_parameters = 3)

test_that("New seed same result as old bain", {
  expect_equal(y$fit$Fit_in[1], y2$fit$Fit_in[1])
})