orfScore {ORFik}R Documentation

Get ORFscore for a GRangesList of ORFs

Description

ORFscore tries to check whether the first frame of the 3 possible frames in an ORF has more reads than second and third frame.

Usage

orfScore(grl, RFP, is.sorted = FALSE)

Arguments

grl

a GRangesList object with ORFs

RFP

ribozomal footprints, given as Galignment object, Granges or GRangesList

is.sorted

logical (F), is grl sorted.

Details

Pseudocode: assume rff - is reads fraction in specific frame

ORFScore = log(rrf1 + rrf2 + rrf3)

For all ORFs where rrf2 or rrf3 is bigger than rff1, negate the resulting value.

ORFScore[rrf1Smaller] <- ORFScore[rrf1Smaller] * -1

As result there is one value per ORF: Positive values say that the first frame have the most reads, negative values say that the first frame does not have the most reads.

Value

a matrix with 4 columns, the orfscore and score of each of the 3 tiles

References

doi: 10.1002/embj.201488411

See Also

Other features: computeFeaturesCage, computeFeatures, disengagementScore, distToCds, distToTSS, entropy, floss, fpkm_calc, fpkm, fractionLength, initiationScore, insideOutsideORF, isInFrame, isOverlapping, kozakSequenceScore, rankOrder, ribosomeReleaseScore, ribosomeStallingScore, subsetCoverage, translationalEff

Examples

ORF <- GRanges(seqnames = "1",
               ranges = IRanges(start = c(1, 10, 20), end = c(5, 15, 25)),
               strand = "+")
names(ORF) <- c("tx1", "tx1", "tx1")
grl <- GRangesList(tx1_1 = ORF)
RFP <- GRanges("1", IRanges(25, 25), "+") # 1 width position based
score(RFP) <- 28 # original width
orfScore(grl, RFP) # negative because more hits on frames 1,2 than 0.

# example with positive result, more hits on frame 0 (in frame of ORF)
RFP <- GRanges("1", IRanges(c(1, 1, 1, 25), width = 1), "+")
score(RFP) <- c(28, 29, 31, 28) # original width
orfScore(grl, RFP)


[Package ORFik version 1.2.1 Index]