RPKM on TSS using DiffBind
1
0
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

library(DiffBind)
DBA <- dba(sampleSheet=region)
DBA$config$mergeOverlap = -1
DBA_Count <- dba.count(DBA, score=DBA_SCORE_RPKM)
rpkm <- dba.peakset(DBA_Count, bRetrieve=TRUE)


ciao Giuseppe

RPKM TSS DiffBind • 206 views
0
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.