R/AcceptanceRate.R

Defines functions AcceptanceRate

Documented in AcceptanceRate

###########################################################################
# AcceptanceRate                                                          #
#                                                                         #
# The purpose of the AcceptanceRate function is to calculate the          #
# acceptance rate of each chain from its samples.                         #
###########################################################################

AcceptanceRate <- function(x)
     {
     if(missing(x)) stop("x is a required argument.")
     if(!is.matrix(x)) x <- as.matrix(x)
     out <- colMeans(x[-nrow(x),] != x[-1,])
     names(out) <- colnames(x)
     return(out)
     }

x <- matrix(rnorm(10*10),10,10)
colnames(x) <- paste("V", 1:10, sep="")
x[2,] <- x[1,]
AcceptanceRate(x)

#End

Try the LaplacesDemon package in your browser

Any scripts or data that you put into this service are public.

LaplacesDemon documentation built on July 1, 2018, 9:02 a.m.