Using GeneSetCollection objects of GSEABase in piano package
1
0
Entering edit mode
@lluis-revilla-sancho
Last seen 39 minutes ago
European Union

I am interested in some function of the piano package. I have read the vignette and the help page, and I assumed that the GSC class needed to use runGSA was the GeneSetCollection provided by the GSEABase package of the project.

However when I read in detail the function and compare with the GeneSetCollection object they aren't the same. I have some GeneSetCollections I was using to perform GSA with GSVA, limma ...

Is there a function or method to convert the GeneSetCollection to the GSC objects needed for piano? The only workaround would be to save in gmt format the GeneSetCollection and then load it again:

toGMT(myGeneSetCollection, "myGSC.gmt")
GSCformat <- loadGSC("myGSC.gmt")

I can understand that some GSA test require the data in other format than GeneSetCollection, but couldn't that transformation be done internally? It would be easier to integrate with existing GSEA packages and tools in Bioconductor.

piano gseabase • 1.2k views
ADD COMMENT
1
Entering edit mode

Hi, thanks for pointing this out. I will aim at including support for the GeneSetCollection class in the next piano update, either directly in runGSA or perhaps by a run through loadGSC.

ADD REPLY
1
Entering edit mode
@lluis-revilla-sancho
Last seen 39 minutes ago
European Union

I found another workaround, to convert GeneSetCollection to GSC

GeneSetCollection2GSC <- function(gsc, addInfo = "none") {
  GSC <- list("gsc" = geneIds(gsc), "addInfo" = addInfo)
  class(GSC) <- "GSC"
  return(GSC)
}
myGSC <- geneSetCollection2GSC(myGeneSetCollection)
myGSC

 

ADD COMMENT

Login before adding your answer.

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