makeTravisCoordsFromTxDb: Make Travis Coordinates from TxDb Object

Usage Arguments Examples

Usage

1
makeTravisCoordsFromTxDb(txdb, maximalAmbiguity = 3, minimalComponentLength = 100, minimalNcRNALength = 300, noBins = 100)

Arguments

txdb
maximalAmbiguity
minimalComponentLength
minimalNcRNALength
noBins

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
##---- Should be DIRECTLY executable !! ----
##-- ==>  Define data, use random,
##--	or do  help(data=index)  for the standard data sets.

## The function is currently defined as
function (txdb, maximalAmbiguity = 3, minimalComponentLength = 100, 
    minimalNcRNALength = 300, noBins = 100) 
{
    require("GenomicFeatures")
    parameter = list()
    parameter$txdb <- txdb
    parameter$maximalAmbiguity <- maximalAmbiguity
    parameter$minimalComponentLength <- minimalComponentLength
    parameter$minimalNcRNALength <- minimalNcRNALength
    parameter$noBins <- noBins
    component <- .extractComponent(parameter)
    print("Building Travis Coordinates. It may take a few minutes ...")
    utr3_bin <- makeTravisCoordsFromGRangesList(component[["utr3"]], 
        parameter$noBins)
    utr5_bin <- makeTravisCoordsFromGRangesList(component[["utr5"]], 
        parameter$noBins)
    cds_bin <- makeTravisCoordsFromGRangesList(component[["cds"]], 
        parameter$noBins)
    ncRNA_bin <- makeTravisCoordsFromGRangesList(component[["ncRNA"]], 
        parameter$noBins)
    print("Travis Coordinates Built ...")
    mcols(utr3_bin) <- data.frame(mcols(utr3_bin), comp = "UTR3", 
        category = "mRNA")
    mcols(utr3_bin)$pos <- mcols(utr3_bin)$pos/3 + 2/3
    mcols(utr5_bin) <- data.frame(mcols(utr5_bin), comp = "UTR5", 
        category = "mRNA")
    mcols(utr5_bin)$pos <- mcols(utr5_bin)$pos/3
    mcols(cds_bin) <- data.frame(mcols(cds_bin), comp = "CDS", 
        category = "mRNA")
    mcols(cds_bin)$pos <- mcols(cds_bin)$pos/3 + 1/3
    mcols(ncRNA_bin) <- data.frame(mcols(ncRNA_bin), comp = "lncRNA", 
        category = "lncRNA")
    TravisCoords <- suppressWarnings(c(utr5_bin, cds_bin, utr3_bin, 
        ncRNA_bin))
    return(TravisCoords)
  }

lzcyzm/Travis documentation built on May 21, 2019, 9:16 a.m.