Wednesday, August 27, 2025
\[ S = \text{Number of observed taxa} \]
\[ E_{\mathrm{Simpson}} = \frac{1}{S \sum_{i=1}^{S} p_i^2} \]
where:
\(p_i\) = relative abundance of taxon \(i\)
\(S\) = number of taxa
Bounded between 0 (low evenness) and 1 (perfectly even)
\[ H' = - \sum_{i=1}^{S} p_i \log(p_i) \]
where:
\(p_i\) = relative abundance of taxon \(i\)
\(S\) = total number of taxa
Combines richness and evenness in one metric
Guidelines for alpha diversity: Cassol, Ibañez, and Bustamante (2025)
library(mia)
tse <- addAlpha(tse)
library(miaViz)
plotBoxplot(tse, col.var = "faith_diversity", x = "patient_status") +
labs(x = "Diagnosis", y = "Faith's diversity")
From OMA online book, Chapter 14: Alpha diversity