inst/unitTests/test_convertToTxFeatures.R

In SGSeq: Splice event prediction and quantification from RNA-seq data

test_convertToTxFeatures <- function()
{
  
    gr <- GRanges(1, IRanges(c(1, 201), c(100, 300)), "+")
    grl <- split(gr, 1)
    current <- convertToTxFeatures(grl)
    target <- TxFeatures(GRanges(1, IRanges(c(1, 100, 201), c(100, 201, 300)),
      "+", type = c("F", "J", "L")), txName = CharacterList(1, 1, 1))
    checkIdentical(target, current)
    
}