tests/testthat/test-readStartCov.R

In RiboProfiling: Ribosome Profiling Data Analysis: from BAM to Data Representation and Interpretation

library(testthat)
library(RiboProfiling)

data(ctrlGAlignments)
aln <- ctrlGAlignments
alnGRanges <- readsToStartOrEnd(aln, what="start")

oneBinRanges <-
    aroundPromoter(
        TxDb.Hsapiens.UCSC.hg19.knownGene::TxDb.Hsapiens.UCSC.hg19.knownGene,
        alnGRanges
    )

listPromoterCov <- readStartCov(
    alnGRanges,
    oneBinRanges,
    matchSize="all",
    fixedInterval=c(-20, 20),
    renameChr="aroundTSS",
    charPerc="perc"
)

test_that("readStartCov output type",{
    expect_is(listPromoterCov,"list")
    expect_is(listPromoterCov[[1]],"GRanges")
}
)

test_that("readStartCov output dim",{
    expect_equal(length(listPromoterCov[[1]]),41)
}
)