[BioC] Bug in featureNames() in R-2.5.0 / BioC 2.0?

Groot, Philip de philip.degroot at wur.nl
Tue May 22 13:50:31 CEST 2007


Hello list,
 
While checking my scripts against the new R/BioC versions, I encountered the following unexpected behaviour. Please consider the following peace of code:
 
library(simpleaffy)
x <- ReadAffy()
 
> x
AffyBatch object
size of arrays=1002x1002 features (14 kb)
cdf=Mouse430_2 (45101 affyids)
number of samples=6
number of genes=1004004
annotation=mouse4302
notes=

> featureNames(x)[1:5]
[1] "1" "2" "3" "4" "5"

> geneNames(x)[1:5]
[1] "1415670_at"   "1415671_at"   "1415672_at"   "1415673_at"   "1415674_a_at"

ADDITIONALLY: when I do the following, the output is correct:
 
x.norm <- rma(x)
 
> featureNames(x.norm)[1:5]
[1] "1415670_at"   "1415671_at"   "1415672_at"   "1415673_at"   "1415674_a_at"
> geneNames(x.norm)[1:5]
[1] "1415670_at"   "1415671_at"   "1415672_at"   "1415673_at"   "1415674_a_at"
Warning message:
'geneNames' is deprecated.
Use 'featureNames' instead.
See help("Deprecated")

 
Because of the "deprecated" message, I changed all geneNames() commands in featureNames() commands and hence found this inconsistency. I guess that my assumption is valid that the commands featureNames and geneNames should always obtain exactly the same results, is it not?
If not, please explain!
 
R/BioC is compiled and runnung on the following linux distribution:
SUSE LINUX Enterprise Server 9 (x86_64)
VERSION = 9
PATCHLEVEL = 3

The gcc version used is:
gcc (GCC) 3.3.3 (SuSE Linux)
Copyright (C) 2003 Free Software Foundation, Inc.
This is free software; see the source for copying conditions.  There is NO
warranty; not even for MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.

Regards,
 
Dr. Philip de Groot
Wageningen University



More information about the Bioconductor mailing list