getGene in biomaRt package
0
0
Entering edit mode
sdurinck ▴ 40
@sdurinck-1480
Last seen 10.2 years ago
Hi James, This is probably due to a change in the VEGA BioMart database and I will check this out. You should be able to use Ensembl though to get gene information from entrezgene identifiers using quotes or not. Try the following: > mart=martConnect() connected to: ensembl_mart_35 > getGene(id=100, type="entrezgene", species="hsapiens", mart=mart) id symbol 1 100 ADA description 1 Adenosine deaminase (EC 3.5.4.4) (Adenosine aminohydrolase). [Source:Uniprot/SWISSPROT;Acc:P00813] band chromosome start end strand martID 1 q13.12 20 42681578 42713790 -1 ENSG00000196839 > getGene(id="100", type="entrezgene", species="hsapiens", mart=mart) id symbol 1 100 ADA description 1 Adenosine deaminase (EC 3.5.4.4) (Adenosine aminohydrolase). [Source:Uniprot/SWISSPROT;Acc:P00813] band chromosome start end strand martID 1 q13.12 20 42681578 42713790 -1 ENSG00000196839 I get an error when I run the example for getGene in the biomaRt package on both windows and linux. It appears that the martTable class expects the id to be character, but for e.g. Entrez Gene the id is numeric. If I substitute id="100", I get a martTable, but with all NA for all slots but the id. >> example(getGene) getGen> if (interactive()) { mart <- martConnect(biomarts = c("ensembl", "vega")) g = getGene(id = "1939_at", array = "hg_u95av2", mart = mart) show(g) g = getGene(id = 100, type = "entrezgene", species = "hsapiens", db = "vega", mar .... [TRUNCATED] connected to: ensembl_mart_35 connected to: vega_mart_35 An object of class "martTable" Slot "id": [1] "1939_at" Slot "table": $symbol [1] "TP53" $description [1] "Cellular tumor antigen p53 (Tumor suppressor p53) (Phosphoprotein p53) (Antigen NY-CO-13). [Source:Uniprot/SWISSPROT;Acc:P04637]" $band [1] "p13.1" $chromosome [1] "17" $start [1] 7512464 $end [1] 7531642 $martID [1] "ENSG00000141510" Error in validObject(.Object) : invalid class "martTable" object: invalid object for slot "id" in class "martTable": got class "numeric", should be or extend class "character" >> sessionInfo() R version 2.2.0, 2005-10-12, i386-pc-mingw32 attached base packages: [1] "tools" "methods" "stats" "graphics" "grDevices" "utils" "datasets" [8] "base" other attached packages: biomaRt XML RMySQL DBI hs133phsug6cdf affy "1.4.0" "0.97-0" "0.5-6" "0.1-9" "6.0.0" "1.8.1" Biobase "1.7.4" Any ideas? Jim James W. MacDonald Affymetrix and cDNA Microarray Core University of Michigan Cancer Center 1500 E. Medical Center Drive 7410 CCGC Ann Arbor MI 48109 734-647-5623 ********************************************************** Electronic Mail is not secure, may not be read every day, and should not be used for urgent or sensitive issues. _______________________________________________ Bioconductor mailing list Bioconductor at stat.math.ethz.ch htt Disclaimer: http://www.kuleuven.be/cwis/email_disclaimer.htm
Microarray Cancer cdf biomaRt Vega Microarray Cancer cdf biomaRt Vega • 1.1k views
ADD COMMENT

Login before adding your answer.

Traffic: 841 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