tests/testthat/test-lefserPlotClad.R

In waldronlab/lefser: R implementation of the LEfSE method for microbiome biomarker discovery

data("zeller14")
z14 <- zeller14[, zeller14$study_condition != "adenoma"]
tn <- get_terminal_nodes(rownames(z14))
z14tn <- z14[tn, ]
z14tn_ra <- relativeAb(z14tn)

res <- lefser(z14tn_ra, classCol = "study_condition")

z14_input <- rowNames2RowData(z14tn_ra)
resClades <- lefserClades(z14_input, classCol = "study_condition")

test_that("lefserClades works", {
    expect_s3_class(resClades, "lefser_df_clades")
})

test_that("lefserPlotClad works", {
    ggt <- lefserPlotClad(df = resClades)
    expect_s3_class(ggt, "ggtree")
})