Hello everyone,
I am learning how to use "biomaRt" following Bioconductor 2018 Workshops.
When I ran the query example, I got a result that is different from the tutorial.
The expected result was:
mart <- useMart("ENSEMBL_MART_ENSEMBL","hsapiens_gene_ensembl")
afyids <- c("1000_at","1001_at","1002_f_at","1007_s_at")
getBM(c("affy_hg_u95av2", "hgnc_symbol"), c("affy_hg_u95av2"), afyids, mart)
#> affy_hg_u95av2 hgnc_symbol
#> 1 1000_at MAPK3
#> 2 1007_s_at DDR1
#> 3 1002_f_at
#> 4 1002_f_at CYP2C19
#> 5 1001_at TIE1
But I got:
affy_hg_u95av2 hgnc_symbol
**1000_at GDPD3**
1007_s_at DDR1
**1000_at MAPK3**
1001_at TIE1
1002_f_at
1002_f_at CYP2C19
The probe "1000_at" can map to GDPD3 and MAPK3, in which the MAPK3 should be right according to the Ensembl website.
Would you please help me to figure out the problem? Thank you very much.
Here is the system information:
Microsoft Windows [Version 10.0.14393] R: 3.5.2 RStudio: Version 1.1.456 biomaRt: 2.38.0
Thank you very much.