R/weighted.entropy.R

Defines functions weighted.entropy

Documented in weighted.entropy

weighted.entropy <-
function(x,w=NULL,alpha=1){
  if(is.null(w)) w<-rep(1/length(x),length(x))
  x <- as.numeric(x)
  n <- length(x)
  fx<- w/sum(w)
  if (is.null(alpha)) 
    alpha <- 1
  if (alpha == 0) { # yields mean logarithmic 
    entropy <- -sum(log(x/weighted.mean(x,w))*fx) 
  }
  else if (alpha == 1) 
    entropy <- sum(x/weighted.mean(x,w)*log(x/weighted.mean(x,w))*fx)
  else {
    entropy <- 1/(alpha * (alpha - 1))*t((x/weighted.mean(x,w))^alpha-1)%*%fx
  }
  return(entropy)
}

Try the acid package in your browser

Any scripts or data that you put into this service are public.

acid documentation built on May 29, 2017, 1:03 p.m.