R/rejectionRate.R

In coda: Output Analysis and Diagnostics for MCMC

Defines functions rejectionRate

Documented in rejectionRate

rejectionRate.mcmc <- function (x) {
  x <- as.matrix(x)
  apply(x[-nrow(x),,drop=FALSE] == x[-1,, drop=FALSE],2,mean)
}

rejectionRate.mcmc.list <- function (x) {
  apply(sapply(x,rejectionRate.mcmc),1,mean)
}

rejectionRate <- function(x) {
  UseMethod("rejectionRate")
}

Try the coda package in your browser