entropy_Zhang <-
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 <- ZhangHzstarEstFun(data)
se <- BootstrapFun(data, B, ZhangHzstarEstFun)
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("Zhang (2012) Hz*")
colnames(out) <- c("Estimator", "Bootstrap s.e.",
paste(conf*100, "% Lower"), paste(conf*100, "% Upper"))
return(out)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.