[BioC] marray: problems loading genepix-files
Jean Yee Hwa Yang
jeany at maths.usyd.edu.au
Sat Mar 17 13:05:04 CET 2007
Hi Heike
I am guessing there is a mismatch between the gal file and the gpr file
you are reading in. Usually the first 5 columns of any gpr file is the
same as the gal file.
Please try the two options and let me know how you go
Option 1
data = read.GenePix(fnames = c("a.gpr", "ar.gpr"))
Option 2
myTargetsa <- read.marrayInfo(file.path(data.dir, "targeta.txt"))
data = read.GenePix(targets = myTargetsa)
Please also show me myTargetsa. If you have further problems, you are
welcome to send me a link where I can download your gpr file and take a
look.
Cheers
Jean
> #################################################################
> > myTargetsa <- read.marrayInfo(file.path(data.dir, "targeta.txt"))
> > galinfo<-read.Galfile("3kArray5.gal")
> >
> data<-read.GenePix(c("a.gpr","ar.gpr"),targets=myTargetsa,layout=galinfo$layout,gnames=galinfo$gnames,)
>
> Reading ... a.gpr
> Error in if (skip > 0) readLines(file, skip) :
> missing value where TRUE/FALSE needed
> #################################################################
> > traceback()
> 4: read.table(f, skip = skip2, header = TRUE, sep = sep, quote = quote,
> check.names = FALSE, as.is = TRUE, comment.char = "", nrows =
> nspots,
> ...)
> 3: read.marrayRaw(fnames = c("a.gpr", "ar.gpr"), name.Gf = "F532 Median",
> name.Gb = "B532 Median", name.Rf = "F635 Median", name.Rb =
> "B635 Median",
> name.W = "Flags", layout = <S4 object of class "marrayLayout">,
> gnames = <S4 object of class "marrayInfo">, targets = <S4 object
> of class "marrayInfo">,
> notes = "GenePix Data", sep = "\t", quote = "\"", DEBUG = FALSE,
> path = NULL, skip = NULL, fill = TRUE)
> 2: do.call("read.marrayRaw", maRaw.args)
> 1: read.GenePix(c("a.gpr", "ar.gpr"), targets = myTargetsa, layout =
> galinfo$layout,
> gnames = galinfo$gnames, )
> #################################################################
More information about the Bioconductor
mailing list