DiffBind Count reads from paired end bamfiles
1
0
Entering edit mode
4.1 years ago
g.persic1991 ▴ 10

Hi guys... I have a (URGENT !!!!) question about how to use dba.count when I have a Paired-End bam files...

For Single-end data processed up to now, after create a consensus peak,

dbob.consensus <- dba.peakset(dbObj, consensus = DBA_CONDITION, minOverlap=0.5)
consensus <- dba.peakset(dbob.consensus, bRetrieve=TRUE,
                             peaks=dbob.consensus$masks$Consensus, minOverlap=1)

I used dba.count with below command

dbObj_01_TMM <- dba.count(dbObj, peaks=consensus, bUseSummarizeOverlaps=TRUE, minOverlap=1)

Now my question is, can I use these code also for BAM PE data???

If not, can you suggest me how to do???

thanks in advance Giuseppe

DiffBind dba.count paired end bamfile • 1.7k views
ADD COMMENT
1
Entering edit mode
4.1 years ago
Rory Stark ★ 2.1k

If you look at the help page for dba.count(), it shows some configuration options that are used when bUseSummarizedOverlaps=TRUE. Specifically, you can set

`myDBA$config$singleEnd <- FALSE`

before calling dba.count() to use paired-end mode.

Note that in the version being released this week (October 28), the default is no bUseSummarizedOverlaps=TRUE, and DiffBind will automatically detect single-end/paired-end mode.

ADD COMMENT

Login before adding your answer.

Traffic: 2560 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