Anyone know why I’m not getting ENSG ids for some of these symbols?
The example below retrieves `NA` for multiple symbols, including AAED1 [whose ENSG is ENSG00000158122][1].
“`
> library(AnnotationHub)
> library(org.Hs.eg.db)
> library(GEOquery)
> temp download.file(getGEO(“GSM4430459″)@header$supplementary_file_1,temp)
> genes unlink(temp)
> ensids = mapIds(org.Hs.eg.db,
keys=genes,
column=”ENSEMBL”,
keytype=”SYMBOL”,
multiVals=”first”)
> ensids[“AAED1”] # here is one of the
AAED1
NA
“`
[1]: grch37.ensembl.org/Homo_sapiens/Gene/Summary?g=ENSG00000158122;r=9:99401859-99417585