Calculates mean read number estimations

Source: R/likelihoods.R

Calculates mean read number estimations

predictExpression(par, pulseData)

Arguments

par

estimated parameters from fitModel
pulseData

a PulseData object.

Value

a named list:

  • predictions, a matrix of the same dimension as of the raw counts

  • llog, a matrix with logarithms of likelihood for the given raw counts.

Examples

# NOT RUN {
# Plot expected values vs the raw counts.
# Let res is the return of fitModel and pulseData is a PulseData object
pr <- predictExpression(pulseData, res)
plot(y = pr$predictions, x = pulseData$counts, xlab = "raw", ylab = "fitted")
# }