[R-sig-eco] Rarefaction curve in phyloseq

[Yogesh Gupta]

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 *

	[[alternative HTML version deleted]]