RPKM on TSS using DiffBind

RPKM on TSS using DiffBind


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






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.

before adding your answer.

Traffic: 988 users visited in the last hour

Read more here: Source link