R/entropy_Grassberger.R

In JohnsonHsieh/ChaoEntropy: Statistical package 'ChaoEntropy'

entropy_Grassberger <-
function(data, B = 200, conf = 0.95) {
  if (is.matrix(data) == TRUE) {
    if (ncol(data) != 1 & nrow(data) != 1)
      stop("Error: The data format is wrong.")
    if (ncol(data) == 1) {
      data <- data[, 1]
    } else {
      data <- data[1, ]
    }
  }
  
  est <- GrassEstFun(data)
  se <- BootstrapFun(data, B, GrassEstFun)
  z <- qnorm(1 - (1 - conf)/2)
  CI <- c(max(est - z * se, 0), est + z * se)
  out <- matrix(c(est, se, CI), nrow = 1)
  rownames(out) <- c("Grassberger (2003)")
  colnames(out) <- c("Estimator", "Bootstrap s.e.",
                     paste(conf*100, "% Lower"), paste(conf*100, "% Upper"))
  return(out)
}