Statistical analysis and graph processing were conducted on R (Version 4.2.3) using RStudio software (Version 2022.07.0). Several packages were used in the analysis: phyloseq (Version 1.42.0) to facilitate the import, storage, handling and analysis of the microbiome data (McMurdie and Holmes, 2013 (link)), as well as to determine alpha diversity. For the latter, both Chao1 and Shannon indexes were calculated. Vegan (Version 2.6-4) package was used for the permutational multivariate analysis of variance (PERMANOVA) using the Bray Curtis metric distance. The mixOmics package was used for sparse partial least squares discriminant analysis (Rohart et al., 2017 (link)) and the microbiomeMarker (Version 1.4) package for the Linear Discriminant Analysis (LDA) effect size (SE). The ggplot2 (Version 3.4.4) package was used for graphical representation of the data.
Free full text: Click here