This is a test version of Biostars. For the public version, visit https://www.biostars.org.
subset sample in phyloseq

I do have three tissue and three treatment conditions in physeq object. My aim is to only want to use the soil sample and want to see the treatment effect on microbial abundance.

How should I remove other tissue and only retain the soil sample in phyloseq object.

sample details:

summary(sample_data(physeq)$Tissue)

Leaf   Root       Soil

And each tissue has given 3 treatment.

summary(sample_data(physeq)$Treatment)

T1     T2     T3

Thanks

r next-gen
• 14,671 views • link
updated by GenoMax GenoMax • written by Bioinfonext Bioinfonext

1 answer

You should read the phyloseq documentation for this, there are a ton of examples for subsetting or pruning data. Also:

https://joey711.github.io/phyloseq/preprocess.html

More specifically, look up the online documentation in R for subset_samples, prune_samples, subset_taxa, and prune_taxa. I also recommend the F1000 paper which is a pretty great walkthrough for phyloseq.

• 0 views • link
written by Chris Fields Chris Fields

Thanks a lot, I tried below code but it is showing error, could please suggest how I can correct this mistake:

> physeq

phyloseq-class experiment-level object

otu_table()   OTU Table:         [ 67405 taxa and 162 samples ]

sample_data() Sample Data:       [ 162 samples by 5 sample variables ]

tax_table()   Taxonomy Table:    [ 67405 taxa by 7 taxonomic ranks ]

phy_tree()    Phylogenetic Tree: [ 67405 tips and 67229 internal nodes ]

> physeq_soil <- subset_samples(physeq, SampleType=="Soil")
Error in eval(e, x, parent.frame()) : object 'SampleType' not found
• 0 views • link
written by Bioinfonext Bioinfonext

Thanks, I got it, if I replace SampleType with Tissue then its works:

physeq_soil <- subset_samples(physeq, Tissue=="Soil")

Thanks

• 0 views • link
written by Bioinfonext Bioinfonext

Hello,

I subset my data (based on only one time point) in phyloseq and wanted to performed a PCoA only with this subset but I still have the all dataset that appear in my PCoA even if the unwanted samples are not any more labelled.

Global subset:

phyloseq-class experiment-level object
otu_table()   OTU Table:         [ 1222 taxa and 40 samples ]
sample_data() Sample Data:       [ 40 samples by 10 sample variables ]
tax_table()   Taxonomy Table:    [ 1222 taxa by 7 taxonomic ranks ]
phy_tree()    Phylogenetic Tree: [ 1222 tips and 1219 internal nodes ]

After subsetting:

phyloseq-class experiment-level object
otu_table()   OTU Table:         [ 1222 taxa and 20 samples ]
sample_data() Sample Data:       [ 20 samples by 10 sample variables ]
tax_table()   Taxonomy Table:    [ 1222 taxa by 7 taxonomic ranks ]
phy_tree()    Phylogenetic Tree: [ 1222 tips and 1219 internal nodes ]

Any idea? Thanks,

• 0 views • link
written by nichecanicheca • nichecanicheca

You only have subset samples not taxa, you need to remove taxa which have 0 counts in your all samples. look at the prune_taxa command.

• 0 views • link
written by Bioinfonext Bioinfonext

Log in to answer this question.

More posts like this