ArrayQuality package bug in agQuality function - marray package
0
0
Entering edit mode
Marcus Davy ▴ 680
@marcus-davy-374
Last seen 10.3 years ago
Hi, there appears to be a bug in the arrayQuality package, specifically the agQuality function which traces back to the function read.Agilent in the marray package. If I specify a path for the location of Agilent microarray files to process which is not the current working directory (path=".", and without specifying 'gnames' information) then I get the following error; >agQuality(fnames=targets$Filename, path=dataDir, resdir = diagPlotDir) [1] "Starting agQuality..." Error in file(con, "r") : cannot open the connection In addition: Warning message: In file(con, "r") : cannot open file '[FILE.txt]': No such file or directory > traceback() 4: file(con, "r") 3: readLines(fullnames[1], n = 100) 2: read.Agilent(fnames = fnames[1], path = path, DEBUG = DEBUG) 1: agQuality(fnames = targets$Filename, path = dataDir, resdir = diagPlotDir) > find("read.Agilent") [1] "package:marray" In 'read.Agilent', the 'fullnames' vector defined on line 15 does not contain fullpath information ('fullfnames') so the readLines function on line 24 fails. If 'gnames' was specified it would probably work. A hack to get it going was to make a symbolic link for the first file to read in the description information from the first file. ln -s /full/path/to/first/file file This temporary workaround works without modifying any R code. If I try putting the full path name 'fnames' in read.Agilent with 'path=NULL' as suggested in help(read.Agilent); > agQuality(fnames=file.path(dataDir, targets$Filename), path=NULL, resdir = diagPlotDir) [1] "Starting agQuality..." Reading ... .//[PATH TO FILE] Error in file(file, "r") : cannot open the connection In addition: Warning message: In file(file, "r") : cannot open file './/[PATH TO FILE]': No such file or directory > > traceback() 6: file(file, "r") 5: read.table(f, skip = skip2, header = TRUE, sep = sep, quote = quote, check.names = FALSE, as.is = TRUE, comment.char = "", nrows = nspots, ...) 4: read.marrayRaw(fnames = "/[PATH TO FILE]", path = ".", name.Gf = "gMedianSignal", name.Gb = "gBGMedianSignal", name.Rf = "rMedianSignal", name.Rb = "rBGMedianSignal", layout = <s4 object="" of="" class="" "marraylayout"="">, gnames = <s4 object="" of="" class="" "marrayinfo"="">, notes = "Agilent Data", skip = 9, sep = "\t", quote = "\"", DEBUG = FALSE) 3: do.call("read.marrayRaw", maRaw.args) 2: read.Agilent(fnames = fnames[1], path = path, DEBUG = DEBUG) 1: agQuality(fnames = file.path(dataDir, targets$Filename), path = NULL, resdir = diagPlotDir) Marcus > sessionInfo() R version 2.7.0 (2008-04-22) powerpc-apple-darwin8.10.1 locale: C attached base packages: [1] grid tools stats graphics grDevices utils datasets [8] methods base other attached packages: [1] arrayQuality_1.19.0 RColorBrewer_1.0-2 gridBase_0.4-3 [4] hexbin_1.14.0 lattice_0.17-6 convert_1.16.0 [7] Biobase_2.0.1 marray_1.19.0 limma_2.14.1 > The contents of this e-mail are privileged and/or confid...{{dropped:5}}
Microarray arrayQuality Microarray arrayQuality • 1.1k views
ADD COMMENT

Login before adding your answer.

Traffic: 625 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6