This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Rarefaction curves for OTU Data

Hi, I'm trying to do a rarefaction curve in R using OTUs but I have no idea how around to do them. How would I procede to do that?

My data is metagenomic 16s sequencing and the files I have are divided by taxonomic level. The genus level for example has the genus on one column and the number of OTUs that align with that genus on the other column. Any idea?

next-gen otus metagenome
• 6,027 views • link
updated by Steven Lakin Steven Lakin • written by ricardo.couto • ricardo.couto

1 answer

Rarefaction is the number of unique OTUs described as a function of the number of units (reads, usually) sampled. So to produce the curve, you have to produce:

number of OTUs that align with that genus on the other column

for each level of subsampling depth you want to investigate.

  1. Subsample your raw data, for example, every 10% from 10 -> 100%
  2. For each of those, find the number of OTUs described at that taxonomic level
  3. Plot the curve with # of unique OTUs on the y axis and the subsampling depth on the x axis

Pipelines like mothur and Qiime have these functions built in for 16S sequences.

• 0 views • link
written by Steven Lakin Steven Lakin

The shuf util can easily be applied to this kind of subsampling. Basically you just increase -n for every round..

• 0 views • link
written by 5heikki 5heikki

Log in to answer this question.

More posts like this