Failed to run with control IgG

[SergioRodLla]

Hi,

First of all thank you for developing SEACR. I'm trying to run the script using a control sample as follows:

bash $seacr AI02_K4me3.normalized.bedgraph controls/AI02_IgG.normalized.bedgraph non stringent /media/sergio/HDD1/cutntag/peakCalling/SEACR/AI02_K4me3.seacr_control.peaks

I'm getting an output error:

Calling enriched regions with control file
Proceeding without normalization of control to experimental bedgraph
Using stringent threshold
Creating experimental AUC file: mié 31 may 2023 16:01:29 CEST
Creating control AUC file: mié 31 may 2023 16:01:29 CEST
Calculating optimal AUC threshold: mié 31 may 2023 16:01:29 CEST
Calculating threshold using non-normalized control: mié 31 may 2023 16:01:29 CEST
Error in while (output == 0) { : missing value where TRUE/FALSE needed
Execution halted

I tried with a different sample and the output message said:

Calling enriched regions with control file
Proceeding without normalization of control to experimental bedgraph
Using stringent threshold
Creating experimental AUC file: mié 31 may 2023 15:59:57 CEST
Creating control AUC file: mié 31 may 2023 15:59:58 CEST
Calculating optimal AUC threshold: mié 31 may 2023 15:59:58 CEST
Calculating threshold using non-normalized control: mié 31 may 2023 15:59:58 CEST
Warning messages:
1: In ((pctremain(x0) + min(pctremain(z)))/2) - pctremain(z) :
longer object length is not a multiple of shorter object length
2: In ((pctremain(x0) + min(pctremain(z)))/2) - pctremain(z) :
longer object length is not a multiple of shorter object length
Creating thresholded feature file: mié 31 may 2023 15:59:59 CEST
Empirical false discovery rate = 0
Merging nearby features and eliminating control-enriched features: mié 31 may 2023 15:59:59 CEST
Removing temporary files: mié 31 may 2023 15:59:59 CEST
Done: mié 31 may 2023 15:59:59 CEST

However, I tried to run the same using the web server and it ran without any problem. I am trying to understand what could be happening.

Thank you,
Sergio