exonmap, xmap and xmapbridge
1
0
Entering edit mode
Paul Hammer ▴ 220
@paul-hammer-2635
Last seen 10.3 years ago
An embedded and charset-unspecified text was scrubbed... Name: not available Url: https://stat.ethz.ch/pipermail/bioconductor/attachments/20080205/ 01ecac15/attachment.pl
• 1.1k views
ADD COMMENT
0
Entering edit mode
@vincent-j-carey-jr-4
Last seen 3 months ago
United States
> dear all, > > I am trying to plot different gene/transcript/exon expressions from exon > affymetrix chips by using the function 'xmapGene' or 'xmapBridge' > (package exonmap). but when I start these functions I get an error > message, that I have to use <s4 object="" of="" class="" "expressionset"=""> for > value 'eset'. I used this code: > > > raw.data = read.exon() > > raw.data at cdfName = "exon.pmcdf" > > affy.rma = rma(raw.data) > Background correcting > Normalizing > Calculating Expression > > xmapDatabase("Human") > done. > > xmapGene("ENSG00000128655", *eset=affy.rma*, > use.bridge=TRUE,gps=list(A=1:3,B=4:6),type="fc",overwrite=FALSE) > Error in xmapGene("ENSG00000128655", *eset = affy.rma*, use.bridge = > TRUE, : > unused argument(s) (*eset = <s4 object="" of="" class="" "expressionset"="">*, > use.bridge = TRUE, gps = list(A = 1:3, B = 4:8), type = "fc", overwrite > = FALSE) > > I thought that the function rma() results in an S4 object of the class > "exprSet", isn't it??? and is "exprSet"the same like "ExpressionSet"? > Is there any chance to get an S4 object of class "ExpressionSet??? Please provide a sessionInfo() result to tell us the versions of R and the various packages you are using. "exprSet" is the name of a legacy class that is not to be used in current applications. There may be software on your system that generates "exprSet" instances but it should be updated. "ExpressionSet" is the name of the class that replaces "exprSet". In the current version of exonmap, xmapGene has a single parameter "x". So your call with numerous parameters shown above (hard to decipher) may not be valid. The information transmitted in this electronic communica...{{dropped:9}}
ADD COMMENT
0
Entering edit mode
Hi Paul, I think you're confusing two functions: Have a look at the man page for xmapGene()... It (and the related functions xmapTranscript(), xmapExon() and xmaProbeset()) simply opens up a browser window on xmap, located at the gene, (or transcript, etc...) specified. It doesn't do anything with the expression data and simply takes one parameter, the gene id. The error you're getting is R complaining that you're calling xmapGene with lots of unused arguments (including eset). The function xmapBridge() is what you want to use to write out expression data for the XMapBridge. It doesn't expect the parameter 'use.bridge', though... This should work, however!: xmapBridge(<my gene="" id="">,eset=affy. rma,gps=list(A=1:3,B=4:6), type="fc",overwrite=FALSE) Crispin On 5/2/08 18:01, "Vincent Carey 525-2265" <stvjc at="" channing.harvard.edu=""> wrote: > >> dear all, >> >> I am trying to plot different gene/transcript/exon expressions from exon >> affymetrix chips by using the function 'xmapGene' or 'xmapBridge' >> (package exonmap). but when I start these functions I get an error >> message, that I have to use <s4 object="" of="" class="" "expressionset"=""> for >> value 'eset'. I used this code: >> >>> raw.data = read.exon() >>> raw.data at cdfName = "exon.pmcdf" >>> affy.rma = rma(raw.data) >> Background correcting >> Normalizing >> Calculating Expression >>> xmapDatabase("Human") >> done. >>> xmapGene("ENSG00000128655", *eset=affy.rma*, >> use.bridge=TRUE,gps=list(A=1:3,B=4:6),type="fc",overwrite=FALSE) >> Error in xmapGene("ENSG00000128655", *eset = affy.rma*, use.bridge = >> TRUE, : >> unused argument(s) (*eset = <s4 object="" of="" class="" "expressionset"="">*, >> use.bridge = TRUE, gps = list(A = 1:3, B = 4:8), type = "fc", overwrite >> = FALSE) >> >> I thought that the function rma() results in an S4 object of the class >> "exprSet", isn't it??? and is "exprSet"the same like "ExpressionSet"? >> Is there any chance to get an S4 object of class "ExpressionSet??? > > Please provide a sessionInfo() result to tell us the versions > of R and the various packages you are using. > > "exprSet" is the name of a legacy class that is not to be used in > current applications. There may be software on your system that > generates "exprSet" instances but it should be updated. > > "ExpressionSet" is the name of the class that replaces "exprSet". > > In the current version of exonmap, xmapGene has a single parameter "x". > So your call with numerous parameters shown above (hard to decipher) > may not be valid. > > The information transmitted in this electronic communica...{{dropped:9}} > > _______________________________________________ > 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 -------------------------------------------------------- This email is confidential and intended solely for the u...{{dropped:13}}
ADD REPLY
0
Entering edit mode
dear all, the function xmapBridge() (package exonmap) worked fine until i did a reboot. now always comes this mysql error message: > xmapBridge("ENSG00000115252",eset=rma.affy,gps=list(A=1:3,B=4:8),type= "fc") Fehler in mysqlExecStatement(conn, statement, ...) : RS-DBI driver: (could not run statement: Unknown database 'homo_sapiens_core_47_36i') when i login to mysql the database with all tables is there. here my session info: > sessionInfo() R version 2.6.0 (2007-10-03) x86_64-unknown-linux-gnu locale: LC_CTYPE=de_DE.UTF-8;LC_NUMERIC=C;LC_TIME=de_DE.UTF-8;LC_COLLATE=de_DE .UTF-8;LC_MONETARY=de_DE.UTF-8;LC_MESSAGES=de_DE.UTF-8;LC_PAPER=de_DE. UTF-8;LC_NAME=C;LC_ADDRESS=C;LC_TELEPHONE=C;LC_MEASUREMENT=de_DE.UTF-8 ;LC_IDENTIFICATION=C attached base packages: [1] splines tools stats graphics grDevices utils datasets [8] methods base other attached packages: [1] exon.pmcdf_1.1 exonmap_1.4.1 plier_1.8.0 [4] RMySQL_0.6-0 DBI_0.2-4 RColorBrewer_1.0-2 [7] simpleaffy_2.14.05 gcrma_2.10.0 matchprobes_1.10.0 [10] genefilter_1.16.0 survival_2.32 affy_1.16.0 [13] preprocessCore_1.0.0 affyio_1.6.1 Biobase_1.16.2 loaded via a namespace (and not attached): [1] annotate_1.16.1 AnnotationDbi_1.0.6 rcompgen_0.1-15 [4] RSQLite_0.6-4 have someone an idea??? thanks paul
ADD REPLY
0
Entering edit mode
Hi Paul, There is a problem with the connection parameters, i.e. the ensembldb port which has been switched from 3306 -> 5306. While the "old" (3306) mysql instance still contains homo_sapiens_core_47_36i the exonmap package already picks up the new 5306 instance which doesn't contain your desired database (see below). Stefan [graef at data graef]$ echo "show databases like 'homo_sapiens_core_4%'" | mysql -hensembldb.ensembl.org -P3306 -uanonymous Database (homo_sapiens_core_4%) homo_sapiens_core_40_36b homo_sapiens_core_41_36c homo_sapiens_core_42_36d homo_sapiens_core_43_36e homo_sapiens_core_44_36f homo_sapiens_core_45_36g homo_sapiens_core_46_36h homo_sapiens_core_47_36i [graef at data graef]$ echo "show databases like 'homo_sapiens_core_4%'" | mysql -hensembldb.ensembl.org -P5306 -uanonymous Database (homo_sapiens_core_4%) homo_sapiens_core_48_36j Paul Hammer wrote: > dear all, > > the function xmapBridge() (package exonmap) worked fine until i did a > reboot. now always comes this mysql error message: > > > > xmapBridge("ENSG00000115252",eset=rma.affy,gps=list(A=1:3,B=4:8),typ e="fc") > Fehler in mysqlExecStatement(conn, statement, ...) : RS-DBI driver: > (could not run statement: Unknown database 'homo_sapiens_core_47_36i') > > when i login to mysql the database with all tables is there. > > here my session info: > > sessionInfo() > R version 2.6.0 (2007-10-03) > x86_64-unknown-linux-gnu > > locale: > LC_CTYPE=de_DE.UTF-8;LC_NUMERIC=C;LC_TIME=de_DE.UTF-8;LC_COLLATE=de_ DE.UTF-8;LC_MONETARY=de_DE.UTF-8;LC_MESSAGES=de_DE.UTF-8;LC_PAPER=de_D E.UTF-8;LC_NAME=C;LC_ADDRESS=C;LC_TELEPHONE=C;LC_MEASUREMENT=de_DE.UTF -8;LC_IDENTIFICATION=C > > > attached base packages: > [1] splines tools stats graphics grDevices utils datasets > [8] methods base > > other attached packages: > [1] exon.pmcdf_1.1 exonmap_1.4.1 plier_1.8.0 > [4] RMySQL_0.6-0 DBI_0.2-4 RColorBrewer_1.0-2 > [7] simpleaffy_2.14.05 gcrma_2.10.0 matchprobes_1.10.0 > [10] genefilter_1.16.0 survival_2.32 affy_1.16.0 > [13] preprocessCore_1.0.0 affyio_1.6.1 Biobase_1.16.2 > > loaded via a namespace (and not attached): > [1] annotate_1.16.1 AnnotationDbi_1.0.6 rcompgen_0.1-15 > [4] RSQLite_0.6-4 > > have someone an idea??? > > thanks > paul > > _______________________________________________ > 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 -- Dr. Stefan Gr?f EMBL-European Bioinformatics Institute Ensembl - Functional Genomics Research Wellcome Trust Genome Campus, Hinxton, Cambridge CB10 1SD, UK phone: +44 (0)1223 494685 fax: +44 (0)1223 494468
ADD REPLY
0
Entering edit mode
Paul Hammer schrieb: > dear all, > > the function xmapBridge() (package exonmap) worked fine until i did a > reboot. now always comes this mysql error message: > > > > xmapBridge("ENSG00000115252",eset=rma.affy,gps=list(A=1:3,B=4:8),typ e="fc") > > Fehler in mysqlExecStatement(conn, statement, ...) : RS-DBI driver: > (could not run statement: Unknown database 'homo_sapiens_core_47_36i') > > when i login to mysql the database with all tables is there. > > here my session info: > > sessionInfo() > R version 2.6.0 (2007-10-03) > x86_64-unknown-linux-gnu > > locale: > LC_CTYPE=de_DE.UTF-8;LC_NUMERIC=C;LC_TIME=de_DE.UTF-8;LC_COLLATE=de_ DE.UTF-8;LC_MONETARY=de_DE.UTF-8;LC_MESSAGES=de_DE.UTF-8;LC_PAPER=de_D E.UTF-8;LC_NAME=C;LC_ADDRESS=C;LC_TELEPHONE=C;LC_MEASUREMENT=de_DE.UTF -8;LC_IDENTIFICATION=C > > > attached base packages: > [1] splines tools stats graphics grDevices utils datasets > [8] methods base > > other attached packages: > [1] exon.pmcdf_1.1 exonmap_1.4.1 plier_1.8.0 > [4] RMySQL_0.6-0 DBI_0.2-4 RColorBrewer_1.0-2 > [7] simpleaffy_2.14.05 gcrma_2.10.0 matchprobes_1.10.0 > [10] genefilter_1.16.0 survival_2.32 affy_1.16.0 > [13] preprocessCore_1.0.0 affyio_1.6.1 Biobase_1.16.2 > > loaded via a namespace (and not attached): > [1] annotate_1.16.1 AnnotationDbi_1.0.6 rcompgen_0.1-15 > [4] RSQLite_0.6-4 > > have someone an idea??? > > thanks > paul > ok i found it by myself. i just restarted the mysql server in the right folder (/etc/init.d/mysqld start)... goodnight from berlin paul
ADD REPLY

Login before adding your answer.

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