RPKM on TSS using DiffBind
Entering edit mode
3 months ago

Hi everyone, I want to extrapolate RPKM value from Diffbind. Investigated regions are TSS (± 2.5kb) but around 3500 TSS are "lost" because merged. I have read that in the new version of DiffBind is possible to use dbObj_region$config$mergeOverlap with negative value to avoid merge but I'm not able to use.

Below my code

region <- read.csv("DiffBind_template_peaks.csv")
DBA <- dba(sampleSheet=region)  
DBA$config$mergeOverlap = -1
DBA_Count <- dba.count(DBA, score=DBA_SCORE_RPKM)  
rpkm <- dba.peakset(DBA_Count, bRetrieve=TRUE)

Thanks a lot in advance

ciao Giuseppe

RPKM TSS DiffBind • 206 views
Entering edit mode
6 weeks ago
Rory Stark ★ 1.2k

The value of DBA$config$mergeOverlap is the number of overlapping bases required for regions to be merged. So if you are starting with regions that overlap and you don't want them to be merged, you should use a larger positive value. For example if you set DBA$config$mergeOverlap <- 100, then regions which overlap by 100 or fewer bases will not be merged.

Negative values for DBA$config$mergeOverlap are used to merge non-overlapping regions that are close to each other. So if you set DBA$config$mergeOverlap <- -100, then regions which are within 100 bases of each other will be merged even if they don't overlap.


Login before adding your answer.

Traffic: 2122 users visited in the last hour
Help About
Access RSS

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

Powered by the version 2.3.6