# Taxa abundance analysis

**URL:** https://forum.qiime2.org/t/taxa-abundance-analysis/845
**Category:** User Support
**Tags:** ancom, taxonomy
**Created:** [July 17, 2017, 7:23pm UTC](https://forum.qiime2.org/t/taxa-abundance-analysis/845 "2017-07-17T19:23:13Z")
**Posts on this page:** 8
**Page:** 1

<div class="post-metadata">

### Author: ![Bing](https://forum.qiime2.org/letter_avatar_proxy/v4/letter/b/839c29/32.png) [@Bing](https://forum.qiime2.org/u/Bing)
#### Post date: [July 17, 2017, 7:23pm UTC](https://forum.qiime2.org/t/taxa-abundance-analysis/845/1 "2017-07-17T19:23:13Z")

</div>

Hi,  
I am curious about the taxa abundance analysis that reports the percentile of abundance. Is there anyway to get the statistical values or pairwise comparisons in QIIME 2?

Thanks,  
BIng

---

<div class="post-metadata">

### Author: ![jairideout](https://forum.qiime2.org/user_avatar/forum.qiime2.org/jairideout/32/9_2.png) [@jairideout](https://forum.qiime2.org/u/jairideout)
#### Post date: [July 19, 2017, 5:31pm UTC](https://forum.qiime2.org/t/taxa-abundance-analysis/845/3 "2017-07-19T17:31:56Z")

</div>

Hi @Bing! I'm not sure what analyses you're referring to -- are you asking about the QIIME 2 taxa barplots generated by `qiime taxa barplot`? If so, we just added an "export" feature to the taxa barplots so that you can export the underlying data that's being plotted as a delimited text file. This export feature will be available in the next QIIME 2 release (2017.7, scheduled for sometime next week).

If that doesn't answer your question, can you provide more details about what you're trying to accomplish? For example, is there a corresponding QIIME 1 script you've used in the past, or can you point me to a paper describing or using this type of analysis?

---

<div class="post-metadata">

### Author: ![Bing](https://forum.qiime2.org/letter_avatar_proxy/v4/letter/b/839c29/32.png) [@Bing](https://forum.qiime2.org/u/Bing)
#### Post date: [July 19, 2017, 10:03pm UTC](https://forum.qiime2.org/t/taxa-abundance-analysis/845/5 "2017-07-19T22:03:12Z")

</div>

Hello @jairideout,

I meant the last step in the tutorial: Differential abundance analysis. After the ANCOM analysis, we can find the different taxa percentile in terms of one variable such as group. My question is that whether we can get specific p-values or pairwise comparisons in the ANCOM analysis in QIIME 2 directly or we have to run it again in R?

Thanks,  
Bing

---

<div class="post-metadata">

### Author: ![jairideout](https://forum.qiime2.org/user_avatar/forum.qiime2.org/jairideout/32/9_2.png) [@jairideout](https://forum.qiime2.org/u/jairideout)
#### Post date: [July 19, 2017, 10:46pm UTC](https://forum.qiime2.org/t/taxa-abundance-analysis/845/7 "2017-07-19T22:46:45Z")

</div>

Thanks for clarifying! My understanding is that ANCOM doesn't report p-values, so I'm not sure that running ANCOM in R would get you that info.

I don't think there's a way to get pairwise comparisons directly in QIIME 2. The [scikit-bio implementation](http://scikit-bio.org/docs/latest/generated/generated/skbio.stats.composition.ancom.html) of ANCOM (which QIIME 2 is using) doesn't support pairwise comparisons either. I'm not certain if the R implementation supports pairwise comparisons, you may want to verify that before exporting and loading your data into R.

@mortonjt do you have any suggestions for pairwise comparisons and/or p-values in ANCOM?

---

<div class="post-metadata">

### Author: ![Bing](https://forum.qiime2.org/letter_avatar_proxy/v4/letter/b/839c29/32.png) [@Bing](https://forum.qiime2.org/u/Bing)
#### Post date: [July 20, 2017, 4:29pm UTC](https://forum.qiime2.org/t/taxa-abundance-analysis/845/9 "2017-07-20T16:29:40Z")

</div>

@jairideout  
so as presented in the tutorial example, only the significant different abundances are shown, right? If not, where can we get that "Reject null hypothesis" table based on the reference [sciket-bio implementation](http://scikit-bio.org/docs/latest/generated/generated/skbio.stats.composition.ancom.html)?

Thanks,  
Bing

---

<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: [July 21, 2017, 1:09am UTC](https://forum.qiime2.org/t/taxa-abundance-analysis/845/10 "2017-07-21T01:09:33Z")

</div>

@jairideout is right - ANCOM doesn't support pvalues or pairwise comparisons. But you can get the W-statistic, which is a proxy for the statistical significance.

---

<div class="post-metadata">

### Author: ![Bing](https://forum.qiime2.org/letter_avatar_proxy/v4/letter/b/839c29/32.png) [@Bing](https://forum.qiime2.org/u/Bing)
#### Post date: [July 21, 2017, 4:34pm UTC](https://forum.qiime2.org/t/taxa-abundance-analysis/845/11 "2017-07-21T16:34:20Z")

</div>

Thank you so much! Bing

---

<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: [August 21, 2017, 10:34pm UTC](https://forum.qiime2.org/t/taxa-abundance-analysis/845/12 "2017-08-21T22:34:30Z")

</div>

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