EnrichmentBrowser (version 1.0.3)

ebrowser: Seamless navigation through enrichment analysis results

Description

This is the all-in-one wrapper function to perform the standard enrichment analysis pipeline implemented in the EnrichmentBrowser package. Given flat gene expression data, the data is read in and subsequently subjected to chosen enrichment analysis methods.

The results from different methods can be combined and investigated in detail in the default browser.

Usage

ebrowser( meth, exprs, pdat, fdat,  gs, grn=NULL, perm=1000, alpha=0.05, beta=1, comb=FALSE, browse=TRUE, nr.show=-1, out.dir=NULL )

Arguments

meth

Enrichment analysis method. Currently, the following enrichment analysis methods are supported: ‘ora’, ‘safe’, ‘gsea’, ‘samgs’, ‘ggea’, ‘spia’, and ‘nea’. See ‘sbea’ and ‘nbea’ for details.

exprs

Expression matrix. A tab separated text file containing *normalized* expression values on a *log* scale. Columns = samples/subjects; rows = features/probes/genes; NO headers, row or column names. Supported data types are log2 counts (microarray single-channel), log2 ratios (microarray two-color), and log2-counts per million (RNA-seq logCPMs). See the limma's user guide http://www.bioconductor.org/packages/release/bioc/vignettes/limma/inst/doc/usersguide.pdf for definition and normalization of the different data types.

pdat

Phenotype data. A tab separated text file containing annotation information for the samples in either *two or three* columns. NO headers, row or column names. The number of rows/samples in this file should match the number of columns/samples of the expression matrix. The 1st colum is reserved for the sample IDs; The 2nd column is reserved for a *BINARY* group assignment. Use '0' and '1' for unaffected (controls) and affected (cases) sample class, respectively. For paired samples or sample blocks a third column is expected that defines the blocks.

fdat

Feature data. A tab separated text file containing annotation information for the features. Exactly *TWO* columns; 1st col = feature IDs; 2nd col = corresponding KEGG gene ID for each feature ID in 1st col; NO headers, row or column names. The number of rows/features in this file should match the number of rows/features of the expression matrix.

gs

Gene sets. Either a list of gene sets (vectors of KEGG gene IDs) or a text file in GMT format storing all gene sets under investigation.

grn

Gene regulatory network. Either an absolute file path to a tabular file or a character matrix with exactly *THREE* cols; 1st col = IDs of regulating genes; 2nd col = corresponding regulated genes; 3rd col = regulation effect; Use '+' and '-' for activation/inhibition.

perm

Number of permutations of the expression matrix to estimate the null distribution. Defaults to 1000.

alpha

Statistical significance level. Defaults to 0.05.

beta

Log2 fold change signifcance level. Defaults to 1 (2-fold).

comb

Logical. Should results be combined if more then one enrichment method is selected? Defaults to FALSE.

browse

Logical. Should results be displayed in the browser for interactive exploration? Defaults to TRUE.

nr.show

Number of gene sets to show. As default all statistical significant gene sets are displayed.

out.dir

Optional output directory all results will be written to.

Value

None, opens the browser to explore results.

See Also

read.eset to read expression data from file; probe.2.gene.eset to transform probe to gene level expression; get.kegg.genesets to retrieve gene set definitions from KEGG; compile.grn.from.kegg to construct a GRN from KEGG pathways; sbea to perform set-based enrichment analysis; nbea to perform network-based enrichment analysis; comb.ea.results to combine results from different methods; ea.browse for exploration of resulting gene sets

Examples

    # expression data from file
    exprs.file <- system.file("extdata/ALL_exprs.tab", package="EnrichmentBrowser")
    pdat.file <- system.file("extdata/ALL_pData.tab", package="EnrichmentBrowser")
    fdat.file <- system.file("extdata/ALL_fData.tab", package="EnrichmentBrowser")
    
    # getting all human KEGG gene sets
    # hsa.gs <- get.kegg.genesets("hsa")
    gs.file <- system.file("extdata/hsa_kegg_gs.gmt", package="EnrichmentBrowser")
    hsa.gs <- parse.genesets.from.GMT(gs.file)

    
    # set-based enrichment analysis
    ebrowser(   meth="ora", 
            exprs=exprs.file, 
            pdat=pdat.file, 
            fdat=fdat.file, 
            gs=hsa.gs)
    

    # compile a gene regulatory network from KEGG pathways
    # hsa.grn <- compile.grn.from.kegg("hsa")
    pwys <- system.file("extdata/hsa_kegg_pwys.zip", package="EnrichmentBrowser")
    hsa.grn <- compile.grn.from.kegg(pwys)
   
    
    # network-based enrichment analysis
    ebrowser(   meth="ggea", 
            exprs=exprs.file, 
            pdat=pdat.file, 
            fdat=fdat.file, 
            gs=hsa.gs,
            grn=hsa.grn)
    

    # combining results
    ebrowser(   meth=c("ora", "ggea"), 
            perm=100,
            exprs=exprs.file, 
            pdat=pdat.file, 
            fdat=fdat.file, 
            gs=hsa.gs,
            grn=hsa.grn)