Using DiffBind with only a counts table?
1
1
Entering edit mode
pgontarz ▴ 10
@pgontarz-20229
Last seen 2.5 years ago

Hello,

I am interested in using diffbind to perform differential expression analysis on some ATAC data. However, all I have is counts table for peaks (no beds or bams) and a meta table with the sample names and conditions they belong to. Is there a way to use diffbind from here and is there any benefit over simply using DESeq2 or edgeR?

Thanks for the help in advance!

Fake example table

peakname condition1r1 condition1r2 condition2r1 condition2r2

peak1 19 21 29 31

peak2 34 30 15 17

diffbind ATACseq ATAC counts differential analysis • 423 views
ADD COMMENT
1
Entering edit mode
Rory Stark ★ 4.1k
@rory-stark-5741
Last seen 3 days ago
CRUK, Cambridge, UK

You can import count data directly into DiffBind using the Counts column in the sample sheet. This is documented in the help page for dba.peakset(). It assume you have already generated a consensus peakset (same intervals for each sample) and counted the reads in each interval, so that you have a file or vector of counts for each sample.

The advantages of using DiffBind for this are mostly that you can use its reporting and plotting functions. The main limitation is that your experimental design need be pretty straightforward, otherwise you are better off using DESeq2 or edgeR.

ADD COMMENT

Login before adding your answer.

Traffic: 231 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6