How to compare taxonomy from a single sample against different groups

Greetings to all,

I am trying to analyze some data from human fecal samples. From the taxonomic tables obtained I need to compare one single sample against several groups of samples (each group correspond with a group of people with certain condition) and asses if the single sample is clustered whith some group in particular or not (or if it is more simmilar/dissimilar to a certain group). And also, which taxons are differentially abundant between the sample and each group.

Until now, I have performed NMDS and PCoa with Anosim to test the clusters.

My question is: wich statistic test or wich tool should I use to asses which taxons are differentially abundant between the sample and each group?
Could ancom serve for this?

(Also I tested Lefse, but it does not accept one single sample as input to compare).

Any help would be appreciated, kind regards

Lucía

Greetings to all,

I would like to know if is possible to use ANCOM for assesing which taxons are differentially abundant between a single sample and a group of samples.

Any help would be appreciated, kind regards

Lucía

Please read the following before posting!

Is this post about a Third-party plugin (i.e., plugins not included in the Core Distribution)? Have you reviewed the QIIME 2 Forum Glossary? Please do not post questions here that have to do with interpretation of results, or general discussion, or technical support for “core” plugins. Posts in this category are not guaranteed a response.

Hi @LuciaGG,

I merged your post with your previous one, because it seems to me they are on the same topic.
I understand it may take long time to get an answer, but please be patients, moderator as well as other user forums are doing their best to help anyone asking in here but it still may take time!

On your specific question, I can try to help but I am not a statistician so I do not have the full picture.
However, as a general observation, it seems a very peculiar experimental design if you have to compare a single sample to one or more sample groups (each group containing n biological replicates).
Could you explain better your experimental design please (or why you end up on having only a single sample? Is anything happening before or during the qiime2 analysis?)
To me, your main problem is to assess how representative your single sample is for is own group! How do you know you are not dealing with an ‘outlier’?
I don’t know how informative would be comparing one sample vs a sample group. I honestly can not tell if there is any statistic comparison tool which can deal with sample size=1, without breaking its beyond assumption.

Hope it helps
Luca

2 Likes