Difference between revisions of "Diversity index"

true <- function(pi, q){
qD <- (sum(pi^q))^(1 / (1 - q))
return(qD)
}
diversity <- function(data, individual = TRUE){
if(individual == FALSE) {pi <- data} else {
     pi <- tapply(data, data, length)
     pi <- pi/sum(pi)
}

out <- data.frame(Name = c("Richness", "Shannon index", "Simpson index", "Inverse Simpson index", 
     "Gini-Simpson index", "Berger-Parker index"), 
     Symbol = c("S or 0D", "H' or log(1D)", "λ or 1/(2D)", "1/λ or 2D", "1-λ or 1-1/(2D)", "∞D"),
     Value = c(true(pi,0), log(true(pi,1)), 1/true(pi,2), true(pi,2), 1-1/true(pi,2), max(pi)))
return(out)
}

library(xtable)
if(is.null(data)) {data <- c(1,4,6,9,3,4,5,6,5,4,3,3,5,5,7,5,5,4,3,4,5,6,8,9,5,4,5,4,3,3,4,9,6,6,4,5,3,2,1,1,2,3,4,3,2,3,4,5,6,7,7)}
out <- diversity(data, individual = individual)
print(xtable(out), type = 'html')