R/givenExpression.R

In segmentSeq: Methods for identifying small RNA loci from high-throughput sequencing data

Defines functions givenExpression

Documented in givenExpression

# modification on git from copied files
givenExpression <- function(cD) {
    nelik = exp(rowSums(log(1 - exp(cD@locLikelihoods)), na.rm = TRUE))
    cD@posteriors <- (cD@posteriors) + log(1 - nelik)
    cD@nullPosts <- matrix(log(nelik), ncol = 1)
    cD
}

Try the segmentSeq package in your browser