R/r.cha.R

In affluenceIndex: Affluence (Richness) Indices

r.cha <- function(x, weight, k, beta){
  n <- length(x)
  if (is.null(weight)) weight <- rep(1, n)
  
  rhow <- k*weighted.median(x, weight)
  ind <- ifelse(x > rhow, 1, 0)
  r.cha <- sum(((1-(rhow/(x[x>0]))^beta)*ind[x>0])*weight[x>0])/sum(weight)
  return(r.cha)
}