SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Rarefaction lisa.brown Bioinformatics 2 02-11-2014 07:53 AM
Species abundance and rarefaction in metagenomics data OliverDeusch Metagenomics 1 12-13-2012 08:17 AM
dexseq simulation Pedrissimo RNA Sequencing 1 06-08-2012 01:47 AM
Rarefaction distance explanation lletourn Metagenomics 2 03-11-2011 08:09 AM
Rarefaction/saturation curve based on NGS data BioGenomics Bioinformatics 0 12-08-2010 06:41 PM

Reply
 
Thread Tools
Old 07-15-2014, 01:30 AM   #1
sphil
Senior Member
 
Location: Stuttgart, Germany

Join Date: Apr 2010
Posts: 192
Default rarefaction curves and simulation

Hiho

i want to create a rarefaction curve using R i.e. vegan {rearecurve/speccaccum}.

This worked out quite well (please see here)

Anyways, as you can see in the figures the curves don't reach saturation. So I know that i have to kind of fit or simulate the missing data but this is exactly the problem. How can I produce a fit or something similar that shows me how large my sample has to be to reach saturation in discovered species?



Thanks for your time and help!



ps. code:

(raremax <- min(rowSums(t(species))))
Srare <- rarefy(t(species), raremax)
plot(specnumber(t(species)), Srare, xlab = "Observed No. of Species", ylab = "Rarefied No. of Species")
abline(0, 1)
rarecurve(t(species), step = 20, sample = raremax, col = "blue", cex = 0.6)

(crosspost on biostars)
sphil is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 07:21 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2019, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO