tests/testthat/test-hyp_hmap.R

hyp_hmap_tests <- function(hyp_obj, return_obj=FALSE) {
    expect_silent(hyp_hmap(hyp_obj, top=5))
    expect_silent(hyp_hmap(hyp_obj, val="pval", top=10))
    expect_silent(hyp_hmap(hyp_obj, val="fdr", top=10))
    p <- hyp_hmap(hyp_obj,  top=10)
    expect_is(p, "visNetwork")
    expect_is(p, "htmlwidget")
    p <- hyp_hmap(hyp_obj,  top=0)
    expect_true(is.null(p))
    if (return_obj) return(hyp_obj)
}

test_that("hyp_hmap() is working", {
    
    testdat <- readRDS(file.path(system.file("extdata", package="hypeR"), "testdat.rds"))
    rgsets_obj <- testdat$rgsets
    
    # Overrepresentation (signature)
    signature <- testdat$signature
    experiment <- testdat$experiment
    
    hypeR(signature, rgsets_obj, test="hypergeometric", background=100) %>%
    hyp_hmap_tests()
    p <- hypeR(experiment, rgsets_obj, test="hypergeometric", background=100) %>%
         hyp_hmap()
    expect_equal(length(p), 3)
    expect_equal(names(p), c("Signature 1", "Signature 2", "Signature 3"))
    expect_is(p[["Signature 3"]], "visNetwork")
    
    # Enrichment (ranked signature)
    signature <- names(testdat$weighted_signature)
    experiment <- lapply(testdat$weighted_experiment, names)

    hypeR(signature, rgsets_obj, test="kstest") %>%
    hyp_hmap_tests()
    p <- hypeR(experiment, rgsets_obj, test="kstest") %>%
         hyp_hmap()
    expect_equal(length(p), 3)
    expect_equal(names(p), c("Signature 1", "Signature 2", "Signature 3"))
    expect_is(p[["Signature 3"]], "visNetwork")
        
    # Enrichment (weighted signature)
    signature <- testdat$weighted_signature
    experiment <- testdat$weighted_experiment
    
    hypeR(signature, rgsets_obj, test="kstest") %>%
    hyp_hmap_tests()
    p <- hypeR(experiment, rgsets_obj, test="kstest") %>%
         hyp_hmap()
    expect_equal(length(p), 3)
    expect_equal(names(p), c("Signature 1", "Signature 2", "Signature 3"))
    expect_is(p[["Signature 3"]], "visNetwork")
})
montilab/hypeR documentation built on Oct. 29, 2023, 12:01 p.m.