[BioC] limma: read.imagene patch

Gordon Smyth smyth at wehi.edu.au
Thu Oct 2 19:23:25 MEST 2003


I have implemented this in limma 1.2.7.

Gordon

At 10:43 PM 1/10/2003, Yves Bastide wrote:
>Gordon Smyth wrote:
>>Dear Yves,
>>Thanks for the suggested code changes. Before I make commit changes to 
>>the code, would you mind showing me an example of some headers which 
>>change length? I know that different versions of Imagene produce 
>>different headers, but read.imagene() like read.maimages() is designed to 
>>read a batch of arrays which should in principle have the same layout, so 
>>I would be uncomfortable about users reading in unrelated arrays. I just 
>>want to see in what ways the header can change within one experiment.
>
>You're right, this patch is way too broad; I don't think the headers for 
>the cy3 and cy5 files of a measure can change.
>
>>Note also that if you have two batches of arrays scanned with different 
>>settings, you can read them into two RGList objects and then combine the 
>>objects using cbind.
>
>This is the case I met, with slightly different quality settings between 
>the first and subsequent scans.  I think it'd be convenient having 
>read.imagene() handle this.
>
>Sketch:
>
>for (i in 1:narrays) {
>   fullname <- ...
>+  headers <- getImageneHeaders(fullname)
>+  if (headers$Field.Dimensions[...] != printer) error(...)
>+  skip <- headers$Begin.Raw.Data
>   obj<-read.table(...)
>
>>Cheers
>>Gordon
>
>Regards,
>
>yves



More information about the Bioconductor mailing list