rnaseq icon indicating copy to clipboard operation
rnaseq copied to clipboard

Error during dupRadar

Open scintilla9 opened this issue 1 year ago • 1 comments

Description of the bug

Hi,

I was trying to run the pipeline of all my samples, however one sample keeped fail during dupRadar. I also attach the log file below. Any suggestion?

Thanks.

Command used and terminal output

ERROR ~ Error executing process > 'NFCORE_RNASEQ:RNASEQ:DUPRADAR (Adh18-S_REP3)'

Caused by:
  Process `NFCORE_RNASEQ:RNASEQ:DUPRADAR (Adh18-S_REP3)` terminated with an error exit status (1)

Command executed:

  dupradar.r \
      Adh18-S_REP3.markdup.sorted.bam \
      Adh18-S_REP3 \
      annot.gtf \
      2 \
      paired \
      1

  cat <<-END_VERSIONS > versions.yml
  "NFCORE_RNASEQ:RNASEQ:DUPRADAR":
      r-base: $(echo $(R --version 2>&1) | sed 's/^.*R version //; s/ .*$//')
      bioconductor-dupradar: $(Rscript -e "library(dupRadar); cat(as.character(packageVersion('dupRadar')))")
  END_VERSIONS

Command exit status:
  1

Command output:
  (empty)

Command error:
  Input bam      (Arg 1): Adh18-S_REP3.markdup.sorted.bam
  Input gtf      (Arg 2): annot.gtf
  Strandness     (Arg 3): reverse
  paired/single  (Arg 4): paired
  Nb threads     (Arg 5): 1
  R package loc. (Arg 6): paired
  Output basename       : Adh18-S_REP3
  Loading required package: dupRadar
  Loading required package: parallel
  Error in if (trim > 0 && n) { : missing value where TRUE/FALSE needed
  Calls: duprateExpDensPlot -> mean -> mean.default
  Execution halted

Relevant files

nextflow.log

System information

No response

scintilla9 avatar Jun 07 '23 09:06 scintilla9