EnrichAB {MAGeCKFlute}R Documentation

Enrichment analysis for Positive and Negative selection genes

Description

Do enrichment analysis for selected genes, in which positive selection and negative selection are termed as GroupA and GroupB

Usage

EnrichAB(data, pvalue = 0.25, enrich_method = "ORT",
  organism = "hsa", limit = c(1, 120), filename = NULL,
  out.dir = ".", width = 6.5, height = 4, ...)

Arguments

data

A data frame.

pvalue

Pvalue cutoff.

enrich_method

One of "ORT"(Over-Representing Test), "GSEA"(Gene Set Enrichment Analysis), and "HGT"(HyperGemetric test).

organism

"hsa" or "mmu".

limit

A two-length vector (default: c(1, 120)), specifying the min and max size of pathways for enrichent analysis.

filename

Suffix of output file name.

out.dir

Path to save plot to (combined with filename).

width

As in ggsave.

height

As in ggsave.

...

Other available parameters in ggsave.

Value

A list containing enrichment results for each group genes. This list contains items four items, keggA, keggB, goA, goB. Four items are all list object, containing subitems of gridPlot and enrichRes. gridPlot is a ggplot object, and enrichRes is a enrichResult instance

Author(s)

Binbin Wang

See Also

EnrichSquare

Examples

data(mle.gene_summary)
# Read beta score from gene summary table in MAGeCK MLE results
dd = ReadBeta(mle.gene_summary)
data=ScatterView(dd, ctrlname = "dmso", treatname = "plx")$data
## Not run: 
  #GO and KEGG enrichment analysis
 enrich_result = EnrichAB(data, pvalue=0.05, organism="hsa")
  print(enrich_result$keggA$gridPlot)
  print(enrich_result$goA$gridPlot)

## End(Not run)


[Package MAGeCKFlute version 1.4.3 Index]