Diversity index

This page is a method. The page identifier is Op_en5455
Moderator:Jouni (see all)
This page is a stub. You may improve it into a full page, and then a rating bar will appear here.
Upload data Show results

Question

How to calculate diversity indices?

Answer

Upload your data to Opasnet Base. Use the function diversity to calculate the most common indices.

Actual function diversity

+ Show code- Hide code 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("0D or S", "log(1D) or H'", "λ 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) }

Example to use function

+ Show code- Hide code library(xtable) 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 = TRUE) print(xtable(out), type = 'html')

Rationale

Diversity indices are thoroughly described in Wikipedia.

Related files

<mfanonymousfilelist></mfanonymousfilelist>