x <- sample(letters[1:4], 100, replace = TRUE, p = c(1, 2, 3, 4))
types <- c("richness", "effective.number", "shannon",
"simpson", "inv.simpson", "gini.simpson", "unb.gini",
"eveness.simpson", "eveness.pielou"
)
sapply(types, function(tp) diversity(x, type = tp))
# hill numbers with increasing order
order <- 0:5
hill.num <- sapply(order, function(q) diversity(x, type = "hill", q = q))
hill.num
plot(order, hill.num, type = "b")
# a matrix of frequencies
spp.freq <- cbind(
sample(letters[1:4], 100, replace = TRUE, p = c(1, 1, 1, 4)),
sample(letters[1:4], 100, replace = TRUE, p = c(4, 1, 1, 1)),
sample(letters[1:4], 100, replace = TRUE, p = c(1, 1, 1, 1))
)
diversity(spp.freq, type = "eff")
Run the code above in your browser using DataLab