tests/testthat/test_read_functions.R

In WeiSong-bio/roryk-bcbioSinglecell: Single-Cell RNA-Seq Utilities

context("Read Functions")



# readCellRanger ===============================================================
test_that("readCellRanger", {
    x <- suppressWarnings(readCellRanger(
        uploadDir = system.file(
            "extdata/cellranger",
            package = "bcbioSingleCell"
        )
    ))
    expect_is(x, "SingleCellExperiment")
    expect_identical(dim(x), c(500L, 500L))
    expect_identical(
        sampleNames(x),
        c(cellranger_1 = "cellranger")
    )
    expect_identical(
        colnames(x)[[1L]],
        "cellranger_1_AAACCTGAGAAGGCCT"
    )
    expect_identical(
        rownames(x)[[1L]],
        "ENSG00000243485"
    )
    expect_identical(
        lapply(metadata(x), class),
        list(
            version = c("package_version", "numeric_version"),
            pipeline = "character",
            level = "character",
            uploadDir = "character",
            sampleDirs = "character",
            sampleMetadataFile = "character",
            sampleData = "data.frame",
            interestingGroups = "character",
            cell2sample = "factor",
            organism = "character",
            genomeBuild = "character",
            ensemblRelease = "integer",
            rowRangesMetadata = c("tbl_df", "tbl", "data.frame"),
            umiType = "character",
            allSamples = "logical",
            call = "call",
            date = "Date",
            wd = "character",
            utilsSessionInfo = "sessionInfo",
            devtoolsSessionInfo = "session_info"
        )
    )
})



# readCellTypeMarkers ==========================================================
test_that("readCellTypeMarkers : Mus musculus", {
    file <- system.file(
        "extdata/cell_type_markers.csv",
        package = "bcbioSingleCell"
    )
    gene2symbol <- makeGene2symbolFromEnsembl("Homo sapiens")
    x <- readCellTypeMarkers(
        file = file,
        gene2symbol = gene2symbol
    )
    group <- dplyr::group_vars(x)
    expect_identical(group, "cellType")
    x <- x[1L, ]
    y <- tibble::tibble(
        "cellType" = "B Cell",
        "geneID" = "ENSG00000177455",
        "geneName" = "CD19"
    ) %>%
        dplyr::group_by(!!rlang::sym("cellType"))
    expect_identical(x, y)
})