getSamplesOf {denoiSeq}

Get posterior samples of a parameter

Description

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

Usage

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.

Value

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')