Hi. i have a question while i’m using the lefse command.
I analyzed my samples from 2 time points(0 time and 1 month).
I got the p-value 0.013 using amova.
As i wanted to know which differences were between 2 time points, i used metastats then lefse.
I knew which OTUs were significantly different between 2 time points form metastats.
But during using lefse, it was said “You did not provide a class”.
I wonder the OTUs significantly different are in which time point and i want to know how i can use the class option and whether i should add the class in which file (design file or metadata file)
You need to provide infomration about the time point in a design file
That being said, you shouldn’t be doing metastats and then lefse. The two test are not independent. You should pick one and use that. Normally, I’d suggest lefse. However, if you have paired data where you are analyzing the same entity at two time points, you really need to perform a paired analysis. You could do that with a wilcoxon test in R to see if there are taxa that have a change that is significantly different from 0.