entropy_ChaoIn <-
function(data, B = 200, conf = 0.95) {
t <- data[1]
y <- data[-1]
est <- ChaoEntropyEstFunIn(y, t)
se <- BootstrapFunIn(y, t, B, ChaoEntropyEstFunIn)
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("Chao entropy")
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.