Hi all,
I'm trying to get probes for particular probset id from my MoGene arrays. From experiment description file (dabg.summary.txt) I can see that there are around 241,500 probset ids for my 12 arrays. When loading .CEL files into bioconductor I see 1,102,500 values for my 12 arrays. Thus I think there should be around 4 probes per 1 probeset on average.
However, when I load an experiment description file into an AnnotatedDataFrame object:
Affy.Expt <- read.AnnotatedDataFrame("dabg.summary.txt", header=TRUE, row.names=1, sep="\t")
 and try to use it as my phenoData when loading .CEL files into Affybatch object :
Affy.Data <- ReadAffy(filenames=colnames(pData(Affy.Expt)), phenoData=Affy.Expt, verbose=TRUE)
I get an error:
Warning message:
In read.affybatch(filenames = l$filenames, phenoData = l$phenoData,  :
  Incompatible phenoData object. Created a new one.
I understand that as a not consistent number of rows between my experiment description file (241,500 probset ids) and number of rows in .CEL files (1,102,500 probes). When it does that it resets the probsets id and starts numbering the rows from 1 to 1,102,500 and thus loosing track of probset ids.

The point is that I need to know which probes belong to which probeset id and have their values stored.
I looked at CDF file but it looks strange and I can't get anything useful from there. I thought that maybe looking into rma algorithm will help me out somehow, but it calls external function, which I don't understand.
Is there some easy way to get that information?

Thank you in advance,
Mark


	[[alternative HTML version deleted]]

