inst/classic-bugs/R/check.R

In nimble: MCMC, Particle Filtering, and Programmable Hierarchical Modeling

library(coda, quietly=TRUE)
x <- read.openbugs(quiet=TRUE)
checkstats <- summary(x)$statistics
if (is.matrix(benchstats)) {
  rnb <- order(rownames(benchstats))
  rnc <- order(rownames(checkstats))
  z <- (benchstats[rnb,1] - checkstats[rnc,1])
  z <- ifelse(z==0, 0, z/benchstats[rnb,2])
} else {
  z <- (benchstats[1] - checkstats[1])
  z <- ifelse(z==0, 0, z/benchstats[2])
}
print(z)
q(status=(any(abs(z) > 0.15)))