data(gCMAPData)
gene.set.collection <- induceCMAPCollection(gCMAPData, "z", higher=2, lower=-2)
## comparison with a single user-provided profile of z-scores
profile <- assayDataElement(gCMAPData, "z")[,1]
gsealm_jg_score(profile, gene.set.collection)
## comparison with of multiple profiles of z-scores to the CMAPCollection
res <- gsealm_jg_score(assayDataElement(gCMAPData, "z"), gene.set.collection)
## first CMAPResult object
res[[1]]
## adjusted p-values from all CMAPResult objects
sapply(res, padj)
## inverted query: CMAPCollection is compared to z-score profiles
gsealm_jg_score(gene.set.collection, assayDataElement(gCMAPData, "z"))[[1]]
Run the code above in your browser using DataLab