Scaling with ranked subsampling curve (SRScurve) | R Documentation |
For each column of the input data, draws a line plot of alpha diversity indices (see metric) at different sample sizes (specified by step) normalized by scaling with ranked subsampling (using SRS). Minimum sample size (cutoff-level) can be evaluated by specifying sample. The function further allows to visualize trade-offs between cutoff-level and alpha diversity and enables direct comparison of SRS and repeated rarefying.
See Beule & Karlovsky (2020) <doi:10.7717/peerj.9593> for details regarding SRS.
SRScurve(data, metric = "richness", step = 50, sample = 0, max.sample.size = 0, rarefy.comparison = FALSE, rarefy.repeats = 10, rarefy.comparison.legend = FALSE, xlab = "sample size", ylab = "richness", label = FALSE, col, lty, ...)
data |
Data frame (species count or OTU table) in which columns are samples and rows are the counts of species or OTUs. Only integers are accepted as data. |
metric |
Character, "richness" (using specnumber) for species richness or "shannon", "simpson" or "invsimpson" (using diversity) for common diversity indices. Default is "richness". |
step |
Numeric, specifying the step used to vary the sample size. Default is 50. |
sample |
Numeric, specifying the cutoff-level to visualize trade-offs between cutoff-level and alpha diversity. |
max.sample.size |
Numeric, specifying the maximum sample size to which SRS curves are drawn. Default is 0 which does not limit the maximum sample size. |
rarefy.comparison |
Logical, if TRUE, median values of rarefy with n repeats (specified by rarefy.repeats) will be drawn for comparison. Default is FALSE. |
rarefy.repeats |
Numeric, specifying the number of repeats used to obtain median values for rarefying. Default is 10. |
rarefy.comparison.legend |
Logical, if TRUE, a legend for the comparison between SRS and rarefy is plotted. Default is FALSE. |
xlab, ylab, label, col, lty, ... |
Graphical parameters. |
See Beule & Karlovsky (2020) <doi:10.7717/peerj.9593> for details regarding scaling with ranked subsampling.
Returns a line plot visualizing the change in alpha diversity indices with changing sample size.
Vitor Heidrich, Petr Karlovsky, Lukas Beule
Beule L, Karlovsky P. 2020. Improved normalization of species count data in ecology by scaling with ranked subsampling (SRS): application to microbial communities. PeerJ 8:e9593
<https://doi.org/10.7717/peerj.9593>
##Samples should be arranged columnwise. ##Input data should not contain any categorial ##data such as taxonomic assignment or barcode sequences. ##An example of the input data can be found below: example_input_data <- matrix(c(sample(1:20, 100, replace = TRUE), sample(1:30, 100, replace = TRUE),sample(1:40, 100, replace = TRUE)), nrow = 100) colnames(example_input_data) <- c("sample_1","sample_2","sample_3") example_input_data <- as.data.frame(example_input_data) example_input_data ##Default settings of SRScurve. SRScurve(example_input_data, metric = "richness", step = 50, ylab = "richness", col = c("#000000", "#E69F00", "#56B4E9")) ##Limit the compution of SRS curves to a sample size of 200. SRScurve(example_input_data, metric = "richness", step = 50, max.sample.size = 200, ylab = "richness", col = c("#000000", "#E69F00", "#56B4E9")) ##SRScurve with comparison of SRS (solid lines) and repeated rarefying (dashed lines). ##Different colors correspond to indiviual samples. Cuttoff-level set to 200. SRScurve(example_input_data, metric = "richness", step = 50, sample = 200, max.sample.size = 200, rarefy.comparison = TRUE, rarefy.repeats = 10, rarefy.comparison.legend = TRUE, ylab = "richness", col = c(rep(c("#000000", "#E69F00", "#56B4E9"),2)), lty = c(1,2))
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.