[BioC] problem read.maimage("Agilent") -limma

Henrik Bengtsson hb at maths.lth.se
Tue Jul 26 17:14:03 CEST 2005


Recently we detected some problems with internal regexpr libraries in R 
v2.1.1. One of the symptoms was that R would crash on Windows, but also 
that the regular expression became corrupt in memory. This was partly 
fixed in the R v2.1.1 patched (2005-07-20). Note that this was 
introduced when the went from R v2.1.0 to v2.1.1, so this might be 
related to your problem.

Cheers

Henrik

Naomi Altman wrote:
> There are "\" and "#" before the offending line.  I could not find any 
> other unusual characters in the offending line.
> 
> --Naomi
> 
> At 09:59 AM 7/26/2005, Sean Davis wrote:
> 
> 
>>On Jul 26, 2005, at 8:13 AM, Gordon K Smyth wrote:
>>
>>
>>>>Date: Mon, 25 Jul 2005 12:22:22 -0400
>>>>From: Naomi Altman <naomi at stat.psu.edu>
>>>>Subject: [BioC] problem read.maimage("Agilent") -limma
>>>>To: bioconductor at stat.math.ethz.ch
>>>>
>>>>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.
>>>
>>>You should have limma 2.0.2.
>>>
>>>
>>>>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.
>>
>>I have had problems with Agilent annotation files containing "special"
>>characters that cause similar "termination" of file reading.  I would
>>look at the annotation for quotation marks, single quotes, # symbols
>>(no idea why this seems to affect things), and backslashes.  I
>>typically write a little perl script to "clean" the files.  I'm not
>>sure why this should vary from one version to the next, though.
>>
>>Sean
>>
> 
> 
> 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
> 
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/bioconductor
> 
>



More information about the Bioconductor mailing list