alpha diversity analysis in qiime2 or in R

Hi everyone.

I was just wondering if it is better to get diversity metrics with Qiime or with the phyloseq package in R. Are there any issues relating to using either or is it the same thing but different methods?

Hello Nengi,

Welcome to the forums! :qiime2:

Both packages should work well for calculating alpha diversity values! Under the hood, Qiime2 uses scikit-bio diversity to calculate different alpha values, and Phyloseq uses Vegan.

Note that there are some differences in the implementation, like the base of the log used to calculate Shannon diversity. R / Vegan / Phyloseq uses natural log (base e) while Python / Scikit Bio / Qiime uses log2 (base 2).

:warning: As I write this, Phyloseq should not be used for calculating phylogenetic BETA diversity values, as there is an outstanding bug with how it calculates UniFrac distances.

:+1: But the alpha values are fine either way.

3 Likes

Right. Thank you very much. I've used qiime2 (core metrics) for both the alpha and beta diversity calculations but I was just curious to know if there was a difference. Thanks again :blush: