Annotation databases--unable to get annotations
1
0
Entering edit mode
Guest User ★ 12k
@guest-user-4897
Last seen 7.4 years ago
Has anyone used the mogene10stprobeset.db annotation package? I have a list of significant probeset ids that I obtained from running an analysis in the limma package. When I try to use mget() to retrieve the corresponding gene names only a couple of probesets get annotated. I've tried with several different gene lists. In a gene list of about 300 genes, only 2 get annotated at the most. The rest of the probesets get a "NA" for the annotation. Here is my R code: ## load library > library("mogene10stprobeset.db") ## map probe ids to gene names... #extract ids from signifcant gene list ids=sig.gene.list[["ID"]] #if not in character form, change to character ids = as.character(ids) #retrieve gene names from annotation package for ids in signficant gene list mget(ids, mogene10stprobesetGENENAME, ifnotfound=NA) I've tried ENTREZID, SYMBOL, and GO, all of them return annotations just for a few probesets. I'm using R version R 2.15.1 for Macs. Any help would be appreciated!! Thanks! -- output of sessionInfo(): > sessionInfo() R version 2.15.1 (2012-06-22) Platform: x86_64-apple-darwin9.8.0/x86_64 (64-bit) locale: [1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8 attached base packages: [1] stats graphics grDevices utils datasets methods base other attached packages: [1] hgu95av2.db_2.7.1 org.Hs.eg.db_2.8.0 limma_3.14.0 mogene10stv1cdf_2.11.0 affy_1.36.0 mogene10stprobeset.db_8.0.1 org.Mm.eg.db_2.8.0 [8] RSQLite_0.11.2 DBI_0.2-5 AnnotationDbi_1.20.0 Biobase_2.18.0 BiocGenerics_0.4.0 BiocInstaller_1.8.2 loaded via a namespace (and not attached): [1] affyio_1.26.0 IRanges_1.16.2 parallel_2.15.1 preprocessCore_1.20.0 stats4_2.15.1 tools_2.15.1 zlibbioc_1.4.0 > -- Sent via the guest posting facility at bioconductor.org.
Annotation GO hgu95av2 probe limma Annotation GO hgu95av2 probe limma • 727 views
ADD COMMENT
0
Entering edit mode
@james-w-macdonald-5106
Last seen 1 hour ago
United States
Hi Manjula, On 10/10/2012 2:31 PM, Manjula Kasoji [guest] wrote: > Has anyone used the mogene10stprobeset.db annotation package? > > I have a list of significant probeset ids that I obtained from running an analysis in the limma package. When I try to use mget() to retrieve the corresponding gene names only a couple of probesets get annotated. I've tried with several different gene lists. In a gene list of about 300 genes, only 2 get annotated at the most. The rest of the probesets get a "NA" for the annotation. > > Here is my R code: That's not all your code, just some. You might also note that you used the affy package to summarize, in which case you have summarized at the transcript level, not the probeset level. So you are using the wrong annotation package. Try the mogene10sttranscriptcluster.db package and see how that works. Best, Jim > > ## load library >> library("mogene10stprobeset.db") > ## map probe ids to gene names... > > #extract ids from signifcant gene list > ids=sig.gene.list[["ID"]] > > #if not in character form, change to character > ids = as.character(ids) > > #retrieve gene names from annotation package for ids in signficant gene list > mget(ids, mogene10stprobesetGENENAME, ifnotfound=NA) > > > I've tried ENTREZID, SYMBOL, and GO, all of them return annotations just for a few probesets. > > I'm using R version R 2.15.1 for Macs. > > Any help would be appreciated!! > > Thanks! > > -- output of sessionInfo(): > >> sessionInfo() > R version 2.15.1 (2012-06-22) > Platform: x86_64-apple-darwin9.8.0/x86_64 (64-bit) > > locale: > [1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8 > > attached base packages: > [1] stats graphics grDevices utils datasets methods base > > other attached packages: > [1] hgu95av2.db_2.7.1 org.Hs.eg.db_2.8.0 limma_3.14.0 mogene10stv1cdf_2.11.0 affy_1.36.0 mogene10stprobeset.db_8.0.1 org.Mm.eg.db_2.8.0 > [8] RSQLite_0.11.2 DBI_0.2-5 AnnotationDbi_1.20.0 Biobase_2.18.0 BiocGenerics_0.4.0 BiocInstaller_1.8.2 > > loaded via a namespace (and not attached): > [1] affyio_1.26.0 IRanges_1.16.2 parallel_2.15.1 preprocessCore_1.20.0 stats4_2.15.1 tools_2.15.1 zlibbioc_1.4.0 > -- > Sent via the guest posting facility at bioconductor.org. > > _______________________________________________ > Bioconductor mailing list > Bioconductor at r-project.org > 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 University of Washington Environmental and Occupational Health Sciences 4225 Roosevelt Way NE, # 100 Seattle WA 98105-6099
ADD COMMENT

Login before adding your answer.

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