# ANCOM Analysis?

**URL:** https://forum.qiime2.org/t/ancom-analysis/5910
**Category:** Other Bioinformatics Tools
**Created:** [September 5, 2018, 4:27pm UTC](https://forum.qiime2.org/t/ancom-analysis/5910 "2018-09-05T16:27:11Z")
**Posts on this page:** 11
**Page:** 1

<div class="post-metadata">

### Author: ![idapantoja](https://forum.qiime2.org/user_avatar/forum.qiime2.org/idapantoja/32/2145_2.png) [@idapantoja](https://forum.qiime2.org/u/idapantoja)
#### Post date: [September 5, 2018, 4:27pm UTC](https://forum.qiime2.org/t/ancom-analysis/5910/1 "2018-09-05T16:27:11Z")

</div>

Hi all,

I did my microbiome analysis already using QIIME 2 v 2017.9 but now I want to do some statistical analysis to confirm the changes I am seeing in my taxa box plot. I was wondering if ANCOM analysis will help on this or if you can suggest another test I can use in this case in QIIME2 or even any other tool like R. Basically, what I want to know if which changes in the organisms abundances are significantly different or not across my samples.  
A would appreciate any suggestion.

Thanks a lot!  
idapantoja

---

<div class="post-metadata">

### Author: ![timanix](https://forum.qiime2.org/user_avatar/forum.qiime2.org/timanix/32/17879_2.png) [@timanix](https://forum.qiime2.org/u/timanix)
#### Post date: [September 6, 2018, 12:59pm UTC](https://forum.qiime2.org/t/ancom-analysis/5910/2 "2018-09-06T12:59:51Z")

</div>

Hi, I think you already answered your question - you can use ANCOM in q2 to reveal which features are significantly changed

---

<div class="post-metadata">

### Author: ![idapantoja](https://forum.qiime2.org/user_avatar/forum.qiime2.org/idapantoja/32/2145_2.png) [@idapantoja](https://forum.qiime2.org/u/idapantoja)
#### Post date: [September 6, 2018, 3:05pm UTC](https://forum.qiime2.org/t/ancom-analysis/5910/3 "2018-09-06T15:05:51Z")

</div>

Thanks so much! Wanted to be sure!

Idapantoja

---

<div class="post-metadata">

### Author: ![idapantoja](https://forum.qiime2.org/user_avatar/forum.qiime2.org/idapantoja/32/2145_2.png) [@idapantoja](https://forum.qiime2.org/u/idapantoja)
#### Post date: [September 6, 2018, 5:50pm UTC](https://forum.qiime2.org/t/ancom-analysis/5910/4 "2018-09-06T17:50:49Z")

</div>

Hey again,  
In the Volcano Plot, what "W" stands for?

Thanks!  
Ida

---

<div class="post-metadata">

### Author: ![timanix](https://forum.qiime2.org/user_avatar/forum.qiime2.org/timanix/32/17879_2.png) [@timanix](https://forum.qiime2.org/u/timanix)
#### Post date: [September 6, 2018, 5:58pm UTC](https://forum.qiime2.org/t/ancom-analysis/5910/5 "2018-09-06T17:58:54Z")

</div>

The bigger meaning, the more difference, if I am not mistaking

---

<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: [September 6, 2018, 7:16pm UTC](https://forum.qiime2.org/t/ancom-analysis/5910/6 "2018-09-06T19:16:25Z")

</div>

Check out the original publication for more details, but here's a great explanation on this forum: [Specify W cutoff for anacom? - #10 by mortonjt](https://forum.qiime2.org/t/specify-w-cutoff-for-anacom/1844/10)

---

<div class="post-metadata">

### Author: ![idapantoja](https://forum.qiime2.org/user_avatar/forum.qiime2.org/idapantoja/32/2145_2.png) [@idapantoja](https://forum.qiime2.org/u/idapantoja)
#### Post date: [September 7, 2018, 8:02pm UTC](https://forum.qiime2.org/t/ancom-analysis/5910/7 "2018-09-07T20:02:55Z")

</div>

Thanks!!

I have another question about the Ancom command. This is what I am trying to run:  
qiime composition ancom   
--i-table MRE-comp-taxonomy\_ResistantStarch\_Only65\_orderedbytimepoints-l7.qza   
--m-metadata-file mapping\_file\_resistant\_starch\_only65\_orderedbytimepoints\_09\_07\_18\_ANCOM   
--m-metadata-category Diet\_Time\_Point   
--o-visualization MRE-l7-ancom-Diet\_Time\_Point.qzv

It worked but in the category I want to be more specific for the analysis and selected only six samples of about 33 samples for the "Diet\_Time\_Point" category. How can I do that? So far, I am getting the analysis for all the samples in that category and the abundance table is big.

I would appreciate your help!

Thanks!  
Idapantoja

---

<div class="post-metadata">

### Author: ![timanix](https://forum.qiime2.org/user_avatar/forum.qiime2.org/timanix/32/17879_2.png) [@timanix](https://forum.qiime2.org/u/timanix)
#### Post date: [September 8, 2018, 6:01pm UTC](https://forum.qiime2.org/t/ancom-analysis/5910/8 "2018-09-08T18:01:06Z")

</div>

There is a tutorial "Filtering data", you can found there how to filter your table according different parameters, including samples.

---

<div class="post-metadata">

### Author: ![idapantoja](https://forum.qiime2.org/user_avatar/forum.qiime2.org/idapantoja/32/2145_2.png) [@idapantoja](https://forum.qiime2.org/u/idapantoja)
#### Post date: [September 10, 2018, 6:24pm UTC](https://forum.qiime2.org/t/ancom-analysis/5910/9 "2018-09-10T18:24:13Z")

</div>

Thanks! I looked it and is working.

What is happening if I don't I Volcano Plot? I am doing something wrong?

Unable to generate volcano plot, please check the ANCOM statistical results (above).

idapantoja

---

<div class="post-metadata">

### Author: ![ebolyen](https://forum.qiime2.org/user_avatar/forum.qiime2.org/ebolyen/32/11_2.png) [@ebolyen](https://forum.qiime2.org/u/ebolyen)
#### Post date: [September 10, 2018, 7:17pm UTC](https://forum.qiime2.org/t/ancom-analysis/5910/10 "2018-09-10T19:17:51Z")

</div>

Hey @idapantoja,

You probably need to update your installation, `2017.9` is pretty old at this point. I do remember something about the volcano plot being conditional, but I don't remember the details, @mortonjt was that the case?

---

<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: [September 10, 2018, 7:45pm UTC](https://forum.qiime2.org/t/ancom-analysis/5910/11 "2018-09-10T19:45:43Z")

</div>

Correct, you may be experiencing an outdated version of bokeh
