simple biomaRt query
2
0
Entering edit mode
@glazko-galina-1653
Last seen 9.6 years ago
An embedded and charset-unspecified text was scrubbed... Name: not available Url: https://stat.ethz.ch/pipermail/bioconductor/attachments/20070928/ 2be50203/attachment.pl
• 634 views
ADD COMMENT
0
Entering edit mode
Steffen ▴ 500
@steffen-2351
Last seen 9.6 years ago
Dear Galina, Where did you get the Ensembl gene identifiers from? It looks like they are not present in Ensembl and this is why getBM returns NULL. Are you sure the identifiers are correct and have the right amount of 0's in them? e.g. ENSG00000008311 does return data. Cheers, Steffen Glazko, Galina wrote: > Dear list, > > > > I am trying to do very simple query using BiomaRt: > > library(biomaRt) > > ensembl=useMart('ensembl') > > ensembl= useDataset('hsapiens_gene_ensembl',mart=ensembl) > > .. > > >> ids<-c("ENSG000000282326","ENSG000000008311","ENSG000000323881") >> > > >> hsdef <- >> > getBM(attributes=c("ensembl_gene_id","ensembl_transcript_id","hgnc_s ymbo > l"),filters="ensembl_gene_id",values=ids, mart=ensembl) > > >> hsdef >> > > NULL > > > > Please let me know, am I doing something wrong? > > > > Best regards > > Galina > > PS > > >> sessionInfo() >> > > R version 2.5.1 (2007-06-27) > > i386-pc-mingw32 > > > > > [[alternative HTML version deleted]] > > _______________________________________________ > Bioconductor mailing list > Bioconductor at stat.math.ethz.ch > https://stat.ethz.ch/mailman/listinfo/bioconductor > Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor > >
ADD COMMENT
0
Entering edit mode
@james-w-macdonald-5106
Last seen 3 hours ago
United States
Hi Galina, Glazko, Galina wrote: > Dear list, > > > > I am trying to do very simple query using BiomaRt: > > library(biomaRt) > > ensembl=useMart('ensembl') > > ensembl= useDataset('hsapiens_gene_ensembl',mart=ensembl) > > .. > >> ids<-c("ENSG000000282326","ENSG000000008311","ENSG000000323881") > >> hsdef <- > getBM(attributes=c("ensembl_gene_id","ensembl_transcript_id","hgnc_s ymbo > l"),filters="ensembl_gene_id",values=ids, mart=ensembl) > >> hsdef > > NULL > > > Please let me know, am I doing something wrong? The only thing you did wrong was to input three Ensembl IDs that don't exist - try entering them individually on Ensembl's website, and you will see that you get no returned value. Best, Jim > > > > Best regards > > Galina > > PS > >> sessionInfo() > > R version 2.5.1 (2007-06-27) > > i386-pc-mingw32 > > > > > [[alternative HTML version deleted]] > > _______________________________________________ > Bioconductor mailing list > Bioconductor at stat.math.ethz.ch > https://stat.ethz.ch/mailman/listinfo/bioconductor > Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor -- James W. MacDonald, M.S. Biostatistician Affymetrix and cDNA Microarray Core University of Michigan Cancer Center 1500 E. Medical Center Drive 7410 CCGC Ann Arbor MI 48109 734-647-5623
ADD COMMENT

Login before adding your answer.

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