FluteMLE {MAGeCKFlute}R Documentation

Downstream analysis based on MAGeCK-MLE result

Description

Integrative analysis pipeline using the gene summary table in MAGeCK MLE results

Usage

FluteMLE(gene_summary, ctrlname, treatname, keytype = "Symbol",
  organism = "hsa", scale_cutoff = 2, top = 10, bottom = 10,
  interestGenes = NA, limit = c(1, 120), pvalueCutoff = 0.25,
  enrich_kegg = "ORT", posControl = NULL, loess = FALSE,
  prefix = "", width = 10, height = 7, outdir = ".",
  view_allpath = FALSE)

Arguments

gene_summary

A data frame, which contains columns of 'Gene', ctrlname.beta and treatname.beta.

ctrlname

A character vector, specifying the names of control samples.

treatname

A character vector, specifying the names of treatment samples.

keytype

"Entrez" or "Symbol".

organism

"hsa" or "mmu".

scale_cutoff

Boolean or numeric, whether scale cutoff to whole genome level, or how many standard deviation will be used as cutoff.

top

An integer, specifying number of top selected genes to be labeled in rank figure.

bottom

An integer, specifying number of bottom selected genes to be labeled in rank figure.

interestGenes

A character vector, specifying interested genes to be labeled in rank figure.

limit

A two-length vector (default: c(1, 120)), specifying the minimal and maximal size of gene sets for enrichent analysis.

pvalueCutoff

A numeric, specifying pvalue cutoff of enrichment analysis, default 1.

enrich_kegg

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

posControl

A character vector, specifying a list of positive control gene symbols.

loess

Boolean, whether include loess normalization in the pipeline.

prefix

A character, indicating the prefix of output file name, which can't contain special characters.

width

The width of summary pdf in inches.

height

The height of summary pdf in inches.

outdir

Output directory on disk.

view_allpath

Boolean, whether output all pathway view figures.

Details

MAGeCK-MLE can be used to analyze screen data from multi-conditioned experiments. MAGeCK-MLE also normalizes the data across multiple samples, making them comparable to each other. The most important ouput of MAGeCK MLE is 'gene_summary' file, which includes the beta scores of multiple conditions and the associated statistics. The 'beta score' for each gene describes how the gene is selected: a positive beta score indicates a positive selection, and a negative beta score indicates a negative selection.

The downstream analysis includes identifying essential, non-essential, and target-associated genes, and performing biological functional category analysis and pathway enrichment analysis of these genes. The function also visualizes genes in the context of pathways to benefit users exploring screening data.

Value

All of the pipeline results is output into the out.dir/prefix_Results, which includes a pdf file and many folders. The pdf file 'prefix_Pipeline_results.pdf' is the summary of pipeline results. For each section in this pipeline, figures and useful data are outputed to corresponding subfolders.

Author(s)

Wubing Zhang

See Also

FluteRRA

Examples

data(mle.gene_summary)
## Not run: 
  # functional analysis for MAGeCK MLE results
  FluteMLE(mle.gene_summary, ctrlname = "dmso", treatname = "plx",
           prefix = "PLX", pvalueCutoff = 0.25, organism = "hsa")

## End(Not run)


[Package MAGeCKFlute version 1.4.3 Index]