rnaseq
rnaseq copied to clipboard
Error during dupRadar
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
System information
No response