# ANCOM result interpretation?

**URL:** https://forum.qiime2.org/t/ancom-result-interpretation/22841
**Category:** User Support
**Tags:** unifrac, ancom, beta
**Created:** [April 21, 2022, 5:31pm UTC](https://forum.qiime2.org/t/ancom-result-interpretation/22841 "2022-04-21T17:31:43Z")
**Posts on this page:** 11
**Page:** 1

<div class="post-metadata">

### Author: ![Yavnika](https://forum.qiime2.org/letter_avatar_proxy/v4/letter/y/e19b73/32.png) [@Yavnika](https://forum.qiime2.org/u/Yavnika)
#### Post date: [April 21, 2022, 5:31pm UTC](https://forum.qiime2.org/t/ancom-result-interpretation/22841/1 "2022-04-21T17:31:43Z")

</div>

Hello,

I'm new to the forum and trying to figure out a few things regarding abundance testing. I have 2 groups (test and control n=8, n=7 respectively) and I see significant differences between the 2 groups in terms of alpha and beta diversity analysis as well (Unweighted p-value 0.003, weighted p-value 0.024 with PERMANOVA) with separation on the emperor plots. My problem is when I try to use ANCOM to look at the differential abundance, I don't get any significant results. I tried filtering at different parameters as well but that didn't work either. When I tried the differential abundance test at a specific taxonomic level, it gives me 2 families with W score of 18,19. Another differential abundance method I tried using was q2dsfdr which gave me multiple taxa that were different but since it only gives out the p-values I can't ascertain which group the taxa belongs to.  
So, I have a few questions regarding the differential abundance testing:

1. Is it possible that although my alpha and beta diversity are significantly different among groups, I just don't have enough power to use ANCOM to actually pin point which specific groups are different ?
2. Can I consider the results given by ANCOM specific taxonomic level analysis to be valid?
3. Are the results from q2 dsfdr valid and if so how can they be processed for me to actually find out which group the different taxa belong to? \*\*Also, I tried using q2dbbact as well but didn't get any results from that either.  
I'm not sure what is actually going on with my data but I would really appreciate some help. Also, if there are any other differential abundance testing methods I can use.

Thanks.

 ![Screen Shot 2022-04-21 at 12.10.38 PM](https://forum-qiime2-org.s3.dualstack.us-west-2.amazonaws.com/original/2X/a/a19b2ec79cea59e747be1eeddc096d3425ef87fd.png)

---

<div class="post-metadata">

### Author: ![lizgehret](https://forum.qiime2.org/user_avatar/forum.qiime2.org/lizgehret/32/17128_2.png) [@lizgehret](https://forum.qiime2.org/u/lizgehret)
#### Post date: [May 17, 2022, 5:57pm UTC](https://forum.qiime2.org/t/ancom-result-interpretation/22841/5 "2022-05-17T17:57:21Z")

</div>



---

<div class="post-metadata">

### Author: ![lizgehret](https://forum.qiime2.org/user_avatar/forum.qiime2.org/lizgehret/32/17128_2.png) [@lizgehret](https://forum.qiime2.org/u/lizgehret)
#### Post date: [May 23, 2022, 7:52pm UTC](https://forum.qiime2.org/t/ancom-result-interpretation/22841/8 "2022-05-23T19:52:14Z")

</div>

Hi @Yavnika,

Welcome back to the :qiime2: forum! Apologies for the delay in response on this.

> [@Yavnika](#):
>
> Is it possible that although my alpha and beta diversity are significantly different among groups, I just don't have enough power to use ANCOM to actually pin point which specific groups are different ?

Yes! This can absolutely be the case. You may just not be seeing significant results at the ASV level, which is not uncommon. You may see more significant results by zooming out, say, to the genus level.

> [@Yavnika](#):
>
> Can I consider the results given by ANCOM specific taxonomic level analysis to be valid?

Yes! 🤓 🧬  
The pattern and differential abundance might only emerge at the genus level.

> [@Yavnika](#):
>
> Are the results from q2 dsfdr valid and if so how can they be processed for me to actually find out which group the different taxa belong to? \*\*Also, I tried using q2dbbact as well but didn't get any results from that either.

I would ask @serenejiang about the validity since they are the plugin developer who created q2-dsfdr; it looks like this plugin hasn't been maintained over the last few years. That being said, you may just not be getting significant results because there actually aren't any for your dataset at this level.

> [@Yavnika](#):
>
> Also, if there are any other differential abundance testing methods I can use.

These are a couple of other community plugins that can be used for differential abundance testing that are currently maintained:

[DEICODE](https://library.qiime2.org/plugins/deicode/19/)  
[q2-aldex2](https://library.qiime2.org/plugins/q2-aldex2/24/)

In addition to these plugins, the following can be used for differential abundance visualizations, which you may find useful as well:

[Qurro](https://library.qiime2.org/plugins/qurro/22/)

Hope this helps!  
Cheers 🦎

---

<div class="post-metadata">

### Author: ![lizgehret](https://forum.qiime2.org/user_avatar/forum.qiime2.org/lizgehret/32/17128_2.png) [@lizgehret](https://forum.qiime2.org/u/lizgehret)
#### Post date: [May 23, 2022, 7:52pm UTC](https://forum.qiime2.org/t/ancom-result-interpretation/22841/9 "2022-05-23T19:52:31Z")

</div>



---

<div class="post-metadata">

### Author: ![Yavnika](https://forum.qiime2.org/letter_avatar_proxy/v4/letter/y/e19b73/32.png) [@Yavnika](https://forum.qiime2.org/u/Yavnika)
#### Post date: [May 24, 2022, 9:51pm UTC](https://forum.qiime2.org/t/ancom-result-interpretation/22841/10 "2022-05-24T21:51:02Z")

</div>

Hi @lizgehret,

Thank you so much for your response. This is really helpful. I tried to move my data to R and use ANCOM-BC for analysis which does show some groups that were differentially abundant but I'm not sure if I was supposed to filter the sequences before using that or not. Another analysis tool I used is LEFse which also showed taxa that were differentially abundant. I am just confused as to which method is more reliable as I do see some common taxa among these different analysis tools but the overall results are different. I would really appreciate if you could please help elaborate on this a little.  
Thank you.

---

<div class="post-metadata">

### Author: ![lizgehret](https://forum.qiime2.org/user_avatar/forum.qiime2.org/lizgehret/32/17128_2.png) [@lizgehret](https://forum.qiime2.org/u/lizgehret)
#### Post date: [May 24, 2022, 9:51pm UTC](https://forum.qiime2.org/t/ancom-result-interpretation/22841/11 "2022-05-24T21:51:54Z")

</div>



---

<div class="post-metadata">

### Author: ![lizgehret](https://forum.qiime2.org/user_avatar/forum.qiime2.org/lizgehret/32/17128_2.png) [@lizgehret](https://forum.qiime2.org/u/lizgehret)
#### Post date: [May 26, 2022, 3:35pm UTC](https://forum.qiime2.org/t/ancom-result-interpretation/22841/13 "2022-05-26T15:35:33Z")

</div>

Hi @Yavnika,

> [@Yavnika](#):
>
> I am just confused as to which method is more reliable as I do see some common taxa among these different analysis tools but the overall results are different.

There will be some variability between tools - the variability will also differ between tools depending on the dataset, and whether or not you are filtering prior to running DA. That being said, @jwdebelius has a paper in pre-print that has shown the following:

> _Limma voom, Wilcoxon, LEfSe, and edgeR tended to find the largest number of significant ASVs compared with other methods._

I would _highly recommend_ reading her paper (linked [here](https://www.nature.com/articles/s41467-022-28034-z)), as it will help clarify your questions/concerns in **much** greater detail than I can. 🤓

---

<div class="post-metadata">

### Author: ![lizgehret](https://forum.qiime2.org/user_avatar/forum.qiime2.org/lizgehret/32/17128_2.png) [@lizgehret](https://forum.qiime2.org/u/lizgehret)
#### Post date: [May 26, 2022, 3:35pm UTC](https://forum.qiime2.org/t/ancom-result-interpretation/22841/14 "2022-05-26T15:35:37Z")

</div>



---

<div class="post-metadata">

### Author: ![jwdebelius](https://forum.qiime2.org/user_avatar/forum.qiime2.org/jwdebelius/32/9655_2.png) [@jwdebelius](https://forum.qiime2.org/u/jwdebelius)
#### Post date: [May 26, 2022, 4:02pm UTC](https://forum.qiime2.org/t/ancom-result-interpretation/22841/15 "2022-05-26T16:02:10Z")

</div>

It's not actually my paper! I nominate papers for the most impactful paper/preprint of the year 🏅, and tI think this one will push the field forward. But yes, highly recommend it

---

<div class="post-metadata">

### Author: ![Yavnika](https://forum.qiime2.org/letter_avatar_proxy/v4/letter/y/e19b73/32.png) [@Yavnika](https://forum.qiime2.org/u/Yavnika)
#### Post date: [May 27, 2022, 4:03pm UTC](https://forum.qiime2.org/t/ancom-result-interpretation/22841/16 "2022-05-27T16:03:14Z")

</div>

Thank you so much, I'll be sure to check out the paper mentioned.

---

<div class="post-metadata">

### Author: ![system](https://forum-qiime2-org.s3.dualstack.us-west-2.amazonaws.com/original/3X/2/1/21af5fe23cb6f4579467c66a9ed94e55274ca7bd.svg) [@system](https://forum.qiime2.org/u/system)
#### Post date: [June 27, 2022, 10:04pm UTC](https://forum.qiime2.org/t/ancom-result-interpretation/22841/17 "2022-06-27T22:04:01Z")

</div>

This topic was automatically closed 31 days after the last reply. New replies are no longer allowed.
