Question: DiffBind: How could I debug estimateDispersionsFit (required by DBA_DESEQ2)
gravatar for zhaolin20042008
17 months 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. 

ADD COMMENTlink written 17 months 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 17 months ago by Gord Brown570
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.2.0
Traffic: 279 users visited in the last hour