R/normalizeFragmentData.r

Defines functions .normalizeFragmentData

.normalizeFragmentData <- function(expData) {

    readCount = length(rawReads(expData))
    relevantFragData = nearCisFragments(expData)

    # simple RPM normalization (reads per million)
    RPMFrags <- relevantFragData
    RPMFrags$leftFragEndReads <- RPMFrags$leftFragEndReads / readCount * 1000000
    RPMFrags$rightFragEndReads <- RPMFrags$rightFragEndReads / readCount * 1000000
    RPMFrags$fragEndReadsAverage <- RPMFrags$fragEndReadsAverage / readCount * 1000000

    return(RPMFrags)
}




setMethod("normalizeFragmentData",
    signature=signature(expData="Data4Cseq"),
    .normalizeFragmentData)

Try the Basic4Cseq package in your browser

Any scripts or data that you put into this service are public.

Basic4Cseq documentation built on Nov. 8, 2020, 6:53 p.m.