Skip to content
Prev 6015 / 7420 Next

Rarefaction curve in phyloseq

Hi,

I am trying to plot rarefaction plot but did not understand how to modify
below code to plot rarefaction based on metadata column instead of sample
name in phyloseq

Like instead of sample name, I want to plot based on season:

 sample_data(ps)
Sample Data:        [34 samples by 4 sample variables]:
                Season Depth Month Year
B0214D1-PL1-D1   Winter     1     2   14
B0214D2-PL1-E1   Winter    10     2   14
B0214D3-PL1-F1   Winter    20     2   14
B0314D1-PL1-G1   Spring     1     3   14
B0314D2-PL1-H1   Spring    10     3   14

rarecurve(t(otu_table(ps)), step=50, cex=0.5)

What does step 50 indicate here?

Thanks
*Yogesh *