data("GlobalPatterns", package = "phyloseq")
data("enterotype", package = "phyloseq")
GP <- GlobalPatterns
data_fungi_2trees <- subset_samples(data_fungi, data_fungi@sam_data$Tree_name %in% c("A10-005", "AD30-abm-X"))
GP_archae <- subset_taxa(GlobalPatterns, GlobalPatterns@tax_table[, 1] == "Archaea")
data_basidio <- subset_taxa(data_fungi, Phylum == "Basidiomycota")
test_that("tsne_pq works with data_fungi_mini dataset", {
skip_on_os("windows")
skip_on_cran()
expect_silent(suppressMessages(res_tsne <- tsne_pq(data_fungi_mini)))
expect_s3_class(res_tsne, "Rtsne")
expect_silent(suppressMessages(res_tsne <- tsne_pq(data_fungi_mini, dims = 3, perplexity = 25)))
})
test_that("plot_tsne_pq works with data_fungi_mini dataset", {
skip_on_os("windows")
skip_on_cran()
expect_silent(suppressMessages(pt <- plot_tsne_pq(data_fungi_mini, fact = "Height", perplexity = 15)))
expect_s3_class(pt, "ggplot")
expect_error(plot_tsne_pq(data_fungi_mini, fact = "HEIgTHT"))
})
test_that("SRS_curve_pq works with data_fungi_mini dataset", {
skip_on_cran()
expect_silent(suppressMessages(sc <- SRS_curve_pq(data_fungi_mini)))
expect_silent(suppressMessages(sc <- SRS_curve_pq(data_fungi_mini, clean_pq = TRUE)))
expect_s3_class(sc, "recordedplot")
expect_silent(suppressMessages(sc <- SRS_curve_pq(data_fungi_mini, metric = "shannon")))
expect_silent(suppressMessages(sc <- SRS_curve_pq(data_fungi_mini, step = 20, rarefy.repeats = 15)))
})
test_that("multiplot works fine", {
skip_on_cran()
res_venn1 <- ggvenn_pq(data_fungi_mini, "Height")
res_venn2 <- ggvenn_pq(data_fungi_mini, "Time")
expect_silent(multiplot(res_venn1, res_venn2))
expect_message(multiplot(res_venn1))
expect_type(multiplot(res_venn1, res_venn2), "NULL")
})
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.