tests/testthat/test_resizeTranscripts.R

In fursham-h/ponder:

context("Test resizeTranscript function")

# setup GRanges
gr1 <- GRanges(seqnames = "chr1", strand = c("+", "+", "+"),
              ranges = IRanges(start = c(1,500,1000), 
                               end = c(100,600,1100)))
gr2 <- GRanges(seqnames = "chr1", strand = c("-", "-", "-"),
               ranges = IRanges(start = c(1,500,1000), 
                                end = c(100,600,1100)))

out1 = resizeGRangesTranscripts(gr1, 20,80)
out2 = resizeGRangesTranscripts(gr2, 20,80)
out3 = resizeGRangesTranscripts(gr1, 110,150)
out4 = resizeGRangesTranscripts(gr2, 110,150)
test_that("Test sample output", {

  expect_equal(start(out1), c(21,500,1000))
  expect_equal(end(out1), c(100,600,1020))
  
  expect_equal(start(out2), c(1000,500,81))
  expect_equal(end(out2), c(1080,600,100))
  
  expect_equal(start(out3), 510)
  expect_equal(end(out3), 551)
  
  expect_equal(start(out4), 550)
  expect_equal(end(out4), 591)
})