search

FredHutch / SEACR

SEACR: Sparse Enrichment Analysis for CUT&RUN
GNU General Public License v2.0
95 stars 41 forks source link

Error (Empirical false discovery rate...) #61

Open kesozm opened 3 years ago

kesozm commented 3 years ago

Hi,

I used several experimental bedgraph files and IgG data bedgraph file (control) as input. Some of the experimental bedgraph files together with the IgG data succeeded. However, one of the experimental files failed, the error was as the following:

The output (if any) follows:

Calling enriched regions with control file Normalizing control to experimental bedgraph Using relaxed threshold Creating experimental AUC file: Fri Jun 11 01:02:19 CDT 2021 Creating control AUC file: Fri Jun 11 01:02:52 CDT 2021 Calculating optimal AUC threshold: Fri Jun 11 01:04:22 CDT 2021 Calculating threshold using normalized control: Fri Jun 11 01:04:22 CDT 2021 Creating thresholded feature file: Fri Jun 11 01:06:11 CDT 2021 Empirical false discovery rate = 0.571428571428571 Merging nearby features and eliminating control-enriched features: Fri Jun 11 01:06:25 CDT 2021 awk: cmd. line:1: fatal: division by zero attempted

What causes this error?

Thank you.

mpmeers commented 3 years ago

Hi,

Apologies for the delay in responding here. There's an identical error reported in Issue #62 that I wrote a patch for--would you mind taking a look at that response to see if that solves your issue?

Mike