R/spcdisc.R

In labdsv: Ordination and Multivariate Analysis for Ecology

Defines functions spcdisc

Documented in spcdisc

spcdisc <- function(x,sort=FALSE)
{
   shannon <- function(y) {
       y <- as.numeric(y)
       frac <- y/sum(y)
       comp <- sum((-1 * frac * log(frac))[frac>0])
       comp <- 1 - (comp/log(length(y)))
       comp
   }
   tmp <- apply(x,1,shannon)
   if (sort) tmp <- tmp[rev(order(tmp))]
   tmp
}

Try the labdsv package in your browser