Good afternoon everyone,
I've been analyzing microbiome data in qiime2, but recently, due to the necessity of performing some adjustments to data analyses, I have migrated to R and started to use qiime2R and phyloseq.
Today, due to major differences in the p-values of 2 permanovas performed on weighted_unifrac data ( both non-significant) one preformed in qiime2, and another preformed in phyloseq (with adonis2), I've realized that the two distance matrixes are different...
I've preformed rarefaction in phyloseq prior to distance calculation (with the same sampling depth used in qiime2)
the major difference that I have observed is 1log in the distance values (other variations might be due to rarefaction)... Does anyone have a clue of what might be happening? are the formulas to calculate weighted unifrac distance different between the 2 pipelines???
code in qiime2:
qiime diversity core-metrics-phylogenetic
(the matrix is one of the outputs of this comand)