Problem with invalid GO term in HyperGResult object - NOT solved
1
0
Entering edit mode
@martin-morgan-1513
Last seen 13 days ago
United States
James MacDonald <jmacdon at="" med.umich.edu=""> writes: > Hi Jenny, > > You want multiassign() in Biobase. > > Something like this should do the trick: > > tmp <- eapply(ath1121501MULTIHIT, > function(ids) if(length(ids) >1) "multiple" else ids) > ath112150ACCNUM <- new.env(hash=TRUE, parent=.GlobalEnv) > multiassign(tmp, envir=ath112150ACCNUM) or l2e ('list to environment'; also from Biobase). > ath112150ACCNUM <- l2e(tmp) Probably in Jim's use it makes sense to have parent=emptyenv(), otherwise you can 'get' symbols that are defined in the global environment rather than the environment you're looking in. > x <- 1 > safer <- new.env(parent=emptyenv()) > unsafe <- new.env(parent=.GlobalEnv) > get("x", safer) Error in get("x", safer) : variable "x" was not found > get("x", unsafe) # look in 'parent', i.e., .GlobalEnv [1] 1 > get("l2e", unsafe) # parents of .GlobalEnv include the search() path function (vals, envir) { if (missing(envir)) envir <- new.env(hash = TRUE, parent = emptyenv(), size = .new_env_size(length(vals))) .Call("listToEnv", vals, envir, PACKAGE = "Biobase") } <environment: namespace:biobase=""> Martin > Best, > > Jim > > > > Jenny Drnevich wrote: >> Hi Herve, >> >> Thanks for looking into this. Unfortunately, your fix didn't quite work: >> >> >>> You can compute the ath1121501ACCNUM yourself with: >>> >>>> ath1121501ACCNUM <- eapply(ath1121501MULTIHIT, >>> function(ids) if (length(ids) > 1) >>> "multiple" else ids) >> >> This changes ath1121501ACCNUM from an environment to a list, and to >> use hyperGTest(), it needs to be an environment. I know very little >> about environments, and the simple as.environment() didn't work, so I >> don't know how to fix it. >> >> Thanks, >> Jenny >> >> >> >>> Let me know if you find any other problem. >>> >>> Thanks! >>> H. >>> >>> >>>> Thanks, >>>> Jenny >>>> >>>>> sessionInfo() >>>> R version 2.6.0 Under development (unstable) (2007-08-28 r42679) >>>> i386-pc-mingw32 >>>> >>>> locale: >>>> LC_COLLATE=English_United States.1252;LC_CTYPE=English_United >>>> States.1252;LC_MONETARY=English_United >>>> States.1252;LC_NUMERIC=C;LC_TIME=English_United States.1252 >>>> >>>> attached base packages: >>>> [1] splines tools stats graphics grDevices utils datasets >>>> [8] methods base >>>> >>>> other attached packages: >>>> [1] limma_2.11.13 ath1121501_1.99.10 GOstats_2.3.17 >>>> [4] Category_2.3.39 genefilter_1.15.11 survival_2.32 >>>> [7] RBGL_1.13.6 annotate_1.15.11 xtable_1.5-1 >>>> [10] GO.db_1.99.3 AnnotationDbi_0.99.1 RSQLite_0.6-2 >>>> [13] DBI_0.2-3 graph_1.15.20 Biobase_1.15.34 >>>> >>>> loaded via a namespace (and not attached): >>>> [1] cluster_1.11.7 simpleaffy_2.13.01 >>>> >>>> >>>> >>>> >>>> >>>> Jenny Drnevich, Ph.D. >>>> >>>> Functional Genomics Bioinformatics Specialist >>>> W.M. Keck Center for Comparative and Functional Genomics >>>> Roy J. Carver Biotechnology Center >>>> University of Illinois, Urbana-Champaign >>>> >>>> 330 ERML >>>> 1201 W. Gregory Dr. >>>> Urbana, IL 61801 >>>> USA >>>> >>>> ph: 217-244-7355 >>>> fax: 217-265-5066 >>>> e-mail: drnevich at uiuc.edu >> >> _______________________________________________ >> 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, MS > Biostatistician > UMCCC cDNA and Affymetrix Core > University of Michigan > 1500 E Medical Center Drive > 7410 CCGC > Ann Arbor MI 48109 > 734-647-5623 > > _______________________________________________ > 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 -- Martin Morgan Bioconductor / Computational Biology http://bioconductor.org
GO Biobase GO Biobase • 862 views
ADD COMMENT
0
Entering edit mode
@james-w-macdonald-5106
Last seen 13 hours ago
United States
Martin Morgan wrote: > Probably in Jim's use it makes sense to have parent=emptyenv(), > otherwise you can 'get' symbols that are defined in the global > environment rather than the environment you're looking in. Yup. My bad. I was thinking .GlobalEnv wasn't right, but I couldn't remember what parent=NULL used to do. Thanks for the heads-up. Best, Jim -- James W. MacDonald, MS Biostatistician UMCCC cDNA and Affymetrix Core University of Michigan 1500 E Medical Center Drive 7410 CCGC Ann Arbor MI 48109 734-647-5623
ADD COMMENT
0
Entering edit mode
Thanks Herve, Jim & Martin! I've got it working now... Cheers, Jenny At 11:43 AM 9/24/2007, James MacDonald wrote: >Martin Morgan wrote: > > Probably in Jim's use it makes sense to have parent=emptyenv(), > > otherwise you can 'get' symbols that are defined in the global > > environment rather than the environment you're looking in. > >Yup. My bad. I was thinking .GlobalEnv wasn't right, but I couldn't >remember what parent=NULL used to do. Thanks for the heads-up. > >Best, > >Jim > > > >-- >James W. MacDonald, MS >Biostatistician >UMCCC cDNA and Affymetrix Core >University of Michigan >1500 E Medical Center Drive >7410 CCGC >Ann Arbor MI 48109 >734-647-5623 > >_______________________________________________ >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 Jenny Drnevich, Ph.D. Functional Genomics Bioinformatics Specialist W.M. Keck Center for Comparative and Functional Genomics Roy J. Carver Biotechnology Center University of Illinois, Urbana-Champaign 330 ERML 1201 W. Gregory Dr. Urbana, IL 61801 USA ph: 217-244-7355 fax: 217-265-5066 e-mail: drnevich at uiuc.edu
ADD REPLY
0
Entering edit mode
Hi Jenny, We've fixed the problem you reported earlier with the ACCNUM map for arabidopsis. The new ag and ath1121501 packages will be online soon (we will post again when they are). Sorry that my temporary workaround was producing a list instead of an environment and thanks to Jim and Martin for providing a quick solution to the problem. Let me add a couple of comments related to this issue: - If you use the ath1121501.db package (the SQLite-based versions of ath1121501), then you'll be able to use the fixed ACCNUM map right now (don't forget to use the latest version of AnnotationDbi though). The most visible difference between our new SQLite-based metadata packages and our classic environment- based metadata packages is that the maps in the formers are S4 objects (AnnDbBimap objects to be precise) instead of environments. Don't be afraid by this: we've put a lot of effort to make those objects look and feel like environments i.e. you can still use length, ls, mget, eapply, get, [[, etc... on them. hyperGTest() should work as expected with ath1121501.db. Please let us know if this is not the case. - In the ath1121501.db package (and soon in the fixed version of the ath1121501 package), the ACCNUM and MULTIHIT maps are complementary: when a probe id is mapped to a single AGI locus ID, then ath1121501MULTIHIT is set to NA: > library(ath1121501.db) > ath1121501ACCNUM[["246779_at"]] [1] "AT5G27520" > ath1121501MULTIHIT[["246779_at"]] [1] NA and when it's mapped to multiple AGI locus IDs, then ath1121501ACCNUM is set to "multiple": > ath1121501ACCNUM[["264451_s_at"]] [1] "multiple" > ath1121501MULTIHIT[["264451_s_at"]] [1] "AT1G10260" "AT1G11265" "AT1G21945" "AT1G48710" "AT1G58140" "AT3G32415" [7] "AT3G59720" "AT3G61330" "AT5G13205" You can easily create a new map that contains the "merged" data from the 2 maps with: > ath1121501ALLACCNUMS <- as(ath1121501ACCNUM, "AnnDbBimap") > ath1121501ALLACCNUMS[["246779_at"]] [1] "AT5G27520" > ath1121501ALLACCNUMS[["264451_s_at"]] [1] "AT1G10260" "AT1G11265" "AT1G21945" "AT1G48710" "AT1G58140" "AT3G32415" [7] "AT3G59720" "AT3G61330" "AT5G13205" This works with the SQLite-based annotations only (and it is very fast). We will add some examples in the man page for those 2 maps to illustrate this "merging" operation. Cheers, H. Jenny Drnevich wrote: > Thanks Herve, Jim & Martin! > > I've got it working now... > > Cheers, > Jenny > > At 11:43 AM 9/24/2007, James MacDonald wrote: >> Martin Morgan wrote: >>> Probably in Jim's use it makes sense to have parent=emptyenv(), >>> otherwise you can 'get' symbols that are defined in the global >>> environment rather than the environment you're looking in. >> Yup. My bad. I was thinking .GlobalEnv wasn't right, but I couldn't >> remember what parent=NULL used to do. Thanks for the heads-up. >> >> Best, >> >> Jim >> >> >> >> -- >> James W. MacDonald, MS >> Biostatistician >> UMCCC cDNA and Affymetrix Core >> University of Michigan >> 1500 E Medical Center Drive >> 7410 CCGC >> Ann Arbor MI 48109 >> 734-647-5623 >> >> _______________________________________________ >> 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 > > Jenny Drnevich, Ph.D. > > Functional Genomics Bioinformatics Specialist > W.M. Keck Center for Comparative and Functional Genomics > Roy J. Carver Biotechnology Center > University of Illinois, Urbana-Champaign > > 330 ERML > 1201 W. Gregory Dr. > Urbana, IL 61801 > USA > > ph: 217-244-7355 > fax: 217-265-5066 > e-mail: drnevich at uiuc.edu > > _______________________________________________ > 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 REPLY
0
Entering edit mode
Source files should be up for these packages now. Windows zip files should follow very soon. Marc Quoting Herve Pages <hpages at="" fhcrc.org="">: > Hi Jenny, > > We've fixed the problem you reported earlier with the ACCNUM map for > arabidopsis. > The new ag and ath1121501 packages will be online soon (we will post again > when > they are). > > Sorry that my temporary workaround was producing a list instead of an > environment > and thanks to Jim and Martin for providing a quick solution to the problem. > > Let me add a couple of comments related to this issue: > > - If you use the ath1121501.db package (the SQLite-based versions of > ath1121501), > then you'll be able to use the fixed ACCNUM map right now (don't forget to > use > the latest version of AnnotationDbi though). The most visible difference > between > our new SQLite-based metadata packages and our classic environment-based > metadata > packages is that the maps in the formers are S4 objects (AnnDbBimap > objects to > be precise) instead of environments. > Don't be afraid by this: we've put a lot of effort to make those objects > look > and feel like environments i.e. you can still use length, ls, mget, > eapply, get, > [[, etc... on them. > hyperGTest() should work as expected with ath1121501.db. Please let us > know if > this is not the case. > > - In the ath1121501.db package (and soon in the fixed version of the > ath1121501 > package), the ACCNUM and MULTIHIT maps are complementary: when a probe id > is > mapped to a single AGI locus ID, then ath1121501MULTIHIT is set to NA: > > > library(ath1121501.db) > > ath1121501ACCNUM[["246779_at"]] > [1] "AT5G27520" > > ath1121501MULTIHIT[["246779_at"]] > [1] NA > > and when it's mapped to multiple AGI locus IDs, then ath1121501ACCNUM is > set > to "multiple": > > > ath1121501ACCNUM[["264451_s_at"]] > [1] "multiple" > > ath1121501MULTIHIT[["264451_s_at"]] > [1] "AT1G10260" "AT1G11265" "AT1G21945" "AT1G48710" "AT1G58140" > "AT3G32415" > [7] "AT3G59720" "AT3G61330" "AT5G13205" > > You can easily create a new map that contains the "merged" data from the > 2 maps with: > > > ath1121501ALLACCNUMS <- as(ath1121501ACCNUM, "AnnDbBimap") > > ath1121501ALLACCNUMS[["246779_at"]] > [1] "AT5G27520" > > ath1121501ALLACCNUMS[["264451_s_at"]] > [1] "AT1G10260" "AT1G11265" "AT1G21945" "AT1G48710" "AT1G58140" > "AT3G32415" > [7] "AT3G59720" "AT3G61330" "AT5G13205" > > This works with the SQLite-based annotations only (and it is very fast). > We will add some examples in the man page for those 2 maps to illustrate > this "merging" operation. > > Cheers, > H. > > > Jenny Drnevich wrote: > > Thanks Herve, Jim & Martin! > > > > I've got it working now... > > > > Cheers, > > Jenny > > > > At 11:43 AM 9/24/2007, James MacDonald wrote: > >> Martin Morgan wrote: > >>> Probably in Jim's use it makes sense to have parent=emptyenv(), > >>> otherwise you can 'get' symbols that are defined in the global > >>> environment rather than the environment you're looking in. > >> Yup. My bad. I was thinking .GlobalEnv wasn't right, but I couldn't > >> remember what parent=NULL used to do. Thanks for the heads-up. > >> > >> Best, > >> > >> Jim > >> > >> > >> > >> -- > >> James W. MacDonald, MS > >> Biostatistician > >> UMCCC cDNA and Affymetrix Core > >> University of Michigan > >> 1500 E Medical Center Drive > >> 7410 CCGC > >> Ann Arbor MI 48109 > >> 734-647-5623 > >> > >> _______________________________________________ > >> 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 > > > > Jenny Drnevich, Ph.D. > > > > Functional Genomics Bioinformatics Specialist > > W.M. Keck Center for Comparative and Functional Genomics > > Roy J. Carver Biotechnology Center > > University of Illinois, Urbana-Champaign > > > > 330 ERML > > 1201 W. Gregory Dr. > > Urbana, IL 61801 > > USA > > > > ph: 217-244-7355 > > fax: 217-265-5066 > > e-mail: drnevich at uiuc.edu > > > > _______________________________________________ > > 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 > > > > _______________________________________________ > 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 REPLY

Login before adding your answer.

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