R/chains2sample.R

In RBesT: R Bayesian Evidence Synthesis Tools

#' Scrambles the order of a mcmc array object for usage as a mcmc
#' sample. It is advisable to set order once per mcmc run, otherwise
#' correlations in the mcmc sample will be lost.
#' @keywords internal
chains2sample <- function(chains, order, drop = TRUE) {
  d <- dim(chains)
  N <- prod(d[2:3])
  if (missing(order)) {
    order <- sample.int(N)
  }
  matrix(chains, nrow = d[1], ncol = N)[, order, drop = drop]
}