getSamplesOf: Get posterior samples of a parameter

In denoiSeq: Differential Expression Analysis Using a Bottom-Up Model

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

denoiSeq documentation built on May 2, 2019, 9:33 a.m.