1 | makeTravisCoordsFromTxDb(txdb, maximalAmbiguity = 3, minimalComponentLength = 100, minimalNcRNALength = 300, noBins = 100)
|
txdb |
|
maximalAmbiguity |
|
minimalComponentLength |
|
minimalNcRNALength |
|
noBins |
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)
}
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.