R/extract_simdat.R

In bayesnec: A Bayesian No-Effect- Concentration (NEC) Algorithm

#'extract_simdat
#'
#' Extracts a range of statistics from a \code{\link{prebayesnecfit}} object.
#'
#' @param modfit An object of class \code{\link{prebayesnecfit}}, as returned
#' by \code{\link{fit_bayesnec}}.
#'
#' @return A \code{\link[base]{list}} of model simulation statistics including
#' iter, thin, warmup, chains, init and nsamples.
#'
#' @noRd
extract_simdat <- function(modfit) {
  x <- modfit$fit$fit@sim
  sample_prior <- attr(modfit$fit$prior, "sample_prior")
  list(iter = x$iter, thin = x$thin, warmup = x$warmup, chains = x$chains,
       init = modfit$fit$fit@inits, sample_prior = sample_prior,
       n_samples = ceiling((x$iter - x$warmup) /
                             x$thin * x$chains))
}