R/simpson_index.R

In compbiomed/animalcules: Interactive microbiome analysis toolkit

#' Get alpha diversity using simpson
#'
#' @param x A list of counts
#' @return A single value
#'
#' @examples
#' simpson_index(seq_len(10))
#'
#' @export

simpson_index <- function(x) {
    # Relative abundances
    p <- x / sum(x)
    
    # Simpson index
    lambda <- sum(p^2)
    
    lambda
}