head.mcmc

head.mcmc.R

head.mcmc <- function(x, n = 6L, ...) {
  mcpar <- attr(x, "mcpar")
  start <- mcpar[1]
  max_end <- mcpar[2]
  thin <- mcpar[3]
  
  niter <- nrow(x)
  end <- start + niter - 1
  attr(x, "mcpar")[2:3] <- c(end, 1)
  out <- window(x, end = min(start + n - 1, end))
  end <- min(start + n * thin - 1, max_end)
  attr(out, "mcpar")[2:3] <- c(end, thin)
  out
}

x$fit@sim$iter * x$fit@sim$chain - (nrow(out) - 1) * x$fit@sim$thin,