[BioC] Limma question: readGenericHeader error message
Gordon K Smyth
smyth at wehi.EDU.AU
Wed Nov 23 04:39:17 CET 2011
Dear Ella,
I have just tried out the read.maimages() call sequence you have used on
an Agilent file that I have, and it works perfectly. So I do not know the
problem you have.
>From the output you give, read.maimages() has successfully read your first
four files, but has failed on the fifth. So it would seem that the
specified columns are not contained in all of your files. The fifth file
specified in your targets data.frame is the culprit.
BTW, what you want to do is a bit more easily achieved by
RG <- read.maimages(targets,source="agilent.median")
although I suspect it will give the same results as you already have.
If you post again, please give sessionInfo() output.
Best wishes
Gordon
> Date: Mon, 21 Nov 2011 17:35:07 -0500
> From: Ella Chang <ella7389 at hotmail.com>
> To: <bioconductor at r-project.org>
> Subject: [BioC] Limma question: readGenericHeader error message
>
>
> Hi all,
> I am trying to load my agilent data to Limma using read.maimages function. Here's what I did:
> library(limma)
> targets<-readTargets("targets.txt")
> RG <- read.maimages(targets, columns = list(G = "gMedianSignal", Gb = "gBGMedianSignal", R = "rMedianSignal", Rb = "rBGMedianSignal"), annotation = c("Row", "Col","FeatureNum", "ControlType","ProbeName"))
> but I got the error message says: Read
> US83000170_252018610364_S01_GE2_105_Dec08_1_3.txt Read
> US83000170_252018610364_S01_GE2_105_Dec08_1_4.txt Read
> US83000170_252018610386_S01_GE2_105_Dec08_1_2.txt Read
> US83000170_252018610386_S01_GE2_105_Dec08_1_4.txt Error in
> readGenericHeader(fullname, columns = columns) : Specified column
> headings not found in fileIn addition: Warning message:In file(file,
> "r") : file("") only supports open = "w+" and open = "w+b": using the
> former
>
> I already checked to made sure that all the columns I specified are
> existed in all my file. Does anyone know how to deal with this?
> thanks a lot for your help
> HC
______________________________________________________________________
The information in this email is confidential and intend...{{dropped:4}}
More information about the Bioconductor
mailing list