[BioC] problem read.maimage("Agilent") -limma
Naomi Altman
naomi at stat.psu.edu
Mon Jul 25 18:22:22 CEST 2005
I am having trouble reading the Agilent arabidopsis 22575 gene array using
read.maimage in Limma under R 2.1.1 (I don't know the limma version, but I
just downloaded using the R packages interface, and also used the update,
so I presume this is the most recent.
Under R 2.0.1, there was no problem reading all the data in the arrays using:
RGf=read.maimages(c("2792.txt","2793.txt","2796.txt","4507.txt","4508.txt","4509.txt"),source="agilent"
)
dim(RGf$R)
22575 6
But under R 2.1.I, I get:
RGf=read.maimages(c("2792.txt","2793.txt","2796.txt","4507.txt","4508.txt","4509.txt"),source="agilent"
)
dim(RGf$R)
12956 6
The last line of RGf$R is all NA.
The problem might be in RGf$genes. When I try to print any row up to the
last one, everything looks normal. Trying to print the last row kills
R. The annotation for this gene appears to be exceptionally long.
Naomi S. Altman 814-865-3791 (voice)
Associate Professor
Bioinformatics Consulting Center
Dept. of Statistics 814-863-7114 (fax)
Penn State University 814-865-1348 (Statistics)
University Park, PA 16802-2111
More information about the Bioconductor
mailing list