[BioC] problems normalizing in limma

Anand C. Patel acpatel at usa.net
Thu Oct 13 18:38:44 CEST 2005


A quick look at the CSV files generated by PerkinElmer's ScanArray  
finds a bottom row that only has in it "END DATA".  Thus, as Dr.  
Smyth predicted:
 > RG2$R[38977,]
slide_13295067 slide_13295072 slide_13295071 slide_13295073  
slide_13295075 slide_13295079 slide_13295080 slide_13295203  
slide_13295068 slide_13295214
             NA             NA             NA              
NA             NA             NA             NA              
NA             NA             NA
slide_13295879 slide_13295887
             NA             NA
 >
Are we the only people using PerkinElmer's ScanArray software?  I did  
my due diligence in terms of searching the list archives prior to  
posting, and didn't see any other posts.  It is frustrating that even  
the "standard" gpr format would not work, and that the csv format  
output is nonstandard.  However, it doesn't seem like it should be  
terribly difficult to modify the genepix routine in read.maimages  
into a "scanarray_gpr" option, much as was done with bluefuse.  Is  
this something I should work through and submit?  Also, besides the  
obvious step of going through the CSV files and erasing the last row,  
is there an easy way to remove row 38977 from the RGlist as a whole?

Much thanks to Dr. Smyth for both limma and his assistance.

Thanks,
Anand C. Patel, MD
Washington University School of Medicine
acpatel at usa.net

 >[BioC] problems normalizing in limma
 >Gordon Smyth smyth at wehi.edu.au
 >Thu Oct 13 12:25:57 CEST 2005
 >
 >Previous message: [BioC] annotations for Codelink arrays
 >Next message: [BioC] Source Code of affy package
 >Messages sorted by: [ date ] [ thread ] [ subject ] [ author ]
 >

 >The output you give shows that all your intensity matrices (R, G,  
Rb and
 >Gb) are one row longer than your annotation data.frame RG2$genes. This
 >cannot be. The layout information agree with the dimension of RG2 
$genes.
 >
 >My guess is that your input gpr files contained a spurious extra  
row which
 >is not real data. If so, you need to fix them to remove this  
spurious row,
 >or simply remove the last row of all the intensity columns. (Check  
first
 >that RG2$R[38977,] is not real data.)
 >
 >Gordon



More information about the Bioconductor mailing list