scde.failure.probability {scde} | R Documentation |
Returns estimated drop-out probability for each cell (row of models
matrix), given either an expression magnitude
scde.failure.probability(models, magnitudes = NULL, counts = NULL)
models |
models determined by |
magnitudes |
a vector ( |
counts |
a vector ( |
a vector or a matrix of drop-out probabilities
data(es.mef.small) cd <- clean.counts(es.mef.small, min.lib.size=1000, min.reads = 1, min.detected = 1) data(o.ifm) # Load precomputed model. Use ?scde.error.models to see how o.ifm was generated o.prior <- scde.expression.prior(models = o.ifm, counts = cd, length.out = 400, show.plot = FALSE) # calculate probability of observing a drop out at a given set of magnitudes in different cells mags <- c(1.0, 1.5, 2.0) p <- scde.failure.probability(o.ifm, magnitudes = mags) # calculate probability of observing the dropout at a magnitude corresponding to the # number of reads actually observed in each cell self.p <- scde.failure.probability(o.ifm, counts = cd)