Question: Cannot open file because of connection error
gravatar for hallie.s
14 months ago by
hallie.s0 wrote:

I am following an NGS ChIP-seq tutorial that requires me to set a working directory and use the function systemArgs. I have successfully created my working directory using this code:

> library(systemPipeRdata)
> genWorkenvir(workflow = "chipseq")
[1] "Generated chipseq directory. Next run in chipseq directory, the R code from *.Rmd (*.Rnw) template interactively. Alternatively, workflows can be exectued with a single command as instructed in the vignette."
> setwd("chipseq")

However, when I try to run the chunk of code that includes the function systemArgs, I get an error reading there is no file and R cannot open the connection. Here is my chunk of code and the error I receive:


args <- systemArgs(sysma = "param/trim.param", mytargets = "targets_chip.txt")
filterFct <- function(fq, cutoff = 20, Nexceptions = 0) {
  qcount <- rowSums(as(quality(fq), "matrix") <= cutoff)
  fq[qcount <= Nexceptions]
preprocessReads(args = args, Fct = "filterFct(fq, cutoff = 20, Nexceptions = 0)", batchsize = 100000)
writeTargetsout(x = args, file = "targets_chip_trim.txt", overwrite = TRUE)


> args <- systemArgs(sysma = "param/trim.param", mytargets = "targets_chip.txt")
cannot open file 'param/trim.param': No such file or directoryError in file(file, "rt") : cannot open the connection

I guess I am asking: how do I fix this error?




ADD COMMENTlink written 14 months ago by hallie.s0
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: 231 users visited in the last hour