camera.EdgeResult: Run CAMERA method using EdgeResult

In bedapub/ribiosGSEA: Gene-Set Enrichment Analysis Tools in Ribios

Run CAMERA method using EdgeResult

Description

Run CAMERA method using EdgeResult

Usage

## S3 method for class 'EdgeResult'
camera(y, gmtList, doParallel = FALSE, ...)

Arguments

y	A EdgeResult object
gmtList	Gene set collections, for example read by readGmt, with namespace.
doParallel	Logical, whether parallel::mclapply should be used. Since at the current setting it makes a job running forever, use TRUE only if you are debugging the code.
...	Not used Note that the EdgeResult object must have a column 'GeneSymbol' in its fData.

Value

A data.frame containing CAMERA results.

Examples

exMat <- matrix(rpois(120, 10), nrow=20, ncol=6)
exGroups <- gl(2,3, labels=c("Group1", "Group2"))
exDesign <- model.matrix(~0+exGroups)
colnames(exDesign) <- levels(exGroups)
exContrast <- matrix(c(-1,1,1,-1), ncol=2, dimnames=list(c("Group1", "Group2"), 
  c("Group2.vs.Group1", "Group1.vs.Group2")))
exDescon <- DesignContrast(exDesign, exContrast, groups=exGroups)
exFdata <- data.frame(GeneSymbol=sprintf("GeneSymbol%d", 1:nrow(exMat)))
exPdata <- data.frame(Name=sprintf("Sample%d", 1:ncol(exMat)),
                     Group=exGroups)
exDgeList <- DGEList(exMat, genes=exFdata, samples=exPdata)
exDgeList <- edgeR::estimateDisp(exDgeList, exDesign)
exEdgeObject <- EdgeObject(exDgeList, exDescon)
exEdgeRes <- ribiosNGS::dgeWithEdgeR(exEdgeObject)
exGmt <- BioQC::GmtList(list(GeneSet1=sprintf("GeneSymbol%d", 1:5),
  GeneSet2=sprintf("GeneSymbol%d", 6:10)))
  
exCameraRes <- camera(exEdgeRes, exGmt)

bedapub/ribiosGSEA documentation built on March 30, 2023, 3:26 p.m.