tests/testthat/test-GeneToSymbol.R

In acidgenomics/r-acidgenomes: Acid Genomics Genome Annotations

test_that("GeneToSymbol", {
    formats <- eval(formals(`GeneToSymbol,DFrame`)[["format"]])
    for (format in formats) {
        object <- gr
        object <- GeneToSymbol(object, format = format)
        expect_s4_class(object, "GeneToSymbol")
        expect_identical(colnames(object), c("geneId", "geneName"))
    }
})

test_that("format: makeUnique", {
    ## Ensure duplicate gene names (symbols) are sanitized using `make.names`,
    ## and that NA gene names are converted to "unannotated".
    object <- GeneToSymbol(
        object = DataFrame(
            "geneId" = c(
                "gene1",
                "gene2",
                "gene3",
                "gene4"
            ),
            "geneName" = c(
                "symbol1",
                "symbol1",
                "symbol2",
                NA_character_
            )
        ),
        format = "makeUnique"
    )
    expected <- DataFrame(
        "geneId" = c(
            "gene1",
            "gene2",
            "gene3",
            "gene4"
        ),
        "geneName" = c(
            "symbol1",
            "symbol1.1",
            "symbol2",
            "unannotated"
        )
    )
    metadata(expected) <- list(
        "date" = Sys.Date(),
        "dupes" = "symbol1",
        "format" = "makeUnique",
        "packageVersion" = .pkgVersion
    )
    expected <- new(Class = "GeneToSymbol", expected)
    expect_identical(object, expected)
    ## Don't allow "geneId" column to contain NA values.
    expect_error(
        GeneToSymbol(
            object = DataFrame(
                "geneId" = c(
                    "gene1",
                    NA_character_
                ),
                "geneName" = c(
                    "symbol1",
                    "symbol2"
                )
            ),
            format = "makeUnique"
        )
    )
    ## Ensure gene identifiers return sorted.
    object <- GeneToSymbol(
        object = DataFrame(
            "geneId" = c(
                "gene2",
                "gene1"
            ),
            "geneName" = c(
                "symbol2",
                "symbol1"
            ),
            row.names = c("B", "A")
        ),
        format = "makeUnique"
    )
    expected <- DataFrame(
        "geneId" = c("gene1", "gene2"),
        "geneName" = c("symbol1", "symbol2"),
        row.names = c("A", "B")
    )
    metadata(expected) <- list(
        "date" = Sys.Date(),
        "format" = "makeUnique",
        "packageVersion" = .pkgVersion
    )
    expected <- new(Class = "GeneToSymbol", expected)
    expect_identical(object, expected)
})

test_that("format: 1:1", {
    object <- GeneToSymbol(
        object = DataFrame(
            "geneId" = c(
                "gene2",
                "gene1",
                "gene4",
                "gene3",
                "gene5",
                NA_character_
            ),
            "geneName" = c(
                "symbol1",
                "symbol1",
                "symbol2",
                "symbol2",
                NA_character_,
                "symbol3"
            ),
            row.names = LETTERS[seq_len(6L)]
        ),
        format = "1:1"
    )
    expected <- DataFrame(
        "geneId" = c(
            "gene1",
            "gene3"
        ),
        "geneName" = c(
            "symbol1",
            "symbol2"
        ),
        row.names = c("B", "D")
    )
    metadata(expected) <- list(
        "date" = Sys.Date(),
        "dropped" = c("E" = 5L, "F" = 6L),
        "format" = "1:1",
        "packageVersion" = .pkgVersion
    )
    expected <- new(Class = "GeneToSymbol", expected)
    expect_identical(object, expected)
})

test_that("format: unmodified", {
    object <- GeneToSymbol(
        object = DataFrame(
            "geneId" = c(
                "gene1",
                "gene2",
                "gene3",
                "gene4",
                NA_character_
            ),
            "geneName" = c(
                "symbol1",
                "symbol1",
                "symbol2",
                NA_character_,
                "symbol3"
            ),
            row.names = LETTERS[seq_len(5L)]
        ),
        format = "unmodified"
    )
    expected <- DataFrame(
        "geneId" = c(
            "gene1",
            "gene2",
            "gene3"
        ),
        "geneName" = c(
            "symbol1",
            "symbol1",
            "symbol2"
        ),
        row.names = LETTERS[seq_len(3L)]
    )
    metadata(expected) <- list(
        "date" = Sys.Date(),
        "dropped" = c("D" = 4L, "E" = 5L),
        "format" = "unmodified",
        "packageVersion" = .pkgVersion
    )
    expected <- new(Class = "GeneToSymbol", expected)
    expect_identical(object, expected)
})

test_that("summary method", {
    object <- gr
    x <- GeneToSymbol(object)
    output <- capture.output(summary(x))
    expect_identical(
        head(output, n = 1L),
        paste("genes:", length(object))
    )
})