[BioC] read.maimages in limma package

john seers (IFR) john.seers at bbsrc.ac.uk
Wed Dec 17 10:36:00 CET 2008


Hi Giusy

If you look at the code for read.maimages you can see which columns
limma is expecting for agilent:

agilent = list(G = "gMeanSignal", 
            Gb = "gBGMedianSignal", R = "rMeanSignal", Rb =
"rBGMedianSignal"),


If your custom arrays do not match these names then you have to specify
them yourself (using the columns parameter of read.maimages).


Regards

John Seers

 
---



Hi!
I am analyzing custom two color agilent array by using limma package and
I am having some problems with the read.maimages
comand:

I launch this comand:

> read.maimages("Targets$Filename, source="agilent")

and the error that I have is the following:

Read Control_2b.txt
Read Control_3b.txt
Error in readGenericHeader(fullname, columns = columns) :
Specified column headings not found in file

It is like the columns in the file that I gave as input (Control_.txt)
are not recognized,
maybe I have this problem because is a custom agilent array?

If this is the problem, I have to specify in my script the following
informations:

1. list with fields R, G, Rb and Gb giving the column names to be used
for red and green foreground and background;
2. character vector of names of other columns to be read containing
spot-specific information;
3. character vector of names of columns containing annotation
information about the probes;

Is there someone that can answer me and that has got some experience in
the 
two color agilent custom array analysis?


Thank you in advance
Giusy

_______________________________________________
Bioconductor mailing list
Bioconductor at stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/bioconductor
Search the archives:
http://news.gmane.org/gmane.science.biology.informatics.conductor



More information about the Bioconductor mailing list