I have a question regarding the metastats function in mothur.
I would like to analyze the differential abundance of the taxa (output of the phylotype command - I’m unable to do the cluster.split due to restricted computational resources) in my groups. For this purpose, I found the metastats function. In the documentation, I found that the function can determine differential abundance for two different groups at a time.
My question is the following:
Can I use the metastats function to do the differential abundance analysis for 4 different groups?
If so, how would you recommend me to do this?
( I found the following article suggesting in their discussion to replace the t-test by a one-way anova)