An R package for Differential Expression Analysis. Given count data from two experimental conditions, denoiSeq helps one determine which transcripts are differentially expressed across the two conditions using Bayesian inference of the parameters of a bottom-up model for PCR amplification developed in "Chromatin conformation governs T cell receptor J beta gene segment usage", by Ndifon et al.
To use the package, one needs to create a readsData
object and invoke the denoiseq
function on it. The results are obtained from the return value of denoiseq using the results
function which then computes the test statistic used in differential analysis.
RD <- new("readsData", counts = ERCC) #creating the readsData object
steps <- 3000 #steps for MCMC
BI <- denoiseq(RD, steps) #invoking denoiseq on the readsData object
rez <- results(BI,steps) #computing the test statistic
This package can be istalled from CRAN using install.packages("denoiSeq") or from github using devtools::install_github("buriom/denoiSeq").