entropy <- function(str) {
  vec   <- strsplit(str, "")[[1]]
  N     <- length(vec)
  p_xi  <- table(vec) / N

  -sum(p_xi * log(p_xi, 2))
}