[BioC] Romer and symbols2indices query

Loren Engrav engrav at u.washington.edu
Wed May 5 17:15:04 CEST 2010


Bingo, thank you
 and romer ran

These missing little tidbits can be brutal

> From: Martin Morgan <mtmorgan at fhcrc.org>
> Date: Wed, 05 May 2010 05:11:02 -0700
> To: Loren Engrav <engrav at u.washington.edu>
> Cc: rbioc <bioconductor at stat.math.ethz.ch>
> Subject: Re: [BioC] Romer and symbols2indices query
> 
> On 05/04/2010 07:57 PM, Loren Engrav wrote:
>> Am back
>> 
>> So I have romer and GSEABase running via previous help thank you, but while
>> running I explore GSEABase
>> 
>> And I have a lesser question for interest
>> 
>> In GSEABase I do
>>  gmtObject <- getGMT("c2all.v2.5.symbols.gmt",
>> collectionType=BroadCollection(category="c2"), geneType=SymbolIdentifier())
> 
> or maybe getBroadSets ?
> 
>>  which finishes without error
>> Then
>>  class(gmtObject) is GeneSetCollection
>> 
>> How do I convert gmtObject to a list of gene sets as required in romer when
>> using
> 
>   gmtl <- geneIds(gmtObject)
>   names(gmtl) <- names(gmtObject)
> 
> ?
> 
> Martin
> 
>> X <- symbols2indices (ListOfGeneSets, MYsymbols)
>> 
>> 
>> 
>> From: Vincent Carey <stvjc at channing.harvard.edu>
>> Date: Tue, 4 May 2010 11:40:39 -0400
>> To: Loren Engrav <engrav at u.washington.edu>
>> Cc: rbioc <bioconductor at stat.math.ethz.ch>
>> Subject: Re: [BioC] Romer and symbols2indices query
>> 
>> Very briefly, the GSEABase package has relevant utilities for gmt file
>> import/export and may be worth considering for these tasks.
>> 
>> On Tue, May 4, 2010 at 10:43 AM, Loren Engrav <engrav at u.washington.edu>
>> wrote:
>>> Thank you, got it
>>> 
>>> Downloading rdata objects saves reading them into an rdata object, cool
>>> 
>>> But for interest, in R/GSA there is
>>> GSA.read.gmt(filename.gmt) to read in a .gmt file
>>> 
>>> Does limma or romer have an equivalent function?
>>> 
>>> 
>>>> From: Matthew Ritchie <mritchie at wehi.EDU.AU>
>>>> Date: Tue, 4 May 2010 14:44:23 +1000 (EST)
>>>> To: Loren Engrav <engrav at u.washington.edu>
>>>> Cc: rbioc <bioconductor at stat.math.ethz.ch>
>>>> Subject: Re: [BioC] Romer and symbols2indices query
>>>> 
>>>> Dear Loren,
>>>> 
>>>> You can find rdata objects of the Broad's MSigDB gene sets at
>>>> 
>>>> http://bioinf.wehi.edu.au/software/MSigDB/index.html
>>>> 
>>>> You are right, the 'symbols' argument in the function symbols2indicies()
>>>> are the gene symbols corresponding to the probes from your microarray
>>>> data.
>>>> 
>>>> For example, to use the human C2 collection, download the rdata file, then
>>>> run the following.
>>>> 
>>>> load("human_c2.rdata")
>>>> c2 = symbols2indices(Hs.gmtl.c2, symbols)
>>>> 
>>>> (this assumes 'symbols' is a vector containing the gene symbols from your
>>>> array data)
>>>> 
>>>> Best wishes,
>>>> 
>>>> Matt
>>>> 
>>>>> Have done GSEA and GSA for set enrichment and am setting out to try romer
>>>>> and have probably "simple" question
>>>>> 
>>>>> To get the Broad set into a list of indices there is
>>>>> symbols2indices(gmtl.official, symbols) but
>>>>> 
>>>>> 1)how do I get the Broad set into gmtl.official? And
>>>>> 2)is symbols a vector of MY probe sets of interest?
>>>>> 
>>>>> I checked gmane and found only one comment about romer
>>>>> Also checked limma reference pdf
>>>>> 
>>>>> Thank you
>>>> 
>>>> 
>>>> ______________________________________________________________________
>>>> The information in this email is confidential and inte...{{dropped:7}}
>>> 
>>> _______________________________________________
>>> 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
> 
> 
> -- 
> Martin Morgan
> Computational Biology / Fred Hutchinson Cancer Research Center
> 1100 Fairview Ave. N.
> PO Box 19024 Seattle, WA 98109
> 
> Location: Arnold Building M1 B861
> Phone: (206) 667-2793



More information about the Bioconductor mailing list