[BioC] R pathview and mouse
Luo Weijun
luo_weijun at yahoo.com
Fri Jan 10 16:32:56 CET 2014
Scott,
Gene symbols and other IDs are case sensitive. For mouse and rat, most OFFICIAL gene symbols are initial upper case only, while all upper cases for human and many other species. For example, some mappable gene symbols (from your list) for mouse would be like Cd74, Mgp and Cd52, while they should be CD74, MGP and CD52 for human.
You may tried gene symbol search at NCBI like http://www.ncbi.nlm.nih.gov/gene/?term=CD52 and look at the results for different species.
Weijun
--------------------------------------------
On Thu, 1/9/14, Scott wrote:
Subject: R pathview and mouse
Date: Thursday, January 9, 2014, 4:41 PM
Hello Weijun,
I am trying to map mouse gene names to entrez ids using the
id2eg in your path view tool. Below shows what I
did. The variable “T” holds the data, and I do a
map, but none of the genes get mapped to entrez gene
ids. However, if I use human (“Hs”) then it
works. Do you know what is wrong?
> biocLite("org.Mm.eg.db")
> T
H2-AB1
HSP105
CD74
GP38
MGP CD52
TMEM176B
IFI27 LAPTM5
TGFBI HSP90AB1
3.880227
3.429302
3.417246
3.164127
2.932691
2.852793
2.665297
2.564207
2.534474
2.366048 2.301527
> T.id.map <- id2eg(ids = names(T), category =
gene.idtype.list[1], org = "Mm")
> T.id.map
SYMBOL
EntrezGene
[1,] "H2-AB1" ""
[2,] "HSP105" ""
[3,] "CD74" ""
[4,] "GP38" ""
[5,] "MGP" ""
[6,] "CD52" ""
[7,] "TMEM176B" ""
[8,] "IFI27" ""
[9,] "LAPTM5" ""
[10,] "TGFBI" ""
[11,] "HSP90AB1" ""
However, if I use human then it works.
> T.id.map <- id2eg(ids = names(T), category =
gene.idtype.list[1], org = "Hs")
> T.id.map
SYMBOL
EntrezGene
[1,] "H2-AB1" ""
[2,] "HSP105" ""
[3,] "CD74"
"972"
[4,] "GP38" ""
[5,] "MGP" "4256"
[6,] "CD52"
"1043"
[7,] "TMEM176B"
"28959"
[8,] "IFI27" "3429"
[9,] "LAPTM5" "7805"
[10,] "TGFBI" "7045"
[11,] "HSP90AB1" "3326"
thanks
Scott
More information about the Bioconductor
mailing list