The P values generated by Maaslin3

Downstream analysis and statistics MaAsLin
1

Hi,

Firstly, thank you so much for developing Maaslin3! We have done some analyses on species abundance using Maaslin3 and we had a few questions regarding the significant outputs:

  1. If we want to draw forest plot using the significant outputs, is it correct to generate the 95% confidence intervals manually by “coef ± 1.96*stderr”?

  2. We have noticed that in some of our outputs, the P values and Q values denote significant results but the CIs contain 0 when we tried to generate CIs (as shown in the picture below). We are wondering whether the P values have been adjusted by other methods or we are using the wrong formula to generate CIs.

    2d3efcc3c248d43d65ee4f81e43bef09
    2d3efcc3c248d43d65ee4f81e43bef091288×370 31 KB

Thanks a lot for your help!

2

Hi,

  1. Yes, that would produce the correct confidence intervals for the relative abundance associations. (Note that these won’t be FDR-corrected, but most people just don’t FDR correct confidence intervals.)
  2. By default, MaAsLin compares coefficients to the median of that metadatum’s coefficients in the hypothesis test in order to test for absolute abundance associations. (The paper has more details on why and how this works.) This means that instead of the usual hypothesis test against 0, it’s testing against the value in the null_hypothesis column. If you don’t care about the absolute abundance associations and just want the relative abundance associations, you can set median_comparison_abundance to FALSE and it’ll hypothesis test against 0. This will give you the expected results with CIs containing 0 iff the association is non-significant.

Will

3

Thank you so much for the explanation!