problem with annotating probeset ids.
0
0
Entering edit mode
9.2 years ago

Hi,

I am trying to get through this tutorial on micrarray data analysis by daniel swan. I am however, getting stuck at the annotation part, when I try to annotate the probeset ids with the gene symbols. Here's what I'm getting

> gene.symbols <- getSYMBOL(probeset.list, "hgu133plus2")
Error in unlist(lookUp(x, data, "SYMBOL")) : 
  error in evaluating the argument 'x' in selecting a method for function 'unlist': Error in .checkKeysAreWellFormed(keys) : 
  keys must be supplied in a character vector with no NAs

Also, the following command gives me a NULL result.

>head(probeset.list$ID)
NULL

So the problem is that there is no 'ID' column right? How do I get around this?

ChIP-Seq R genome • 3.4k views
ADD COMMENT
1
Entering edit mode

can you do a head() & class() on probeset.list and show what does it look like and what is its class?

ADD REPLY
1
Entering edit mode

Ok a quick Google search took me here. So you need to do:

gene.symbols <- getSYMBOL(rownames(probeset.list), "hgu133plus2")
ADD REPLY
0
Entering edit mode

Thanks a ton. That sorted out my problem :)

ADD REPLY
0
Entering edit mode

Which tutorial? Can you link to it please.

ADD REPLY

Login before adding your answer.

Traffic: 1362 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6