Question: ArrayQuality package bug in agQuality function - marray package
0
gravatar for Marcus Davy
11.0 years ago by
Marcus Davy680
Marcus Davy680 wrote:
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 • 580 views
ADD COMMENTlink written 11.0 years ago by Marcus Davy680
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: 145 users visited in the last hour