ObservedEstFunIn <- function(y, t) { y <- y[y > 0] U <- sum(y) pi <- (y / t) / (U / t) est <- -sum(pi * log(pi)) return(est) }
Add the following code to your website.
REMOVE THIS Copy to clipboard
For more information on customizing the embed code, read Embedding Snippets.