Question: DiffBind: How could I debug estimateDispersionsFit (required by DBA_DESEQ2)
gravatar for zhaolin20042008
2.7 years ago by
U of Michigan
zhaolin200420080 wrote:


I got an error when I used DiffBind to treat a human ATAC-seq library data set. (see is the data set like DNase-seq, the enrichment of open chromatin region are variable. 

Error in estimateDispersionsFit(object, fitType = fitType, quiet = quiet) :
  all gene-wise dispersion estimates are within 2 orders of magnitude
from the minimum value, and so the standard curve fitting techniques will not work.
You can instead use the gene-wise estimates as final estimates:
dds <- estimateDispersionsGeneEst(dds)
dispersions(dds) <- mcols(dds)$dispGeneEst

I knew the reason but I failed to debug using trace () at DBA_DESEQ2 because it is not a function.Could anyone suggest a solution? Thanks a lot. 

diffbind atac-seq • 639 views
ADD COMMENTlink written 2.7 years ago by zhaolin200420080


If you want to debug estimateDispersionsFit you can debug it directly:

> debug(DESeq2::estimateDispersionsFit)

> x = dba.analyze(...)

Then use the debug commands to step through it, when R stops in that function.

ADD REPLYlink written 2.7 years ago by Gord Brown590
Please log in to add an answer.


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