Dear MaAsLin2 developpers and users,
Is there any way to control or at least test the effect of sequencing depth to make sure observed differences are not resulting from underlining differences in sequencing depth between groups/samples.
I typically do that using anova in order to check if sequencing depth explain differences in alpha diversity e.g., alpha_metric ~ Group * LibrairySize and i can confirm if Librairy Size has an effect alone or depending ou the Group I am testing.
I have added the number of reads I subbmitted to humann2 in my metadata and when running MaAsLin2 I detect and see strong significant correlations, but is there any syntax to check for Sequencing Depth as an interaction factor to other factor tested through MaAsLin2?
Thanks a ton.