I’m working on a project that uses MaAsLin2 for multivariable modeling of feature (taxa)-metadata association. My ASVs are agglomerated at the genus level, and I have about ~200 of them. I’ve run several different models (adjusting for different sets of potential confounders) and my results are null: no features associated with my exposure of interest.
I have also run my data through the QIIME2 PICRUSt2 plugin and have a feature table of pathways as well as mapfiles. Can I run a similar MaAsLin2 pipeline for these pathway features, in a similar manner to the standard taxa pipeline?
- Edit: I just realized that the documentation specifies that we can use genes and/or pathways with MaAsLin2 as well. - Nevermind!