tests/testthat/test-ExpressionSet.R

In StoreyLab/biobroom: Turn Bioconductor objects into tidy data frames

context("ExpressionSet")

test_that("ExpressionSet tidier works as expected", {
    dds <- makeExampleDataSet("ExpressionSet")
    td <- tidy(dds)
    tdp <- tidy(dds, addPheno=TRUE)

    ## Check that both objects have the bare minimum columns with expected
    ## types
    req.cols <- c(gene='character', sample='character', value='numeric')
    for (rc in names(req.cols)) {
        expect_is(td[[rc]], req.cols[[rc]], info=sprintf('td (%s)', rc))
        expect_is(tdp[[rc]], req.cols[[rc]], info=sprintf('tdp (%s)', rc))
    }

    ## tdp has the extra sample data
    pd <- pData(dds)
    pd.cols <- setNames(sapply(pd, class), names(pd))
    for (pc in names(pd.cols)) {
        expect_is(tdp[[pc]], pd.cols[[pc]], info=sprintf("tdp pheno (%s)", pc))
    }
})