Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

SEACR run failed for control mode #86

Open
sfpacman opened this issue Jun 3, 2022 · 1 comment
Open

SEACR run failed for control mode #86

sfpacman opened this issue Jun 3, 2022 · 1 comment

Comments

@sfpacman
Copy link

sfpacman commented Jun 3, 2022

I wasn't unable to run SEACR (top 0.01 run just fine) when I use control bedgraphs with the following command:
bash SEACR_1.3.sh $sample_bed $control_bed non stringent ${sample}_seacr_ctrl.peak

Using stringent threshold
Creating experimental AUC file: Fri Jun  3 03:17:29 PDT 2022
Creating control AUC file: Fri Jun  3 03:18:20 PDT 2022
Calculating optimal AUC threshold: Fri Jun  3 03:18:30 PDT 2022
Calculating threshold using non-normalized control: Fri Jun  3 03:18:30 PDT 2022
Creating thresholded feature file: Fri Jun  3 03:19:41 PDT 2022
Empirical false discovery rate = 7.62799780318346e-06
Merging nearby features and eliminating control-enriched features: Fri Jun  3 03:19:42 PDT 2022
Error: Unable to determine type for file oS9kQIRuTq0Ya.auc.threshold.bed

The bed file in question either seems to be either corrupted or something in binary.
image

@mpmeers
Copy link
Collaborator

mpmeers commented Jun 27, 2022

Hello,

I'm guessing this is an input issue because there are no compression steps in the SEACR code--is it possible you're inputting a control bigwig instead of a bedgraph? It's possible that wouldn't be caught.

Mike

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants