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
library(DiffBind)
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
• 104 views
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.
Traffic: 988 users visited in the last hour
Read more here: Source link