Question: R's DiffBind. Error in pv.DBAplotVolcano(DBA, contrast = contrast, method = method,: object 'sigSites' not found
gravatar for msimmer92
11 months ago by
msimmer92220 wrote:

Hello all, I am using R´s DiffBind to find differentially bound sites on my ChIP-seq data. I have ¨healthy¨ condition, and two other ¨unhealthy¨ conditions that I want to compare. When I compare one unhealthy to the healthy, it all works fine with my script. I get nice volcano and MA plots (which is what I want, in this case):

results=dba.count(results, minOverlap=1)
counts=dba.peakset(results, bRetrieve=T, DataType=DBA_DATA_FRAME)
results=dba.contrast(results, results$masks$G, results$masks$M, unique(Condition)[2], unique(Condition)[3], 
minMembers=min(c(table(Condition))), categories=DBA_CONDITION) 
results=dba.analyze(results, method=DBA_DESEQ2)

When I compare the two unhealthy ones, I encounter the following error:

> dba.plotVolcano(results)
Error in pv.DBAplotVolcano(DBA, contrast = contrast, method = method,  : 
  object 'sigSites' not found
In addition: Warning message:
No sites above threshold 

Does someone know what can be wrong in this case? I couldn´t find any thread about this DiffBind problem. Thank you in advance!

Computer and R Specifications:

platform x86_64-apple-darwin15.6.0
arch x86_64
os darwin15.6.0
system x86_64, darwin15.6.0
major 3
minor 5.2
year 2018
month 12
day 20
svn rev 75870
language R
version.string R version 3.5.2 (2018-12-20) nickname Eggshell Igloo

chip-seq diffbind R • 491 views
ADD COMMENTlink modified 11 months ago by Friederike5.2k • written 11 months ago by msimmer92220
gravatar for Friederike
11 months ago by
United States
Friederike5.2k wrote:

It seems that there simply aren't any sites that are different between your two unhealthy samples.

You can check this, e.g.:


You can probably still get a volcano plot by changing the threshold, i.e. dba.plotVolcano(results, UsePval=1)

ADD COMMENTlink written 11 months ago by Friederike5.2k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 585 users visited in the last hour