# Load Paracou data (number of trees per species in two 1-ha plot of a tropical forest)
data(Paracou618)
# Ps is the vector of probabilities
Ps <- Paracou618.MC$Ps
# Calculate Simpson entropy
Simpson(Ps)
# Ns is the vector of abundances of the metacommunity
Ns <- Paracou618.MC$Ns
# Calculate an unbiased estimator of Simpson's index of diversity
bcSimpson(Ns)
Run the code above in your browser using DataLab