make_gsea_overlap_igraph: this function takes a bioConductor EList (eg....

Description Usage Arguments Value Examples

View source: R/make_gsea_overlap_igraph.R

Description

this function takes a bioConductor EList (eg. voom-transformed RNASeq counts) GSEA output matrix from limma's camera() and a user-specified geneSet and create an igraph graph where the nodes are gene lists and the edge weight is how many expressed genes are common between the gene lists Phu T. Van, FHCRC 2017, w/ substantial help and input from C.Murie & V.Voillet

Usage

1
make_gsea_overlap_igraph(expressionList, cameraMat, geneSets, verbose = FALSE)

Arguments

expressionList

a matrix of gene expression values

cameraMat

matrix contain output from CAMERA

geneSets

a list of list of geneIds

verbose

Boolean whether to print out progress

Value

table of genes' mean expression values for each group

Examples

1
2
3
4
5
6
7
8
## Not run: 
vDat <- voom(exprs(eDat), design=designMat, plot=FALSE, lib.size=libNorm)
res <- camera(vDat, setsIndices, design=designMat, contrast=cons[i], sort=TRUE)
geneSets <- geneIds(getGmt(gmtFile))
n <- make_gsea_igraph(vDat, res, geneSets)
plot(n, edge.width=E(n)$overlap*.1, vertex.color="white", vertex.label.cex=0.5)

## End(Not run)

ptvan/Pmisc documentation built on Nov. 19, 2020, 10:27 p.m.