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

PHANTOMPEAKQUALTOOLS throws stack overflow exception #370

Closed
krokicki opened this issue Nov 22, 2023 · 1 comment
Closed

PHANTOMPEAKQUALTOOLS throws stack overflow exception #370

krokicki opened this issue Nov 22, 2023 · 1 comment
Assignees
Labels
bug Something isn't working

Comments

@krokicki
Copy link

Description of the bug

PHANTOMPEAKQUALTOOLS throws a stack overflow exception on certain types of input data. This issue suggests a fix by increasing the default --max-ppsize. This is already implemented in MULTIQC_CUSTOM_PHANTOMPEAKQUALTOOLS and should be added to PHANTOMPEAKQUALTOOLS as well.

Command used and terminal output

nextflow run main.nf -profile janelia --input samples.csv --read_length 50 --outdir ./results --fasta internal.fasta --gff internal.gff
ERROR ~ Error executing process > 'NFCORE_CHIPSEQ:CHIPSEQ:PHANTOMPEAKQUALTOOLS (neg_rep1)'

Caused by:
  Process `NFCORE_CHIPSEQ:CHIPSEQ:PHANTOMPEAKQUALTOOLS (neg_rep1)` terminated with an error exit status (1)

Command executed:

  RUN_SPP=`which run_spp.R`
  Rscript  -e "library(caTools); source(\"$RUN_SPP\")" -c="neg_rep1.mLb.clN.sorted.bam" -savp="neg_rep1.spp.pdf" -savd="neg_rep1.spp.Rdata" -out="neg_rep1.spp.out" -p=6

  cat <<-END_VERSIONS > versions.yml
  "NFCORE_CHIPSEQ:CHIPSEQ:PHANTOMPEAKQUALTOOLS":
      phantompeakqualtools: 1.2.2
  END_VERSIONS

Command exit status:
  1

Command output:
  ################
  ChIP data: neg_rep1.mLb.clN.sorted.bam
  Control data: NA
  strandshift(min): -500
  strandshift(step): 5
  strandshift(max) 1500
  user-defined peak shift NA
  exclusion(min): 10
  exclusion(max): NaN
  num parallel nodes: 6
  FDR threshold: 0.01
  NumPeaks Threshold: NA
  Output Directory: .
  narrowPeak output file name: NA
  regionPeak output file name: NA
  Rdata filename: neg_rep1.spp.Rdata
  plot pdf filename: neg_rep1.spp.pdf
  result filename: neg_rep1.spp.out
  Overwrite files?: FALSE

  Reading ChIP tagAlign/BAM file neg_rep1.mLb.clN.sorted.bam
  opened /tmp/RtmpEulNd9/neg_rep1.mLb.clN.sorted.tagAlign2961e49a6f
  done. read 326050 fragments

Command error:
  INFO:    Environment variable SINGULARITYENV_NXF_DEBUG is set, but APPTAINERENV_NXF_DEBUG is preferred
  WARNING: Skipping mount /var/apptainer/mnt/session/etc/resolv.conf [files]: /etc/resolv.conf doesn't exist in container
  Loading required package: Rcpp
  Error: protect(): protection stack overflow
  Execution halted

Relevant files

No response

System information

No response

@krokicki krokicki added the bug Something isn't working label Nov 22, 2023
@krokicki krokicki mentioned this issue Nov 22, 2023
9 tasks
@krokicki
Copy link
Author

Fixed by #371.

@krokicki krokicki self-assigned this Nov 27, 2023
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
bug Something isn't working
Projects
None yet
Development

No branches or pull requests

1 participant