Problems with the abundance table by group

Hello everyone! I'm trying to create a table with relative abundance of phyla and genre by groups of patient I have in my dataset. The groups are HPV, CT, and healthy. The problem is when I subset my phyloseq object "ps1" by group I get the same results. Like there is no difference in the number of phyla in ps1 and e.g. the healthy group. I'm going to leave the code and the output below. The healthy group is composed of
The code:
tab <- table(tax_table(ps1)[, "Phylum"], exclude = NULL)

tab

#subset phyloseq object
Healthy_group <- subset_samples (ps1, Group=="Healthy")

tab1 <- table(tax_table(Healthy_group)[, "Phylum"], exclude = NULL)
tab1

I tried checking Healthy_group and ps1:

image

They have the same taxa number. Btw it's the same result for the other 2 groups.
Can someone please help me out? I don't know what I'm doing wrong.
Thank you!

1 Like

Remember put a reproducible example for better help all the community.

Im saw that you select only the Phylum column in both cases.

Yes, I wanted to compare phylum abundance between different groups. But I get the same result everytime. I'll try to post a reproducible example. Thank you.

Please, attach data so we can see the structure of it. You can use dput function if the dataset is not too big

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

If you have a query related to it or one of the replies, start a new topic and refer back with a link.