The R Project SVN R-packages

Rev

Rev 1938 | Blame | Compare with Previous | Last modification | View Log | Download | RSS feed 

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")
}