tests/testthat/test_extend.R

# Create example granges
bsgenome <- BSgenome.Mmusculus.UCSC.mm10::BSgenome.Mmusculus.UCSC.mm10
granges <-  GenomicRanges::GRanges('chr1', 
                    '100-200',
                    strand   = '-',
                    seqinfo  = seqinfo(bsgenome))
chrlength <- GenomeInfoDb::seqlengths(bsgenome)[['chr1']]

# Test
context('extend')

test_that('extend works', {
    expect_equal(
        start(extend(granges, -5, 5)),  
        start(granges) - 5)
    expect_equal(  
        end(extend(granges, -5, 5)),
        end(granges) + 5)
})

test_that('extend warns for coordinates < 1', {
    expect_warning(extend(granges,  -500))
})

test_that('extend warns for coordinates > chrlength', {
    expect_warning(extend(granges, 1, 1 + chrlength, plot = FALSE))
})

Try the multicrispr package in your browser

Any scripts or data that you put into this service are public.

multicrispr documentation built on Nov. 8, 2020, 5:10 p.m.