Hello,

I tried to find gene names that match affyIDs, and I tried two ways as below, but neither worked. I got all NAs for the Symbol column. I have no idea what happened. Thank you very much for helping me out!

(1) library(affy)
mydata <- ReadAffy()
eset <- rma(mydata)
library(annotate)
library(hgu95av2.db)
ID <- featureNames(eset)
Symbol <- getSYMBOL(ID,"hgu95av2.db")
fData(eset) <- data.frame(ID=ID,Symbol=Symbol)

All NAs for Symbol in the topTable

(2) fit$genes$Symbol <- getSYMBOL(fit$genes$ID, "hgu95av2.db")

>fit
$genes
       ID        Symbol
1 7892501     NA
2 7892502     NA
3 7892503     NA
4 7892504     NA
5 7892505     NA
32316 more rows ...

Thanks,
Xiayu

	[[alternative HTML version deleted]]

