ATACseqQC icon indicating copy to clipboard operation
ATACseqQC copied to clipboard

estimateLibComplexity error

Open noorisotoudeh opened this issue 4 years ago • 0 comments

Hi ATACseqQC team, thanks for your nice package for atac-seq quality control. I get an error while i am doing estimateLibComplexity(readsDupFreq("test.bam")) it gives me the following error:


Error in f(t * t.scale): could not find function "f" Traceback:

  1. estimateLibComplexity(readsDupFreq(poorR1))
  2. suppressWarnings({ . result = ds.rSAC.bootstrap(histFile, r = 1, times = times) . })
  3. withCallingHandlers(expr, warning = function(w) if (inherits(w, . classes)) tryInvokeRestart("muffleWarning"))
  4. ds.rSAC.bootstrap(histFile, r = 1, times = times)
  5. f.rSACs[times]

when I try it on estimateLibComplexity(readsDupFreq("GL1.bam")) it works well but if I work with GL2.bam and GL3.bam from package extdata files similar error comes up.

I am using ATACseqQC Version: 1.12.5

could you please let me know what is difference between them and how can i resolve it?

Thanks,

Noori

noorisotoudeh avatar Oct 16 '20 20:10 noorisotoudeh