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", pathway_limit = c(3, 50), pvalueCutoff = 0.25,
      adjust = "BH", 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".
`pathway_limit` | A two-length vector (default: c(3, 50)), specifying the minimal and maximal size of gene sets for enrichent analysis.
`pvalueCutoff` | A numeric, specifying pvalue cutoff of enrichment analysis, default 1.
`adjust` | One of "holm", "hochberg", "hommel", "bonferroni", "BH", "BY", "fdr", "none".
`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)
    ## Not run:
        # Run the FluteRRA pipeline
        FluteRRA(rra.gene_summary, prefix="GSC", organism="hsa")

    ## End(Not run)

[Package MAGeCKFlute version 1.2.3 Index]