tests/test_seeds.R

cat("Testing for non duplicated results in seeded simulations")

library(bkmrhat)
set.seed(111)
dat <- bkmr::SimData(n = 50, M = 4)
y <- dat$y
Z <- dat$Z
X <- dat$X

future::plan(strategy = future::sequential)
fitkm.list <- kmbayes_parallel(nchains=2, y = y, Z = Z, X = X, iter = 10,
                               verbose = FALSE, varsel = TRUE)

rr = as.mcmc.list(fitkm.list)
tt = all.equal(as.data.frame(rr[[1]]), as.data.frame(rr[[2]]))
stopifnot(mode(tt)=="character")


fitkm.list[[1]] = fitkm.list[[2]]
badpred = SamplePred_parallel(fitkm.list, method="exact", Znew = dat$Z[1,], Xnew = dat$X[1,])
ff = tapply(badpred$znew1, badpred$chain, mean)
tt = all.equal(ff[1], ff[2])
stopifnot(mode(tt)=="character")

Try the bkmrhat package in your browser

Any scripts or data that you put into this service are public.

bkmrhat documentation built on March 29, 2022, 9:08 a.m.