[BioC] Bug in probeNames()
Groot, Philip de
philip.degroot at wur.nl
Thu Jul 26 21:23:14 CEST 2007
Hello list,
I ran into the following particular problem using the Custom CDF-files from Bioconductor 2.0. Please consider the following code:
library(simpleaffy)
x <- ReadAffy(cdfname="mm430mmentrezg")
require(sprintf("%sprobe", cleancdfname(cdfName(x), addcdf=F)), character.only=T)
ProbeNames_CustomCDF <- as.character(as.data.frame(mm430mmentrezgprobe)[,4])
grep("66736_at", ProbeNames_CustomCDF, fixed=T)
173151 173152 173153 173154 173155 173156 173157
However, when I do it the following way different indices are returned:
grep("66736_at", probeNames(x), fixed=T)
172057 172058 172059 172060 172061 172062 172063
I do not understand the reason for this! To get the order of the basepairs (for this particular probe) I need to use the upper grep-command (I extract the basepairs from the probes library), but to get the CORRECT pm() and mm() intensities I need to use the lower grep-command. Why is this? For the Affymetrix annotation libraries everything works fine (as I would expect).
Regards,
Dr Philip de Groot
Wageningen University
More information about the Bioconductor
mailing list