maximum overlapping internal junction score
internal_junc_score( juncGR, scoreColname = "score", sampleColname = "sample_id", minScore = 0, verbose = FALSE, ... )
juncGR | GRanges containing splice junctions, with numeric
column |
---|---|
scoreColname, sampleColname | colnames in |
minScore | optional |
verbose | logical indicating whether to print verbose output. |
... | additional arguments are ignored. |
numeric vector named by names(juncGR)
whose values
are the maximum score of internal overlapping junction ends.
This function is intended for internal use, and calculates the maximum score for junction GRanges for the special case where:
a junction range overlaps the start or end of another junction
the start or end of the other junction is internal, that is it does not overlap the same start or end.
overlaps are constrained to the same "sample_id"
stored
in values(juncGR)$sample_id
.
The goal is to determine the highest score contained inside each junction region, so that the junction arc height can be defined higher in order to minimize junction arc overlaps.
Other jam RNA-seq functions:
assignGRLexonNames()
,
closestExonToJunctions()
,
combineGRcoverage()
,
defineDetectedTx()
,
detectedTxInfo()
,
exoncov2polygon()
,
flattenExonsBy()
,
getGRcoverageFromBw()
,
groups2contrasts()
,
makeTx2geneFromGtf()
,
make_ref2compressed()
,
prepareSashimi()
,
runDiffSplice()
,
sortSamples()
,
spliceGR2junctionDF()