# ANCOM clr results confusion

**URL:** https://forum.qiime2.org/t/ancom-clr-results-confusion/8794
**Category:** User Support
**Tags:** ancom
**Created:** [March 15, 2019, 7:42pm UTC](https://forum.qiime2.org/t/ancom-clr-results-confusion/8794 "2019-03-15T19:42:20Z")
**Posts on this page:** 12
**Page:** 1

<div class="post-metadata">

### Author: ![Zach\_Burcham](https://forum.qiime2.org/user_avatar/forum.qiime2.org/zach_burcham/32/18784_2.png) [@Zach\_Burcham](https://forum.qiime2.org/u/Zach_Burcham)
#### Post date: [March 15, 2019, 7:42pm UTC](https://forum.qiime2.org/t/ancom-clr-results-confusion/8794/1 "2019-03-15T19:42:20Z")

</div>

Hello,

I am having some confusion with my ANCOM results. From my understanding, the clr is like the fold-change and the 100% abundance shows the counts in each group. I would expect if I have group A and B then anytime I have more counts in A than B I would obtain either positive or negative clr's consistently. But in my data set, I have two significant taxa that are both higher in group A than B, but one of them gives a negative clr while the other gives a positive clr. How is this possible? Do I report these taxa as being more abundant in group A based on the percentile abundance or do I report one as more abundant and the other less since the clr's are positive and negative even though the percentile abundace says both are more abundant?

Thanks,  
Zach

---

<div class="post-metadata">

### Author: ![Nicholas\_Bokulich](https://forum.qiime2.org/user_avatar/forum.qiime2.org/nicholas_bokulich/32/19937_2.png) [@Nicholas\_Bokulich](https://forum.qiime2.org/u/Nicholas_Bokulich)
#### Post date: [March 18, 2019, 12:31pm UTC](https://forum.qiime2.org/t/ancom-clr-results-confusion/8794/2 "2019-03-18T12:31:53Z")

</div>



---

<div class="post-metadata">

### Author: ![mortonjt](https://forum.qiime2.org/user_avatar/forum.qiime2.org/mortonjt/32/3488_2.png) [@mortonjt](https://forum.qiime2.org/u/mortonjt)
#### Post date: [March 18, 2019, 2:39pm UTC](https://forum.qiime2.org/t/ancom-clr-results-confusion/8794/3 "2019-03-18T14:39:23Z")

</div>

Hi @Zach_Burcham - having a plot showing the results could help us better answer this question.

clr coordinate do not necessarily imply increase / decrease (since it is not possible do this in the first place). Instead, these coordinates represent log-fold change relative to the _average_ microbe.

So it is possible that both taxa are higher, but one of them is increased greater than the average microbe and the other is increased less than the average microbe (provided that the average microbe is also increased between these two conditions).

---

<div class="post-metadata">

### Author: ![Zach\_Burcham](https://forum.qiime2.org/user_avatar/forum.qiime2.org/zach_burcham/32/18784_2.png) [@Zach\_Burcham](https://forum.qiime2.org/u/Zach_Burcham)
#### Post date: [March 18, 2019, 2:54pm UTC](https://forum.qiime2.org/t/ancom-clr-results-confusion/8794/4 "2019-03-18T14:54:46Z")

</div>

[ancom-adult\_or\_youth-L6.qzv](https://cdck-file-uploads-global.s3.dualstack.us-west-2.amazonaws.com/flex002/uploads/qiime21/original/2X/3/3825b79f56881f7955ab438b8447b7325a8ba1e9.qzv) (465.9 KB)

Hi @mortonjt - Thank you for your response. I have attached a .qzv file from my data to better show what I mean. For example, Treponema and Lautropia both show higher counts in the Adult samples and the CLR is practically opposite. So from my understanding of the "average" microbe approach, CLR isn't a way to determine which sample group contained higher relative abundances. So is it appropriate to say that the relative abundances of Treponema and Lautropia are higher in Adult samples than Youth when reporting these results?

---

<div class="post-metadata">

### Author: ![mortonjt](https://forum.qiime2.org/user_avatar/forum.qiime2.org/mortonjt/32/3488_2.png) [@mortonjt](https://forum.qiime2.org/u/mortonjt)
#### Post date: [March 18, 2019, 3:11pm UTC](https://forum.qiime2.org/t/ancom-clr-results-confusion/8794/5 "2019-03-18T15:11:17Z")

</div>

Correct - no one can actually do differential abundance without absolute abundances: see below (the clr transform is also described there)

[Analysis and correction of compositional bias in sparse sequencing count data | BMC Genomics | Full Text](https://bmcgenomics.biomedcentral.com/articles/10.1186/s12864-018-5160-5)

You can report Treponema and Lautropia have differing relative abundances according to ANCOM. I wouldn't put too much weight on the clr transform in the plot -- that is just to help visualize the effect size.

---

<div class="post-metadata">

### Author: ![Zach\_Burcham](https://forum.qiime2.org/user_avatar/forum.qiime2.org/zach_burcham/32/18784_2.png) [@Zach\_Burcham](https://forum.qiime2.org/u/Zach_Burcham)
#### Post date: [March 18, 2019, 3:13pm UTC](https://forum.qiime2.org/t/ancom-clr-results-confusion/8794/6 "2019-03-18T15:13:37Z")

</div>

Great, thanks for the help!

---

<div class="post-metadata">

### Author: ![Nicholas\_Bokulich](https://forum.qiime2.org/user_avatar/forum.qiime2.org/nicholas_bokulich/32/19937_2.png) [@Nicholas\_Bokulich](https://forum.qiime2.org/u/Nicholas_Bokulich)
#### Post date: [March 18, 2019, 4:15pm UTC](https://forum.qiime2.org/t/ancom-clr-results-confusion/8794/7 "2019-03-18T16:15:03Z")

</div>



---

<div class="post-metadata">

### Author: ![Zach\_Burcham](https://forum.qiime2.org/user_avatar/forum.qiime2.org/zach_burcham/32/18784_2.png) [@Zach\_Burcham](https://forum.qiime2.org/u/Zach_Burcham)
#### Post date: [April 15, 2019, 6:34pm UTC](https://forum.qiime2.org/t/ancom-clr-results-confusion/8794/8 "2019-04-15T18:34:17Z")

</div>

Hi @mortonjt, one more question. When looking at my feature table, it shows higher counts of Neisseria in my youth group samples compared to adults, but when I run ANCOM I get that there is a significant difference between my age groups with Neisseria being higher in my adult samples. Why does there seem to be discrepancy here? The other genera ANCOM is telling me are significantly different correspond to the feature table except one.

---

<div class="post-metadata">

### Author: ![thermokarst](https://forum.qiime2.org/letter_avatar_proxy/v4/letter/t/8e7dd6/32.png) [@thermokarst](https://forum.qiime2.org/u/thermokarst)
#### Post date: [April 16, 2019, 2:07pm UTC](https://forum.qiime2.org/t/ancom-clr-results-confusion/8794/9 "2019-04-16T14:07:37Z")

</div>



---

<div class="post-metadata">

### Author: ![mortonjt](https://forum.qiime2.org/user_avatar/forum.qiime2.org/mortonjt/32/3488_2.png) [@mortonjt](https://forum.qiime2.org/u/mortonjt)
#### Post date: [April 16, 2019, 2:11pm UTC](https://forum.qiime2.org/t/ancom-clr-results-confusion/8794/10 "2019-04-16T14:11:49Z")

</div>

Really not crazy about using the terms `lower` or `higher` in the context of relative abundances, since relative abundances can't tell you which microbes actually increased or decreased.

I'm not an expert on ANCOM - best consult with the corresponding author on the [ANCOM paper](https://www.ncbi.nlm.nih.gov/pubmed/26028277)

---

<div class="post-metadata">

### Author: ![Nicholas\_Bokulich](https://forum.qiime2.org/user_avatar/forum.qiime2.org/nicholas_bokulich/32/19937_2.png) [@Nicholas\_Bokulich](https://forum.qiime2.org/u/Nicholas_Bokulich)
#### Post date: [April 16, 2019, 2:49pm UTC](https://forum.qiime2.org/t/ancom-clr-results-confusion/8794/11 "2019-04-16T14:49:03Z")

</div>



---

<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: [May 18, 2019, 4:39am UTC](https://forum.qiime2.org/t/ancom-clr-results-confusion/8794/17 "2019-05-18T04:39:08Z")

</div>

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