Description Usage Arguments Value Author(s) Examples
use DEseq2 to test differential peaks
1 |
counts |
output of countTable |
... |
parameters could be passed to DESeqDataSet |
an object of GRanges
Jianhong Ou
1 2 3 4 5 6 7 | path <- system.file("extdata", package = "diffPeaks", mustWork = TRUE)
bamfiles <- dir(path, "bam$")
peaks <- dir(path, "bed$")
p <- mergePeaks(file.path(path, peaks))
colData <- DataFrame(samples=bamfiles, condition=sub(".rep..bam", "", bamfiles))
cnt <- countTable(p, file.path(path, bamfiles), colData)
DEpeaks(cnt, design= ~condition)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.