[BioC] Affymetrix probeset ids to gene symbols

Peter Robinson peter.robinson at t-online.de
Thu Jul 3 22:31:51 CEST 2008


Dear all,

I have a list of affymetrix probeset ids from another program and would 
like to use annaffy to extract the corresponding gene names. I am still 
something of a novice at R and am probably doing something silly, but 
found no answer in the package vignette. My script:


library(annaffy)

dat <- read.table('sign.txt.cdt',header=T)
psets<-dat[,3]
symbols<-aafSymbol(as.character(psets),"moe430b.db")
s<-as.character(symbols)

I was surprisied that so few of the probeset ids got identified by this 
script. WHat am I doing wrong?

THanks Peter
s<-as.character(symbols)
 > s
  [1] "character(0)"       "character(0)"       "character(0)"
  [4] "character(0)"       "character(0)"       "character(0)"
  [7] "character(0)"       "character(0)"       "character(0)"
 [10] "character(0)"       "character(0)"       "character(0)"
 [13] "character(0)"       "character(0)"       "Egr3"
 [16] "character(0)"       "character(0)"       "character(0)"
 [19] "character(0)"       "character(0)"       "character(0)"
 [22] "character(0)"       "character(0)"       "character(0)"
 [25] "Irak2"              "character(0)"       "Coq10b"
 [28] "character(0)"       "BC063749"           "character(0)"
 [31] "4631422O05Rik"      "character(0)"       "Coq10b"
 [34] "character(0)"       "character(0)"       "AI452195"
 [37] "character(0)"       "character(0)"       "character(0)"
 [40] "Mobkl2a"            "character(0)"       "character(0)"

(...snip....)



More information about the Bioconductor mailing list