Question: File not found error when using dupRadar::analyzeDuprates
0
gravatar for b.nota
9 months ago by
b.nota320
Netherlands
b.nota320 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:

library(dupRadar)

bamDuprm <- markDuplicates(dupremover="picard",
               bam="Sample01Aligned.sortedByCoord.out.bam",
               path="~/picard-2.17.6",
               rminput=FALSE)

# 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:

library(Rsubread)
featureCounts()
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/libblas.so.3.6.0
LAPACK: /usr/lib/lapack/liblapack.so.3.6.0

locale:
 [1] LC_CTYPE=en_US.UTF-8       LC_NUMERIC=C
 [3] LC_TIME=nl_NL.UTF-8        LC_COLLATE=en_US.UTF-8
 [5] LC_MONETARY=nl_NL.UTF-8    LC_MESSAGES=en_US.UTF-8
 [7] LC_PAPER=nl_NL.UTF-8       LC_NAME=C
 [9] LC_ADDRESS=C               LC_TELEPHONE=C
[11] LC_MEASUREMENT=nl_NL.UTF-8 LC_IDENTIFICATION=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 • 263 views
ADD COMMENTlink modified 9 months ago by Gordon Smyth37k • written 9 months ago by b.nota320

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 9 months ago • written 9 months ago by Gordon Smyth37k

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

ADD REPLYlink written 9 months ago by b.nota320
Answer: C: File not found error when using dupRadar::analyzeDuprates
3
gravatar for Wei Shi
9 months ago by
Wei Shi3.1k
Australia
Wei Shi3.1k 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 9 months ago by Wei Shi3.1k

Thank you, I will check the newest version.

ADD REPLYlink written 9 months ago by b.nota320
Please log in to add an answer.

Help
Access

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