context("poolCoverage")
normal.coverage.file <- system.file("extdata", "example_normal.txt.gz",
package = "PureCN")
normal2.coverage.file <- system.file("extdata", "example_normal2.txt.gz",
package = "PureCN")
normal.coverage.files <- c(normal.coverage.file, normal2.coverage.file)
test_that("Example coverage is averaged", {
coverage <- lapply(normal.coverage.files, readCoverageFile)
pool <- poolCoverage(coverage)
expect_equal(coverage[[1]]$average.coverage + coverage[[2]]$average.coverage,
pool$average.coverage)
expect_equal(coverage[[1]]$coverage + coverage[[2]]$coverage,
pool$coverage)
pool2 <- poolCoverage(coverage, w = c(0.5, 0.5))
expect_equal((coverage[[1]]$coverage + coverage[[2]]$coverage) / 2,
pool2$coverage)
})
test_that("Exceptions happend with wrong input", {
coverage <- lapply(normal.coverage.files, readCoverageFile)
expect_error(poolCoverage(coverage, w = seq(3)), "different lengths")
})
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.