This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Seeks good way to convert GB_LIST to SYMBOL

Take GPL5188 as example. This GPL map probe id to GB_LIST(genebank).

> gpl <- getGEO("GPL5188", AnnotGPL=FALSE, destdir=".")
> as_tibble(Table(gpl)) %>% dplyr::select(ID, GB_LIST) %>% head(n = 3)
# A tibble: 3 x 2
       ID GB_LIST                     
    <dbl> <chr>                       
1 2315101 NR_024005,NR_024004,AK093685
2 2315102 NR_024004,NR_024005,AK093685
3 2315103 ""

I want to map probe id to SYMBOL, so I use ACCNUM and SYMBOL field in org.Hs.eg.db.

> gblist <- c("NR_024004", "AK093685", "NR_024005", "AB007962", "BC063535")
> gbMap <- AnnotationDbi::select(org.Hs.eg.db, keys=gblist, keytype="ACCNUM", columns=c("ACCNUM", "SYMBOL"))
'select()' returned 1:1 mapping between keys and columns
> head(gbMap)
     ACCNUM    SYMBOL
1 NR_024004   DDX11L2
2  AK093685   DDX11L2
3 NR_024005   DDX11L2
4  AB007962 LINC00869
5  BC063535 LINC00623

This GPL5188 has multiple id in GB_LIST separate by ",", I only take the first one to map to symbol. But I lost most of my probe after doing that. Maybe some symbol are actually map to 2/3 one in GB_LIST. Any other way to convert such GB_LIST to SYMBOL? Thanks.

geo

0 answers

No answers yet.

Log in to answer this question.