Question: File not found error when using dupRadar::analyzeDuprates
gravatar for b.nota
15 months ago by
b.nota340 wrote:

I am trying to use dupRadar on bam files, but get an Rsubread error. However, Rsubread is properly installed and no errors there. Is this a bug? Here my code:


bamDuprm <- markDuplicates(dupremover="picard",

# The call parameters:
bamDuprm <- "Sample01Aligned.sortedByCoord.out_duprm.bam"    # the duplicate marked bam file
gtf <- "~/Reference/mm10/Mus_musculus/Mus_musculus.GRCm38.92.gtf"  # the gene model
stranded <- 0       # '0' (unstranded), '1' (stranded) and '2' (reversely stranded)
paired   <- FALSE   # is the library paired end?
threads  <- 20       # number of threads to be used

# Duplication rate analysis
dm <- analyzeDuprates(bamDuprm,gtf,stranded,paired,threads)

Error in file(file, "rt") : cannot open the connection
In addition: Warning message:
In file(file, "rt") :
  cannot open file './.Rsubread_featureCounts_pid17093': No such file or directory

But Rsubread is installed:

NCBI RefSeq annotation for mm10 (build 38.1) is used.
Error in paste(files, collapse = ";") :
  argument "files" is missing, with no default

> sessionInfo()
R version 3.4.4 (2018-03-15)
Platform: x86_64-pc-linux-gnu (64-bit)
Running under: Ubuntu 16.04.5 LTS

Matrix products: default
BLAS: /usr/lib/libblas/
LAPACK: /usr/lib/lapack/

 [1] LC_CTYPE=en_US.UTF-8       LC_NUMERIC=C
 [3] LC_TIME=nl_NL.UTF-8        LC_COLLATE=en_US.UTF-8
 [7] LC_PAPER=nl_NL.UTF-8       LC_NAME=C
 [9] LC_ADDRESS=C               LC_TELEPHONE=C

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base

other attached packages:
[1] Rsubread_1.26.1 dupRadar_1.6.0

loaded via a namespace (and not attached):
[1] compiler_3.4.4 tools_3.4.4

Anyone seen this before? Any solutions? Thanks.

rsubread dupradar • 357 views
ADD COMMENTlink modified 15 months ago by Gordon Smyth39k • written 15 months ago by b.nota340

I thought originally that Rsubread doesn't make the file that is mentioned in the error message, but it turns out that it does. See Wei Shi's answer below.

If you install R 3.5.1 and the latest version of Bioconductor, the code should run fine.

ADD REPLYlink modified 15 months ago • written 15 months ago by Gordon Smyth39k

Thank you, I think it is the tilde (~) in my path name. I thought R would understand that.

ADD REPLYlink written 15 months ago by b.nota340
Answer: C: File not found error when using dupRadar::analyzeDuprates
gravatar for Wei Shi
15 months ago by
Wei Shi3.2k
Wei Shi3.2k wrote:

This file path issue has been fixed in the recent Rsubread releases. Please update your Rsubread and R to their latest. The latest Rsubread version is 1.30.5.

ADD COMMENTlink written 15 months ago by Wei Shi3.2k

Thank you, I will check the newest version.

ADD REPLYlink written 15 months ago by b.nota340
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 385 users visited in the last hour