tests/testthat/test-local.R

In selinj/chromswitch: An R package to detect chromatin state switches from epigenomic data

context("Local processing functions")

test_that("retrievePeaks finds peaks in the query region", {
    skip_on_os("windows")

    metadata <- data.frame(Sample = c("A", "B"),
                           TestMark = c("../test_data/basic.bed",
                                        "../test_data/basic.bed"),
                           stringsAsFactors = FALSE)

    pks <- list(A = GRanges(seqnames = rep("chr1", 3),
                            ranges = IRanges(start = c(100, 150, 500),
                                             end = c(200, 250, 600))),
                B = GRanges(seqnames = rep("chr1", 3),
                            ranges = IRanges(start = c(100, 150, 500),
                                             end = c(200, 250, 600))))

    pks_in_region <- list(A = GRanges(seqnames = rep("chr1", 2),
                                      ranges = IRanges(start = c(100, 150),
                                                       end = c(200, 250))),
                          B = GRanges(seqnames = rep("chr1", 2),
                                      ranges = IRanges(start = c(100, 150),
                                                       end = c(200, 250))))

    pks_in_region$A <- `elementMetadata<-`(x = pks_in_region$A,
                        value = data.frame(length = as.integer(c(101, 101))))
    pks_in_region$B <- `elementMetadata<-`(x = pks_in_region$B,
                        value = data.frame(length = as.integer(c(101, 101))))

    region <- GRanges(seqnames = "chr1",
                      ranges = IRanges(start = 100, end = 300))

    expect_equal(retrievePeaks(pks, metadata, region),
                 LocalPeaks(region, pks_in_region, c("A", "B")))

    region2 <- GRanges(seqnames = "chr2",
                      ranges = IRanges(start = 100, end = 300))

    # No peaks
    expect_warning(retrievePeaks(pks, metadata, region2))

})


test_that("reducePeaks reduces nearby peaks within a gap", {
    skip_on_os("windows")

    pks <- list(A = GRanges(seqnames = rep("chr1", 3),
                            ranges = IRanges(start = c(100, 210, 500),
                                             end = c(200, 250, 600))),
                B = GRanges(seqnames = rep("chr1", 3),
                            ranges = IRanges(start = c(100, 210, 500),
                                             end = c(200, 250, 600))))

    region <- GRanges(seqnames = "chr1",
                      ranges = IRanges(start = 100, end = 1000))

    lpk <- LocalPeaks(region, pks, c("A", "B"))

    pks_red <- list(A = GRanges(seqnames = rep("chr1", 2),
                            ranges = IRanges(start = c(100, 500),
                                             end = c(250, 600))),
                B = GRanges(seqnames = rep("chr1", 2),
                            ranges = IRanges(start = c(100, 500),
                                             end = c(250, 600))))

    expect_equal(reducePeaks(lpk, 80), LocalPeaks(region, pks_red, c("A", "B")))
    expect_equal(reducePeaks(lpk, 5), lpk)
    expect_error(reducePeaks(lpk, -100), "must be a positive integer")
    expect_error(reducePeaks(lpk, 0), "must be a positive integer")

})