This is a test version of Biostars. For the public version, visit https://www.biostars.org.
dba.count reduces binding matrix, why?
db <- dba(samplesheet)

This gives me:

db
9 Samples, 138044 sites in matrix (193226 total)

then I do db.cnt <- dba.count(db)

db.cnt
9 Samples, 55893 sites in matrix

I wonder why dba.count is removing more than half of the intervals?

diffbind3 bioconductor diffbind atac-seq chip-seq
• 1,934 views • link
updated by Ram Ram • written by halffedelf • halffedelf

I figured filter and filterFun was causing removal of such large number of intervals. I had the score parameter set to DBA_SCORE_RPKM. I am wondering, in this case, if "filter" was applied using DBA_SCORE_RPKM or raw read counts? Because by default, filter = 5 and filterFun = max, so does it mean peaks are filtered for 5 read counts, or 5 RPKM (since my score=DBA_SCORE_RPKM), in at least 1 sample?

• 0 views • link
written by halffedelf • halffedelf

I removed score=DBA_SCORE_RPKM and now it retains 98% of the intervals. I think it is a very important thing to mention in the documentation. Filter = 5 will filter very differently between raw reads and RPKM. RPKM>5 is a very very stringent filtering, whereas raw reads>5 is not so much.

• 0 views • link
written by halffedelf • halffedelf

1 answer

What version is this?

The current default should be max RPKM of 1.

• 0 views • link
written by Rory Stark Rory Stark

Log in to answer this question.

More posts like this