FluteRRA {MAGeCKFlute}R Documentation

Downstream analysis based on MAGeCK-RRA result

Description

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

Usage

FluteRRA(gene_summary, sgrna_summary, lfcCutoff = c(-1, 1),
  organism = "hsa", limit = c(1, 120), pvalueCutoff = 0.25,
  prefix = "Test", width = 12, height = 6, outdir = ".")

Arguments

gene_summary

A file path or a data frame of gene summary data.

sgrna_summary

A file path or a data frame of sgRNA summary data.

lfcCutoff

A two-length vector (default: c(-1, 1)), specifying the logFC cutoff for negative selection and positive selection.

organism

"hsa" or "mmu".

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.

prefix

A character, indicating the prefix of output file name.

width

The width of summary pdf in inches.

height

The height of summary pdf in inches.

outdir

Output directory on disk.

Details

MAGeCK RRA allows for the comparison between two experimental conditions. It can identify genes and sgRNAs are significantly selected between the two conditions. The most important output of MAGeCK RRA is the file 'gene_summary.txt'. MAGeCK RRA will output both the negative score and positive score for each gene. A smaller score indicates higher gene importance. MAGeCK RRA will also output the statistical value for the scores of each gene. Genes that are significantly positively and negatively selected can be identified based on the p-value or FDR.

The downstream analysis of this function includes identifying positive and negative selection genes, and performing biological functional category analysis and pathway enrichment analysis of these genes.

Value

All of the pipeline results is output into the out.dir/prefix_Results, which includes a pdf file and a folder named 'RRA'.

Author(s)

Wubing Zhang

See Also

FluteMLE

Examples

data("rra.gene_summary")
data("rra.sgrna_summary")
## Not run: 
    # Run the FluteRRA pipeline
    FluteRRA(rra.gene_summary, rra.sgrna_summary, prefix="RRA", organism="hsa")

## End(Not run)



[Package MAGeCKFlute version 1.4.3 Index]