getSamplesOf function

Get posterior samples of a parameter

Extracts posterior samples of individual parameters contained in the output slot of the readsData object returned by denoiSeq.

getSamplesOf(RDobject, parm, steps, condition = "A")

Arguments

• RDobject: A readsData object with a filled output slot.
• parm: A parameter name string i.e p, f or gene name.
• steps: An integer representing number of iterations used while calling denoiseq.
• condition: A character (either A or B) representing the two experimental conditions.

Returns

A vector of parameter samples, of length equal to steps.

Examples

#pre-filtering to remove lowly expressed genes
ERCC <- ERCC[rowSums(ERCC)>0, ]
RD <- new('readsData', counts = ERCC)
steps <- 30
#30 steps are just for illustration here. Atleast 5000 steps are adequate.
BI <- denoiseq(RD, steps)
samples <- getSamplesOf(BI, 'ERCC-00051', steps)
plot(samples, type='l', main = 'History plot of ERCC-00051')