gs1 <- GeneSet(setName="set1", setIdentifier="101")
gs2 <- GeneSet(setName="set2", setIdentifier="102")
## construct from indivdiual elements...
gsc <- GeneSetCollection(gs1, gs2)
## or from a list
gsc <- GeneSetCollection(list(gs1, gs2))
## set names must be unique
try(GeneSetCollection(gs1, gs1))
data(sample.ExpressionSet)
gsc <- GeneSetCollection(sample.ExpressionSet[200:250],
setType = GOCollection())
## Not run:
# ## from KEGG identifiers, for example
# library(KEGG.db)
# lst <- head(as.list(KEGGEXTID2PATHID))
# gsc <- GeneSetCollection(mapply(function(geneIds, keggId) {
# GeneSet(geneIds, geneIdType=EntrezIdentifier(),
# collectionType=KEGGCollection(keggId),
# setName=keggId)
# }, lst, names(lst)))
# ## End(Not run)
Run the code above in your browser using DataLab